新增B18：CRM-HL高升力外形模拟算例

B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU/README.txt

@@ -0,0 +1,2 @@
+本算例在定常状态下采用RANS方法模拟CRM-HL高升力外形标模（B型网格）在大攻角状态下的气动特性。
+网格下载地址：http://hiliftpw.larc.nasa.gov/Workshop4/grids_downloads.html

B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU/bin/cfd_para_subsonic.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       = 50000;
+int intervalStepFlow  = 1000;
+int intervalStepPlot  = 1000;
+int intervalStepForce = 50；
+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.2;
+double attackd       = 17.05;
+double angleSlide    = 0.00;
+
+int inflowParaType = 0;
+double refReNumber = 7.8369E05;
+double refDimensionalTemperature = 289.44;
+
+//int inflowParaType = 1;
+//double height = 0.001;
+
+double gridScaleFactor = 0.0254;
+
+double forceReferenceLengthSpanWise = 29.38145;    // unit of meter.
+double forceReferenceLength = 7.00532;            // unit of meter.
+double forceReferenceArea = 191.8488;              // unit of meter^2.
+double TorqueRefX = 33.67786;                      // unit of meter. 
+double TorqueRefY = 4.51993;                      // 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                              *
+#*******************************************************************
+// 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 = 100.0;
+
+int nLUSGSSweeps = 5;
+
+#########################################################################
+#                       Multi-Grid parameters                           #
+#########################################################################
+// nMGLevel: The number of Multi-Grid level.
+//           = 1 -- single-level.
+//           > 1 -- multi-level.
+// flowInitStep: Flow initialization step, 0 - 500 is suggested.
+//               Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
+//               Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
+
+int nMGLevel = 1;
+int flowInitStep = 100;
+
+#########################################################################
+#                       File In or Out                                  #
+#########################################################################
+// gridfile: The partitioned Grid file path, using relative path,
+//           which is relative to the working directory.
+// IMPORTANT WARNING: The file index should be ignored,
+//                     e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
+//                     Please use 'rae2822_hybrid2d__4.fts' here!
+// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
+
+string gridfile = "./grid/CRM-HL_40-37__Level-B__1024.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 = 3;
+int visualVariables[] = [ 1, 8, 12];
+
+// 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;

B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU/bin/grid_para.hypara

@@ -0,0 +1,31 @@
+#########################################################################
+#                           Grid data type                              #
+#########################################################################
+// 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.
+// 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.
+//             9 -- Gridgen type of structured grid, *.dat/*.grd.
+int gridtype   = 0;
+int axisup     = 2;
+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/CRM-HL_40-37_Nominal_Unstr_LevelB_Q1.cgns";
+string out_gfile  = "./grid/CRM-HL_40-37_Nominal_Unstr_LevelB_Q1.fts";

B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU/bin/key.hypara

@@ -0,0 +1,58 @@
+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/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    = 5;
+//string parafilename = "./bin/overset_grid_view.hypara";
+
+//int    nsimutask    = 13;
+//string parafilename = "./bin/lbm_para.hypara";
+
+//int    nsimutask    = 14;
+//string parafilename = "./bin/integrative_solver.hypara";
+
+//int    nsimutask    = 15;
+//string parafilename = "./bin/gridcheck_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 = "";

B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU/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 = 0;
+int maxproc   = 1024;
+
+string original_grid_file  = "./grid/CRM-HL_40-37_Nominal_Unstr_LevelB_Q1.fts";
+string partition_grid_file = "./grid/CRM-HL_40-37_Nominal_Unstr_LevelB_Q1__1024.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;

README.md

@@ -71,6 +71,7 @@ B 非结构
 | 二维非结构NLR7301翼型低速绕流 | B15_TwoD_NLR7301_SA_Unstruct_4CPU |
 | 三维非结构空心圆柱裙 | B16_ThreeD_Hollow_Cylinder_Flare_Laminar_Unstruct_16CPU |
 | 二维非结构S-K低速平板转捩模拟 | B17_TwoD_Plate_S-KSR_SST_Unstruct_4CPU |
+| 三维非结构CRM-HL高升力外形模拟 | B18-ThreeD_CRM-HL_Level-B_SA_Unstruct_1024CPU |
 
 C 混合