Inviscid compressible flow in two space dimensions. A geometry according to the picture below is defined and the problem consists of computing numerically the 2D steady-state flow including compression regions, shocks and expansion regions. The flow domain consists of a lower ramp and an upper straight wall, giving a converging 2D channel. A supersonic flow is applied at the inlet.

Overall view of computational domain
Figure 1: Overall view of computational domain, showing shape and extent of outer boundary


Literature review


Suggested search topics

  • Oblique shock
  • Prandtl-Meyer expansion
  • Shock-expansion theory
  • Slip line
  • Compression ramp
  • Shock reflection
  • Mach reflection


Specifications


The smooth ramp section starts at \((x,y)=(x1,0.0)\) and reaches \(y=h\) at the axial coordinate \(x=x_2\). For the lower wall and axial coordinates greater than \(x_2\), \(y\) is constant (\(y=h\)).


	x1 = 1.0 m
	x2 = 4.0 m
	x3 = 7.0 m
	y1 = 2.0 m
	h  = 0.6 m

For an inlet Mach number of 2.0, the inlet and outlet flow conditions are as follows:

inlet
\(\rho\) 0.267124 \(kg/m^3\)
\(p\) 12780 \(Pa\)
\(u\) 517.61 \(m/s\)
outlet
\(\rho\) 0.4045 \(kg/m^3\)
\(p\) 22840 \(Pa\)
\(u\) 455.481 \(m/s\)


Task 1


Do an inviscid simulation for the smooth compression ramp as specified above



Task 2


Investigate the effect of viscosity for the compression ramp flow by running a viscous simulation using the Spallart-Almaras turbulence model.



Expected results and presentation


The flow field will contain compression regions, oblique shocks, and expansion fans. Although contour plots should be used for the visualization of these flow features it might be difficult to do a qualitative comparison of different simulations (simulations made using different meshes and numerical settings) and analytical results. Therefore, when comparing shock strength and location of flow features, it is recommended to extract data along axial lines at different \(y\) coordinates and compare the data in terms of \(xy\)-plots.



Grid generation guidelines


There should not be any large jump in cell sizes anywhere. The changes in cell size must be smooth otherwise you might run into problems with convergence.

The smooth ramp contour can be generated using the spline functionality in the STAR-CCM+ sketch tool. Start by making two points using the Create point tool.

point tool

Adjust the coordinates according to the specifications for \(x_1\) and \(x_2\) above. Then right click on each of the two points and click Apply Fixation Constraint. This will make sure that the spline generation tool does not modify the location of these coordinates later.

ramp end points

Now, use the Create Spline tool to generate the gradual ramp.

spline tool

Select the \(x_1\) coordinate as the first point on the spline curve, then select two control points between the \(x_1\) and \(x_2\) (these points can be moved later). Select the \(x_2\) coordinate as the end point pf the spline curve. Press esc on your keyboard to exit the spline tool. Move the two control points such that the ramp looks smooth.

ramp profile

Now, connect all other corner points with lines and the sketch is done.

ramp domain


CFD guidelines


Some general guidelines for the simulation:

  • Don't forget to set the reference pressure to zero:
Continua -> Fluid -> Reference Values -> Reference Pressure
  • In your project, you should investigate the performance of different numerical schemes. The numerical scheme is modified under Continua -> Fluid -> Models -> Coupled Flow. Note, in case you have not renamed the physics entry under Continua, the default name is Physics 1.
numerical scheme
  • Some of the cases, or at least some of the operation conditions for some cases, will have convergence problems due to flow instabilities. In case you get that type of solver behavior, try to find another convergence criterion by for example measuring forces or massflow. You may also consider updating the mesh as it is crucial to resolve important flow features to reach convergence.
  • Some of the cases should be simulated using a viscous solver. Choose Turbulent instead of Inviscid when you select models and then select an appropriate turbulence model.
turbulence models
  • In case you would like to start your simulation from scratch, the flow field can be reset as indicated in the picture below (don't forget to initialize the flow field after the reset).
reset flow
  • Time stepping approach (implicit/explicit) is selected under Continua -> Fluid -> Models -> Coupled Flow
reset flow
  • CFL number and other solver settings are found under Solver -> Coupled Implicit (or Solver -> Coupled Explicit)
reset flow
reset flow