|
[Sponsors] |
supersonic (Laval) nozzle, unphysical Mach number |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 13, 2014, 05:27 |
supersonic (Laval) nozzle, unphysical Mach number
|
#1 |
New Member
Judith Richter
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
Hello,
I am fairly new to ANSYS CFX and am struggling with the simulation of a Ma = 1,7 covergent-divergent nozzle. My area ratio is 1.337 which leads to an isentropic exit Mach number of 1,7. From my understanding the exit Mach number is indepent of the pressure ratio (as long as sonic condition is reached at the nozzle throat). However, my simulation obviously depends on the exit pressure: If I decrease the exit pressure, the Mach number increses (over 1.7). It looks like the nozzle flow is even over-expanded (see the shock waves in the picture). I tried many boundary conditions (subsonic, supersonic, opening) but only the subsonic BC converges: IN subsonic total pressure = 2 bar turbulence: medium static temperature = 380 K OUT subsonic average static pressure: 0.3 bar WALL free slip Does anybody have an idea, what my mistake could be? |
|
August 13, 2014, 13:48 |
|
#2 |
New Member
Judith Richter
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
Here is the Mach number plot of the nozzle: The red color indicates the areas where the Mach number exceeds 1.7. Hope that helps understanding my problem.
Thank you! |
|
August 13, 2014, 19:53 |
|
#3 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144 |
I do not understand your problem. Your results look believeable to me - you are close to your expected Mach 1.7 and superimposed on that you have some Mach wave reflections. This is all as expected.
|
|
August 14, 2014, 04:31 |
|
#4 |
New Member
Judith Richter
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
That's true, but I don't understand where the Mach waves come from. The nozzle contour comes from the method of characteristics and the simulation is inviscid...
|
|
August 14, 2014, 07:43 |
|
#5 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144 |
You have a sharp transition from the expansion section to the straight duct. This will create a Mach wave - and that is what you see. This is not a viscous effect, it is a compressibility effect. So no surprises there.
I assume you are comparing the simulation to a 1D model or an analytical solution of supersonic flow. These approaches ignore flow details and assume there are no sharp transition. |
|
August 14, 2014, 09:54 |
|
#6 | |
New Member
Judith Richter
Join Date: Jul 2014
Posts: 8
Rep Power: 12 |
Quote:
I discussed this topic with my colleague and he suggested that the disturbance comes from the grid. However, my grid is very fine (>500.000 nodes) and I cannot imagine that this would have just a big influence... |
||
August 14, 2014, 18:15 |
|
#7 |
Senior Member
Edmund Singer P.E.
Join Date: Aug 2010
Location: Minneapolis, MN
Posts: 511
Rep Power: 21 |
You are below nozzle design pressure on your outlet. You should expect to see what that series of mach waves.
If you dont want to see those, increase your back pressure. |
|
August 14, 2014, 19:28 |
|
#8 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144 |
If you think the Mach wave artefact is from your grid then simply changing the grid (significantly finer preferably, but coarser will probably work too) will show whether you are right. If the pattern changes then yes, it is a grid artefact.
|
|
Tags |
cfx, laval, mach number, nozzle, supersonic |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] snappyHexMesh sticking point | natty_king | OpenFOAM Meshing & Mesh Conversion | 11 | February 20, 2024 10:12 |
AMI speed performance | danny123 | OpenFOAM | 21 | October 24, 2020 05:13 |
[blockMesh] --> foam fatal error: | lillo763 | OpenFOAM Meshing & Mesh Conversion | 0 | March 5, 2014 11:27 |
How obtain a Mach number value from a case run with Xifoam solver?? | lfgmarc | OpenFOAM Programming & Development | 0 | June 15, 2011 04:00 |
Laval Nozzle Mach Contours | Sohail Ahmed | Main CFD Forum | 0 | May 19, 2004 06:42 |