Supersonic Expansion (Prandtl–Meyer Fan)

verification

Velocity-magnitude field through the expansion fan
Velocity-magnitude field on the quasi-2D plane; the flow accelerates through the centered fan that springs from the convex corner and follows the deflected wall.

Overview

When a supersonic stream flows around a convex corner it expands: the flow turns away from itself through a continuous fan of Mach waves centered on the corner, accelerating to a higher Mach number while its pressure, density and temperature fall. Unlike a shock, the expansion is isentropic and spread over a finite angular fan, so the process is reversible and the post-expansion state follows from the Prandtl–Meyer function together with the isentropic flow relations. For a given upstream Mach number M1 and turning angle, the downstream Mach number M2 is fixed by ν(M2)=ν(M1)+Δθ, and the pressure and temperature ratios then follow from the isentropic relations evaluated at M1 and M2.

The solution is closed-form and free of any modeling assumption, so the case verifies that the solver captures a smooth supersonic expansion, the compression complement of the oblique-shock wedge, at the correct strength.

Problem Setup

A Mach-2 uniform stream flows along a flat wall that turns away by 15 at a convex corner; the opposite wall is parallel and the domain is sized so the expansion fan leaves through the downstream boundary without reflecting. The flow is modeled as inviscid (Euler), compressible and steady, with an ideal gas (γ=1.4). The walls are slip boundaries, the inlet imposes the supersonic freestream, and the outlet is a supersonic pressure outlet that extrapolates the interior state. A second-order FDS scheme with no limiter is used.

The grid is a structured convex-corner mesh built so the turning angle is exactly 15. The flow downstream of the fan is uniform and aligned with the deflected wall, where the post-expansion state is sampled.

Quantities of Interest

The verification quantities are the post-expansion flow state:

Each is compared to the exact Prandtl–Meyer and isentropic result for M1=2, a 15 turn and γ=1.4.

Sources

Ames Research Staff, "Equations, Tables, and Charts for Compressible Flow," NACA Report 1135, 1953.: https://ntrs.nasa.gov/citations/19930091059
Authoritative tabulation of the Prandtl–Meyer function and the isentropic flow relations used as the analytic reference.

Results

Post-expansion state vs the Prandtl–Meyer solution

SourceM2p2/p1T2/T1
Luminary 2.60.393050.76688
Analytic 2.59840.393070.76583

Post-expansion Mach number M2 and the isentropic p2/p1, T2/T1 sampled in the uniform region above the deflected wall.

Solver configuration

features exercised
Verification2DInviscidIdeal gasSteadyEnergy equation