rho in basicPsiThermo

megacrout · January 30, 2012, 12:47

Hi,

1) Why does basicRhoThermo possess a "non-const access to the local density field", i.e.

Quote:

virtual volScalarField& rho();

but not basicPsiThermo?

2) I need to access rho as a non-tmp object using basicPsiThermo (need Janaf and Sutherland). How do I have modify this utility in order to make it work?

Thx.

Tibo

Chris Lucas · January 31, 2012, 04:22

Hi,

About question no. 1)

Because basicRhoThermo has an internal rho volScalarField and basicPsiThermo doesn't. In basicPsiThermo rho is calculated by psi*p.

About question no. 2).

I'm not sure about janaf, but I can use sutherland in combination with basicPsiThermo. The way janaf works (not 100% sure), it should work as well in combination with basicPsiThermo. None of the two models need the density field.

Best Regards,
Christian

megacrout · February 2, 2012, 04:25

Hey Chris,

Thanks for your help.

1) I obviously noticed that rho is not computed the same, but I was confused noticing that this influences the way it is defined and the way it can be used. (And still am a bit).

2) I guess my question was not detailed enough. I know that Janaf AND Sutherland models work in combination with basicPsiThermo (s. e.g. reactingFoam) and that is precisely why I want to use this thermo package.
But I need to use and modify the value of thermo.rho() in the pEqn.H file and thought I needed to change the way it is defined to do so. I found an alternative solution though, as explained in the original thread about that topic: http://www.cfd-online.com/Forums/ope...-operator.html (see post #15).

Regards.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
rhoSimplecFoam, solving rho	j-blindi	OpenFOAM Running, Solving & CFD	4	October 14, 2011 05:54
Low Mach number Compressible jet flow using LES	ankgupta8um	OpenFOAM Running, Solving & CFD	7	January 15, 2011 14:38
rho euation in rhoPisoFoam	nileshjrane	OpenFOAM Running, Solving & CFD	4	October 5, 2010 08:44
Error while running rhoPisoFoam..	nileshjrane	OpenFOAM Running, Solving & CFD	8	August 26, 2010 13:50
Error log	vw.cfd	OpenFOAM	6	August 7, 2009 06:44

January 31, 2012, 04:22		#2
Chris Lucas Senior Member Christian Lucas Join Date: Aug 2009 Location: Braunschweig, Germany Posts: 202 Rep Power: 18	Hi, About question no. 1) Because basicRhoThermo has an internal rho volScalarField and basicPsiThermo doesn't. In basicPsiThermo rho is calculated by psi*p. About question no. 2). I'm not sure about janaf, but I can use sutherland in combination with basicPsiThermo. The way janaf works (not 100% sure), it should work as well in combination with basicPsiThermo. None of the two models need the density field. Best Regards, Christian

February 2, 2012, 04:25		#3
megacrout Member Tibo Join Date: Jun 2011 Posts: 68 Rep Power: 15	Hey Chris, Thanks for your help. 1) I obviously noticed that rho is not computed the same, but I was confused noticing that this influences the way it is defined and the way it can be used. (And still am a bit). 2) I guess my question was not detailed enough. I know that Janaf AND Sutherland models work in combination with basicPsiThermo (s. e.g. reactingFoam) and that is precisely why I want to use this thermo package. But I need to use and modify the value of thermo.rho() in the pEqn.H file and thought I needed to change the way it is defined to do so. I found an alternative solution though, as explained in the original thread about that topic: http://www.cfd-online.com/Forums/ope...-operator.html (see post #15). Regards.