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

OpenFOAM supersonic Nozzle rhoCentralFoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 13, 2023, 06:38
Lightbulb OpenFOAM supersonic Nozzle rhoCentralFoam
  #1
New Member
 
Baptiste
Join Date: Mar 2023
Posts: 25
Rep Power: 3
K-al-Eps-o is on a distinguished road
Hi.

I am simulating the flow inside and at the exit of a supersonic nozzle (convergent followed by a long straight tube - no divergent) using rhoCentralFoam.

My inlet BC is totalPressure = 5.00bar and the oulet BC is waveTransmissive = 1.01bar.
To take advantage of the geometry, I am solving the flow with axisymmetric BC. My domain is then a simple wedge - of less than 5°.
rhoCentralFoam is used as a density-based solver and the turbulent model used here is k-omega SST.

The issue I have is that the simulation evolves a lot between the first iterations, the iterations in the middle and the last iteration. Of course, I know rhoCentralFoam is a transient solver, but my system should be steady, when the right regime is reached. The solution I get seems to give me pretty good results according to experimental results at step 2000 but the solution at final step 38600 is far from the experimental results.

I attached figures showing the pressure gradient at step 2000 and step 38600 as well as the "numerical schlieren" and my experimental schlieren results.

Does anyone have a clue ?

Best
Attached Images
File Type: png 2000.png (26.3 KB, 57 views)
File Type: png 38000.png (18.0 KB, 50 views)
File Type: jpg Capture.jpg (34.0 KB, 49 views)
K-al-Eps-o is offline   Reply With Quote

Old   June 13, 2023, 06:40
Default
  #2
New Member
 
Baptiste
Join Date: Mar 2023
Posts: 25
Rep Power: 3
K-al-Eps-o is on a distinguished road
Quote:
I attached figures showing the pressure gradient at step 2000 and step 38600 as well as the "numerical schlieren" and my experimental schlieren results.
*pressure field
not "pressure gradient"
K-al-Eps-o is offline   Reply With Quote

Old   June 21, 2023, 05:10
Default Any thoughts ?
  #3
New Member
 
Baptiste
Join Date: Mar 2023
Posts: 25
Rep Power: 3
K-al-Eps-o is on a distinguished road
Has anyone any idea by any chance ?
K-al-Eps-o is offline   Reply With Quote

Old   June 22, 2023, 13:09
Default
  #4
Senior Member
 
Join Date: Oct 2011
Posts: 242
Rep Power: 17
naffrancois is on a distinguished road
I would try first without any viscosity, Euler equations.

Also I don't know precisely what is behind the outlet bc you use, but if flow is mostly supersonic there I would give a shot to a supersonic bc, which is basically an extrapolation of all interior variables, e.g. all waves going out of the domain. If your outlet is not truly transmissive it will reflect waves and spoil the shape of the jet. Alternatively you can try to increase the size of the "receiver", e.g. place the outlet further away from the jet and use bigger cells close to the outlet to damp spurious waves.
naffrancois is offline   Reply With Quote

Old   August 8, 2023, 06:10
Default
  #5
New Member
 
Baptiste
Join Date: Mar 2023
Posts: 25
Rep Power: 3
K-al-Eps-o is on a distinguished road
Thanks for your insight, I will follow your recommandations.

What do you mean by supersonic bc? Do you mean 'supersonicFreestream'?
K-al-Eps-o is offline   Reply With Quote

Old   August 8, 2023, 08:38
Default
  #6
Senior Member
 
Joern Beilke
Join Date: Mar 2009
Location: Dresden
Posts: 516
Rep Power: 20
JBeilke is on a distinguished road
What do you mean with "iteration"? It is a transient solver. So we are talking about timesteps, time-step-size, countant-numbers ....
JBeilke is offline   Reply With Quote

Old   August 8, 2023, 08:52
Default
  #7
New Member
 
Baptiste
Join Date: Mar 2023
Posts: 25
Rep Power: 3
K-al-Eps-o is on a distinguished road
Yeah yeah, sorry. I meant time-step.
I was mistaken because my solution should be steady (solver should converge towards a steady solution) but the solver is transient.
Thanks
K-al-Eps-o is offline   Reply With Quote

Old   August 8, 2023, 09:13
Default
  #8
Senior Member
 
Joern Beilke
Join Date: Mar 2009
Location: Dresden
Posts: 516
Rep Power: 20
JBeilke is on a distinguished road
Here you have an example with a "subsonicSupersonicPressureOutlet". It was used for the following paper:


https://ieeexplore.ieee.org/document/8636345

We used rhoPimpleCentralFOAM from:

https://github.com/unicfdlab/hybridCentralSolvers

and the libcompressibleTools (contains this BC) from:

https://github.com/unicfdlab/libcompressibleTools


It was done with OF4.0. I did not check with a new version.
Attached Files
File Type: gz Ladenburg_rhoPimpleCentralFOAM.tar.gz (3.7 KB, 11 views)
JBeilke is offline   Reply With Quote

Old   August 8, 2023, 10:27
Default
  #9
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,748
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
Would you verify that the solution at 38000 is still supersonic near the dump? Do you have a train of curved bow shocks or did the solution revert back to a subsonic solution? If subsonic then we need to steer the solution towards the supersonic case. If you have bow shocks then we need to think whether it is just an issue of discretization or the BCs needing to be adjusted to get the overexpanded jet which will give you the shock diamonds you are looking for (i.e. increase or lower the back pressure).

It might be easier to visualize if you plot density or Mach number instead of the pressure. Obviously plotting the Mach number would directly answer my question. Pressure is fine too but you have to really zoom in to see the gradient and then do the logic in your head.
LuckyTran is offline   Reply With Quote

Old   August 8, 2023, 10:56
Default
  #10
Senior Member
 
Joern Beilke
Join Date: Mar 2009
Location: Dresden
Posts: 516
Rep Power: 20
JBeilke is on a distinguished road
The problem might be the waveTransmissive BC. There is a parameter "lInf" which specifies a virtual length, which is added to the outlet, where the pressure has its specified value. If you don't explicitly set the value of "lInf" to something meaningful, it is assumed to be very large (infinity ??) and your solution is very stable but wrong.
JBeilke is offline   Reply With Quote

Reply

Tags
openfoam 2212, rhocentrafoam, supersonic nozzle


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Frequently Asked Questions about Installing OpenFOAM wyldckat OpenFOAM Installation 3 November 14, 2023 12:58
LES of a nozzle with rhoCentralFoam vcattoir OpenFOAM Running, Solving & CFD 7 February 22, 2017 21:37
De laval Supersonic Nozzle Exit Mach number in CFX bharath CFX 7 December 19, 2014 03:44
Suggestion for a new sub-forum at OpenFOAM's Forum wyldckat Site Help, Feedback & Discussions 20 October 28, 2014 10:04
reversed flow in over-expanded supersonic nozzle imnull FLUENT 0 March 28, 2013 11:48


All times are GMT -4. The time now is 12:46.