p45328170
/
PHengLEI-TestCases
forked from PHengLEI/PHengLEI-TestCases
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

更新F05-F06（参数，文档，流程，网格）；新增D05、A30、A31算例

This commit is contained in:
hechao 2022-06-29 23:05:35 +08:00
parent fda314f0ec
commit edab4af668
107 changed files with 486 additions and 1777063 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

197
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/bin/cfd_para_subsonic.hypara
View File

@ -1,197 +0,0 @@
#########################################################################
# 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 = 60000;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
int intervalStepForce = 50;
int intervalStepRes = 50;
// 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 = 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";
double MUSCLCoefXk = 0.333333;
#*******************************************************************
# 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 tscheme = 8;
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;

52
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/bin/key.hypara
View File

@ -1,52 +0,0 @@
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 = "";

24
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/bin/partition.hypara
View File

@ -1,24 +0,0 @@
// 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 = 1;
int maxproc = 4;
string original_grid_file = "./grid/30p30n_str.fts";
string partition_grid_file = "./grid/30p30n_str__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;

BIN
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/grid/30p30n_str.grd
View File

Binary file not shown.

161
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/grid/30p30n_str.inp
View File

@ -1,161 +0,0 @@
1
8
81 617
A
21
1 41 1 1 4
1 1 65 1 -1
353 353 -65 -1 7
1 1 65 129 -1
281 281 -1 -65 6
1 1 189 129 -1
121 121 -61 -1 8
1 1 221 189 -1
121 121 -93 -61 8
1 1 221 245 -1
121 121 -93 -117 8
1 1 309 245 -1
41 41 65 1 2
1 1 309 341 -1
41 41 -65 -97 2
1 1 -481 -553 -1
41 41 -65 -137 5
1 1 -481 -449 -1
41 41 -65 -33 5
1 1 -373 -449 -1
121 121 -33 -109 4
1 1 -341 -373 -1
121 121 -1 -33 4
41 1 617 617 4
1 1 -553 -617 -1
1 1 -65 -1 7
81 81 553 617 4
81 41 617 617 4
81 81 341 553 4
81 81 245 341 4
81 81 65 245 4
81 81 65 1 4
41 81 1 1 4
41 97
B
6
41 41 65 1 -1
1 1 309 245 1
-1 -41 1 1 -1
81 121 117 117 8
1 1 65 1 2
41 41 -65 -97 -1
1 1 309 341 1
1 1 97 65 2
1 41 97 97 -1
81 121 1 1 4
81 49
C
6
81 49 49 49 -1
81 49 1 1 4
81 81 1 49 2
81 49 1 1 -1
81 49 117 117 8
-1 -49 49 49 -1
1 49 1 1 4
1 1 49 1 2
-49 -1 1 1 -1
49 1 117 117 8
121 141
D
14
121 121 -1 -33 -1
1 1 -341 -373 1
81 121 1 1 -1
1 41 97 97 2
81 49 1 1 -1
81 49 49 49 3
1 1 33 1 2
1 49 1 1 -1
-1 -49 49 49 3
121 121 -33 -109 -1
1 1 -373 -449 1
1 1 77 33 2
1 1 93 77 2
1 1 109 93 2
1 1 125 109 2
1 1 -125 -141 -1
1 1 -17 -1 8
49 1 141 141 -1
113 161 65 65 6
-121 -49 141 141 -1
41 113 65 65 6
121 121 141 109 -1
41 41 -1 -33 5
41 137
E
8
41 41 -1 -33 -1
121 121 141 109 4
1 1 33 1 2
1 41 1 1 -1
1 41 65 65 6
41 41 -65 -33 -1
1 1 -481 -449 1
1 1 65 33 2
41 41 -65 -137 -1
1 1 -481 -553 1
1 1 137 65 2
41 1 137 137 -1
1 41 65 65 7
281 65
F
8
1 1 65 1 2
1 161 1 1 -1
73 233 65 65 7
113 161 65 65 -1
49 1 141 141 4
41 113 65 65 -1
-121 -49 141 141 4
1 41 65 65 -1
1 41 1 1 5
161 281 65 65 -1
1 121 1 1 8
281 281 -1 -65 -1
1 1 65 129 1
161 281 1 1 -1
233 353 65 65 7
353 65
G
7
1 1 -65 -1 -1
1 1 -553 -617 1
1 41 65 65 -1
41 1 137 137 5
73 41 65 65 2
73 233 65 65 -1
1 161 1 1 6
233 353 65 65 -1
161 281 1 1 6
353 353 -65 -1 -1
1 1 65 1 1
1 353 1 1 4
121 117
H
11
1 1 61 17 2
1 1 -17 -1 -1
1 1 -125 -141 4
1 121 1 1 -1
161 281 65 65 6
121 121 -61 -1 -1
1 1 189 129 1
1 1 93 61 2
121 121 -93 -61 -1
1 1 221 189 1
1 1 117 93 2
81 49 117 117 -1
81 49 1 1 3
49 1 117 117 -1
-49 -1 1 1 3
81 121 117 117 -1
-1 -41 1 1 2
121 121 -93 -117 -1
1 1 221 245 1

BIN
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/二维结构30p30n翼型低速绕流(Matrix LUSGS)算例说明文档.pdf
View File

Binary file not shown.

197
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/bin/cfd_para_subsonic.hypara
View File

@ -1,197 +0,0 @@
#########################################################################
# 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.
// 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 = 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";
double MUSCLCoefXk =0.333333;
#*******************************************************************
# 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 tscheme = 8;
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;

31
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/bin/grid_para.hypara
View File

@ -1,31 +0,0 @@
#########################################################################
# 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 = 1;
int axisup = 1;
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";

BIN
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/grid/Mesh4_137_97.grd
View File

Binary file not shown.

10
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/grid/Mesh4_137_97.inp
View File

@ -1,10 +0,0 @@
1
1
137 97
A
5
1 25 1 1 3
25 137 1 1 2
137 137 1 97 6
1 137 97 97 4
1 1 1 97 4

BIN
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/二维结构低速湍流平板(Matrix LUSGS)算例说明文档.pdf
View File

Binary file not shown.

30
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/bin/grid_para.hypara
View File

@ -1,30 +0,0 @@
#########################################################################
# 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.
int gridtype = 1;
int axisup = 1;
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";

BIN
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/grid/rae2822_vis2d.grd
View File

Binary file not shown.

13
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/grid/rae2822_vis2d.inp
View File

@ -1,13 +0,0 @@
1
1
369 65
A
6
-33 -1 1 1 -1
-337 -369 1 1 1
33 337 1 1 2
-337 -369 1 1 -1
-33 -1 1 1 1
369 369 1 65 4
1 369 65 65 4
1 1 1 65 4

BIN
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/二维结构Rae2822翼型跨声速绕流(Matrix LUSGS)算例说明文档.pdf
View File

Binary file not shown.

196
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/bin/cfd_para_subsonic.hypara
View File

@ -1,196 +0,0 @@
#########################################################################
# 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 = 60000;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
int intervalStepForce = 50;
int intervalStepRes = 50;
// 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 = 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 = 6;
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";
double MUSCLCoefXk = 0.333333;
#*******************************************************************
# 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;

31
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/bin/grid_para.hypara
View File

@ -1,31 +0,0 @@
#########################################################################
# 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 = 1;
int axisup = 1;
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";

52
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/bin/key.hypara
View File

@ -1,52 +0,0 @@
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 = "";

24
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/bin/partition.hypara
View File

@ -1,24 +0,0 @@
// 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 = 1;
int maxproc = 4;
string original_grid_file = "./grid/30p30n_str.fts";
string partition_grid_file = "./grid/30p30n_str__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;

BIN
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/grid/30p30n_str.grd
View File

Binary file not shown.

161
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/grid/30p30n_str.inp
View File

@ -1,161 +0,0 @@
1
8
81 617
A
21
1 41 1 1 4
1 1 65 1 -1
353 353 -65 -1 7
1 1 65 129 -1
281 281 -1 -65 6
1 1 189 129 -1
121 121 -61 -1 8
1 1 221 189 -1
121 121 -93 -61 8
1 1 221 245 -1
121 121 -93 -117 8
1 1 309 245 -1
41 41 65 1 2
1 1 309 341 -1
41 41 -65 -97 2
1 1 -481 -553 -1
41 41 -65 -137 5
1 1 -481 -449 -1
41 41 -65 -33 5
1 1 -373 -449 -1
121 121 -33 -109 4
1 1 -341 -373 -1
121 121 -1 -33 4
41 1 617 617 4
1 1 -553 -617 -1
1 1 -65 -1 7
81 81 553 617 4
81 41 617 617 4
81 81 341 553 4
81 81 245 341 4
81 81 65 245 4
81 81 65 1 4
41 81 1 1 4
41 97
B
6
41 41 65 1 -1
1 1 309 245 1
-1 -41 1 1 -1
81 121 117 117 8
1 1 65 1 2
41 41 -65 -97 -1
1 1 309 341 1
1 1 97 65 2
1 41 97 97 -1
81 121 1 1 4
81 49
C
6
81 49 49 49 -1
81 49 1 1 4
81 81 1 49 2
81 49 1 1 -1
81 49 117 117 8
-1 -49 49 49 -1
1 49 1 1 4
1 1 49 1 2
-49 -1 1 1 -1
49 1 117 117 8
121 141
D
14
121 121 -1 -33 -1
1 1 -341 -373 1
81 121 1 1 -1
1 41 97 97 2
81 49 1 1 -1
81 49 49 49 3
1 1 33 1 2
1 49 1 1 -1
-1 -49 49 49 3
121 121 -33 -109 -1
1 1 -373 -449 1
1 1 77 33 2
1 1 93 77 2
1 1 109 93 2
1 1 125 109 2
1 1 -125 -141 -1
1 1 -17 -1 8
49 1 141 141 -1
113 161 65 65 6
-121 -49 141 141 -1
41 113 65 65 6
121 121 141 109 -1
41 41 -1 -33 5
41 137
E
8
41 41 -1 -33 -1
121 121 141 109 4
1 1 33 1 2
1 41 1 1 -1
1 41 65 65 6
41 41 -65 -33 -1
1 1 -481 -449 1
1 1 65 33 2
41 41 -65 -137 -1
1 1 -481 -553 1
1 1 137 65 2
41 1 137 137 -1
1 41 65 65 7
281 65
F
8
1 1 65 1 2
1 161 1 1 -1
73 233 65 65 7
113 161 65 65 -1
49 1 141 141 4
41 113 65 65 -1
-121 -49 141 141 4
1 41 65 65 -1
1 41 1 1 5
161 281 65 65 -1
1 121 1 1 8
281 281 -1 -65 -1
1 1 65 129 1
161 281 1 1 -1
233 353 65 65 7
353 65
G
7
1 1 -65 -1 -1
1 1 -553 -617 1
1 41 65 65 -1
41 1 137 137 5
73 41 65 65 2
73 233 65 65 -1
1 161 1 1 6
233 353 65 65 -1
161 281 1 1 6
353 353 -65 -1 -1
1 1 65 1 1
1 353 1 1 4
121 117
H
11
1 1 61 17 2
1 1 -17 -1 -1
1 1 -125 -141 4
1 121 1 1 -1
161 281 65 65 6
121 121 -61 -1 -1
1 1 189 129 1
1 1 93 61 2
121 121 -93 -61 -1
1 1 221 189 1
1 1 117 93 2
81 49 117 117 -1
81 49 1 1 3
49 1 117 117 -1
-49 -1 1 1 3
81 121 117 117 -1
-1 -41 1 1 2
121 121 -93 -117 -1
1 1 221 245 1

BIN
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/二维结构30p30n翼型低速绕流(熵修正)算例说明文档.pdf
View File

Binary file not shown.

52
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/bin/key.hypara
View File

@ -1,52 +0,0 @@
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 = "";

BIN
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/~WRL1422.tmp
View File

Binary file not shown.

BIN
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/二维结构圆柱超声速绕流(熵修正)算例说明文档.pdf
View File

Binary file not shown.

1396
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/bin/cfd_para.hypara
View File

File diff suppressed because it is too large Load Diff

BIN
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/grid/rae2822_vis2d.grd
View File

Binary file not shown.

13
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/grid/rae2822_vis2d.inp
View File

@ -1,13 +0,0 @@
1
1
369 65
A
6
-33 -1 1 1 -1
-337 -369 1 1 1
33 337 1 1 2
-337 -369 1 1 -1
-33 -1 1 1 1
369 369 1 65 4
1 369 65 65 4
1 1 1 65 4

BIN
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/二维结构Rae2822翼型跨声速绕流(熵修正)算例说明文档.pdf
View File

Binary file not shown.

0
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/bin/cfd_para.hypara → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/bin/cfd_para.hypara
View File

51
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/bin/cfd_para_supersonic.hypara → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/bin/cfd_para_hypersonic.hypara
View File

@ -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 = 5000;
int maxSimuStep = 10000;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
@ -26,8 +26,11 @@ int intervalStepRes = 10;
// 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:
@ -40,15 +43,16 @@ int intervalStepRes = 10;
// forceReferenceLength, forceReferenceLengthSpanWise, forceReferenceArea: Reference length, SpanWise length and area, independent of grid unit.
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
double refMachNumber = 3.0;
double refMachNumber = 10.00;
double attackd = 0.00;
double angleSlide = 0.00;
double wallTemperature = -1.0;
int inflowParaType = 0;
double refReNumber = 1.0e5;
double refReNumber = 1.835e5;
double refDimensionalTemperature = 293;
double freestream_vibration_temperature = 300.0;
//int inflowParaType = 1;
//double height = 0.001;
@ -57,6 +61,16 @@ double refDimensionalTemperature = 293;
//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)).
//The velocity, temperature and density are fixed.
//int inflowParaType = 4;
//double refDimensionalVelocity = 1000.0;
//double refDimensionalDensity = 1.0e3;
//The velocity, temperature and pressure are fixed.
//int inflowParaType = 5;
//double refDimensionalVelocity = 1000.0;
//double refDimensionalPressure = 1.0e5;
double gridScaleFactor = 1.0;
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
@ -69,7 +83,7 @@ double TorqueRefZ = 0.0; // unit of meter.
#########################################################################
# Physical models #
#########################################################################
// viscousType : Viscous model.
// viscousType: Viscous model.
// 0 -- Euler.
// 1 -- Lamilar.
// 3 -- 1eq turbulent.
@ -83,11 +97,11 @@ double TorqueRefZ = 0.0; // unit of meter.
// 2 -- DDES.
// 3 -- IDDES.
int viscousType = 0;
string viscousName = "Euler";
//int viscousType = 0;
//string viscousName = "Euler";
//int viscousType = 1;
//string viscousName = "laminar";
int viscousType = 1;
string viscousName = "laminar";
//int viscousType = 3;
//string viscousName = "1eq-sa";
@ -109,7 +123,7 @@ double roeEntropyScale = 1.0;
// Using this when solve structered grid or hybrid.
// -- "vanleer", "steger", "ausmpw", "ausmpw+".
// str_limiter_name: Limiter of struct grid.
// -- "minmod", "3rd_minmod_smooth".
// -- "minmod", "3rd_minmod_smooth", "vanalbada".
string inviscidSchemeName = "roe";
string str_limiter_name = "minvan";
@ -145,14 +159,11 @@ double venkatCoeff = 0.5;
int iunsteady = 0;
double CFLStart = 0.1;
double CFLEnd = 20.0;
double CFLEnd = 20.0;
int CFLVaryStep = 100;
double ktmax = 1.0e10;
double ktmax = 1.0e10;
double LUSGSTolerance = 0.1;
#########################################################################
# File In or Out #
#########################################################################
@ -163,7 +174,7 @@ double LUSGSTolerance = 0.1;
// Please use 'rae2822_hybrid2d__4.fts' here!
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
string gridfile = "./grid/cylinder.fts";
string gridfile = "./grid/cylinder.fts";
int plotFieldType = 0;
// ----------------- Advanced Parameters, DO NOT care it ----------------
@ -173,11 +184,13 @@ int plotFieldType = 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),
// -- 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 nVisualVariables = 8;
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// limitVariables: Limit model (It is useful only if limitVector is 0).
@ -190,6 +203,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 = 0;
int reconmeth = 0;
int limitVariables = 0;
int limitVector = 1;
int limitVector = 1;

0
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/bin/grid_para.hypara → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/bin/grid_para.hypara
View File

4
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/bin/key.hypara → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/bin/key.hypara
View File

@ -18,10 +18,10 @@ int ndim = 2;
int nparafile = 1;
int nsimutask = 0;
string parafilename = "./bin/cfd_para_subsonic.hypara";
//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_hypersonic.hypara";
//string parafilename = "./bin/cfd_para_incompressible.hypara";
//int nsimutask = 1;

0
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/grid/cylinder.grd → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/grid/cylinder.grd
View File

0
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/grid/cylinder.inp → A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/grid/cylinder.inp
View File

BIN
A30_TwoD_Cylinder_Inv_Ma10d0_Entropyfix6_Struct_1CPU/二维结构圆柱高超声速绕流(熵修正)算例说明文档.pdf Normal file
View File

Binary file not shown.

0
A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU/bin/cfd_para.hypara → A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/bin/cfd_para.hypara
View File

47
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/bin/cfd_para_transonic.hypara → A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/bin/cfd_para_transonic.hypara
View File

@ -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 = 10000;
int maxSimuStep = 15000;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
int intervalStepForce = 50;
int intervalStepRes = 50;
int intervalStepFlow = 1000;
int intervalStepPlot = 1000;
int intervalStepForce = 100;
int intervalStepRes = 10;
#########################################################################
# Inflow Parameter #
@ -39,13 +39,13 @@ int intervalStepRes = 50;
// 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.73;
double attackd = 2.31;
double refMachNumber = 0.785;
double attackd = 2.61;
double angleSlide = 0.00;
int inflowParaType = 0;
double refReNumber = 6.5e6;
double refDimensionalTemperature = 288.15;
double refReNumber = 1.703665e7;
double refDimensionalTemperature = 298.15;
//int inflowParaType = 1;
//double height = 0.001;
@ -57,9 +57,9 @@ double refDimensionalTemperature = 288.15;
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 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.
@ -86,11 +86,11 @@ double TorqueRefZ = 0.0; // unit of meter.
//int viscousType = 1;
//string viscousName = "laminar";
int viscousType = 3;
string viscousName = "1eq-sa";
//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;
@ -111,7 +111,7 @@ double roeEntropyScale = 1.0;
string inviscidSchemeName = "roe";
string str_limiter_name = "smooth";
double MUSCLCoefXk = 0.333333;
double MUSCLCoefXk = 0.333333;
#*******************************************************************
# UnStruct Solver *
#*******************************************************************
@ -142,11 +142,12 @@ double venkatCoeff = 5.0;
// 1 -- is recommended for structured solver.
// 1-3 -- is recommended for unstructured solver.
int tscheme = 8;
int iunsteady = 0;
double CFLEnd = 10.0;
double CFLEnd = 30.0;
int tscheme = 8;
int nLUSGSSweeps = 1;
#########################################################################
@ -172,7 +173,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/rae2822_vis2d__4.fts";
string gridfile = "./grid/chn-t-structure__256.fts";
int plotFieldType = 0;
// ----------------- Advanced Parameters, DO NOT care it ----------------
@ -186,7 +187,7 @@ int plotFieldType = 0;
// 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).
@ -199,6 +200,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;

33
A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/bin/grid_para.hypara Normal file
View File

@ -0,0 +1,33 @@
#########################################################################
# 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 -- 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 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/chn-t-structure.cgns";
string out_gfile = "./grid/chn-t-structure.fts";

2
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/bin/key.hypara → A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/bin/key.hypara
View File

@ -14,7 +14,7 @@ string defaultParaFile = "./bin/cfd_para.hypara";
// 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 ndim = 3;
int nparafile = 1;
int nsimutask = 0;

22
A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/bin/partition.hypara Normal file
View File

@ -0,0 +1,22 @@
// 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 pgridtype = 1;
int maxproc = 256;
string original_grid_file = "./grid/chn-t-structure.fts";
string partition_grid_file = "./grid/chn-t-structure__256.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;

1
A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/readme.txt Normal file
View File

@ -0,0 +1 @@
如需该算例网格，请通过邮箱phenglei@126.com联系风雷团队获取。

BIN
A31_ThreeD_CHNT_SST_MatrixLUSGS_Struct_256CPU/三维结构CHNT跨声速流动(Matrix LUSGS)算例说明文档.pdf Normal file
View File

Binary file not shown.

0
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/bin/cfd_para.hypara → D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/bin/cfd_para.hypara
View File

97
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/bin/cfd_para_transonic.hypara → D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/bin/cfd_para_hypersonic.hypara
View File

@ -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 = 10000;
int maxSimuStep = 50000;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
int intervalStepForce = 50;
int intervalStepRes = 50;
int intervalStepFlow = 2000;
int intervalStepPlot = 2000;
int intervalStepForce = 2000;
int intervalStepRes = 10;
#########################################################################
# Inflow Parameter #
@ -20,6 +20,13 @@ int intervalStepRes = 50;
// refMachNumber: Mach number.
// attackd: Angle of attack.
// angleSlide: Angle of sideslip.
// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
// dump_Q: Dump out thermal flux Q of solid wall.
// 0 -- no dump out.
// 1 -- dump out wall Q only.
// 2 -- dump out wall Q & the typical position Q of ball.
// 3 -- dump out wall Q & the typical position Q of cone.
// 4 -- dump out wall Q & the typical position Q of double sphere.
// inflowParaType: The type of inflow parameters.
// 0 -- the nondimensional conditions.
// 1 -- the flight conditions.
@ -27,6 +34,7 @@ int intervalStepRes = 50;
// 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:
@ -36,16 +44,20 @@ int intervalStepRes = 50;
// 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.
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
double refMachNumber = 0.73;
double attackd = 2.31;
double refMachNumber = 8.15;
double attackd = 0.00;
double angleSlide = 0.00;
double wallTemperature = 288.0;
int dump_Q = 1;
int inflowParaType = 0;
double refReNumber = 6.5e6;
double refDimensionalTemperature = 288.15;
double refReNumber = 1.67e7;
double refDimensionalTemperature = 56.0;
double freestream_vibration_temperature = 10000.0;
//int inflowParaType = 1;
//double height = 0.001;
@ -54,7 +66,7 @@ double refDimensionalTemperature = 288.15;
//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 gridScaleFactor = 0.01;
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
double forceReferenceLength = 1.0; // unit of meter.
@ -83,18 +95,18 @@ double TorqueRefZ = 0.0; // unit of meter.
//int viscousType = 0;
//string viscousName = "Euler";
//int viscousType = 1;
//string viscousName = "laminar";
int viscousType = 1;
string viscousName = "laminar";
int viscousType = 3;
string viscousName = "1eq-sa";
//int viscousType = 3;
//string viscousName = "1eq-sa";
//int viscousType = 4;
//string viscousName = "2eq-kw-menter-sst";
int DESType = 0;
int roeEntropyFixMethod = 6;
int roeEntropyFixMethod = 3;
double roeEntropyScale = 1.0;
#########################################################################
# Spatial Discretisation #
@ -104,14 +116,13 @@ 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".
// str_limiter_name: Limiter of struct grid.
// -- "3rdsmooth", "smooth".
// -- "nolim", no limiter.
// -- "minmod", "3rd_minmod_smooth".
string inviscidSchemeName = "steger";
string str_limiter_name = "wenn3_prm211";
string inviscidSchemeName = "roe";
string str_limiter_name = "smooth";
double MUSCLCoefXk = 0.333333;
#*******************************************************************
# UnStruct Solver *
#*******************************************************************
@ -126,9 +137,9 @@ double MUSCLCoefXk = 0.333333;
// 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";
double venkatCoeff = 5.0;
double venkatCoeff = 0.5;
#########################################################################
# Temporal Discretisation #
@ -138,28 +149,17 @@ 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.
// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
int iunsteady = 0;
double CFLEnd = 10.0;
double CFLEnd = 2.0;
int CFLVaryStep = 100;
int nLUSGSSweeps = 1;
int nLUSGSSweeps = 4;
double LUSGSTolerance = 1.0e-20;
#########################################################################
# 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;
double ktmax = 1.0e10;
#########################################################################
# File In or Out #
@ -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!
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
int plotFieldType = 1;
string gridfile = "./grid/rae2822_vis2d__4.fts";
int plotFieldType = 0;
string gridfile = "./grid/dbl_yzg__48.fts";
int plotFieldType = 1;
// ----------------- Advanced Parameters, DO NOT care it ----------------
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@ -181,11 +181,11 @@ int plotFieldType = 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), iblank(81).
// -- 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 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,7 @@ 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 = 0;
int limitVariables = 0;
int limitVector = 0;
int limitVector = 1;

6
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/bin/grid_para.hypara → D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/bin/grid_para.hypara
View File

@ -19,12 +19,12 @@
// 8 -- GMSH, *.msh.
int gridtype = 1;
int axisup = 1;
int from_gtype = 3;
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/rae2822_vis2d.grd";
string out_gfile = "./grid/rae2822_vis2d.fts";
string from_gfile = "./grid/dbl_yzg.cgns";
string out_gfile = "./grid/dbl_yzg.fts";

6
F11_ThreeD_M6_GridParallelRefine_Unstruct_4CPU/bin/key.hypara → D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/bin/key.hypara
View File

@ -21,11 +21,11 @@ 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_hypersonic.hypara";
//string parafilename = "./bin/cfd_para_incompressible.hypara";
int nsimutask = 1;
string parafilename = "./bin/grid_para.hypara";
//int nsimutask = 1;
//string parafilename = "./bin/grid_para.hypara";
//int nsimutask = 2;
//string parafilename = "./bin/cfd_para.hypara";

6
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/bin/partition.hypara → D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/bin/partition.hypara
View File

@ -9,10 +9,10 @@
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
int pgridtype = 1;
int maxproc = 4;
int maxproc = 48;
string original_grid_file = "./grid/rae2822_vis2d.fts";
string partition_grid_file = "./grid/rae2822_vis2d__4.fts";
string original_grid_file = "./grid/dbl_yzg.fts";
string partition_grid_file = "./grid/dbl_yzg__48.fts";
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
// 1 -- single level.

BIN
D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/grid/dbl_yzg.cgns Normal file
View File

Binary file not shown.

BIN
D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU/三维结构双椭球高精度WENN格式算例说明文档.pdf Normal file
View File

Binary file not shown.

65
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/boundary_condition_ref.hypara Normal file
View File

@ -0,0 +1,65 @@
# 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 = 5;
string bcName = "missile";
{
string bodyName = "body1";
int bcType = 2;
}
string bcName = "wing";
{
string bodyName = "body0";
int bcType = 2;
}
string bcName = "symmetry";
{
int bcType = 3;
}
string bcName = "farfield";
{
int bcType = 4;
}
string bcName = "overset";
{
int bcType = 1000;
}
# 'bcType' is defined as following:
# -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

0
A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU/bin/cfd_para.hypara → F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/cfd_para.hypara
View File

12
F11_ThreeD_M6_GridParallelRefine_Unstruct_4CPU/bin/grid_para.hypara → F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/grid_para.hypara
View File

@ -17,9 +17,8 @@
// 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 = 1;
int axisup = 2;
int from_gtype = 2;
#########################################################################
@ -27,5 +26,10 @@ int from_gtype = 2;
#########################################################################
// 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_unstr.cgns";
string out_gfile = "./grid/M6_unstr.fts";
int numberOfGridFile = 4;
string from_gfile = "./grid/Missile.cgns";
string from_gfile1 = "./grid/MissileIn.cgns";
string from_gfile2 = "./grid/Wing.cgns";
string from_gfile3 = "./grid/WingIn.cgns";

11
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/key.hypara → F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/key.hypara
View File

@ -23,8 +23,9 @@ int nparafile = 1;
//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;
//string parafilename = "./bin/grid_para.hypara";
@ -35,11 +36,11 @@ int nparafile = 1;
//int nsimutask = 3;
//string parafilename = "./bin/partition.hypara";
//int nsimutask = 4;
//string parafilename = "./bin/repository.hypara";
//int nsimutask = 1;
//string parafilename = "./bin/grid_deform_para.hypara";
//int nsimutask = 5;
//string parafilename = "./bin/overset_grid_view.hypara";
//int nsimutask = 1;
//string parafilename = "./bin/grid_refine_para.hypara";
int nsimutask = 6;
string parafilename = "./bin/overset_config.hypara";

12
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/overset_config.hypara Normal file
View File

@ -0,0 +1,12 @@
int numberOfGridGroups = 2;
string gridfile = "./grid/Wing__8.fts";
string gridfile1 = "./grid/Missile__8.fts";
int codeOfOversetGrid = 1; // 0: 无重叠 1: 重叠;
int symetryOrNot = 1; // 是否只进行半场计算
int readInAuxiliaryInnerGrid = 1; // 是否需要辅助网格（inner）
string auxiliaryInnerGrid0 = "./grid/WingIn.fts";
string auxiliaryInnerGrid1 = "./grid/MissileIn.fts";
int twoOrderInterpolationOrNot = 1; // 1-采用二阶插值（多一排插值点）；0-一阶插值
int keyEnlargeOfActiveNodes = 3; // 活跃区域扩展次数
int outTecplotOverset = 1;

12
F11_ThreeD_M6_GridParallelRefine_Unstruct_4CPU/bin/partition.hypara → F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/bin/partition.hypara
View File

@ -8,11 +8,15 @@
// 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 = 4;
int numberOfGridFile = 2;
int pgridtype = 0;
int pgridtype1 = 0;
int maxproc = 8;
int maxproc1 = 8;
string original_grid_file = "./grid/Missile.fts";
string original_grid_file1 = "./grid/Wing.fts";
string original_grid_file = "./grid/M6_unstr.fts";
string partition_grid_file = "./grid/M6_unstr__4.fts";
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
// 1 -- single level.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/grid/Missile.cgns Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/grid/MissileIn.cgns Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/grid/Wing.cgns Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/grid/WingIn.cgns Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/~$构三维AEDC弹体分离_流场计算_算例说明文档.doc Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/~WRL1112.tmp Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/~WRL3023.tmp Normal file
View File

Binary file not shown.

BIN
F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU/非结构三维AEDC吊舱分离_重叠装配算例说明文档.pdf Normal file
View File

Binary file not shown.

58
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/boundary_condition.hypara
View File

@ -1,58 +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 = 5;
string bcName = "1";
{
int bcType = 2;
}
string bcName = "4";
{
int bcType = 1000;
}
string bcName = "0";
{
int bcType = 2;
}
string bcName = "3";
{
int bcType = 3;
}
string bcName = "2";
{
int bcType = 4;
}
# '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

1220
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/cfd_para.hypara
View File

File diff suppressed because it is too large Load Diff

41
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/grid_para.hypara
View File

@ -1,41 +0,0 @@
#########################################################################
# 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.
int gridtype = 0;
int axisup = 1;
int from_gtype = 5;
int dumpOldGrid = 0;
#########################################################################
# 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/Missile.cas";
//string out_gfile = "./grid/Missile.fts";
//string from_gfile = "./grid/MissileIn.cas";
//string out_gfile = "./grid/MissileIn.fts";
//string from_gfile = "./grid/Wing.cas";
//string out_gfile = "./grid/Wing.fts";
string from_gfile = "./grid/WingIn.cas";
string out_gfile = "./grid/WingIn.fts";

25
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/overset_config.hypara
View File

@ -1,25 +0,0 @@
int parallelStrategy = 1;
int numberOfGridGroups = 2;
string gridfile = "./grid/Wing__6.fts";
string gridfile1 = "./grid/Missile__2.fts";
int codeOfOversetGrid = 1; // 0: 无重叠 1: 重叠;
int codeOfOversetSlipGrid = 0; // 0: 无滑移 1: 滑移;
int readOversetFileOrNot = 0; // 0: 适时装配 1: 读ovs文件
int symetryOrNot = 1; //是否只进行半场计算
int readInAuxiliaryInnerGrid = 1; //是否需要辅助网格（inner）
int readInAuxiliaryOuterGrid = 0; //是否需要辅助网格（outer）
int readInSklFileOrNot = 0; //读入skl文件有助于点属性的准确通讯，但是比较耗时
string auxiliaryInnerGrid0 = "./grid/WingIn.fts";
string auxiliaryInnerGrid1 = "./grid/MissileIn.fts";
string oversetGridFileName = "./grid/overlap.ovs";
double walldistMainZone = 1.0
double toleranceForOversetSearch = 1e-3;
double toleranceForOversetBox = 1e-3;
int twoOrderInterpolationOrNot = 1; //1-采用二阶插值（多一排插值点）；0-一阶插值
int keyEnlargeOfActiveNodes = 1; //活跃区域扩展次数
int outTecplotOverset = 1; //
// kinetic
int numberOfMovingBodies = 2;
int morphing_0 = 0; //0-不变形；1-二维fish
int morphing_1 = 0; //0-不变形；1-二维fish
int morphing_2 = 0; //0-不变形；1-二维fish

24
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/partition.hypara
View File

@ -1,24 +0,0 @@
// 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 = 6;
//string original_grid_file = "./grid/Missile.fts";
//string partition_grid_file = "./grid/Missile__2.fts";
string original_grid_file = "./grid/Wing.fts";
string partition_grid_file = "./grid/Wing__6.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;

213140
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/grid/Missile.cas
View File

File diff suppressed because it is too large Load Diff

44663
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/grid/MissileIn.cas
View File

File diff suppressed because it is too large Load Diff

554423
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/grid/Wing.cas
View File

File diff suppressed because it is too large Load Diff

72784
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/grid/WingIn.cas
View File

File diff suppressed because it is too large Load Diff

BIN
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/三维AEDC弹翼分离非结构重叠网格装配算例说明文档.pdf
View File

Binary file not shown.

3
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/补充说明.txt
View File

@ -1,3 +0,0 @@
针对本算例，风雷代码cmake构建的时候请不要勾选USE_OVERSET选项！！！
请不要勾选USE_OVERSET选项！！！
请不要勾选USE_OVERSET选项！！！

57
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/boundary_condition.hypara
View File

@ -1,57 +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 = 5;
string bcName = "0";
{
int bcType = 2;
}
string bcName = "1";
{
int bcType = 2;
}
string bcName = "4";
{
int bcType = 1000;
}
string bcName = "3";
{
int bcType = 3;
}
string bcName = "2";
{
int bcType = 4;
}
# '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

1220
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/cfd_para.hypara
View File

File diff suppressed because it is too large Load Diff

41
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/grid_para.hypara
View File

@ -1,41 +0,0 @@
#########################################################################
# 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.
int gridtype = 0;
int axisup = 1;
int from_gtype = 5;
int dumpOldGrid = 0;
#########################################################################
# 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/Missile.cas";
//string out_gfile = "./grid/Missile.fts";
//string from_gfile = "./grid/MissileIn.cas";
//string out_gfile = "./grid/MissileIn.fts";
//string from_gfile = "./grid/Wing.cas";
//string out_gfile = "./grid/Wing.fts";
string from_gfile = "./grid/WingIn.cas";
string out_gfile = "./grid/WingIn.fts";

57
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/key.hypara
View File

@ -1,57 +0,0 @@
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 = 3;
int nparafile = 2;
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";
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 = 4;
//string parafilename = "./bin/repository.hypara";
//int nsimutask = 5;
//string parafilename = "./bin/overset_grid_view.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 = "";

25
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/overset_config.hypara
View File

@ -1,25 +0,0 @@
int parallelStrategy = 1;
int numberOfGridGroups = 2;
string gridfile = "./grid/Wing__6.fts";
string gridfile1 = "./grid/Missile__2.fts";
int codeOfOversetGrid = 1; // 0: 无重叠 1: 重叠;
int codeOfOversetSlipGrid = 0; // 0: 无滑移 1: 滑移;
int readOversetFileOrNot = 0; // 0: 适时装配 1: 读ovs文件
int symetryOrNot = 1; //是否只进行半场计算
int readInAuxiliaryInnerGrid = 1; //是否需要辅助网格（inner）
int readInAuxiliaryOuterGrid = 0; //是否需要辅助网格（outer）
int readInSklFileOrNot = 0; //读入skl文件有助于点属性的准确通讯，但是比较耗时
string auxiliaryInnerGrid0 = "./grid/WingIn.fts";
string auxiliaryInnerGrid1 = "./grid/MissileIn.fts";
string oversetGridFileName = "./grid/overlap.ovs";
double walldistMainZone = 1.0
double toleranceForOversetSearch = 1e-3;
double toleranceForOversetBox = 1e-3;
int twoOrderInterpolationOrNot = 1; //1-采用二阶插值（多一排插值点）；0-一阶插值
int keyEnlargeOfActiveNodes = 1; //活跃区域扩展次数
int outTecplotOverset = 1; //
// kinetic
int numberOfMovingBodies = 2;
int morphing_0 = 0; //0-不变形；1-二维fish
int morphing_1 = 0; //0-不变形；1-二维fish
int morphing_2 = 0; //0-不变形；1-二维fish

24
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/partition.hypara
View File

@ -1,24 +0,0 @@
// 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 = 6;
//string original_grid_file = "./grid/Missile.fts";
//string partition_grid_file = "./grid/Missile__2.fts";
string original_grid_file = "./grid/Wing.fts";
string partition_grid_file = "./grid/Wing__6.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;

213140
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/grid/Missile.cas
View File

File diff suppressed because it is too large Load Diff

44663
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/grid/MissileIn.cas
View File

File diff suppressed because it is too large Load Diff

554423
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/grid/Wing.cas
View File

File diff suppressed because it is too large Load Diff

72784
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/grid/WingIn.cas
View File

File diff suppressed because it is too large Load Diff

BIN
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/三维AEDC弹翼分离非结构重叠网格流场计算算例说明文档.pdf
View File

Binary file not shown.

3
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/补充说明.txt
View File

@ -1,3 +0,0 @@
针对本算例，风雷代码cmake构建的时候请不要勾选USE_OVERSET选项！！！
请不要勾选USE_OVERSET选项！！！
请不要勾选USE_OVERSET选项！！！

65
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/boundary_condition_ref.hypara Normal file
View File

@ -0,0 +1,65 @@
# 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 = 5;
string bcName = "missile";
{
string bodyName = "body1";
int bcType = 2;
}
string bcName = "wing";
{
string bodyName = "body0";
int bcType = 2;
}
string bcName = "symmetry";
{
int bcType = 3;
}
string bcName = "farfield";
{
int bcType = 4;
}
string bcName = "overset";
{
int bcType = 1000;
}
# 'bcType' is defined as following:
# -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

0
A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU/bin/cfd_para.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/cfd_para.hypara
View File

54
F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU/bin/cfd_para_transonic.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/cfd_para_transonic.hypara
View File

@ -7,11 +7,11 @@
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual 'res.dat' saved.
int maxSimuStep = 50;
int maxSimuStep = 10000;
int intervalStepFlow = 10;
int intervalStepPlot = 10;
int intervalStepForce = 10;
int intervalStepFlow = 500;
int intervalStepPlot = 500;
int intervalStepForce = 100;
int intervalStepRes = 10;
#########################################################################
@ -43,12 +43,12 @@ double refMachNumber = 0.95;
double attackd = 0.00;
double angleSlide = 0.00;
int inflowParaType = 0;
double refReNumber = 7.87E6;
double refDimensionalTemperature = 288.15;
//int inflowParaType = 0;
//double refReNumber = 7.87E6;
//double refDimensionalTemperature = 288.15;
//int inflowParaType = 1;
//double height = 0.001;
int inflowParaType = 1;
double height = 8.0;
//int inflowParaType = 2;
//double refDimensionalTemperature = 6051.024; // The total temperature, T*(1+(refGama-1)*M*M/2).
@ -80,8 +80,8 @@ double TorqueRefZ = 0.0; // unit of meter.
// 2 -- DDES.
// 3 -- IDDES.
//int viscousType = 0;
//string viscousName = "Euler";
int viscousType = 0;
string viscousName = "Euler";
//int viscousType = 1;
//string viscousName = "laminar";
@ -89,12 +89,12 @@ 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;
int roeEntropyFixMethod = 3;
int roeEntropyFixMethod = 6;
double roeEntropyScale = 0.1;
#########################################################################
# Spatial Discretisation #
@ -128,7 +128,7 @@ string str_limiter_name = "smooth";
string uns_scheme_name = "roe";
string uns_limiter_name = "vencat";
double venkatCoeff = 0.5;
double venkatCoeff = 1.0;
#########################################################################
# Temporal Discretisation #
@ -142,11 +142,18 @@ double venkatCoeff = 0.5;
// 1 -- is recommended for structured solver.
// 1-3 -- is recommended for unstructured solver.
int iunsteady = 0;
int iunsteady = 0;
double physicalTimeStepDimensional = 0.003;
int min_sub_iter = 20;
int max_sub_iter = 30;
double tol_sub_iter = 1.0e-5;
int aleStartStrategy = 1;
int ifStartFromSteadyResults = 1;
int nLUSGSSweeps = 2;
double CFLEnd = 0.5;
int nLUSGSSweeps = 1;
int ifLocalTimeStep = 0;
double CFLEnd = 10.0;
int CFLVaryStep = 200;
#########################################################################
# Multi-Grid parameters #
@ -171,7 +178,6 @@ int flowInitStep = 0;
// Please use 'rae2822_hybrid2d__4.fts' here!
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
string gridfile = "./grid/Wing__6.fts";
int plotFieldType = 1;
// ----------------- Advanced Parameters, DO NOT care it ----------------
@ -181,11 +187,11 @@ 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).
// 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 +204,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;

61
F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU/bin/cfd_para_transonic.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/cfd_para_transonic_unsteady.hypara
View File

@ -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 = 50;
int maxSimuStep = 100;
int intervalStepFlow = 10;
int intervalStepPlot = 10;
int intervalStepForce = 10;
int intervalStepRes = 10;
int intervalStepFlow = 5;
int intervalStepPlot = 5;
int intervalStepForce = 1;
int intervalStepRes = 1;
#########################################################################
# Inflow Parameter #
@ -43,12 +43,12 @@ double refMachNumber = 0.95;
double attackd = 0.00;
double angleSlide = 0.00;
int inflowParaType = 0;
double refReNumber = 7.87E6;
double refDimensionalTemperature = 288.15;
//int inflowParaType = 0;
//double refReNumber = 7.87E6;
//double refDimensionalTemperature = 288.15;
//int inflowParaType = 1;
//double height = 0.001;
int inflowParaType = 1;
double height = 8.0;
//int inflowParaType = 2;
//double refDimensionalTemperature = 6051.024; // The total temperature, T*(1+(refGama-1)*M*M/2).
@ -80,8 +80,8 @@ double TorqueRefZ = 0.0; // unit of meter.
// 2 -- DDES.
// 3 -- IDDES.
//int viscousType = 0;
//string viscousName = "Euler";
int viscousType = 0;
string viscousName = "Euler";
//int viscousType = 1;
//string viscousName = "laminar";
@ -89,12 +89,12 @@ 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;
int roeEntropyFixMethod = 3;
int roeEntropyFixMethod = 6;
double roeEntropyScale = 0.1;
#########################################################################
# Spatial Discretisation #
@ -128,7 +128,7 @@ string str_limiter_name = "smooth";
string uns_scheme_name = "roe";
string uns_limiter_name = "vencat";
double venkatCoeff = 0.5;
double venkatCoeff = 1.0;
#########################################################################
# Temporal Discretisation #
@ -142,11 +142,18 @@ double venkatCoeff = 0.5;
// 1 -- is recommended for structured solver.
// 1-3 -- is recommended for unstructured solver.
int iunsteady = 0;
int iunsteady = 1;
double physicalTimeStepDimensional = 0.003;
int min_sub_iter = 20;
int max_sub_iter = 30;
double tol_sub_iter = 1.0e-5;
int aleStartStrategy = 1;
int ifStartFromSteadyResults = 1;
int nLUSGSSweeps = 2;
double CFLEnd = 0.5;
int nLUSGSSweeps = 1;
int ifLocalTimeStep = 0;
double CFLEnd = 10.0;
int CFLVaryStep = 20;
#########################################################################
# Multi-Grid parameters #
@ -169,10 +176,10 @@ int flowInitStep = 0;
// 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.
// isPlotVolumeField: If dump out the whole field results to tecplot or not, 0 / 1.
string gridfile = "./grid/Wing__6.fts";
int plotFieldType = 1;
int plotFieldType = 1;
int visualfileType = 0;
// ----------------- Advanced Parameters, DO NOT care it ----------------
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
@ -181,11 +188,11 @@ 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).
// 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 +205,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;

16
A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU/bin/grid_para.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/grid_para.hypara
View File

@ -17,15 +17,19 @@
// 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 = 1;
int axisup = 1;
int from_gtype = 3;
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/30p30n_str.grd";
string out_gfile = "./grid/30p30n_str.fts";
int numberOfGridFile = 4;
string from_gfile = "./grid/Missile.cgns";
string from_gfile1 = "./grid/MissileIn.cgns";
string from_gfile2 = "./grid/Wing.cgns";
string from_gfile3 = "./grid/WingIn.cgns";

10
A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU/bin/key.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/key.hypara
View File

@ -14,8 +14,9 @@ string defaultParaFile = "./bin/cfd_para.hypara";
// 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;
// 4 -- Knowledge repository / examples of PHengLEI-API.
int ndim = 3;
int nparafile = 3;
int nsimutask = 0;
//string parafilename = "./bin/cfd_para_subsonic.hypara";
@ -23,6 +24,8 @@ 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";
string parafilename2 = "./bin/kinetic_para.hypara";
//int nsimutask = 1;
//string parafilename = "./bin/grid_para.hypara";
@ -39,6 +42,9 @@ string parafilename = "./bin/cfd_para_transonic.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";

47
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/kinetic_para.hypara Normal file
View File

@ -0,0 +1,47 @@
int codeOfAleModel = 1;
#########################################################################
# motive information #
#########################################################################
int numberOfMovingBodies = 2;
double referenceLength = 1.0;
double referenceVelocity = 292.715;
double referenceDensity = 0.526445;
int methodForKineticEquation = 1;
############################## body0 ##############################
//部件的质量
double mass_0 = 0;
//部件的质量矩阵 Ixx Iyy Izz Ixy Ixz Iyz
double massMatrix_0[] = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0;
//部件的初始六自由度位置信息 xc yc zc
double massCenter_0[] = 0.0, 0.0, 0.0;
//部件的初始六自由度位置信息 angleX angleY angleZ
double attitudeAngle_0[] = 0.0 , 0.0, 0.0;
//部件的初始六自由度运动信息 vc vy vz
double massCenterVelocity_0[] = 0.0, 0.0, 0.0;
//部件的初始六自由度运动信息 omigX omigY omigZ
double angularVelocity_0[] = 0.0, 0.0, 0.0;
//部件所属的物体
int fartherIndex_0 = -1;
//部件运动方式
int RBDMethod_0 = 0;
############################## body1 ##############################
//部件的质量
double mass_1 = 907.184;
//部件的质量矩阵 Ixx Iyy Izz Ixy Ixz Iyz
double massMatrix_1[] = 27.116, 488.094, 488.094, 0.0, 0.0, 0.0;
//部件的初始六自由度位置信息 xc yc zc
double massCenter_1[] = 4.2615, -0.89916, 3.302;
//部件的初始六自由度位置信息 angleX angleY angleZ
double attitudeAngle_1[] = 0.0 , 0.0, 0.0;
//部件的初始六自由度运动信息 vc vy vz
double massCenterVelocity_1[] = 0.0, 0.0, 0.0;
//部件的初始六自由度运动信息 omigX omigY omigZ
double angularVelocity_1[] = 0.0, 0.0, 0.0;
//部件所属的物体
int fartherIndex_1 = -1;
//部件运动方式
int RBDMethod_1 = 1;

12
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/overset_config.hypara Normal file
View File

@ -0,0 +1,12 @@
int numberOfGridGroups = 2;
string gridfile = "./grid/Wing__8.fts";
string gridfile1 = "./grid/Missile__8.fts";
int codeOfOversetGrid = 1; // 0: 无重叠 1: 重叠;
int symetryOrNot = 1; // 是否只进行半场计算
int readInAuxiliaryInnerGrid = 1; // 是否需要辅助网格（inner）
string auxiliaryInnerGrid0 = "./grid/WingIn.fts";
string auxiliaryInnerGrid1 = "./grid/MissileIn.fts";
int twoOrderInterpolationOrNot = 1; // 1-采用二阶插值（多一排插值点）；0-一阶插值
int keyEnlargeOfActiveNodes = 3; // 活跃区域扩展次数
int oversetInterpolationMethod = 1;

12
A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU/bin/partition.hypara → F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/bin/partition.hypara
View File

@ -8,11 +8,15 @@
// 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;
int numberOfGridFile = 2;
int pgridtype = 0;
int pgridtype1 = 0;
int maxproc = 8;
int maxproc1 = 8;
string original_grid_file = "./grid/Missile.fts";
string original_grid_file1 = "./grid/Wing.fts";
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.

BIN
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/grid/Missile.cgns Normal file
View File

Binary file not shown.

BIN
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/grid/MissileIn.cgns Normal file
View File

Binary file not shown.

BIN
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/grid/Wing.cgns Normal file
View File

Binary file not shown.

BIN
F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU/grid/WingIn.cgns Normal file
View File

Binary file not shown.

Some files were not shown because too many files have changed in this diff Show More