[AmirBaqa1987]

Hi everyone,
i want to develop a solver by adding energy equation
to momentum equation(like adding temperature to icoFoam solver) but the difference is that the variables are coupled in these two equations for instance, viscosity in momentum equation varies with temperature in energy equation.
what should i do?

any help is appreciated.

Arjang

[fabian_roesler]

Hi

stick to this howTo:

http://openfoamwiki.net/index.php/How_to_add_temperature_to_icoFoam

1. You can add a new scalarField nu for the kinematic viscosity to the createFields.H like this:

volScalarField nu (     IOobject     (          "nu",          runTime.timeName(),          mesh,          IOobject::NO_READ,          IOobject::AUTO_WRITE      ),      nu0*(scalar(1)-nuk*(T-T0) );

2. Read the constants of your viscosity function. For example nu0, nuk and T0 for a linear function

nu = nu0*(scalar(1)-nuk*(T-T0));

dimensionedScalar nu0 (      transportProperties.lookup("nu0") );
dimensionedScalar nuk (      transportProperties.lookup("nuk") );
...

3. Update your viscosity function after you solved for temperature in the TEqn.H

nu = nu0*(scalar(1)-nuk*(T-T0));

This should be all.

Regards

Fabian

[AmirBaqa1987]

Hallo Fabian

Vielen Dank.