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.

1<PARAMETERS> 
2  <P> TimeStep            = 0.0000001                  </P> 
3  <P> FinTime             = 1.0                        </P> 
4  <P> NumSteps            = FinTime/TimeStep           </P> 
5  <P> IO_CheckSteps       = 5000                       </P> 
6  <P> IO_InfoSteps        = 1                          </P> 
7  <P> Gamma               = 1.4                        </P> 
8  <P> pInf                = 101325                     </P> 
9  <P> rhoInf              = 1.225                      </P> 
10  <P> GasConstant         = 287.058                    </P> 
11  <P> TInf                = pInf/(287.058*rhoInf)      </P> 
12  <P> Twall               = pInf/(287.058*rhoInf)+15.0 </P> 
13  <P> uInf                = 147.4                      </P> 
14  <P> vInf                = 0.0                        </P> 
15  <P> wInf                = 0.0                        </P> 
16  <P> mu                  = 1e-5                       </P> 
17  <P> Pr                  = 0.72                       </P> 
18  <P> thermalConductivity = 0.02                       </P> 
19  <P> IO_Timer_Level      = 3                          </P> 
20</PARAMETERS>

Time Integration Scheme

Under this section it is possible to set the time integration scheme information.

1<TIMEINTEGRATIONSCHEME> 
2  <METHOD> RungeKutta </METHOD> 
3  <VARIANT> SSP </VARIANT> 
4  <ORDER> 3 </ORDER> 
5</TIMEINTEGRATIONSCHEME>

Solver info

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

1<SOLVERINFO> 
2  <I PROPERTY="EQTYPE"                        VALUE="NavierStokesImplicitCFE" /> 
3  <I PROPERTY="Projection"                    VALUE="DisContinuous"           /> 
4  <I PROPERTY="TimeIntegrationMethod"         VALUE="DIRKOrder2"              /> 
5  <I PROPERTY="AdvectioType"                  VALUE="WeakDG"                  /> 
6  <I PROPERTY="DiffusionType"                 VALUE="InteriorPenalty"         /> 
7  <I PROPERTY="UpwindType"                    VALUE="Roe"                     /> 
8  <I PROPERTY="ProblemType"                   VALUE="General"                 /> 
9  <I PROPERTY="ViscosityType"                 VALUE="Constant"                /> 
10  <I PROPERTY="EquationOfState"               VALUE="IdealGas"                /> 
11  <I PROPERTY="Driver"                        VALUE="Standard"                /> 
12</SOLVERINFO>

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:

1<VARIABLES> 
2  <V ID="0"> rho  </V> 
3  <V ID="1"> rhou </V> 
4  <V ID="2"> E    </V> 
5</VARIABLES>

For a 2D problem we have

1<VARIABLES> 
2  <V ID="0"> rho  </V> 
3  <V ID="1"> rhou </V> 
4  <V ID="2"> rhov </V> 
5  <V ID="3"> E    </V> 
6</VARIABLES>

For a 3D problem we have:

1<VARIABLES> 
2  <V ID="0"> rho  </V> 
3  <V ID="1"> rhou </V> 
4  <V ID="2"> rhov </V> 
5  <V ID="3"> rhow </V> 
6  <V ID="4"> E    </V> 
7</VARIABLES>

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:

1<BOUNDARYREGIONS> 
2  <B ID="0"> C[100] </B> 
3</BOUNDARYREGIONS>

Finally we can specify the boundary conditions on the regions specified under BOUNDARYREGIONS. Note that the no-slip, isothermal, wall boundary condition requires Twall to be specified. In the following are some examples for a 2D problem:

Initial conditions and exact solution

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

1<FUNCTION NAME="InitialConditions"> 
2  <E VAR="rho"    VALUE="rhoInf"/> 
3  <E VAR="rhou"   VALUE="rhoInf*uInf"   /> 
4  <E VAR="rhov"   VALUE="rhoInf*vInf"   /> 
5  <E VAR="E" 
6  VALUE="pInf/(Gamma-1)+0.5*rhoInf*(uInf*uInf+vInf*vInf)"/> 
7</FUNCTION> 
8 
9<FUNCTION NAME="ExactSolution"> 
10  <E VAR="rho"    VALUE="rhoInf"          /> 
11  <E VAR="rhou"   VALUE="rhoInf*uInf"   /> 
12  <E VAR="rhov"   VALUE="rhoInf*vInf"   /> 
13  <E VAR="E" 
14  VALUE="pInf/(Gamma-1)+0.5*rhoInf*(uInf*uInf+vInf*vInf)"/> 
15</FUNCTION>