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

Unable to validate results with sonicFoam and rhoCentralFoam

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

Like Tree2Likes
  • 1 Post By schuyler
  • 1 Post By schuyler

Reply
 
LinkBack Thread Tools Display Modes
Old   June 9, 2016, 16:15
Default Unable to validate results with sonicFoam and rhoCentralFoam
  #1
Senior Member
 
Robert
Join Date: Sep 2010
Posts: 158
Rep Power: 8
lordvon is on a distinguished road
What I thought was going to be a simulation of trivial physics, was in fact something the supersonic solvers in OpenFOAM could not handle. I was trying to validate an axisymmetric simulation of a converging-diverging nozzle with the data found here:

http://www.grc.nasa.gov/WWW/wind/valid/cdv/cdv.html

The pressure ratio of 0.89 (subsonic, peaking at 0.8+ Mach) simply could not match; peaks were much too low by about 10%, when it should be almost exactly matching. I tried different numerical schemes, grid resolutions, to no avail. The pressure ratio of 0.16 (supersonic, no shocks in throat) actually predicted a shock right at the exit, causing the predicted pressures to be much too high and the Mach numbers much too low.

So you can observe this for yourself, the boundary condition directories that I used can be found in 'case_sonic/0' and 'case_rhocentral/0' at the following github page:

https://github.com/rlee32/upcoming/tree/master/CDNozzle

Let me know what you think...


[Moderator note: For future reference, the original title was: PSA: sonicFoam and rhoCentralFoam are bunk...]

Last edited by wyldckat; September 9, 2016 at 18:13. Reason: see "Moderator note:"
lordvon is offline   Reply With Quote

Old   June 18, 2016, 16:57
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,860
Blog Entries: 39
Rep Power: 108
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick note: I haven't checked the case set-up you have (and this topic is beyond my expertises ), but there have been others who have stumbled on this issue. A few examples:
wyldckat is offline   Reply With Quote

Old   July 14, 2016, 14:00
Default
  #3
New Member
 
W. Schuyler Hinman
Join Date: Apr 2013
Location: Calgary, Alberta, Canada
Posts: 26
Rep Power: 6
schuyler is on a distinguished road
My experience with rhoCentralFoam has been very positive! I believe the problem you are experiencing is due to your setup. I think the main issue is that the totalPressure boundary condition requires a definition of psi for compressible flow.

Since this is a compressible flow the definition of psi should be
Code:
thermo:psi
calling psi from the thermophysical models set for the compressible flow.

With regards to the supersonic flow with an unwanted shock, you should not specify an outlet pressure for this case. You CAN. But remember that the supersonic flow solution for this problem is unique (contrasting to flow with shock, or subsonic where the solution is not unique). It does not require this information. All you are doing is making the simulation more rigid. Unless you input the exact answer that the solver is looking for (and this will depend slightly on schemes, resolution etc) you will have some difficulties. The shock formed because of this probably. I suggest using zeroGradient for the outlet pressure. Once you get the answer, how closely the exit pressure matches the analytical, expected result, will be an indication of the quality of the solution.

I started a blog on CFD about 6 months ago. I just made a post regarding this test case:

https://curiosityfluids.com/2016/07/...hocentralfoam/

I am still building up the content. But now this case is covered! Hopefully this helps.

Last edited by schuyler; July 14, 2016 at 14:43. Reason: thermo:psi giving smiley face
schuyler is offline   Reply With Quote

Old   July 15, 2016, 22:36
Default
  #4
Senior Member
 
Robert
Join Date: Sep 2010
Posts: 158
Rep Power: 8
lordvon is on a distinguished road
Great job schuyler! I wish I could change the title now, but it seems I cannot.
lordvon is offline   Reply With Quote

Old   September 1, 2016, 10:03
Default
  #5
Member
 
Bruno Blais
Join Date: Sep 2013
Location: Canada
Posts: 63
Rep Power: 5
blais.bruno is on a distinguished road
Dear Schuyler,
Do you have the mesh and files for the tests cases you discuss on your blog post?

Best regards,
Bruno

Quote:
Originally Posted by schuyler View Post
My experience with rhoCentralFoam has been very positive! I believe the problem you are experiencing is due to your setup. I think the main issue is that the totalPressure boundary condition requires a definition of psi for compressible flow.

Since this is a compressible flow the definition of psi should be
Code:
thermo:psi
calling psi from the thermophysical models set for the compressible flow.

With regards to the supersonic flow with an unwanted shock, you should not specify an outlet pressure for this case. You CAN. But remember that the supersonic flow solution for this problem is unique (contrasting to flow with shock, or subsonic where the solution is not unique). It does not require this information. All you are doing is making the simulation more rigid. Unless you input the exact answer that the solver is looking for (and this will depend slightly on schemes, resolution etc) you will have some difficulties. The shock formed because of this probably. I suggest using zeroGradient for the outlet pressure. Once you get the answer, how closely the exit pressure matches the analytical, expected result, will be an indication of the quality of the solution.

I started a blog on CFD about 6 months ago. I just made a post regarding this test case:

https://curiosityfluids.com/2016/07/...hocentralfoam/

I am still building up the content. But now this case is covered! Hopefully this helps.
blais.bruno is offline   Reply With Quote

Old   September 2, 2016, 18:11
Default
  #6
New Member
 
W. Schuyler Hinman
Join Date: Apr 2013
Location: Calgary, Alberta, Canada
Posts: 26
Rep Power: 6
schuyler is on a distinguished road
Hi Bruno,

I will update the blog this weekend and upload the files! I usually do... dont really know why I forgot this time.

Schuyler

Quote:
Originally Posted by blais.bruno View Post
Dear Schuyler,
Do you have the mesh and files for the tests cases you discuss on your blog post?

Best regards,
Bruno
blais.bruno likes this.
schuyler is offline   Reply With Quote

Old   September 6, 2016, 17:38
Default
  #7
New Member
 
W. Schuyler Hinman
Join Date: Apr 2013
Location: Calgary, Alberta, Canada
Posts: 26
Rep Power: 6
schuyler is on a distinguished road
The blog post is updated now with a downloadable link.

Quote:
Originally Posted by blais.bruno View Post
Dear Schuyler,
Do you have the mesh and files for the tests cases you discuss on your blog post?

Best regards,
Bruno
blais.bruno likes this.
schuyler is offline   Reply With Quote

Old   December 5, 2016, 15:16
Default
  #8
Member
 
Bruno Blais
Join Date: Sep 2013
Location: Canada
Posts: 63
Rep Power: 5
blais.bruno is on a distinguished road
Hello!
Thank you very much for putting the files online.
However I have a question. I am trying to reproduce the subsonic and the sonic (well the case where there is a shock wave within the nozzle) and I am getting confusing results. How do you set the pressure outlet boundary? do you set it as a totalPressure boundary condition (specifying thermosi and etc.) or as a regular pressure BC with fixedValue?

I think I am failing to understand the maths and physics behind the totalPressure BC.
Sorry, I am not too familiar with compressible flows....


Quote:
Originally Posted by schuyler View Post
The blog post is updated now with a downloadable link.
blais.bruno is offline   Reply With Quote

Old   December 5, 2016, 15:25
Default
  #9
New Member
 
W. Schuyler Hinman
Join Date: Apr 2013
Location: Calgary, Alberta, Canada
Posts: 26
Rep Power: 6
schuyler is on a distinguished road
Quote:
Originally Posted by blais.bruno View Post
Hello!
Thank you very much for putting the files online.
However I have a question. I am trying to reproduce the subsonic and the sonic (well the case where there is a shock wave within the nozzle) and I am getting confusing results. How do you set the pressure outlet boundary? do you set it as a totalPressure boundary condition (specifying thermosi and etc.) or as a regular pressure BC with fixedValue?

I think I am failing to understand the maths and physics behind the totalPressure BC.
Sorry, I am not too familiar with compressible flows....
For the outlet boundary, I set the pressure as a fixedValue. You can see how the outlet pressure dictates the location of a normal shock in a nozzle here:

https://curiosityfluids.com/2016/03/...l-shock-waves/

also, an alternative case to the one discussed in this feed is posted here:

https://curiosityfluids.com/2016/04/...hocentralfoam/

Schuyler
schuyler 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
Which equations are solved by rhoCentralFoam, rhoSimplecFoam and sonicFoam? gamma_user OpenFOAM Verification & Validation 1 April 19, 2015 12:13
sonicFoam vs rhoCentralFoam forwardStep tutorial hk318i OpenFOAM Running, Solving & CFD 2 March 12, 2015 12:15
rhoCentralFoam / sonicFoam VSass OpenFOAM Running, Solving & CFD 1 June 26, 2013 12:39
is sonicFoam better in convergence than rhoPimplefoam or rhoCentralFoam? immortality OpenFOAM Running, Solving & CFD 1 February 11, 2013 23:26


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