Blog : Noms de champs dans OpenFOAM(R)

Récemment, un problème d’unités est survenu :

--> FOAM FATAL ERROR: (openfoam-2312)
Different dimensions for '(a = b)'
     dimensions : [0 2 -1 0 0 0 0] != [1 -1 -1 0 0 0 0]

Manifestement, il y avait un problème lié à la division par la masse volumique. L’erreur est localisée dans :

#5  Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::multiphaseInter::multiphaseSystem> > >::correctNut()

Il semblait donc que le problème venait des unités de nut. Pourtant, après une analyse plus approfondie, le problème provenait du champ alphat :

template<class BasicTurbulenceModel>
void Foam::EddyDiffusivity<BasicTurbulenceModel>::correctNut()
{
    // Read Prt if provided
    Prt_ = dimensionedScalar("Prt", dimless, 1.0, this->coeffDict());
    alphat_ = this->rho_*this->nut()/Prt_;
    alphat_.correctBoundaryConditions();
}

Cependant, si l’on considère nut comme la diffusivité de la quantité de mouvement et alphat comme la diffusivité thermique, on peut par extension qualifier alphat de “nut thermique”.