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

sonicFoam #include[s] rhoEqn.H...

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

Like Tree2Likes
  • 1 Post By l_r_mcglashan
  • 1 Post By JasonWang3

Reply
 
LinkBack Thread Tools Display Modes
Old   January 29, 2010, 11:21
Default sonicFoam #include[s] rhoEqn.H...
  #1
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 9
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
...but the rhoEqn.H cannot be found in the solver directory (version 1.6). I was curious when I was using the solver before, it gave results for everything, but the rho was always calculated as 0. Does anybody know the purpose of this? Is this intentional? Can I incorporate some other rhoEqn.H easily?
thekay is offline   Reply With Quote

Old   January 29, 2010, 12:55
Default
  #2
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 15
l_r_mcglashan will become famous soon enough
It's in the folder OpenFOAM-1.6.x/src/finiteVolume/cfdTools/compressible.
Lisandro Maders likes this.
__________________
Laurence R. McGlashan :: Website
l_r_mcglashan is offline   Reply With Quote

Old   February 15, 2010, 10:14
Default
  #3
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 9
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
Revisiting this thread after some time:

I am running a case at the moment, using sonicFoam. The case is pretty much the same as /home/username/OpenFOAM/OpenFOAM-1.6/tutorials/compressible/sonicFoam/ras/prism/ but instead of "prismWall", I have "wall" which is the wall of my domain. While running, in each time step, I get:

Code:
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Obviously I don't have a rho file in 0/ but even if I do, I get the same output for rho.

My questions are:

  1. Can I get sonicFoam to solve for rho as well?
  2. Why is the above output produced? In no tutorial it solves for rho, so I am curious if it was just left there from another solver.
Thank you.
thekay is offline   Reply With Quote

Old   April 26, 2010, 13:57
Default
  #4
Senior Member
 
KGN
Join Date: Oct 2009
Location: Chennai, India
Posts: 102
Rep Power: 9
mecbe2002 is on a distinguished road
Hi,

Any solution for this problem.
mecbe2002 is offline   Reply With Quote

Old   October 30, 2013, 05:17
Default
  #5
New Member
 
Richard Mathie
Join Date: Oct 2012
Posts: 3
Rep Power: 6
RicMat is on a distinguished road
Probably you need to supply an equation of state in your case dir? that is dont you need a thermoPhysicalProperties dict? I cant say for shure as I am looking at engineFoam and rhoPimpleFOAM at the moment, but that seems to be what's happening there. rhoEqn.h is just the mass conservation equation.
RicMat is offline   Reply With Quote

Old   March 18, 2015, 12:58
Default
  #6
Member
 
Xinguang Wang
Join Date: Feb 2015
Posts: 36
Rep Power: 3
JasonWang3 is on a distinguished road
Quote:
Originally Posted by thekay View Post
Revisiting this thread after some time:

I am running a case at the moment, using sonicFoam. The case is pretty much the same as /home/username/OpenFOAM/OpenFOAM-1.6/tutorials/compressible/sonicFoam/ras/prism/ but instead of "prismWall", I have "wall" which is the wall of my domain. While running, in each time step, I get:

Code:
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
Obviously I don't have a rho file in 0/ but even if I do, I get the same output for rho.

My questions are:

  1. Can I get sonicFoam to solve for rho as well?
  2. Why is the above output produced? In no tutorial it solves for rho, so I am curious if it was just left there from another solver.
Thank you.
I met the same situation. sonicFoam is a pressure based solver, and the element of density don't need to solve alone. So that's maybe the reason of rho=0.
But still don't know how to get the rho values of the fields, two options:
1. state equation
2. phi=rho*velocity*A
JasonWang3 is offline   Reply With Quote

Old   October 2, 2015, 10:18
Default Solved?
  #7
Member
 
Lisandro Maders
Join Date: Feb 2013
Posts: 67
Rep Power: 5
Lisandro Maders is on a distinguished road
Do you guys solved this issue? Why is rho equation not being calculated?
Lisandro Maders is offline   Reply With Quote

Old   October 5, 2015, 06:27
Default
  #8
Member
 
Xinguang Wang
Join Date: Feb 2015
Posts: 36
Rep Power: 3
JasonWang3 is on a distinguished road
Hi Maders
Try to add these extra sentence at the end of controlDict, and the rho will output in the calculation folder.
Quote:
functions
{
rhofunc
{
type writeRegisteredObject;
functionObjectLibs ("libIOFunctionObjects.so");
outputControl outputTime;
outputInterval 1;
objectNames
(
"rho"
);
}
}
kcp18 likes this.
JasonWang3 is offline   Reply With Quote

Old   October 5, 2015, 07:13
Default Worked!
  #9
Member
 
Lisandro Maders
Join Date: Feb 2013
Posts: 67
Rep Power: 5
Lisandro Maders is on a distinguished road
Quote:
Originally Posted by JasonWang3 View Post
Hi Maders
Try to add these extra sentence at the end of controlDict, and the rho will output in the calculation folder.
Hi Wang,

that worked! I have already done it by inserting write() functions at runTime, but I found your way more straightforward.. Thanks!
Lisandro Maders is offline   Reply With Quote

Old   October 5, 2015, 07:18
Default Actually the answer is..
  #10
Member
 
Lisandro Maders
Join Date: Feb 2013
Posts: 67
Rep Power: 5
Lisandro Maders is on a distinguished road
Actually the main answer iin this topic was about the rho calculation.. That seemed not to be doing by the "0" values of iterations and residuals. Here is the aswer, which is very logic by the way: rho euation in rhoPisoFoam

Regards, Lisandro
Lisandro Maders is offline   Reply With Quote

Old   June 8, 2016, 17:45
Default Additional note
  #11
New Member
 
Kristopher Pierson
Join Date: Jun 2016
Posts: 2
Rep Power: 0
kcp18 is on a distinguished road
I would just like to add to this that by adding the rho function I was able to use the R utility to compute Reynolds stress when using the sonicFoam solver. This was very helpful if you are using some other turbulence model at the start of the solution with the intent to switch to LRR model after the flow field has been initialized.
kcp18 is offline   Reply With Quote

Reply

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
Setting a case with sonicFoam achilleslaststand OpenFOAM Pre-Processing 0 December 4, 2009 12:29
define walls in sonicFoam shangzung OpenFOAM 0 October 21, 2009 08:18
Different response from 1.5 and 1.6 version of sonicFoam solver nishant_hull OpenFOAM Running, Solving & CFD 1 October 12, 2009 13:45
energy in sonicFoam joern OpenFOAM Running, Solving & CFD 0 June 16, 2009 15:53
SonicFoam forwardStepTutorial doesnbt complete the run alberto OpenFOAM Bugs 1 June 10, 2007 15:35


All times are GMT -4. The time now is 18:19.