9.3 Session file configuration

In the following we describe the session file configuration. Specifically we consider the sections under the tag <CONDITIONS> in the session (.xml) file.

Parameters

Under this section it is possible to set the parameters of the simulation.

• TimeStep is the time-step we want to use;
• FinTime is the final physical time at which we want our simulation to stop;
• NumSteps is the equivalent of FinTime but instead of specifying the physical final time we specify the number of time-steps;
• IO_CheckSteps sets the number of steps between successive checkpoint files;
• IO_InfoSteps sets the number of steps between successive info stats are printed to screen;
• Gamma ratio of the specific heats. Default value = 1.4;
• pInf farfield pressure (i.e. p_∞). Default value = 101325 Pa;
• rhoInf farfield density (i.e. ρ_∞). Default value = 1.225 Kg∕m³;
• TInf farfield temperature (i.e. T_∞). Default value = 288.15 K;
• Twall temperature at the wall when isothermal boundary conditions are employed (i.e. T_w). Default value = 300.15K;
• uint farfield X-component of the velocity (i.e. u_∞). Default value = 0.1 m∕s;
• vInf farfield Y -component of the velocity (i.e. v_∞). Default value = 0.0 m∕s;
• wInf farfield Z-component of the velocity (i.e. w_∞). Default value = 0.0 m∕s;
• mu dynamic viscosity (i.e. μ_∞). Default value = 1.78e-05 Pas;
• Pr Prandtl number. Default value = 0.72;
• thermalConductivity thermal conductivity (i.e. κ_∞). Default value = 0.0257 W∕(Km);

Solver info

Under this section it is possible to set the solver information.

• EQType is the tag which specify the equations we want solve:
  • NavierStokesCFE (Compressible Navier-Stokes equations);
  • EulerCFE (Compressible Euler equations).
• Projection is the type of projection we want to use:
  • DisContinuous.
    Note that the Continuous projection is not supported in the Compressible Flow Solver.
• AdvectionType is the advection operator we want to use:
  • WeakDG (classical DG in weak form);
  • FRDG (Flux-Reconstruction recovering nodal DG scheme);
  • FRSD (Flux-Reconstruction recovering a spectral difference (SD) scheme);
  • FRHU (Flux-Reconstruction recovering Huynh (G2) scheme);
  • FRcmin (Flux-Reconstruction with c = c_min);
  • FRcinf (Flux-Reconstruction with c = ∞).
• DiffusionType is the diffusion operator we want to use:
  • WeakDG (classical DG in weak form);
  • FRDG (Flux-Reconstruction recovering nodal DG scheme);
  • FRSD (Flux-Reconstruction recovering a spectral difference (SD) scheme);
  • FRHU (Flux-Reconstruction recovering Huynh (G2) scheme);
  • FRcmin (Flux-Reconstruction with c = c_min);
  • FRcinf (Flux-Reconstruction with c = ∞).
• TimeIntegrationMethod is the time-integration scheme we want to use. Note that only an explicit discretisation is supported:
  • ForwardEuler;
  • RungeKutta2_ImprovedEuler;
  • ClassicalRungeKutta4.
• UpwindType is the numerical interface flux (i.e. Riemann solver) we want to use for the advection operator:
  • AUSM0;
  • AUSM1;
  • AUSM2;
  • AUSM3;
  • Average;
  • ExactToro;
  • HLL;
  • HLLC;
  • LaxFriedrichs;
  • Roe.
• ProblemType is the problem type we want to solve. This tag is supported for solving ad hoc problems such as the isentropic vortex or the Ringleb flow.
  • General;
  • IsentropicVortex;
  • RinglebFlow;
    
• ViscosityType is the viscosity type we want to use:
  • Constant (Constant viscosity);
  • Variable (Variable viscosity through the Sutherland’s law.);

Boundary conditions

In this section we can specify the boundary conditions for our problem. First we need to define the variables under the section VARIABLES. For a 1D problem we have:

For a 2D problem we have

For a 3D problem we have:

After having defined the variables depending on the dimensions of the problem we want to solve it is necessary to specify the boundary regions on which we want to define the boundary conditions:

Finally we can specify the boundary conditions on the regions specified under BOUNDARYREGIONS. In the following some examples for a 2D problem:

• Slip wall boundary conditions:
  1<BOUNDARYCONDITIONS> 
  2  <REGION REF="0"> 
  3    <D VAR="rho"  USERDEFINEDTYPE="Wall" VALUE="0" /> 
  4    <D VAR="rhou" USERDEFINEDTYPE="Wall" VALUE="0" /> 
  5    <D VAR="rhov" USERDEFINEDTYPE="Wall" VALUE="0" /> 
  6    <D VAR="E"    USERDEFINEDTYPE="Wall" VALUE="0" /> 
  7  </REGION> 
  8</BOUNDARYCONDITIONS>
• No-slip wall boundary conditions:
  1<BOUNDARYCONDITIONS> 
  2  <REGION REF="0"> 
  3    <D VAR="rho"  USERDEFINEDTYPE="WallViscous" VALUE="0" /> 
  4    <D VAR="rhou" USERDEFINEDTYPE="WallViscous" VALUE="0" /> 
  5    <D VAR="rhov" USERDEFINEDTYPE="WallViscous" VALUE="0" /> 
  6    <D VAR="E"    USERDEFINEDTYPE="WallViscous" VALUE="0" /> 
  7  </REGION> 
  8</BOUNDARYCONDITIONS>
• Adiabatic wall boundary conditions:
  1<BOUNDARYCONDITIONS> 
  2  <REGION REF="0"> 
  3    <D VAR="rho"  USERDEFINEDTYPE="WallAdiabatic" VALUE="0" /> 
  4    <D VAR="rhou" USERDEFINEDTYPE="WallAdiabatic" VALUE="0" /> 
  5    <D VAR="rhov" USERDEFINEDTYPE="WallAdiabatic" VALUE="0" /> 
  6    <D VAR="E"    USERDEFINEDTYPE="WallAdiabatic" VALUE="0" /> 
  7  </REGION> 
  8</BOUNDARYCONDITIONS>
• Farfield boundary conditions (including inviscid characteristic boundary conditions):
  1<BOUNDARYCONDITIONS> 
  2  <REGION REF="0"> 
  3    <D VAR="rho"  VALUE="rhoInf" /> 
  4    <D VAR="rhou" VALUE="rhoInf*uInf" /> 
  5    <D VAR="rhov" VALUE="rhoInf*vInf" /> 
  6    <D VAR="E" 
  7    VALUE="pInf/(Gamma-1)+0.5*rhoInf*(uInf*uInf+vInf*vInf+wInf*wInf)"/> 
  8  </REGION> 
  9</BOUNDARYCONDITIONS>
• Pressure outflow boundary conditions:
  1<BOUNDARYCONDITIONS> 
  2  <REGION REF="0"> 
  3    <D VAR="rho"  USERDEFINEDTYPE="PressureOutflow" VALUE="0" /> 
  4    <D VAR="rhou" USERDEFINEDTYPE="PressureOutflow" VALUE="0" /> 
  5    <D VAR="rhov" USERDEFINEDTYPE="PressureOutflow" VALUE="0" /> 
  6    <D VAR="E"    USERDEFINEDTYPE="PressureOutflow" VALUE="0" /> 
  7  </REGION> 
  8</BOUNDARYCONDITIONS>

Initial conditions and exact solution

Under the two following sections it is possible to define the initial conditions and the exact solution (if existent).