PHengLEI-TestCases/TwoD_30p30n_SST_Struct/bin/cfd_para_subsonic.hypara

190 lines
8.8 KiB
Plaintext
Raw Permalink Normal View History

2021-03-22 16:44:54 +08:00
#########################################################################
# 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 = 50000;
int intervalStepFlow = 1000;
int intervalStepPlot = 1000;
int intervalStepForce = 100;
int intervalStepRes = 10;
// ifLowSpeedPrecon: precondition to accelerate convergence for low speed flow.
// ONLY available for unstructured solver!
// 0 -- no precondition process (default, mach > 0.2).
// 1 -- carry out precondition process ( mach number < 0.2 ).
int ifLowSpeedPrecon = 0;
#########################################################################
# Inflow Parameter #
#########################################################################
// refMachNumber: Mach number per meter.
// attackd: Angle of attack.
// angleSlide: Angle of sideslip.
// 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.
// inflowParaType: the type of inflow parameters.
// 0 - the nondimensional conditions.
// 1 - the flight conditions.
// height: Fly height, unit of km.
// gridScaleFactor: The customizable unit of the grid, default value 1.0 metre.
// forceRefenenceLength, forceRefenenceArea: Reference length and area, independent to grid scale.
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent to grid scale.
double refMachNumber = 0.2;
2021-04-15 09:34:10 +08:00
double attackd = 19.0;
double angleSlide = 0.00;
2021-03-22 16:44:54 +08:00
2021-04-15 09:34:10 +08:00
int inflowParaType = 0;
double refReNumber = 9.0e6;
double refDimensionalTemperature = 288.0;
2021-03-22 16:44:54 +08:00
2021-04-15 09:34:10 +08:00
//int inflowParaType = 1;
//double height = 0.001;
2021-03-22 16:44:54 +08:00
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 *
#*******************************************************************
// iviscous: Viscous model
// 0 - Euler
// 1 - Lamilar
// 3 - 1eq turbulent
// 4 - 2eq turbulent
// viscousName : Laminar or tubulent model
// - "1eq-sa", when iviscous = 3
// - "2eq-kw-menter-sst", when iviscous = 4
// DESType : Type of DES
// 0 - RANS (default);
// 1 - DES;
// 2 - DDES;
// 3 - IDDES;
//int iviscous = 0;
2021-04-15 09:34:10 +08:00
//string viscousName = "Euler";
2021-03-22 16:44:54 +08:00
//int iviscous = 1;
2021-04-15 09:34:10 +08:00
//string viscousName = "laminar";
2021-03-22 16:44:54 +08:00
//int iviscous = 3;
2021-04-15 09:34:10 +08:00
//string viscousName = "1eq-sa";
2021-03-22 16:44:54 +08:00
int viscousType = 4;
2021-04-15 09:34:10 +08:00
string viscousName = "2eq-kw-menter-sst";
2021-03-22 16:44:54 +08:00
2021-04-15 09:34:10 +08:00
int DESType = 0;
2021-03-22 16:44:54 +08:00
int roeEntropyFixMethod = 2;
2021-04-15 09:34:10 +08:00
double roeEntropyScale = 1.0;
2021-03-22 16:44:54 +08:00
#########################################################################
# Spatial Discretisation #
#########################################################################
#*******************************************************************
# Struct Solver *
#*******************************************************************
// str_limiter_name: Limiter of struct grid
// - "3rdsmooth", "smooth"
// - "nolim" - no limiter
string str_limiter_name = "smooth";
#*******************************************************************
# UnStruct Solver *
#*******************************************************************
// uns_limiter_name: Limiter of Unstruct grid
// - "vencat"
// - "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_limiter_name = "vencat";
double venkatCoeff = 50.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.
2021-04-15 09:34:10 +08:00
int iunsteady = 0;
2021-03-22 16:44:54 +08:00
double CFLEnd = 10.0;
2021-04-15 09:34:10 +08:00
int nLUSGSSweeps = 1;
2021-03-22 16:44:54 +08:00
#########################################################################
# 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.
2021-04-15 09:34:10 +08:00
int nMGLevel = 1;
int flowInitStep = 100;
2021-03-22 16:44:54 +08:00
#########################################################################
# File In or Out #
#########################################################################
// gridfile: The partitioned Grid file path, using relative path,
// which is relative to the working directory.
// IMPORTANT WARNNING: 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.
2021-04-15 09:34:10 +08:00
string gridfile = "./grid/30p30n_str__4.fts";
int isPlotVolumeField = 0;
2021-03-22 16:44:54 +08:00
// ---------------- 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!!!.
// Arriables order must from small to larger.
2021-04-15 09:34:10 +08:00
int nVisualVariables = 8;
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
2021-03-22 16:44:54 +08:00
// limitVariables: limit model.
// 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.
2021-04-15 09:34:10 +08:00
int reconmeth = 1;
int limitVariables = 0;
int limitVector = 0;