PHengLEI-TestCases/TwoD_Rae2822_SA_Unstruct/bin/cfd_para_transonic.hypara

#########################################################################
#                     General Control Parameter                         #
#########################################################################
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual 'res.dat' saved.

int maxSimuStep       = 20000;

int intervalStepFlow  = 1000;
int intervalStepPlot  = 1000;
int intervalStepForce = 100;
int intervalStepRes   = 10;

#########################################################################
#                       Inflow Parameter                                #
#########################################################################
// refMachNumber: Mach number.
// attackd: Angle of attack.
// angleSlide: Angle of sideslip.
// inflowParaType: The type of inflow parameters.
//                 0 -- the nondimensional conditions.
//                 1 -- the flight conditions.
//                 2 -- the experiment conditions.
//                 3 -- the subsonic boundary conditions.
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
// height: Fly height, unit of km.
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
//                    1 dm = 0.1  m.
//                    1 cm = 0.01 m.
//                    1 mm = 0.001m.
//                    1 inch             = 0.0254m.
//                    1 foot = 12 inches = 0.3048m.
//                    1 yard = 3  feet   = 0.9144m.
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.

double refMachNumber = 0.73;
double attackd       = 2.79;
double angleSlide    = 0.00;

int inflowParaType = 0;
double refReNumber = 6.5e6;
double refDimensionalTemperature = 288.15;

//int inflowParaType = 1;
//double height = 0.001;

//int inflowParaType = 2;
//double refDimensionalTemperature = 6051.024;    // The total temperature, T*(1+(refGama-1)*M*M/2).
//double refDimensionalPressure = 4.299696E09;    // The total pressure, p*(T0/T)^(refGama/(refGama-1)).

double gridScaleFactor = 1.0;

double forceRefenenceLengthSpanWise = 1.0;    // unit of meter.
double forceRefenenceLength = 1.0;            // unit of meter.
double forceRefenenceArea = 1.0;              // unit of meter^2.
double TorqueRefX = 0.0;                      // unit of meter. 
double TorqueRefY = 0.0;                      // unit of meter.
double TorqueRefZ = 0.0;                      // unit of meter.

#########################################################################
#                       Physical models                                 #
#########################################################################
// viscousType : Viscous model.
//           0 -- Euler.
//           1 -- Lamilar.
//           3 -- 1eq turbulent.
//           4 -- 2eq turbulent.
// viscousName: Laminar or tubulent model.
//              -- "1eq-sa", when viscousType = 3.
//              -- "2eq-kw-menter-sst", when viscousType = 4.
// DESType: Type of DES.
//          0 -- RANS.(default)
//          1 -- DES.
//          2 -- DDES.
//          3 -- IDDES.

//int    viscousType    = 0;
//string viscousName = "Euler";

//int    viscousType    = 1;
//string viscousName = "laminar";

int    viscousType    = 3;
string viscousName = "1eq-sa";

//int    viscousType    = 4;
//string viscousName = "2eq-kw-menter-sst";

int DESType = 0;

int roeEntropyFixMethod = 3;
double roeEntropyScale  = 1.0;
#########################################################################
#                       Spatial Discretisation                          #
#########################################################################
#*******************************************************************
#                       Struct Solver                              *
#*******************************************************************
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
//                     Using this when solve structered grid or hybrid.
//                     -- "roe", "vanleer", "ausm+up", "ausmpw".
// str_limiter_name: Limiter of struct grid.
//                   -- "3rdsmooth", "smooth".
//                   -- "nolim", no limiter.

string inviscidSchemeName = "roe";
string str_limiter_name   = "smooth";

#*******************************************************************
#                       UnStruct Solver                            *
#*******************************************************************
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
//                  Using this when solve Unstructered grid or hybrid.
//                  -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
//                  -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
// uns_limiter_name: Limiter of Unstruct grid.
//                   -- "vencat", "barth".
//                   -- "1st", meaning accuracy of first-order.
//                   --  "nolim", no limiter.
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
//              The smaller the value, the more robust it is.

string uns_scheme_name  = "roe";
string uns_limiter_name = "nolim";
double venkatCoeff      = 5.0;

#########################################################################
#                       Temporal Discretisation                         #
#########################################################################
// iunsteady: Steady or unsteady.
//            0 -- steady.
//            1 -- unsteay.
// CFLEnd: The CFL number, [0.1, 100].
//         The bigger the value, the convergence faster but lower robustness.
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
//               1 -- is recommended for structured solver.
//             1-3 -- is recommended for unstructured solver.

int iunsteady = 0;

double CFLEnd = 20.0;

int nLUSGSSweeps = 5;

#########################################################################
#                       Multi-Grid parameters                           #
#########################################################################
// nMGLevel: The number of Multi-Grid level.
//           = 1 -- single-level.
//           > 1 -- multi-level.
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
//               Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
//               Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.

int nMGLevel = 1;
int flowInitStep = 100;

#########################################################################
#                       File In or Out                                  #
#########################################################################
// gridfile: The partitioned Grid file path, using relative path,
//           which is relative to the working directory.
// IMPORTANT WARNING: The file index should be ignored,
//                     e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
//                     Please use 'rae2822_hybrid2d__4.fts' here!
// isPlotVolumeField: If dump out the whole field results to tecplot or not, 0 / 1.

string gridfile = "./grid/rae2822_hybrid2d_fine.fts";
int isPlotVolumeField = 0;

// ----------------- Advanced Parameters, DO NOT care it ----------------
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
// visualVariables:  Variable types dumped, listed as following:
//                   -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
//                   -- viscosityLaminar(7), viscosityTurbulent(8),
//                   -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
//                   -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
//                   -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21).
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
// Variables order must from small to big.

int nVisualVariables = 8;
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];

// limitVariables: Limit model (It is useful only if limitVector is 0).
//                 0 -- limit only for pressure and denstiny, then get the min value.
//                 1 -- limit for every variables, then get the min value.
// limitVector:
//              0 -- Each variable use the same limiter coefficient.
//              1 -- Each variable use the respective limiter coefficients.
// reconmeth:
//            0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
//            1 -- Q+, Q- use the min limiter coefficients of left and right cell.

int reconmeth = 1;
int limitVariables = 0;
int limitVector = 0;