F06算例说明文档更新

2）D01-D02、F01-F02文件参数更新

A01_TwoD_Plate_Laminar_Struct_1CPU/bin/cfd_para_subsonic.hypara

@@ -0,0 +1,197 @@
+#########################################################################
+#                     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       = 150000;
+
+int intervalStepFlow  = 1000;
+int intervalStepPlot  = 1000;
+int intervalStepForce = 100;
+int intervalStepRes   = 10;
+
+// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
+//                   0 -- no precondition process. (default, mach > 0.3)
+//                   1 -- carry out precondition process. (mach number <= 0.3)
+int ifLowSpeedPrecon = 0;
+
+#########################################################################
+#                       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.1;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 2.0e5;
+double refDimensionalTemperature = 288.15;
+
+//int inflowParaType = 1;
+//double height = 0.001;
+
+double gridScaleFactor = 1.0;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// str_limiter_name: Limiter of struct grid.
+//                   -- "3rdsmooth", "smooth".
+//                   -- "nolim", no limiter.
+
+string str_limiter_name   = "nolim";
+
+#*******************************************************************
+#                       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.
+
+int iunsteady = 0;
+
+double CFLEnd = 100.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+		
+string gridfile   = "./grid/flat_laminar_73_81.fts";
+int plotFieldType = 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;
+
+double MUSCLCoefXk = 0.333333;

A01_TwoD_Plate_Laminar_Struct_1CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0]; 
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/flat_laminar_73_81.grd";
+string out_gfile  = "./grid/flat_laminar_73_81.fts";

A01_TwoD_Plate_Laminar_Struct_1CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A01_TwoD_Plate_Laminar_Struct_1CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/flat_laminar_73_81.fts";
+string partition_grid_file = "./grid/flat_laminar_73_81.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A01_TwoD_Plate_Laminar_Struct_1CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A02_TwoD_Plate_Laminar_Ma5_Struct_1CPU/bin/cfd_para_hypersonic.hypara

@@ -0,0 +1,196 @@
+#########################################################################
+#                     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;
+
+#########################################################################
+#                       Inflow Parameter                                #
+#########################################################################
+// refMachNumber: Mach number.
+// attackd: Angle of attack.
+// angleSlide: Angle of sideslip.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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.
+// freestream_vibration_temperature: Dimensional freestream vibration temperature.
+// 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 = 5.00;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 327.36;
+
+int inflowParaType = 0;
+double refReNumber = 1.5e7;
+double refDimensionalTemperature = 300.33;
+double freestream_vibration_temperature = 10000.0;
+
+//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 = 0.001;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 2.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 = 2;
+double roeEntropyScale  = 0.0002;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// inviscidSchemeName: Spatial discretisation scheme of struct grid.
+//                     Using this when solve structered grid or hybrid.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+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  = "vanleer";
+string uns_limiter_name = "vencat";
+double venkatCoeff      = 0.5;
+
+#########################################################################
+#                       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.
+// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
+
+int iunsteady = 0;
+
+double CFLEnd      = 100;
+int    CFLVaryStep = 500;
+double ktmax       = 1.0e10;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/plate_137_97.fts";
+int plotFieldType = 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),
+//                   -- vibration temperature(Tv, 33), electron temperature(Te, 34), vibrational energy(Ev, 35), electric energy(Ee, 36),
+//                   -- number density of electron(Ne, 37), dimensioanl density(rho, 38), dimensioanl pressure(p, 39), dimensioanl temperature(T, 40).
+// 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      = 0;
+int limitVariables = 0;
+int limitVector    = 1;
+
+double MUSCLCoefXk = 0.333333;

A02_TwoD_Plate_Laminar_Ma5_Struct_1CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order. 
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/plate_137_97.grd";
+string out_gfile  = "./grid/plate_137_97.fts";

A02_TwoD_Plate_Laminar_Ma5_Struct_1CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A02_TwoD_Plate_Laminar_Ma5_Struct_1CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/plate_137_97.fts";
+string partition_grid_file = "./grid/plate_137_97.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A02_TwoD_Plate_Laminar_Ma5_Struct_1CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A03_TwoD_plate_SST_LowMach_Struct/bin/cfd_para_subsonic.hypara

@@ -0,0 +1,198 @@
+#########################################################################
+#                     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       = 100000;
+
+int intervalStepFlow  = 5000;
+int intervalStepPlot  = 5000;
+int intervalStepForce = 100;
+int intervalStepRes   = 100;
+
+// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
+//                   0 -- no precondition process. (default, mach > 0.3)
+//                   1 -- carry out precondition process. (mach number <= 0.3)
+int ifLowSpeedPrecon = 0;
+
+#########################################################################
+#                       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.2;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 5.0e6;
+double refDimensionalTemperature = 288.15;
+
+//int inflowParaType = 1;
+//double height = 0.001;
+
+double gridScaleFactor = 0.001;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 2.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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// str_limiter_name: Limiter of struct grid.
+//                   -- "3rdsmooth", "smooth".
+//                   -- "nolim", no limiter.
+
+string str_limiter_name   = "nolim";
+
+#*******************************************************************
+#                       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.
+
+int iunsteady = 0;
+
+double CFLEnd = 500.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+		
+string gridfile   = "./grid/Mesh4_137_97.fts";
+int plotFieldType = 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;
+
+double MUSCLCoefXk = 0.333333;
+

A03_TwoD_plate_SST_LowMach_Struct/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/Mesh4_137_97.grd";
+string out_gfile  = "./grid/Mesh4_137_97.fts";

A03_TwoD_plate_SST_LowMach_Struct/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A03_TwoD_plate_SST_LowMach_Struct/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/Mesh4_137_97.fts";
+string partition_grid_file = "./grid/Mesh4_137_97.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A03_TwoD_plate_SST_LowMach_Struct/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A04_TwoD_Plate_SST_Ma5_Struct_1CPU/bin/cfd_para_supersonic.hypara

@@ -0,0 +1,192 @@
+#########################################################################
+#                     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  = 2000;
+int intervalStepPlot  = 2000;
+int intervalStepForce = 100;
+int intervalStepRes   = 10;
+
+#########################################################################
+#                       Inflow Parameter                                #
+#########################################################################
+// refMachNumber: Mach number.
+// attackd: Angle of attack.
+// angleSlide: Angle of sideslip.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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 = 5.00;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 327.36;
+
+int inflowParaType = 0;
+double refReNumber = 1.5e7;
+double refDimensionalTemperature = 300.33;
+
+//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 = 0.001;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 2.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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// inviscidSchemeName: Spatial discretisation scheme of struct grid.
+//                     Using this when solve structered grid or hybrid.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+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  = "vanleer";
+string uns_limiter_name = "vencat";
+double venkatCoeff      = 0.5;
+
+#########################################################################
+#                       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.
+// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
+
+int iunsteady = 0;
+
+double CFLEnd = 50.0;
+
+double ktmax  = 1.0e10;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/Mesh4_137_97.fts";
+int plotFieldType = 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      = 0;
+int limitVariables = 0;
+int limitVector    = 1;
+
+double MUSCLCoefXk = 0.333333;

A04_TwoD_Plate_SST_Ma5_Struct_1CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order. 
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/Mesh4_137_97.grd";
+string out_gfile  = "./grid/Mesh4_137_97.fts";

A04_TwoD_Plate_SST_Ma5_Struct_1CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A04_TwoD_Plate_SST_Ma5_Struct_1CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/Mesh4_137_97.fts";
+string partition_grid_file = "./grid/Mesh4_137_97.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A04_TwoD_Plate_SST_Ma5_Struct_1CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A05_TwoD_Cylinder_Laminar_Ma8d03_Struct/bin/cfd_para_hypersonic.hypara

@@ -0,0 +1,195 @@
+#########################################################################
+#                     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       = 30000;
+
+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.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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.
+// freestream_vibration_temperature: Dimensional freestream vibration temperature.
+// 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 = 8.03;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 294.0 ;
+
+int    inflowParaType = 0;
+double refReNumber    = 3.67e5;
+double refDimensionalTemperature = 125.0;
+
+//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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 0.0001;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// inviscidSchemeName: Spatial discretisation scheme of struct grid.
+//                     Using this when solve structered grid or hybrid.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+string inviscidSchemeName = "roe";
+string str_limiter_name   = "vanleer";
+
+#*******************************************************************
+#                       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  = "vanleer";
+string uns_limiter_name = "vencat";
+double venkatCoeff      = 0.5;
+
+#########################################################################
+#                       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.
+// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
+
+int iunsteady = 0;
+
+double CFLEnd = 3.0;
+
+double ktmax  = 1.0e10;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/input.fts";
+int plotFieldType = 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),
+//                   -- vibration temperature(Tv, 33), electron temperature(Te, 34), vibrational energy(Ev, 35), electric energy(Ee, 36),
+//                   -- number density of electron(Ne, 37), dimensioanl density(rho, 38), dimensioanl pressure(p, 39), dimensioanl temperature(T, 40).
+// 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      = 0;
+int limitVariables = 0;
+int limitVector    = 1;
+
+double MUSCLCoefXk = 0.333333;

A05_TwoD_Cylinder_Laminar_Ma8d03_Struct/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order. 
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string   from_gfile = "./grid/input.grd";
+string   out_gfile  = "./grid/input.fts";

A05_TwoD_Cylinder_Laminar_Ma8d03_Struct/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A05_TwoD_Cylinder_Laminar_Ma8d03_Struct/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/input.fts";
+string partition_grid_file = "./grid/input.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A05_TwoD_Cylinder_Laminar_Ma8d03_Struct/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A06_ThreeD_NACA0012_SA_Struct_4CPU/bin/cfd_para_transonic.hypara

@@ -7,12 +7,12 @@
 // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
 // intervalStepRes: The step intervals for residual 'res.dat' saved.
 
-int maxSimuStep       = 20000;
+int maxSimuStep       = 50000;
 
 int intervalStepFlow  = 1000;
 int intervalStepPlot  = 1000;
 int intervalStepForce = 100;
-int intervalStepRes   = 50;
+int intervalStepRes   = 10;
 
 #########################################################################
 #                       Inflow Parameter                                #
@@ -39,13 +39,13 @@ int intervalStepRes   = 50;
 // forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
 // TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
 
-double refMachNumber = 7.11;
-double attackd       = 0.00;
+double refMachNumber = 0.799;
+double attackd       = 2.26;
 double angleSlide    = 0.00;
 
 int inflowParaType = 0;
-double refReNumber = 5.706e6;
-double refDimensionalTemperature = 80;
+double refReNumber = 9.0e6;
+double refDimensionalTemperature = 288.15;
 
 //int inflowParaType = 1;
 //double height = 0.001;
@@ -54,11 +54,11 @@ double refDimensionalTemperature = 80;
 //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 = 0.01;
+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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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.
@@ -66,14 +66,14 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 #########################################################################
 #                       Physical models                                 #
 #########################################################################
-// iviscous: Viscous model.
+// viscousType : 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.
+//              -- "1eq-sa", when viscousType = 3.
+//              -- "2eq-kw-menter-sst", when viscousType = 4.
 // DESType: Type of DES.
 //          0 -- RANS.(default)
 //          1 -- DES.
@@ -110,7 +110,7 @@ double roeEntropyScale  = 1.0;
 //                   -- "nolim", no limiter.
 
 string inviscidSchemeName = "roe";
-string str_limiter_name   = "minvan";
+string str_limiter_name   = "smooth";
 
 #*******************************************************************
 #                       UnStruct Solver                            *
@@ -139,12 +139,12 @@ double venkatCoeff      = 5.0;
 // 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.
-//               0 -- is recommended for structured solver.
+//               1 -- is recommended for structured solver.
 //             1-3 -- is recommended for unstructured solver.
 
 int iunsteady = 0;
 
-double CFLEnd = 10.0;
+double CFLEnd = 30.0;
 
 int nLUSGSSweeps = 1;
 
@@ -169,10 +169,10 @@ int flowInitStep = 100;
 // 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/3D_shock_wave__4.fts";
-int isPlotVolumeField = 1;
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+		
+string gridfile   = "./grid/NACA0012__4.fts";
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@@ -181,11 +181,11 @@ int isPlotVolumeField = 1;
 //                   -- 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).
+//                   -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
 // Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
 // Variables order must from small to big.
 
-int nVisualVariables = 8;
+int nVisualVariables  = 8;
 int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 
 // limitVariables: Limit model (It is useful only if limitVector is 0).
@@ -198,7 +198,6 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 //            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 reconmeth      = 1;
 int limitVariables = 0;
-int limitVector = 0;
-double MUSCLCoefXk = 0.333333;
+int limitVector    = 0;

A06_ThreeD_NACA0012_SA_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0]; 
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/NACA0012.grd";
+string out_gfile  = "./grid/NACA0012.fts";

A06_ThreeD_NACA0012_SA_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A06_ThreeD_NACA0012_SA_Struct_4CPU/bin/partition.hypara

@@ -11,8 +11,8 @@
 int pgridtype = 1;
 int maxproc   = 4;
 
-string original_grid_file  = "./grid/3D_shock_wave.fts";
-string partition_grid_file = "./grid/3D_shock_wave__4.fts";
+string original_grid_file  = "./grid/NACA0012.fts";
+string partition_grid_file = "./grid/NACA0012__4.fts";
 
 // numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
 //                    1 -- single level.

A06_ThreeD_NACA0012_SA_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A07_TwoD_Rae2822_SST_Struct_4CPU/bin/cfd_para_transonic.hypara

@@ -56,9 +56,9 @@ double refDimensionalTemperature = 288.15;
 
 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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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.
@@ -66,14 +66,14 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 #########################################################################
 #                       Physical models                                 #
 #########################################################################
-// iviscous: Viscous model.
+// viscousType : 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.
+//              -- "1eq-sa", when viscousType = 3.
+//              -- "2eq-kw-menter-sst", when viscousType = 4.
 // DESType: Type of DES.
 //          0 -- RANS.(default)
 //          1 -- DES.
@@ -89,8 +89,8 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 //int    viscousType    = 3;
 //string viscousName = "1eq-sa";
 
-int    viscousType  = 4;
-string viscousName        = "2eq-kw-menter-sst";
+int    viscousType    = 4;
+string viscousName = "2eq-kw-menter-sst";
 
 int DESType = 0;
 
@@ -139,7 +139,7 @@ double venkatCoeff      = 5.0;
 // 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.
-//               0 -- is recommended for structured solver.
+//               1 -- is recommended for structured solver.
 //             1-3 -- is recommended for unstructured solver.
 
 int iunsteady = 0;
@@ -169,10 +169,10 @@ int flowInitStep = 100;
 // 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.
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
-string gridfile = "./grid/rae2822_vis2d__4.fts";
-int isPlotVolumeField = 0;
+string gridfile   = "./grid/rae2822_vis2d__4.fts";
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@@ -181,11 +181,11 @@ int isPlotVolumeField = 0;
 //                   -- 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).
+//                   -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
 // Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
 // Variables order must from small to big.
 
-int nVisualVariables = 8;
+int nVisualVariables  = 8;
 int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 
 // limitVariables: Limit model (It is useful only if limitVector is 0).
@@ -198,6 +198,6 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 //            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 reconmeth      = 1;
 int limitVariables = 0;
-int limitVector = 0;
+int limitVector    = 0;

A07_TwoD_Rae2822_SST_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0]; 
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/rae2822_vis2d.grd";
+string out_gfile  = "./grid/rae2822_vis2d.fts";

A07_TwoD_Rae2822_SST_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A07_TwoD_Rae2822_SST_Struct_4CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/rae2822_vis2d.fts";
+string partition_grid_file = "./grid/rae2822_vis2d__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A07_TwoD_Rae2822_SST_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A08_TwoD_30p30n_SST_Struct_4CPU/bin/cfd_para_subsonic.hypara

@@ -0,0 +1,195 @@
+#########################################################################
+#                     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 process to accelerate convergence for low speed flow.
+//                   0 -- no precondition process. (default, mach > 0.3)
+//                   1 -- carry out precondition process. (mach number <= 0.3)
+int ifLowSpeedPrecon = 0;
+
+#########################################################################
+#                       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.2;
+double attackd       = 19.0;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 9.0e6;
+double refDimensionalTemperature = 288.0;
+
+//int inflowParaType = 1;
+//double height = 0.001;
+
+double gridScaleFactor = 1.0;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       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.
+
+int iunsteady = 0;
+
+double CFLEnd = 10.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/30p30n_str__4.fts";
+int plotFieldType = 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;

A08_TwoD_30p30n_SST_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/30p30n_str.grd";
+string out_gfile  = "./grid/30p30n_str.fts";

A08_TwoD_30p30n_SST_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A08_TwoD_30p30n_SST_Struct_4CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/30p30n_str.fts";
+string partition_grid_file = "./grid/30p30n_str__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A08_TwoD_30p30n_SST_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A09_ThreeD_M6_SST_Struct_MG2_4CPU/bin/cfd_para_transonic.hypara

@@ -0,0 +1,205 @@
+#########################################################################
+#                     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       = 10000;
+
+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.8395;
+double attackd       = 3.06;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 1.171e7;
+double refDimensionalTemperature = 288;
+
+//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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+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 = "vencat";
+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 = 10.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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 = 2;
+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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/m6_str__4.fts";
+int plotFieldType = 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), iblank(81).
+// 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;
+
+double MUSCLCoefXk = 0.333333;

A09_ThreeD_M6_SST_Struct_MG2_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes     = 1;
+int axisRotateOrder[]    = [1];
+double axisRotateAngles[]= [90.0];
+int from_gtype = 2;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/m6_str.cgns";
+string out_gfile  = "./grid/m6_str.fts";

A09_ThreeD_M6_SST_Struct_MG2_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A09_ThreeD_M6_SST_Struct_MG2_4CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/m6_str.fts";
+string partition_grid_file = "./grid/m6_str__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 2;

A09_ThreeD_M6_SST_Struct_MG2_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A10_ThreeD_CHNT_SST_Struct_16CPU/bin/cfd_para_transonic.hypara

@@ -56,9 +56,9 @@ double refDimensionalTemperature = 298.15;
 
 double gridScaleFactor = 1.0;
 
-double forceRefenenceLengthSpanWise = 1.0;    // unit of meter.
-double forceRefenenceLength = 0.1937;         // unit of meter.
-double forceRefenenceArea = 0.12892;          // unit of meter^2.
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 0.1937;         // unit of meter.
+double forceReferenceArea = 0.12892;          // unit of meter^2.
 double TorqueRefX = 0.66087;                  // unit of meter. 
 double TorqueRefY = 0.0;                      // unit of meter.
 double TorqueRefZ = 0.0;                      // unit of meter.
@@ -80,16 +80,16 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 //          2 -- DDES.
 //          3 -- IDDES.
 
-//int    viscousType    = 0;
+//int    viscousType = 0;
 //string viscousName = "Euler";
 
-//int    viscousType    = 1;
+//int    viscousType = 1;
 //string viscousName = "laminar";
 
-//int    viscousType    = 3;
+//int    viscousType = 3;
 //string viscousName = "1eq-sa";
 
-int    viscousType    = 4;
+int    viscousType = 4;
 string viscousName = "2eq-kw-menter-sst";
 
 int DESType = 0;
@@ -169,10 +169,10 @@ int flowInitStep = 100;
 // 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.
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
-string gridfile = "./grid/stm-dkbm-c-all__16.fts";
-int isPlotVolumeField = 0;
+string gridfile   = "./grid/stm-dkbm-c-all__16.fts";
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@@ -181,11 +181,11 @@ int isPlotVolumeField = 0;
 //                   -- 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).
+//                   -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
 // Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
 // Variables order must from small to big.
 
-int nVisualVariables = 8;
+int nVisualVariables  = 8;
 int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 
 // limitVariables: Limit model (It is useful only if limitVector is 0).
@@ -198,8 +198,8 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
 //            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 reconmeth      = 1;
 int limitVariables = 0;
-int limitVector = 0;
+int limitVector    = 0;
 
 double MUSCLCoefXk = 0.333333;

A10_ThreeD_CHNT_SST_Struct_16CPU/bin/grid_para.hypara

@@ -4,9 +4,12 @@
 // gridtype: Grid type for generation, conversion, reconstruction, merging.
 //           0 -- Unstructured grid.
 //           1 -- Structured grid.
-// axisup:   Type of Cartisien coordinates system, used in grid conversion.
-//           1 -- Y upward. (default)
-//           2 -- Z upward.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
 // from_gtype: Type of grid data type in grid conversion process.
 //            -1 -- MULTI_TYPE.
 //             1 -- PHengLEI, *.fts.
@@ -18,7 +21,9 @@
 //             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
 //             8 -- GMSH, *.msh.
 int gridtype   = 1;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 3;
 
 #########################################################################

A10_ThreeD_CHNT_SST_Struct_16CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A10_ThreeD_CHNT_SST_Struct_16CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A11_ThreeD_Sphere_Laminar_Ma10_Struct/bin/cfd_para_hypersonic.hypara

@@ -0,0 +1,193 @@
+#########################################################################
+#                     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       = 25000;
+
+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.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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.
+// freestream_vibration_temperature: Dimensional freestream vibration temperature.
+// 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 = 10.00;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 294.0;
+
+int inflowParaType = 0;
+double refReNumber = 1.0e5;
+double refDimensionalTemperature = 79.0;
+
+//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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 0.0001;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// inviscidSchemeName: Spatial discretisation scheme of struct grid.
+//                     Using this when solve structered grid or hybrid.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+string inviscidSchemeName = "steger";
+string str_limiter_name   = "minmod";
+
+#*******************************************************************
+#                       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  = "vanleer";
+string uns_limiter_name = "vencat";
+double venkatCoeff      = 0.5;
+
+#########################################################################
+#                       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.
+// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
+
+int iunsteady = 0;
+
+double CFLEnd   = 10.0;
+int CFLVaryStep = 100;
+double ktmax    = 1.0e10;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/3dball.fts";
+int plotFieldType = 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),
+//                   -- vibration temperature(Tv, 33), electron temperature(Te, 34), vibrational energy(Ev, 35), electric energy(Ee, 36),
+//                   -- number density of electron(Ne, 37), dimensioanl density(rho, 38), dimensioanl pressure(p, 39), dimensioanl temperature(T, 40).
+// 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      = 0;
+int limitVariables = 0;
+int limitVector    = 1;

A11_ThreeD_Sphere_Laminar_Ma10_Struct/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/3dball.grd";
+string out_gfile  = "./grid/3dball.fts";

A11_ThreeD_Sphere_Laminar_Ma10_Struct/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A11_ThreeD_Sphere_Laminar_Ma10_Struct/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 1;
+
+string original_grid_file  = "./grid/3dball.fts";
+string partition_grid_file = "./grid/3dball.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A11_ThreeD_Sphere_Laminar_Ma10_Struct/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A12_ThreeD_DoubleEllipse_Laminar_Struct_4CPU/bin/cfd_para_hypersonic.hypara

@@ -0,0 +1,194 @@
+#########################################################################
+#                     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       = 30000;
+
+int intervalStepFlow  = 2000;
+int intervalStepPlot  = 2000;
+int intervalStepForce = 100;
+int intervalStepRes   = 10;
+
+#########################################################################
+#                       Inflow Parameter                                #
+#########################################################################
+// refMachNumber: Mach number.
+// attackd: Angle of attack.
+// angleSlide: Angle of sideslip.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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.
+// freestream_vibration_temperature: Dimensional freestream vibration temperature.
+// 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 = 8.15;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 296.0;
+
+int inflowParaType = 0;
+double refReNumber = 1.6e7;
+double refDimensionalTemperature = 56.0;
+double freestream_vibration_temperature = 10000.0;
+
+//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 = 0.01;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+string inviscidSchemeName = "steger";
+string str_limiter_name   = "minmod";
+
+#*******************************************************************
+#                       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  = "vanleer";
+string uns_limiter_name = "vencat";
+double venkatCoeff      = 0.5;
+
+#########################################################################
+#                       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.
+// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
+
+int iunsteady = 0;
+
+double CFLEnd = 1.0;
+
+double ktmax  = 1.0e10;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/DoubleEllipse__4.fts";
+int plotFieldType = 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),
+//                   -- vibration temperature(Tv, 33), electron temperature(Te, 34), vibrational energy(Ev, 35), electric energy(Ee, 36),
+//                   -- number density of electron(Ne, 37), dimensioanl density(rho, 38), dimensioanl pressure(p, 39), dimensioanl temperature(T, 40).
+// 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      = 0;
+int limitVariables = 0;
+int limitVector    = 1;

A12_ThreeD_DoubleEllipse_Laminar_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/DoubleEllipse.grd";
+string out_gfile  = "./grid/DoubleEllipse.fts";

A12_ThreeD_DoubleEllipse_Laminar_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A12_ThreeD_DoubleEllipse_Laminar_Struct_4CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/DoubleEllipse.fts";
+string partition_grid_file = "./grid/DoubleEllipse__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A12_ThreeD_DoubleEllipse_Laminar_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/bin/boundary_condition_ref.hypara

@@ -1,8 +1,8 @@
-# nBoundaryConditons : number of global boundary conditions.
+# nBoundaryConditions : number of global boundary conditions.
 # bcName             : Boundary Condition Name.
 # bcType(in PHengLEI): Boundary Condition Type.
 
-# Account of how to set boundaryconditon.
+# Account of how to set boundarycondition.
 # string bcName = "Farfield";          
 # {                                      
 #   int bcType = 4;                      
@@ -13,7 +13,7 @@
 #   double angleSlide = 0;               
 # }                                      
 
-int nBoundaryConditons = 4;
+int nBoundaryConditions = 4;
 string bcName = "SOLID_SURFACE";
 {
   int bcType = 2;

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/bin/cfd_para_subsonic.hypara

@@ -57,9 +57,9 @@ double refDimensionalTemperature = 298.33;
 
 double gridScaleFactor = 0.0127;
 
-double forceRefenenceLengthSpanWise = 1.0;    // unit of meter.
-double forceRefenenceLength = 1.0;            // unit of meter.
-double forceRefenenceArea = 1.0;              // unit of meter^2.
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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.
@@ -67,14 +67,14 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 #########################################################################
 #                       Physical models                                 #
 #########################################################################
-// iviscous: Viscous model.
+// viscousType : 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.
+//              -- "1eq-sa", when viscousType = 3.
+//              -- "2eq-kw-menter-sst", when viscousType = 4.
 // DESType: Type of DES.
 //          0 -- RANS.(default)
 //          1 -- DES.
@@ -107,7 +107,7 @@ double roeEntropyScale  = 1.0;
 //                   -- "3rdsmooth", "smooth".
 //                   -- "nolim", no limiter.
 
-string str_limiter_name   = "nolim";
+string str_limiter_name = "nolim";
 
 #*******************************************************************
 #                       UnStruct Solver                            *
@@ -161,10 +161,10 @@ int flowInitStep = 100;
 // 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.
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
 string gridfile = "./grid/2D_BackwardStep__4.fts";
-int isPlotVolumeField = 0;
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,35 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 3;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/2D_BackwardStep.grd";
+string out_gfile  = "./grid/2D_BackwardStep.fts";

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/bin/partition.hypara

@@ -0,0 +1,21 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/2D_BackwardStep.fts";
+string partition_grid_file = "./grid/2D_BackwardStep__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;

A13_TwoD_BackwardStep_PressureOutlet_SA_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A14_TwoD_Plate_TotalPressTempBC_SA_Struct_4CPU/bin/boundary_condition_ref.hypara

@@ -1,8 +1,8 @@
-# nBoundaryConditons : number of global boundary conditions.
+# nBoundaryConditions : number of global boundary conditions.
 # bcName             : Boundary Condition Name.
 # bcType(in PHengLEI): Boundary Condition Type.
 
-# Account of how to set boundaryconditon.
+# Account of how to set boundarycondition.
 # string bcName = "Farfield";          
 # {                                      
 #   int bcType = 4;                      
@@ -13,7 +13,7 @@
 #   double angleSlide = 0;               
 # }                                      
 
-int nBoundaryConditons = 5;
+int nBoundaryConditions = 5;
 string bcName = "SOLID_SURFACE";
 {
   int bcType = 2;
@@ -26,14 +26,14 @@ string bcName = "FARFIELD";
 {
   int bcType = 4;
 }
-string bcName = "Inflow";
+string bcName = "INFLOW";
 {
   int bcType = 52;
   double totalPressure = 118309.784;
   double totalTemperature = 302.4;
   double direction_inlet[] = 1, 0, 0;
 }
-string bcName = "Outflow";
+string bcName = "OUTFLOW";
 {
   int bcType = 62;
   double staticPressure = 115056.0;

A14_TwoD_Plate_TotalPressTempBC_SA_Struct_4CPU/bin/cfd_para_subsonic.hypara

@@ -55,14 +55,14 @@ double refDimensionalTemperature = 300.00;
 //double height = 0.001;
 
 //int inflowParaType = 2;
-//double refDimensionalTemperature = 302.4;    // The total temperature, T*(1+(refGama-1)*M*M/2).
+//double refDimensionalTemperature = 302.4;      // The total temperature, T*(1+(refGama-1)*M*M/2).
 //double refDimensionalPressure = 118309.784;    // 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 forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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.
@@ -70,30 +70,30 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 #########################################################################
 #                       Physical models                                 #
 #########################################################################
-// viscousType: Viscous model.
+// viscousType : 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.
+//              -- "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;
+//int    viscousType = 0;
 //string viscousName = "Euler";
 
-//int    viscousType    = 1;
+//int    viscousType = 1;
 //string viscousName = "laminar";
 
-int    viscousType    = 3;
+int    viscousType = 3;
 string viscousName = "1eq-sa";
 
-//int    viscousType    = 4;
+//int    viscousType = 4;
 //string viscousName = "2eq-kw-menter-sst";
 
 int DESType = 0;
@@ -110,7 +110,7 @@ double roeEntropyScale  = 1.0;
 //                   -- "3rdsmooth", "smooth".
 //                   -- "nolim", no limiter.
 
-string str_limiter_name   = "smooth";
+string str_limiter_name = "smooth";
 
 #*******************************************************************
 #                       UnStruct Solver                            *
@@ -134,7 +134,7 @@ double venkatCoeff      = 50.0;
 // 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.
-//               0 -- is recommended for structured solver.
+//               1 -- is recommended for structured solver.
 //             1-3 -- is recommended for unstructured solver.
 
 int iunsteady = 0;
@@ -164,10 +164,10 @@ int flowInitStep = 100;
 // 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.
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
 string gridfile = "./grid/mesh0_545x385_plate__4.fts";
-int isPlotVolumeField = 0;
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.

A14_TwoD_Plate_TotalPressTempBC_SA_Struct_4CPU/bin/grid_para.hypara

@@ -4,9 +4,12 @@
 // gridtype: Grid type for generation, conversion, reconstruction, merging.
 //           0 -- Unstructured grid.
 //           1 -- Structured grid.
-// axisup:   Type of Cartisien coordinates system, used in grid conversion.
-//           1 -- Y upward. (default)
-//           2 -- Z upward.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
 // from_gtype: Type of grid data type in grid conversion process.
 //            -1 -- MULTI_TYPE.
 //             1 -- PHengLEI, *.fts.
@@ -18,7 +21,9 @@
 //             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
 //             8 -- GMSH, *.msh.
 int gridtype   = 1;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 3;
 
 #########################################################################

A14_TwoD_Plate_TotalPressTempBC_SA_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 2;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+//string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A14_TwoD_Plate_TotalPressTempBC_SA_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A15_ThreeD_ShockWave_PeriodicBoundary_SA_Struct_4CPU/bin/cfd_para_hypersonic.hypara

@@ -0,0 +1,204 @@
+#########################################################################
+#                     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   = 50;
+
+#########################################################################
+#                       Inflow Parameter                                #
+#########################################################################
+// refMachNumber: Mach number.
+// attackd: Angle of attack.
+// angleSlide: Angle of sideslip.
+// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
+// 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 = 7.11;
+double attackd       = 0.00;
+double angleSlide    = 0.00;
+
+double wallTemperature = 311.0;
+int inflowParaType = 0;
+double refReNumber = 5.706e6;
+double refDimensionalTemperature = 80;
+
+//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 = 0.01;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       Spatial Discretisation                          #
+#########################################################################
+#*******************************************************************
+#                       Struct Solver                              *
+#*******************************************************************
+// inviscidSchemeName: Spatial discretisation scheme of struct grid.
+//                     Using this when solve structered grid or hybrid.
+//                     -- "vanleer", "steger", "ausmpw".
+// str_limiter_name: Limiter of struct grid.
+//                   -- "minmod", "3rd_minmod_smooth".
+
+string inviscidSchemeName = "roe";
+string str_limiter_name   = "minvan";
+
+#*******************************************************************
+#                       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 = "vencat";
+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 = 10.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile = "./grid/3D_shock_wave__4.fts";
+int plotFieldType = 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),
+//                   -- vibration temperature(Tv, 33), electron temperature(Te, 34), vibrational energy(Ev, 35), electric energy(Ee, 36),
+//                   -- number density of electron(Ne, 37), dimensioanl density(rho, 38), dimensioanl pressure(p, 39), dimensioanl temperature(T, 40).
+// 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;
+double MUSCLCoefXk = 0.333333;
+
+// ----------------- Periodic Parameters --------------------------------
+// Notice: Rotational periodicity only support rotation along the X axis!
+// periodicType: Which periodic boundary is used.
+//               0 -- without Periodic Boundary.
+//               1 -- Translational periodicity.
+//               2 -- Rotational periodicity.
+int periodicType = 2; 
+double translationLength[] = [0.0,0.0,0.0];
+double rotationAngle = 1.0;

A15_ThreeD_ShockWave_PeriodicBoundary_SA_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,45 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// gridtype: Grid type for generation, conversion, reconstruction, merging.
+//           0 -- Unstructured grid.
+//           1 -- Structured grid.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
+// axisRotateOrder : axis rotating order.
+//             1 -- X-axis.
+//             2 -- Y-axis.
+//             3 -- Z-axis.
+// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
+// from_gtype: Type of grid data type in grid conversion process.
+//            -1 -- MULTI_TYPE.
+//             1 -- PHengLEI, *.fts.
+//             2 -- CGNS, *.cgns.
+//             3 -- Plot3D type of structured grid, *.dat/*.grd.
+//             4 -- Fieldview type of unstructured grid, *.dat/*.inp.
+//             5 -- Fluent, *.cas/*.msh.
+//             6 -- Ustar, mgrid.in.
+//             7 -- Hybrid, include both of unstructured and structured grid, *.fts.
+//             8 -- GMSH, *.msh.
+int gridtype   = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 2;
+
+#########################################################################
+#                           File path                                   #
+#########################################################################
+// from_gfile: path of original data file for unstructure grid convert from.
+// out_gfile:  path of target file for grid convert to, *.fts type of file usually.
+string from_gfile = "./grid/3D_shock_wave.cgns";
+string out_gfile  = "./grid/3D_shock_wave.fts";
+
+// ----------------- Periodic Parameters --------------------------------
+// Notice: Rotational periodicity only support rotation along the X axis!
+// periodicType: Which periodic boundary is used.
+//               0 -- without Periodic Boundary.
+//               1 -- Translational periodicity.
+//               2 -- Rotational periodicity.
+int periodicType = 2; 
+double translationLength[] = [0.0,0.0,0.0];
+double rotationAngle = 1.0;

A15_ThreeD_ShockWave_PeriodicBoundary_SA_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,52 @@
+string title = "PHengLEI Main Parameter Control File";
+
+// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
+string defaultParaFile = "./bin/cfd_para.hypara";
+
+// ndim: Dimensional of the grid, 2 or 3.
+// nparafile: the number of parameter files.
+// nsimutask: simulation task type.
+//            0 -- CFD Solver of NS or Turbulation.
+//            1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
+//                 Grid conversion: from other format to PHengLEI format (.fts).
+//                 Grid reconstruction: such as grid adaptation.
+//                 Grid merging: merge two blocks into one block.
+//                 Grid repairing: repair the original grid in order to remove the negative volume cells.
+//            2 -- Wall distance computation for turb-solver.
+//            3 -- Grid partition.
+int ndim      = 3;
+int nparafile = 1;
+
+int    nsimutask    = 0;
+//string parafilename = "./bin/cfd_para_subsonic.hypara";
+//string parafilename = "./bin/cfd_para_transonic.hypara";
+//string parafilename = "./bin/cfd_para_supersonic.hypara";
+string parafilename = "./bin/cfd_para_hypersonic.hypara";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_para.hypara";
+
+//int    nsimutask    = 2;
+//string parafilename = "./bin/cfd_para.hypara";
+
+//int    nsimutask    = 3;
+//string parafilename = "./bin/partition.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_deform_para.hypara";
+
+//int    nsimutask    = 1;
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

A15_ThreeD_ShockWave_PeriodicBoundary_SA_Struct_4CPU/bin/partition.hypara

@@ -0,0 +1,31 @@
+// pgridtype: The grid type. 
+//            0 -- unstruct grid.
+//            1 -- struct grid.
+// maxproc:   The number of partition zones that want to be divided into,
+//            which is equal to the number of CPU processors you want.
+//            Usually, 50~100 thousands structured cells per CPU-Core is suggested.
+//            30~70  thousands unstructured cells per CPU-Core is suggested.
+// original_grid_file:  Original grid file that want to be divided(PHengLEI type, *.fts).
+// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
+
+int pgridtype = 1;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/3D_shock_wave.fts";
+string partition_grid_file = "./grid/3D_shock_wave__4.fts";
+
+// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
+//                    1 -- single level.
+//                    2 -- 2 level. 
+//                    N -- N level,..., et al.
+int numberOfMultigrid = 1;
+
+// ----------------- Periodic Parameters --------------------------------
+// Notice: Rotational periodicity only support rotation along the X axis!
+// periodicType: Which periodic boundary is used.
+//               0 -- without Periodic Boundary.
+//               1 -- Translational periodicity.
+//               2 -- Rotational periodicity.
+int periodicType = 2; 
+double translationLength[] = [0.0,0.0,0.0];
+double rotationAngle = 1.0;

A15_ThreeD_ShockWave_PeriodicBoundary_SA_Struct_4CPU/grid/网格地址.txt

@@ -0,0 +1,5 @@
+红山开源风雷算例库原始网格获取百度网盘链接：
+链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA 
+提取码：w47m
+
+注：plot3D格式网格需同时下载.grd和.inp文件

A16_TwoD_NLR7301_SST_Struct_1CPU/bin/cfd_para_subsonic.hypara

@@ -0,0 +1,195 @@
+#########################################################################
+#                     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       = 80000;
+
+int intervalStepFlow  = 1000;
+int intervalStepPlot  = 1000;
+int intervalStepForce = 100;
+int intervalStepRes   = 10;
+
+// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
+//                   0 -- no precondition process. (default, mach > 0.3)
+//                   1 -- carry out precondition process. (mach number <= 0.3)
+int ifLowSpeedPrecon = 0;
+
+#########################################################################
+#                       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.
+// forceReferenceLength, forceReferenceLengthSpanWise, forceReferenceArea: Reference length, SpanWise length and area, independent of grid unit.
+// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
+
+double refMachNumber = 0.185;
+double attackd       = 6.00;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 2.51e6;
+double refDimensionalTemperature = 288.15;
+
+//int inflowParaType = 1;
+//double height = 0.001;
+
+double gridScaleFactor = 1.0;
+
+double forceReferenceLengthSpanWise = 1.0;    // unit of meter.
+double forceReferenceLength = 1.0;            // unit of meter.
+double forceReferenceArea = 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 = 2;
+double roeEntropyScale  = 1.0;
+#########################################################################
+#                       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.
+
+int iunsteady = 0;
+
+double CFLEnd = 10.0;
+
+int nLUSGSSweeps = 1;
+
+#########################################################################
+#                       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!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile   = "./grid/2D_NLR7301_Str.fts";
+int plotFieldType = 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;