SupersonicWavesReflectionsDuct SimulationExample
SupersonicWavesReflectionsDuct SimulationExample
Reflections of 2D Supersonic
Waves
Results
• The static pressure contours show clearly how the shock reflects off the duct walls.
• We can also identify where the expansion fan forms. Looking at the areas of smooth variations of
pressure contours color, we can get an idea on how it reflects and widens along the duct.
• This alternating pattern of shock and expansion reflections continue throughout the duct.
Reflected Expansion
Expansion Fan
2
Results
3
Results
• We can also observe that the shock interacts with the boundary layer.
• At the shock nears the wall, it impinges on the boundary layer and forced it to separate.
• This fascinating phenomenon of the physical shock-induced boundary layer separation would be
overlooked if we used the inviscid assumption.
4
Results
• The pathlines show how the flow suddenly changes direction past the shock and its reflections, while it
smoothly turns as it crosses the expansion fan and its reflections.
• Interestingly, the fluid experiences a smoother change of direction where the expansion fan and the
shock waves interact.
• The reflections in the duct cause the fluid to have an undulatory motion.
5
Appendix
Reflections of 2D Supersonic Waves
Boundary Conditions
• The velocity inlet and the pressure outlet
boundary conditions require the temperature
input to provide well-posed set of the boundary
conditions.
• The velocity inlet require a static value of the
temperature, while the pressure outlet requires
backflow total temperature.
• No reversed flow develops in this problem, thus
the temperature setting at the outlet is not
applied.
• Still, it is a good practice to apply realistic values
of the backflow total temperature at the outlet.
• In order to estimate the total temperature, we
use the isentropic relations as follows:
𝑇 = 150 𝐾 𝑀=3 𝛾 = 1.4
𝛾−1 2
𝑇0 = 𝑇 1 + 𝑀 = 420 𝐾
2
7
Mesh Adaption
• The gradients across the shock are large and require higher mesh resolution
for resolving them. Mesh Adaption can be used to resolve the resolution
across the shocks to obtain thinner, sharper shocks.
• Mesh adaption can be applied using these two steps:
‐ Cell Register: create a cell register using the gradient derivative option for the static Initial Mesh
pressure in order to identify the cells in the shock region.
‐ Refinement: Fluent refines the region identified by the cell register created in the
previous step.
Refinement