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’ .

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)

 (a) Inflow

 (b) Outflow

 (d) Outer wall

 (c) Inner wall 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

 G. Chiocchia, Exact solutions to transonic and supersonic flows. Technical Report AR-211, AGARD, 1985.