F系列、G系列算例参数及文档更新；新增F01及F02结构走非结构装配流程算例、新增G08二维储气罐泄漏算例

F01_TwoD_30p30n_OversetGrid_SA_Struct_4CPU/README.txt

@@ -0,0 +1,2 @@
+本算例结构网格重叠装配过程在读网格阶段会在结构网格上重新构建非结构网格信息，然后按照非结构的方式进行重叠网格装配过程
+用户如需使用结构重叠，则需使用2171版本代码及相应算例库版本的F01结构重叠算例

F01_TwoD_30p30n_OversetGrid_SA_Struct_4CPU/bin/grid_para.hypara

@@ -0,0 +1,41 @@
+#########################################################################
+#                           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 -- HyperFLOW( 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.
+int numberOfGridFile = 4;
+string   from_gfile  = "./grid/30p30n-background.cgns";
+string   from_gfile1 = "./grid/30p30n-main.cgns"; 
+string   from_gfile2 = "./grid/30p30n-flap.cgns"; 
+string   from_gfile3 = "./grid/30p30n-slat.cgns"; 
+
+
+

F01_TwoD_30p30n_OversetGrid_SA_Struct_4CPU/bin/key.hypara

@@ -0,0 +1,57 @@
+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.
+//            4 -- Knowledge repository / examples of PHengLEI-API.
+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";
+//string parafilename1  = "./bin/overset_config.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    = 6;
+string parafilename = "./bin/overset_config.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 = "";

F01_TwoD_30p30n_OversetGrid_SA_Struct_4CPU/bin/overset_config.hypara

@@ -0,0 +1,35 @@
+// ----------------- overlap configuration ------------------------------
+// numberOfGridGroups: The number of grid groups.
+// codeOfOversetGrid: Overlapping(overset) grid or not.
+//             0 -- NON-overlapping grid.
+//             1 -- Overlapping grid.
+// readOversetFileOrNo: Whether to read overset-file(.ovs) that has been generated.
+//             0 -- no.
+//             1 -- yes.
+// symetryOrNot: If there exist symetryplanes(XY plane, the coordinate of Z direction is 0) in 
+                 the current overset grid(only for three dimension).
+//             0 -- no.
+//             1 -- yes.
+// readInAuxiliaryInnerGrid: Whether to read auxiliary inner grid
+//             0 -- no.
+//             1 -- yes.
+// twoOrderInterpolationOrNot: The number of interpolated cell layers in the overlapping boundary.
+//             0 -- one layer.
+//             1 -- two layers.
+// keyEnlargeOfActiveNodes: The number of enlarged overset-boundary layers in the buffer region.
+// outTecplotOverset: Whether to dump out the flowfield data after the progress of overset configure.
+//             0 -- no.
+//             1 -- yes.
+
+int    numberOfGridGroups         = 4;
+string gridfile                   = "./grid/30p30n-background__4.fts";
+string gridfile1                  = "./grid/30p30n-main__4.fts";
+string gridfile2                  = "./grid/30p30n-slat__4.fts";
+string gridfile3                  = "./grid/30p30n-flap__4.fts";
+int    codeOfOversetGrid          = 1; 
+int    twoOrderInterpolationOrNot = 1;
+int    keyEnlargeOfActiveNodes    = 0;
+double gridScaleFactor            = 0.001;
+int    outTecplotOverset          = 1;      
+
+

F01_TwoD_30p30n_OversetGrid_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(HyperFLOW/PHengLEI type, *.fts).
+// partition_grid_file    : target partition grid file(HyperFLOW/PHengLEI type, *.fts).
+
+int numberOfGridFile = 4;
+int pgridtype  = 1;
+int pgridtype1 = 1;
+int pgridtype2 = 1;
+int pgridtype3 = 1;
+int maxproc  = 4;
+int maxproc1 = 4;
+int maxproc2 = 4;
+int maxproc3 = 4;
+
+string original_grid_file  = "./grid/30p30n-background.fts";
+string original_grid_file1 = "./grid/30p30n-main.fts"; 
+string original_grid_file2 = "./grid/30p30n-flap.fts"; 
+string original_grid_file3 = "./grid/30p30n-slat.fts"; 
+
+// Number of multi-grid levels, ONLY used for structured grid.
+//               1 -- single level, 2 -- 2 level, N -- N level, ..., et al.
+int   numberOfMultigrid = 1;
+
+
+

F01_TwoD_30p30n_OversetGrid_SA_Struct_4CPU/grid/网格地址.txt

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

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/README.txt

@@ -0,0 +1,2 @@
+本算例结构网格重叠装配过程在读网格阶段会在结构网格上重新构建非结构网格信息，然后按照非结构的方式进行重叠网格装配过程
+用户如需使用结构重叠，则需使用2171版本代码及相应算例库版本的F01结构重叠算例

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/boundary_condition.hypara

@@ -1,64 +0,0 @@
-# nBoundaryConditons : number of global boundary conditions.
-# bcName             : Boundary Condition Name.
-# bcType(in PHengLEI): Boundary Condition Type.
-
-# Account of how to set boundaryconditon.
-# string bcName = "Farfield";          
-# {                                      
-#   int bcType = 4;                      
-#   int inflowParaType = 1;              
-#   double attackd = 0;                  
-#   double refReNumber = 6.5e6;          
-#   double refMachNumber = 3.5;          
-#   double angleSlide = 0;               
-# }                                      
-
-int nBoundaryConditons = 8;
-string bcName = "SOLID_SURFACE";
-{
-  int bcType = 2;
-}
-string bcName = "Wall_11";
-{
-  int bcType = 11;
-}
-string bcName = "Wall_12";
-{
-  int bcType = 12;
-}
-string bcName = "Wall_13";
-{
-  int bcType = 13;
-}
-string bcName = "Wall_8";
-{
-  int bcType = 8;
-}
-string bcName = "SYMMETRY";
-{
-  int bcType = 3;
-}
-string bcName = "FARFIELD";
-{
-  int bcType = 4;
-}
-string bcName = "POLE";
-{
-  int bcType = 7;
-}
-
-# 'bcType' is defined as following:
-#    99: PERIODIC                  
-#    -2: WAKE                      
-#    -1: INTERFACE                 
-#    0 : NO_BOUNDARY_CONDITION     
-#    1 : EXTRAPOLATION             
-#    2 : SOLID_SURFACE             
-#    3 : SYMMETRY                  
-#    4 : FARFIELD                  
-#    5 : INFLOW                    
-#    6 : OUTFLOW                   
-#    52: PRESSURE_INLET            
-#    62: PRESSURE_OUTLET           
-#    61: OUTFLOW_CONFINED          
-#    7 : POLE                      

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/grid_para.hypara

@@ -4,12 +4,15 @@
 // gridtype: Grid type for generation, conversion, reconstruction, merging.
 //               0 -- Unstructured grid.
 //               1 -- Structured grid.
-// axisup:   Type of Cartisien coordinates system, used in grid conversion.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//         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 -- MULTI_TYPE
-//             1 -- PHengLEI, *.fts.
+//               1 -- HyperFLOW( PHengLEI ), *.fts.
 //               2 -- CGNS, *.cgns.
 //               3 -- Plot3D type of structured grid, *.dat/*.grd.
 //               4 -- Fieldview type of unstructured grid, *.dat/*.inp.
@@ -18,13 +21,22 @@
 //               7 -- Hybrid, include both of unstructured and structured grid, *.fts.
 //               8 -- GMSH, *.msh.
 int gridtype         = 1;
-int axisup     = 1;
+int nAxisRotateTimes      = 1;
-int from_gtype = 3;
+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/updatedgrid.grd";
+int numberOfGridFile = 4;
-string out_gfile  = "./grid/updatedgrid.fts";
+string from_gfile    = "./grid/wing.cgns";
+string from_gfile1   = "./grid/missile1.cgns";
+string from_gfile2   = "./grid/missile2.cgns";
+string from_gfile3   = "./grid/missile3.cgns";
+
+
+
+

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/key.hypara

@@ -18,34 +18,36 @@ string defaultParaFile = "./bin/cfd_para.hypara";
 int ndim      = 3;
 int nparafile = 1;
 
-int    nsimutask    = 0;
+//int    nsimutask    = 0;
 //string parafilename = "./bin/cfd_para_subsonic.hypara";
-string parafilename = "./bin/cfd_para_transonic.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";
+//string parafilename = "./bin/cfd_para_incompressible.hypara";
+//string parafilename1  = "./bin/overset_config.hypara";
+//string parafilename2  = "./bin/kinetic_para.hypara";
 
-//int    nsimutask    = 1;
+int    nsimutask    = 1;
-//string parafilename = "./bin/grid_para.hypara";
+string parafilename = "./bin/grid_para.hypara";
 
 //int    nsimutask    = 2;
 //string parafilename = "./bin/cfd_para.hypara";
 
-//int    nsimutask    = 3;
+int    nsimutask    = 3;
-//string parafilename = "./bin/partition.hypara";
+string parafilename = "./bin/partition.hypara";
 
-//int    nsimutask    = 4;
+//int    nsimutask    = 1;
-//string parafilename = "./bin/repository.hypara";
+//string parafilename = "./bin/grid_deform_para.hypara";
 
-//int    nsimutask    = 5;
+//int    nsimutask    = 1;
-//string parafilename = "./bin/overset_grid_view.hypara";
+//string parafilename = "./bin/grid_refine_para.hypara";
+
+int    nsimutask    = 6;
+string parafilename = "./bin/overset_config.hypara";
 
 //int    nsimutask    = 14;
 //string parafilename = "./bin/integrative_solver.hypara";
 
-//int    nsimutask    = 22;
-//string parafilename = "./bin/grid_para.hypara";
-
 //int    nsimutask    = 99;
 //string parafilename = "./bin/post_processing.hypara";
 

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/overset.txt

@@ -1,17 +0,0 @@
-int taskSelector       = 1;
-
-int codeOfLargeScale   = 0;
-int numberOfMultigrid  = 1;
-int numberOfProcessors = 8;
-
-string originalGridFile     = "./grid/duodan.grd";
-string originalBoundaryFile = "./grid/duodan.inp";
-
-int numberOfGridGroups = 4;
-int zoneSpan[numberOfGridGroups] = 34,34,34,32;
-
-int solidBcSize = 4;
-int solidColorList[solidBcSize] = 2,11,12,13;
-
-int outerBcSize = 1;
-int outerColorList[outerBcSize] = 8;

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/overset_config.hypara

@@ -0,0 +1,36 @@
+// ----------------- overlap configuration ------------------------------
+// numberOfGridGroups: The number of grid groups.
+// codeOfOversetGrid: Overlapping(overset) grid or not.
+//             0 -- NON-overlapping grid.
+//             1 -- Overlapping grid.
+// symetryOrNot: If there exist symetryplanes(XY plane, the coordinate of Z direction is 0) in 
+                 the current overset grid(only for three dimension).
+//             0 -- no.
+//             1 -- yes.
+// readInAuxiliaryInnerGrid: Whether to read auxiliary inner grid
+//             0 -- no.
+//             1 -- yes.
+// twoOrderInterpolationOrNot: The number of interpolated cell layers in the overlapping boundary.
+//             0 -- one layer.
+//             1 -- two layers.
+// keyEnlargeOfActiveNodes: The number of enlarged overset-boundary layers in the buffer region.
+// outTecplotOverset: Whether to dump out the flowfield data after the progress of overset configure.
+//             0 -- no.
+//             1 -- yes.
+
+int    numberOfGridGroups = 4;
+string gridfile           = "./grid/wing__8.fts"; 
+string gridfile1          = "./grid/missile1__8.fts";
+string gridfile2          = "./grid/missile2__8.fts";
+string gridfile3          = "./grid/missile3__8.fts";
+int    codeOfOversetGrid        = 1; 
+int    symetryOrNot             = 0; 
+double walldistMainZone          = 1.0;
+double toleranceForOversetSearch = 1.0e-5;
+double toleranceForOversetBox    = 1.0e-5;
+
+int twoOrderInterpolationOrNot = 0;  
+int keyEnlargeOfActiveNodes    = 0;
+int oversetInterpolationMethod = 1;  
+int outTecplotOverset          = 1;
+int outPutOversetVisualization = 0;

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/bin/partition.hypara

@@ -0,0 +1,30 @@
+// 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 numberOfGridFile = 4;
+int pgridtype        = 1;
+int pgridtype1       = 1;
+int pgridtype2       = 1;
+int pgridtype3       = 1;
+int maxproc          = 8;
+int maxproc1         = 8;
+int maxproc2         = 8;
+int maxproc3         = 8;
+
+string original_grid_file  = "./grid/wing.fts";
+string original_grid_file1 = "./grid/missile1.fts";
+string original_grid_file2 = "./grid/missile2.fts";
+string original_grid_file3 = "./grid/missile3.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;

F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU/补充说明.txt

@@ -1,2 +0,0 @@
-针对本算例，风雷代码cmake构建的时候需要勾选USE_OVERSET选项！！！
-该算例当前只支持2171版本

F03_TwoD_Dual0012_OversetConfig_Unstruct_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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

F03_TwoD_Dual0012_OversetConfig_Unstruct_4CPU/bin/overset_config.hypara

@@ -29,6 +29,7 @@ int    codeOfOversetGrid          = 1;
 int    readOversetFileOrNot       = 0;      
 int    symetryOrNot               = 0;      
 int    readInAuxiliaryInnerGrid   = 1;    
+double gridScaleFactor            = 0.001;        
 string auxiliaryInnerGrid0        = "./grid/aux-upper.fts";
 string auxiliaryInnerGrid1        = "./grid/aux-lower.fts";
 string oversetGridFileName        = "./grid/overlap.ovs";

F04_TwoD_Dual0012_OversetGrid_InvisCal_Unstruct_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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 2;
+int nAxisRotateTimes      = 1;
+int axisRotateOrder[]     = [1];
+double axisRotateAngles[] = [90.0];
 int from_gtype = 2;
 
 #########################################################################

F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/overset_config.hypara

@@ -24,8 +24,8 @@ string gridfile1          = "./grid/Missile__8.fts";
 int    codeOfOversetGrid        = 1; 
 int    symetryOrNot             = 1; 
 int    readInAuxiliaryInnerGrid = 1; 
-string auxiliaryInnerGrid0      = ".
+string auxiliaryInnerGrid0      = "./grid/WingIn.fts";
-string auxiliaryInnerGrid1      = ".
+string auxiliaryInnerGrid1      = "./grid/MissileIn.fts";
 
 int twoOrderInterpolationOrNot = 1;  
 int keyEnlargeOfActiveNodes    = 3;  

F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/cfd_para_transonic_unsteady.hypara

@@ -179,7 +179,6 @@ int flowInitStep = 0;
 // isPlotVolumeField: If dump out the whole field results to tecplot or not, 0 / 1.
 
 int plotFieldType  = 1;
-int visualfileType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.

F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 2;
+int nAxisRotateTimes      = 1;
+int axisRotateOrder[]     = [1];
+double axisRotateAngles[] = [90.0];
 int from_gtype = 2;
 
 #########################################################################

F07-TwoD_NACA0012_PitchingMovement_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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//               1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//               2 -- Z upward.
+//             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 -- HyperFLOW( 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 = 2;
 
 #########################################################################

F08-ThreeD_Finner_Laminar_Ma2d5_Struct_16CPU/bin/boundary_condition.hypara

@@ -1,69 +0,0 @@
-# nBoundaryConditions : Number of global boundary conditions.
-# bcName             : Boundary condition name.
-# bcType(in PHengLEI): Boundary condition type.
-
-# How to set boundary condition, for example: 
-# string bcName = "Wall";                   
-# {                                           
-#   int bcType = 2;                           
-#   int viscousType = 1;                      
-#   double wallTemperature = -1.0;            
-#   double uWall = 0.0;                       
-#   double vWall = 0.0;                       
-#   double wWall = 0.0;                       
-# }                                           
-# string bcName = "Inflow";                 
-# {                                           
-#   int bcType = 5;                           
-#   int inflowParaType = 0;                   
-#   double refMachNumber = 0.73;              
-#   double attackd = 2.79;                    
-#   double angleSlide = 0.0;                  
-#   double refReNumber = 6.5e6;               
-#   double refDimensionalTemperature = 288.15;
-# }                                           
-
-# For more information, see examples/bin/boundary_condition.hypara file!!!
-
-int nBoundaryConditions = 6;
-string bcName = "BCWall";
-{
-  string bodyName = "body";
-  int bcType = 2;
-}
-string bcName = "BCFarfield";
-{
-  int bcType = 4;
-}
-string bcName = "BCInflow";
-{
-  int bcType = 5;
-}
-string bcName = "BCOutflow";
-{
-  int bcType = 6;
-}
-string bcName = "BCDegenerateLine";
-{
-  int bcType = 71;
-}
-string bcName = "BCDegenerateLine";
-{
-  int bcType = 73;
-}
-
-# 'bcType' is defined as following:
-#    99: PERIODIC                  
-#    -2: WAKE                      
-#    -1: INTERFACE                 
-#    0 : NO_BOUNDARY_CONDITION     
-#    1 : EXTRAPOLATION             
-#    2 : SOLID_SURFACE             
-#    3 : SYMMETRY                  
-#    4 : FARFIELD                  
-#    5 : INFLOW                    
-#    6 : OUTFLOW                   
-#    52: PRESSURE_INLET            
-#    62: PRESSURE_OUTLET           
-#    61: OUTFLOW_CONFINED          
-#    7 : POLE                      

F08-ThreeD_Finner_Laminar_Ma2d5_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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//               1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//               2 -- Z upward.
+//             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 -- HyperFLOW( 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 = 2;
 
 #########################################################################

F09-ThreeD_Brid_GridDeformationSPRING_Unstruct_1CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

F10-ThreeD_Brid_GridDeformationRBF_Unstruct_1CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

F11_ThreeD_M6_GridParallelRefine_Unstruct_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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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.
@@ -19,7 +22,9 @@
 //             8 -- GMSH, *.msh.
 //             9 -- Gridgen type of structured grid, *.dat/*.grd.
 int gridtype   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

G01_TwoD_Cylinder_Laminar_OneTemperMode_FCW_Struct/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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;
 
 #########################################################################

G02_ThreeD_Electre_Laminar_OneTemperMode_NCW_Struct/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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;
 
 #########################################################################

G05_TwoD_Cylinder_N2-O2_Premixing-Unstruct_1CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 #########################################################################
 #                           File path                                   #

G06_ThreeD_Sphere_N2-O2_Premixing-Unstruct_1CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 #########################################################################
 #                           File path                                   #

G07_TwoD_Pipe_Air-CH4_Mixing-Unstruct_1CPU/bin/cfd_para_subsonic.hypara

@@ -7,7 +7,7 @@
 // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
 // intervalStepRes: The step intervals for residual 'res.dat' saved.
 
-int maxSimuStep       = 30000;
+int maxSimuStep       = 10000;
 int intervalStepFlow  = 1000;
 int intervalStepPlot  = 1000;
 int intervalStepForce = 100;
@@ -71,7 +71,7 @@ double freestream_vibration_temperature = 300.0;
 
 //The velocity, temperature and pressure are fixed.
 int inflowParaType = 5;
-double refDimensionalVelocity = 9.2;
+double refDimensionalVelocity = 340.5;
 double refDimensionalPressure = 101325.0;
 
 double gridScaleFactor = 0.001;
@@ -146,7 +146,7 @@ string str_limiter_name   = "minmod";
 // venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
 //              The smaller the value, the more robust it is.
 
-string uns_scheme_name  = "steger";
+string uns_scheme_name  = "roe";
 string uns_limiter_name = "vencat";
 double venkatCoeff      = 0.5;
 
@@ -163,7 +163,7 @@ double venkatCoeff      = 0.5;
 //             1-3 -- is recommended for unstructured solver.
 
 int iunsteady = 0;
-double CFLEnd = 50.0;
+double CFLEnd = 10.0;
 
 int nLUSGSSweeps = 1;
 
@@ -190,7 +190,7 @@ int flowInitStep = 100;
 //                     Please use 'rae2822_hybrid2d__4.fts' here!
 // plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
-string gridfile = "./grid/multi-species.fts";
+string gridfile = "../grid/multi-species.fts";
 int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------

G07_TwoD_Pipe_Air-CH4_Mixing-Unstruct_1CPU/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.
+// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
-//           1 -- Y upward. (default)
+// axisRotateOrder : axis rotating order.
-//           2 -- Z upward.
+//             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   = 0;
-int axisup     = 1;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
 int from_gtype = 2;
 
 #########################################################################

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/bin/boundary_condition_ref.hypara

@@ -0,0 +1,35 @@
+int nBoundaryConditions = 4;
+string bcName = "farinlet";
+{
+	int bcType = 5;
+	int inflowParaType = 5;
+	double angleSlide = 0;
+	double attackd = 0;
+	int directionMethod = 0;
+	double initMassFraction[] = [1.0,0.0,0.0,0.0];
+	double refDimensionalPressure = 101325;
+	double refDimensionalTemperature = 308.15;
+	double refDimensionalVelocity = 2;
+	string speciesName = "Air,CH4,H2S,CO2";
+}
+string bcName = "outlet";
+{
+	int bcType = 6;
+}
+string bcName = "wall";
+{
+	int bcType = 2;
+	int dumpHingeMoment = 0;
+	int twall_control_select = 0;
+	double wallTemperature = -1.0;
+}
+string bcName = "ycinlet";
+{
+	int bcType = 52;
+	int directionMethod = 1;
+	double direction_inlet[] = [0,1,0];
+	double initMassFraction[] = [0.0,0.75,0.16,0.09];
+	string speciesName = "Air,CH4,H2S,CO2";
+	double totalPressure = 8.3e6;
+	double totalTemperature = 308.15;
+}

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/bin/cfd_para_hypersonic.hypara

@@ -7,11 +7,12 @@
 // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
 // intervalStepRes: The step intervals for residual 'res.dat' saved.
 
-int maxSimuStep       = 1;
+int maxSimuStep       = 5000;
 
-int intervalStepFlow  = 1;
+int intervalStepFlow  = 100;
-int intervalStepPlot  = 1;
+int intervalStepPlot  = 10;
-int intervalStepForce = 1;
+int intervalStepSample =  1;
+int intervalStepForce = 100;
 int intervalStepRes   = 1;
 
 #########################################################################
@@ -20,13 +21,17 @@ int intervalStepRes   = 1;
 // 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.
+//                 4 -- the condition that the velocity, temperature and density are given.
+//                 5 -- the condition that the velocity, temperature and pressure are given.
 // 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:
@@ -36,23 +41,21 @@ int intervalStepRes   = 1;
 //                    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.
+// 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.95;
+double refMachNumber = 0.20;
-double attackd       = 2.79;
+double attackd       = 0.00;
 double angleSlide    = 0.00;
 
-int inflowParaType = 0;
+double wallTemperature = -1.0;
-double refReNumber = 6.5e6;
+double freestream_vibration_temperature = 300.0;
-double refDimensionalTemperature = 288.15;
 
-//int inflowParaType = 1;
+//The velocity, temperature and pressure are fixed.
-//double height = 0.001;
+int inflowParaType = 5;
-
+double refDimensionalVelocity = 2.0;
-//int inflowParaType = 2;
+double refDimensionalPressure = 101325;
-//double refDimensionalTemperature = 6051.024;    // The total temperature, T*(1+(refGama-1)*M*M/2).
+double refDimensionalTemperature = 308.15;
-//double refDimensionalPressure = 4.299696E09;    // The total pressure, p*(T0/T)^(refGama/(refGama-1)).
 
 double gridScaleFactor = 1.0;
 
@@ -62,11 +65,12 @@ 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.
+double knudsenLength = 1.0;		              // unit of meter.
 
 #########################################################################
 #                       Physical models                                 #
 #########################################################################
-// viscousType : Viscous model.
+// viscousType: Viscous model.
 //           0 -- Euler.
 //           1 -- Lamilar.
 //           3 -- 1eq turbulent.
@@ -83,11 +87,11 @@ double TorqueRefZ = 0.0;                      // unit of meter.
 //int    viscousType = 0;
 //string viscousName = "Euler";
 
-//int    viscousType    = 1;
+int    viscousType = 1;
-//string viscousName = "laminar";
+string viscousName = "laminar";
 
-int    viscousType    = 3;
+//int    viscousType = 3;
-string viscousName = "1eq-sa";
+//string viscousName = "1eq-sa";
 
 //int    viscousType = 4;
 //string viscousName = "2eq-kw-menter-sst";
@@ -104,13 +108,12 @@ double roeEntropyScale  = 1.0;
 #*******************************************************************
 // inviscidSchemeName: Spatial discretisation scheme of struct grid.
 //                     Using this when solve structered grid or hybrid.
-//                     -- "roe", "vanleer", "ausm+up", "ausmpw".
+//                     -- "vanleer", "steger", "ausmpw", "ausmpw+".
 // str_limiter_name: Limiter of struct grid.
-//                   -- "3rdsmooth", "smooth".
+//                   -- "minmod", "3rd_minmod_smooth", "vanalbada".
-//                   -- "nolim", no limiter.
 
-string inviscidSchemeName = "roe";
+string inviscidSchemeName = "steger";
-string str_limiter_name   = "smooth";
+string str_limiter_name   = "minmod";
 
 #*******************************************************************
 #                       UnStruct Solver                            *
@@ -126,9 +129,9 @@ string str_limiter_name   = "smooth";
 // 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_scheme_name  = "vanleer";
-string uns_limiter_name = "vencat";
+string uns_limiter_name = "1st";
-double venkatCoeff      = 5.0;
+double venkatCoeff      = 0.5;
 
 #########################################################################
 #                       Temporal Discretisation                         #
@@ -138,28 +141,14 @@ double venkatCoeff      = 5.0;
 //            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;
+int iunsteady = 1;
+double physicalTimeStep =  0.2;
+int min_sub_iter =  50;
+int max_sub_iter =  50;
+double tol_sub_iter =  0.001;
 
-double CFLEnd = 10.0;
+double CFLEnd = 50.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                                  #
@@ -171,8 +160,8 @@ int flowInitStep = 100;
 //                     Please use 'rae2822_hybrid2d__4.fts' here!
 // plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
 
-string gridfile = "./grid/updatedgrid.fts";
+string gridfile = "./grid/yanchongCase14__4.fts";
-int plotFieldType = 1;
+int plotFieldType = 0;
 
 // ----------------- Advanced Parameters, DO NOT care it ----------------
 // nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@@ -181,12 +170,30 @@ int plotFieldType = 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), iblank(81).
+//                   -- 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).
+//                   -- Knudsen number(kn, 60), Damkohler number(Da, 61), vibrational nonequilibrium number(Vi, 62)
 // Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
 // Variables order must from small to big.
 
-int nVisualVariables = 1;
+int nVisualVariables = 13;
-int visualVariables[] = [81];
+int visualVariables[] = [0,1,2,3,4,5,6,38,39,40,53,54,55];
+
+int ifSetDataMonitor =  1;
+
+//-----------the optional parameters list for the wall boundary condition----------------
+// nVisualWallVariables: The number of visual variables on wall.
+// visualWallVariables : dumped variable types, listed as following:
+// -coefficient of pressure(cp, 0), -coefficient of friction(cf, 1), yplus(2), -non-dimensional heat flux(Q_NonDim, 3), -dimensional heat flux(Q_Dim, 4),
+// -pressure on wall(pw, 5), -temperature on wall(Tw, 6), -density on wall(rhow, 7), -heat flux of translational-rotational temperature term(Qtr, 8),
+// -heat flux of species diffusion term(Qs, 9), -heat flux of vibrational temperature term(Qv, 10), -heat flux of electron temperature term(Qe, 11),
+// -species mass fractions(Ns, 12), -x component of wall velocity(Vx, 13), -y component of wall velocity(Vy, 14), -z component of wall velocity(Vz, 15)
+// -slip translational-rotational temperature(Tts, 16), -slip vibrational temperature(Tvs, 17), -slip electron temperature(Tes, 18), -absolute wall velocity(Vs, 19)
+// -Stanton number(St, 20), -coefficient of heat rate(Ch, 21), -temperature jump(deltaT, 22), -Grid Reynolds number on wall(Re_w, 23) , -Knudsen number(Kn_wall, 24)
+
+int nVisualWallVariables = 6;
+int visualWallVariables[] = [0, 1, 2, 3, 4, 5];
 
 // limitVariables: Limit model (It is useful only if limitVector is 0).
 //                 0 -- limit only for pressure and denstiny, then get the min value.
@@ -201,3 +208,18 @@ int visualVariables[] = [81];
 int reconmeth = 1;
 int limitVariables = 0;
 int limitVector = 0;
+
+
+int nchem =  1;
+int nchemsrc =  0;
+int nchemrad =  0;
+int ntmodel =  1;
+int nIdealState =  1;
+
+string gasfile  =  "Gas-Mixture";
+string speciesName  =  "Air,CH4,H2S,CO2";
+string initMassFraction  =  "1.0,0.0,0.0,0.0";
+int nFraction =  1;
+int nLeakageMonitor =  1;
+double totalLeakageVolume =  1000.0;
+double monitorThresholdValue =  0.05;

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/bin/grid_para.hypara

@@ -0,0 +1,36 @@
+#########################################################################
+#                           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.
+//             9 -- Gridgen type of structured grid, *.dat/*.grd.
+int gridtype   = 0;
+int nAxisRotateTimes      = 0;
+int axisRotateOrder[]     = [1, 2, 3];
+double axisRotateAngles[] = [0.0, 0.0, 0.0];
+int from_gtype = 5;
+
+#########################################################################
+#                           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/yanchongCase14.cas";
+string out_gfile  = "./grid/yanchongCase14.fts";

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/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      = 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    = 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    = 99;
+//string parafilename = "./bin/post_processing.hypara";
+
+// ---------------- Advanced Parameters, DO NOT care it ----------------
+int numberOfGridProcessor = 0;
+// ATP read
+//@string parafilename1 = ""
+//@string parafilename2 = "";

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/bin/partition.hypara

@@ -0,0 +1,24 @@
+// 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).
+// numberOfMultifile: The number of partition grid files that want to be dumped out.
+
+int pgridtype = 0;
+int maxproc   = 4;
+
+string original_grid_file  = "./grid/yanchongCase14.fts";
+string partition_grid_file = "./grid/yanchongCase14__4.fts";
+
+int numberOfMultifile = 1;
+
+// 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;

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/bin/probes_XYZ.dat

@@ -0,0 +1,5 @@
+4
+28.678 55.643 0
+80.785 116.821 0
+125.335 158.357 0
+162.272 308.038 0

G08_TwoD_TankLeak_NG_Multispecies_Unstruct_Unsteady_4CPU/grid/网格地址.txt

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