[chicky]

Hello everyone, I am new to cfd especially to openfoam. I try to implement a solver Tsimplefoam wich is simplefoam with temperature fields.

When I discover the file simpleFoam.C I see that the transportProperties was not called. Does it mean that simpleFoam doesn't use nu ? I find it very strange because Navier-Stoke steady-state still have the viscosity parameter.

[Swagga5aur]

Hello Nicky,
Well the simpleFoam solver does include the nu property from the transportProperties, however, not everything of the solver is defined in the simpleFoam.C folder or createFields. If you consider the options file in the make folder, -I$(LIB_SRC)/transportModels/incompressible/singlePhaseTransportModel \ is included and if you look at the singlePhaseTransportModel.H the following piece of code secures the loading of the nu viscosity:

    IOdictionary
    (
        IOobject
        (
            "transportProperties",
            U.time().constant(),
            U.db(),
            IOobject::MUST_READ_IF_MODIFIED,
            IOobject::NO_WRITE
        )
    ),
    viscosityModelPtr_(viscosityModel::New("nu", *this, U, phi))

which is quite similar to the code in the createFields dictionary.

If you want a varying nu model as a function of the temperature you could implement it yourself.

I have attached a TemperatureSimpleFoam solver and a viscosity model, which varies with temperature depending on the specified coefficients.

Hope its of any use and let me know if you have any questions.

[chicky]

Thank you very much for your reply. I understand better how it works. Implementing nu as a function of temperature is a great idea to train my self with openFOAM. I will try to do it with your file as a model.

Thank you it is great

[fede32]

I have found different results in simulations using in transportProperties:

nu              [0 2 -1 0 0 0 0] 4.43e-5

and

nu              nu[0 2 -1 0 0 0 0] 4.43e-5

why could be that?

[Swagga5aur]

Could you perhaps elaborate upon this or attach your case and describe how you have tested this?

Its hard to determine what the reason for this may be without additional information.

Regards Lasse

[fede32]

Originally Posted by Swagga5aur

Could you perhaps elaborate upon this or attach your case and describe how you have tested this?

Its hard to determine what the reason for this may be without additional information.

Regards Lasse

Sorry, I found that the error was to write nu glued to the bracket, this way:
nu nu[...]
It seems openfoam couldn''t read that and so apply a default value.
Thanks!

[Swagga5aur]

Ah good didn't even notice

Glad that it works now!