OpenFOAM® versão 2.0.0: mudanças na biblioteca thermophysicalModels

No início do mês passado a OpenCFD lançou a versão 2.0 do OpenFOAM. A própria OpenCFD disponibilizou uma lista contendo as principais mudanças na nova versão. A proposta deste post é ressaltar aquela que eu considero a principal mudança na biblioteca thermophysicalModels: a nova sintaxe do arquivo thermophysicalProperties.

A biblioteca thermophysicalModels é responsável pela definição dos modelos para o cálculo de propriedades termofísicas. Propriedades como a massa específica, a viscosidade e o calor específico, por exemplo, dependem da  temperatura e da pressão e, por este motivo, o usuário precisa especificar os modelos para o cálculo destas propriedades durante a simulação do seu caso.

Essa é função do arquivo thermophysicalProperties na pasta constant dos casos e a biblioteca thermophysicalModels estabelece as regras de entrada de dados neste arquivo. A principal mudança na versão 2.0 esta justamente na forma de entrada de dados no aquivo thermophysicalProperties. Portanto, se você usa algum solver contendo este arquivo, fique atento, pois o seu caso criado em uma versão 1.X não irá funcionar na versão 2.0! Será necessário atualizar a forma de entrada de dados do arquivo thermophysicalProperties!


A boa notícia é que essa mudança transformou a leitura de dados muito mais palatável para um ser humano! Antes o arquivo thermophysicalProperties se resumia a um conjunto de números e nomes de modelos que, a principio, não fazia o menor sentido para um primeiro leitor! Agora fica mais fácil identificar os parâmetros relacionados a cada modelo e o impacto para um usuário novato é menor! Mas como o mundo não é perfeito, o usuário ainda precisa consultar o User's Guide para saber o que cada nome de modelo especificado no thermoType significa na prática.


Somente para ilustrar a diferença na forma de entrada de dados, um trecho do arquivo thermophysicalProperties  do solver  XiFoam tinha a seguinte estrutura:

reactants       reactants 24.8095 29.4649 200 5000 1000 3.28069 0.00195035 -6.53483e-07 1.00239e-10 -5.64653e-15 -1609.55 4.41496 3.47696 0.000367499 1.84866e-06 -9.8993e-10 -3.10214e-14 -1570.81 3.76075 1.67212e-06 170.672;

Agora, na versão 2.0, esta mesma informação é entrada da seguinte forma:

reactants
{
    specie
    {
        nMoles          24.8095;
        molWeight       29.4649;
    }
    thermodynamics
    {
        Tlow            200;
        Thigh           5000;
        Tcommon         1000;
        highCpCoeffs    ( 3.28069 0.00195035 -6.53483e-07 1.00239e-10 -5.64653e-15 -1609.55 4.41496 );
        lowCpCoeffs     ( 3.47696 0.000367499 1.84866e-06 -9.8993e-10 -3.10214e-14 -1570.81 3.76075 );
    }
    transport
    {
        As              1.67212e-06;
        Ts              170.672;
    }
}

Muito melhor, não!? Agora resta a expectativa que a sintaxe mais palatável seja estendida também para a definição dos modelos no thermoType! Tomara que não demore muito!

Até a próxima.

Comments

Popular posts from this blog

Portal CFD Brasil

Gnuplot: Tipos de pontos e linhas

Avaliação da qualidade da malha