CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   how to develope a solver (http://www.cfd-online.com/Forums/openfoam/117226-how-develope-solver.html)

AmirBaqa1987 May 5, 2013 04:01

how to develope a solver
 
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 May 5, 2013 05:13

new volScalarField nu
 
Hi

stick to this howTo:
HTML Code:

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:


Code:

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

Code:

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


Code:

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



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

Code:

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


This should be all.

Regards

Fabian

AmirBaqa1987 May 6, 2013 02:57

appreciation
 
Hallo Fabian

Vielen Dank.


All times are GMT -4. The time now is 20:17.