CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

rho in basicPsiThermo

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 30, 2012, 12:47
Default rho in basicPsiThermo
  #1
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 6
megacrout is on a distinguished road
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
megacrout is offline   Reply With Quote

Old   January 31, 2012, 04:22
Default
  #2
Senior Member
 
Christian Lucas
Join Date: Aug 2009
Location: Braunschweig, Germany
Posts: 198
Rep Power: 7
Chris Lucas is on a distinguished road
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
Chris Lucas is offline   Reply With Quote

Old   February 2, 2012, 04:25
Default
  #3
Member
 
Tibo
Join Date: Jun 2011
Posts: 68
Rep Power: 6
megacrout is on a distinguished road
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: error: no match for ´operator+=´ (see post #15).

Regards.
megacrout is offline   Reply With Quote

Reply

Tags
basicpsithermo, basicrhothermo, janaf, sutherland

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
rhoSimplecFoam, solving rho j-blindi OpenFOAM Running, Solving & CFD 4 October 14, 2011 04: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 07:44
Error while running rhoPisoFoam.. nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50
Error log vw.cfd OpenFOAM 6 August 7, 2009 05:44


All times are GMT -4. The time now is 14:54.