Problem C1.2. Ringleb Problem
Overview
This problem tests the spatial accuracy of high-order methods. The flow is transonic and smooth. The geometry is also smooth, and high-order curved boundary representation appears to be critical. The exact solution is known via ¡®hodograph transformation¡¯ [1].
Governing Equations
The governing equations are the 2D Euler
equations with 
.
Geometry
Let k be a streamline parameter, i.e., k = constant
on each streamline. The two stream lines for the two wall boundaries are 
for the inner wall, and 
for the outer wall. Let 
be the velocity magnitude. For each fixed 
, 
, the variable varies between 
and . For each , define the
speed of sound , density , pressure , and a quantity
denoted by 
by
|
(1) |
For each pair 
, set
|
(2) |
|
|
|
|
Fig. 1.2 Ringleb geometry; thick curves: walls; thin curves: inflow and outflow boundaries.
With again , 
, and 
, the four boundaries are:
(a) inflow, 
, , and 
; (b) outflow, 
, , and ; (c) inner wall, 
and 
; and outer wall, 
and . See Fig. 1.
Exact Solution
The exact solution is given
by (1) and (2). The flow is irrotational and
isentropic. It reaches a supersonic speed of Mach number 1.5 at location of the inner wall. The sonic line for the
problem is shown in Fig. 2. Entropy should be a constant in the flow field.
Fig. 2. Plot of velocity magnitude and streamlines for the region .
Requirements
1. Start the simulation from the exact solution, and monitor convergence using the L2 norm of the density residual. Compute the work units required to achieve steady state.
2.
Similar to the smooth bump problem, use the non-dimensional 
entropy errors
as accuracy indicator. Perform the calculations for at least three different
meshes (quadruple the number of cells each time); use methods with different p to assess the performance of schemes
of various orders.
3. To assess the efficiency of the spatial discretization, calculate the residual (right hand side) on the densest mesh 1000 times (equivalent to 1000 Euler forward time steps). Provides the results in work units.
4. Submit three sets of data to the workshop contact (above)
o Errors vs mesh sizes for different p
o Errors and work units for different h (
) and p
o Work units to complete the residual evaluation 1000 times for various p.
References
[1] G. Chiocchia, Exact solutions to transonic and supersonic flows. Technical Report AR-211, AGARD, 1985.