PLEASE EXPLAIN: why reactingFoam is not outputting density and viscosity

sahmed · December 7, 2016, 15:35

Dear all
I am using reactingFoam for one of my flame simulations. But I dont know why it is not spitting out properties like density and viscosity as an output. The density should be a function of P and T from ideal gas law and viscosity as a function of temperature from Sutherland law (used in this version of the solver). Could anyone PLEASE explain the reason why it is not spitting out in the output and how can I tell the solver to spit it out.
Your kind help is highly appreciated.
Thanks

Verse · December 7, 2016, 15:48

Check your createFields.H file, and for density, there should be something like:

Code:

volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::AUTO_READ,
        IOobject::NO_WRITE
    ),
    thermo.rho()
);

Change the NO_WRITE to AUTO_WRITE. That should solve the density problem. Can't comment too much about viscosity myself.

sahmed · December 7, 2016, 17:31

Thanks a lot Dominic for your prompt reply.
Anyone has any idea of the viscosity to spit out? Can I use the same procedure as Dominic mentioned to change in createFields.H??

December 7, 2016, 15:35	PLEASE EXPLAIN: why reactingFoam is not outputting density and viscosity	#1
sahmed Member Sheikh Ahmed Join Date: Dec 2015 Location: South Carolina, USA Posts: 88 Rep Power: 11	Dear all I am using reactingFoam for one of my flame simulations. But I dont know why it is not spitting out properties like density and viscosity as an output. The density should be a function of P and T from ideal gas law and viscosity as a function of temperature from Sutherland law (used in this version of the solver). Could anyone PLEASE explain the reason why it is not spitting out in the output and how can I tell the solver to spit it out. Your kind help is highly appreciated. Thanks

December 7, 2016, 15:48		#2
Verse New Member Dominic Join Date: May 2016 Posts: 27 Rep Power: 10	Check your createFields.H file, and for density, there should be something like: Code: volScalarField rho ( IOobject ( "rho", runTime.timeName(), mesh, IOobject::AUTO_READ, IOobject::NO_WRITE ), thermo.rho() ); Change the NO_WRITE to AUTO_WRITE. That should solve the density problem. Can't comment too much about viscosity myself.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Density correction function for turbulent viscosity	Bando	STAR-CCM+	1	July 9, 2017 03:44
How to find density without having dynamic viscosity?	ooo	OpenFOAM Running, Solving & CFD	6	January 28, 2014 13:46
Variable density and viscosity	ravibhadauria	FLUENT	0	September 15, 2011 18:09
icoFoam >> density >> kinematic viscosity for air 20°C	mgolbs	OpenFOAM Pre-Processing	11	February 15, 2010 11:09
REAL GAS UDF	brian	FLUENT	6	September 11, 2006 09:23

December 7, 2016, 17:31		#3
sahmed Member Sheikh Ahmed Join Date: Dec 2015 Location: South Carolina, USA Posts: 88 Rep Power: 11	Thanks a lot Dominic for your prompt reply. Anyone has any idea of the viscosity to spit out? Can I use the same procedure as Dominic mentioned to change in createFields.H??