small density change for isothermal flows

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

 April 6, 2013, 10:02 small density change for isothermal flows #1 Senior Member   Join Date: Nov 2012 Posts: 168 Rep Power: 7 Hi All, When I simulate the isothermal flows, I used the rhoPimeplFoam. I also made some modification on that solver (1) comment the statements for solving the enthalpy equation but retain therm.correct() (2) use Low Ma number assumption, pressure for the EoS is fixed to be 101325 Pa. However, when I use this solver to simulate the isothermal flows, the density shows very small change. I attached an example: an isothermal air jet, uj=7m/s, temperature is always 294K in the initial field and during the simulation. However, as I said, the density in the domain varies between 1.197 and 1.188. This is the density for air. In Openfoam, the density is updated through rho=p*psi, where psi=1/(R*T). Along this line, it seems that the density change should be from R because in this case p and T are fixed. But I do not know why R is changed since I use pureMixture , which should be also fixed. The complete options for thermalphysical models I used are: thermoType hPsiThermo>>>>; Another possible explanation is that this small variation of density is completely caused by numerical error, which should be tolerable. I am a little confused about this problem. Thank you in advance if someone can give me comments on this. Thanks.

April 6, 2013, 10:05
#2
Senior Member

Join Date: Nov 2012
Posts: 168
Rep Power: 7
Attach a contour of density here. Thanks.

Quote:
 Originally Posted by hz283 Hi All, When I simulate the isothermal flows, I used the rhoPimeplFoam. I also made some modification on that solver (1) comment the statements for solving the enthalpy equation but retain therm.correct() (2) use Low Ma number assumption, pressure for the EoS is fixed to be 101325 Pa. However, when I use this solver to simulate the isothermal flows, the density shows very small change. I attached an example: an isothermal air jet, uj=7m/s, temperature is always 294K in the initial field and during the simulation. However, as I said, the density in the domain varies between 1.197 and 1.188. This is the density for air. In Openfoam, the density is updated through rho=p*psi, where psi=1/(R*T). Along this line, it seems that the density change should be from R because in this case p and T are fixed. But I do not know why R is changed since I use pureMixture , which should be also fixed. The complete options for thermalphysical models I used are: thermoType hPsiThermo>>>>; Another possible explanation is that this small variation of density is completely caused by numerical error, which should be tolerable. I am a little confused about this problem. Thank you in advance if someone can give me comments on this. Thanks.
Attached Images
 air_jets.jpg (20.6 KB, 7 views)

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Confused student. Main CFD Forum 27 March 18, 2017 13:25 thiez03 CFX 4 April 5, 2012 17:04 kel85uk OpenFOAM 1 July 22, 2011 04:01 Astrid CFX 5 April 26, 2002 05:47 Elliot Schwartz Main CFD Forum 6 May 28, 1999 02:11

All times are GMT -4. The time now is 03:15.