forked from PHengLEI/PHengLEI-TestCases
添加A24-A29,F11新算例;修改A13、A15参数、更新README
This commit is contained in:
parent
edab4af668
commit
35fd719b7a
|
@ -1,48 +0,0 @@
|
||||||
# nBoundaryConditions : number of global boundary conditions.
|
|
||||||
# bcName : Boundary Condition Name.
|
|
||||||
# bcType(in PHengLEI): Boundary Condition Type.
|
|
||||||
|
|
||||||
# Account of how to set boundarycondition.
|
|
||||||
# string bcName = "Farfield";
|
|
||||||
# {
|
|
||||||
# int bcType = 4;
|
|
||||||
# int inflowParaType = 1;
|
|
||||||
# double attackd = 0;
|
|
||||||
# double refReNumber = 6.5e6;
|
|
||||||
# double refMachNumber = 3.5;
|
|
||||||
# double angleSlide = 0;
|
|
||||||
# }
|
|
||||||
|
|
||||||
int nBoundaryConditions = 4;
|
|
||||||
string bcName = "SOLID_SURFACE";
|
|
||||||
{
|
|
||||||
int bcType = 2;
|
|
||||||
}
|
|
||||||
string bcName = "SYMMETRY";
|
|
||||||
{
|
|
||||||
int bcType = 3;
|
|
||||||
}
|
|
||||||
string bcName = "FARFIELD";
|
|
||||||
{
|
|
||||||
int bcType = 4;
|
|
||||||
}
|
|
||||||
string bcName = "OUTFLOW";
|
|
||||||
{
|
|
||||||
int bcType = 6;
|
|
||||||
}
|
|
||||||
|
|
||||||
# '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
|
|
|
@ -1,44 +0,0 @@
|
||||||
# nBoundaryConditions : number of global boundary conditions.
|
|
||||||
# bcName : Boundary Condition Name.
|
|
||||||
# bcType(in PHengLEI): Boundary Condition Type.
|
|
||||||
|
|
||||||
# Account of how to set boundarycondition.
|
|
||||||
# string bcName = "Farfield";
|
|
||||||
# {
|
|
||||||
# int bcType = 4;
|
|
||||||
# int inflowParaType = 1;
|
|
||||||
# double attackd = 0;
|
|
||||||
# double refReNumber = 6.5e6;
|
|
||||||
# double refMachNumber = 3.5;
|
|
||||||
# double angleSlide = 0;
|
|
||||||
# }
|
|
||||||
|
|
||||||
int nBoundaryConditions = 3;
|
|
||||||
string bcName = "Wall_viscous";
|
|
||||||
{
|
|
||||||
int bcType = 2;
|
|
||||||
}
|
|
||||||
string bcName = "Farfield";
|
|
||||||
{
|
|
||||||
int bcType = 4;
|
|
||||||
}
|
|
||||||
string bcName = "Outflow";
|
|
||||||
{
|
|
||||||
int bcType = 6;
|
|
||||||
}
|
|
||||||
|
|
||||||
# '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
|
|
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,197 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 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;
|
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// axisup: Type of Cartisien coordinates system, used in grid conversion.
|
||||||
|
// 1 -- Y upward. (default)
|
||||||
|
// 2 -- Z upward.
|
||||||
|
// from_gtype: Type of grid data type in grid conversion process.
|
||||||
|
// -1 -- MULTI_TYPE.
|
||||||
|
// 1 -- PHengLEI, *.fts.
|
||||||
|
// 2 -- CGNS, *.cgns.
|
||||||
|
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||||
|
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||||
|
// 5 -- Fluent, *.cas/*.msh.
|
||||||
|
// 6 -- Ustar, mgrid.in.
|
||||||
|
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||||
|
// 8 -- GMSH, *.msh.
|
||||||
|
// 9 -- Gridgen type of structured grid, *.dat/*.grd.
|
||||||
|
int gridtype = 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";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
|
@ -0,0 +1,24 @@
|
||||||
|
// pgridtype: The grid type.
|
||||||
|
// 0 -- unstruct grid.
|
||||||
|
// 1 -- struct grid.
|
||||||
|
// maxproc: The number of partition zones that want to be divided into,
|
||||||
|
// which is equal to the number of CPU processors you want.
|
||||||
|
// Usually, 50~100 thousands structured cells per CPU-Core is suggested.
|
||||||
|
// 30~70 thousands unstructured cells per CPU-Core is suggested.
|
||||||
|
// original_grid_file: Original grid file that want to be divided(PHengLEI type, *.fts).
|
||||||
|
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||||
|
// numberOfMultifile: The number of partition grid files that want to be dumped out.
|
||||||
|
|
||||||
|
int pgridtype = 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;
|
Binary file not shown.
|
@ -0,0 +1,161 @@
|
||||||
|
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
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,197 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 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;
|
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// axisup: Type of Cartisien coordinates system, used in grid conversion.
|
||||||
|
// 1 -- Y upward. (default)
|
||||||
|
// 2 -- Z upward.
|
||||||
|
// from_gtype: Type of grid data type in grid conversion process.
|
||||||
|
// -1 -- MULTI_TYPE.
|
||||||
|
// 1 -- PHengLEI, *.fts.
|
||||||
|
// 2 -- CGNS, *.cgns.
|
||||||
|
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||||
|
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||||
|
// 5 -- Fluent, *.cas/*.msh.
|
||||||
|
// 6 -- Ustar, mgrid.in.
|
||||||
|
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||||
|
// 8 -- GMSH, *.msh.
|
||||||
|
// 9 -- Gridgen type of structured grid, *.dat/*.grd.
|
||||||
|
int gridtype = 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";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
Binary file not shown.
|
@ -0,0 +1,10 @@
|
||||||
|
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
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,204 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 10000;
|
||||||
|
|
||||||
|
int intervalStepFlow = 500;
|
||||||
|
int intervalStepPlot = 500;
|
||||||
|
int intervalStepForce = 50;
|
||||||
|
int intervalStepRes = 50;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# 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.73;
|
||||||
|
double attackd = 2.31;
|
||||||
|
double angleSlide = 0.00;
|
||||||
|
|
||||||
|
int inflowParaType = 0;
|
||||||
|
double refReNumber = 6.5e6;
|
||||||
|
double refDimensionalTemperature = 288.15;
|
||||||
|
|
||||||
|
//int inflowParaType = 1;
|
||||||
|
//double height = 0.001;
|
||||||
|
|
||||||
|
//int inflowParaType = 2;
|
||||||
|
//double refDimensionalTemperature = 6051.024; // The total temperature, T*(1+(refGama-1)*M*M/2).
|
||||||
|
//double refDimensionalPressure = 4.299696E09; // The total pressure, p*(T0/T)^(refGama/(refGama-1)).
|
||||||
|
|
||||||
|
double gridScaleFactor = 1.0;
|
||||||
|
|
||||||
|
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||||
|
double forceReferenceLength = 1.0; // unit of meter.
|
||||||
|
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||||
|
double TorqueRefX = 0.0; // unit of meter.
|
||||||
|
double TorqueRefY = 0.0; // unit of meter.
|
||||||
|
double TorqueRefZ = 0.0; // unit of meter.
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Physical models #
|
||||||
|
#########################################################################
|
||||||
|
// viscousType : Viscous model.
|
||||||
|
// 0 -- Euler.
|
||||||
|
// 1 -- Lamilar.
|
||||||
|
// 3 -- 1eq turbulent.
|
||||||
|
// 4 -- 2eq turbulent.
|
||||||
|
// viscousName: Laminar or tubulent model.
|
||||||
|
// -- "1eq-sa", when viscousType = 3.
|
||||||
|
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||||
|
// DESType: Type of DES.
|
||||||
|
// 0 -- RANS.(default)
|
||||||
|
// 1 -- DES.
|
||||||
|
// 2 -- DDES.
|
||||||
|
// 3 -- IDDES.
|
||||||
|
|
||||||
|
//int viscousType = 0;
|
||||||
|
//string viscousName = "Euler";
|
||||||
|
|
||||||
|
//int viscousType = 1;
|
||||||
|
//string viscousName = "laminar";
|
||||||
|
|
||||||
|
int viscousType = 3;
|
||||||
|
string viscousName = "1eq-sa";
|
||||||
|
|
||||||
|
//int viscousType = 4;
|
||||||
|
//string viscousName = "2eq-kw-menter-sst";
|
||||||
|
|
||||||
|
int DESType = 0;
|
||||||
|
|
||||||
|
int roeEntropyFixMethod = 2;
|
||||||
|
double roeEntropyScale = 1.0;
|
||||||
|
#########################################################################
|
||||||
|
# Spatial Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
#*******************************************************************
|
||||||
|
# Struct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||||
|
// Using this when solve structered grid or hybrid.
|
||||||
|
// -- "roe", "vanleer", "ausm+up", "ausmpw".
|
||||||
|
// str_limiter_name: Limiter of struct grid.
|
||||||
|
// -- "3rdsmooth", "smooth".
|
||||||
|
// -- "nolim", no limiter.
|
||||||
|
|
||||||
|
string inviscidSchemeName = "roe";
|
||||||
|
string str_limiter_name = "smooth";
|
||||||
|
double MUSCLCoefXk = 0.333333;
|
||||||
|
#*******************************************************************
|
||||||
|
# UnStruct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||||
|
// Using this when solve Unstructered grid or hybrid.
|
||||||
|
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||||
|
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||||
|
// uns_limiter_name: Limiter of Unstruct grid.
|
||||||
|
// -- "vencat", "barth".
|
||||||
|
// -- "1st", meaning accuracy of first-order.
|
||||||
|
// -- "nolim", no limiter.
|
||||||
|
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||||
|
// The smaller the value, the more robust it is.
|
||||||
|
|
||||||
|
string uns_scheme_name = "roe";
|
||||||
|
string uns_limiter_name = "vencat";
|
||||||
|
double venkatCoeff = 5.0;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Temporal Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
// iunsteady: Steady or unsteady.
|
||||||
|
// 0 -- steady.
|
||||||
|
// 1 -- unsteay.
|
||||||
|
// CFLEnd: The CFL number, [0.1, 100].
|
||||||
|
// The bigger the value, the convergence faster but lower robustness.
|
||||||
|
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
|
||||||
|
// 1 -- is recommended for structured solver.
|
||||||
|
// 1-3 -- is recommended for unstructured solver.
|
||||||
|
|
||||||
|
int iunsteady = 0;
|
||||||
|
|
||||||
|
double CFLEnd = 10.0;
|
||||||
|
|
||||||
|
int 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/rae2822_vis2d__4.fts";
|
||||||
|
int plotFieldType = 0;
|
||||||
|
|
||||||
|
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
|
||||||
|
// visualVariables: Variable types dumped, listed as following:
|
||||||
|
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
|
||||||
|
// -- viscosityLaminar(7), viscosityTurbulent(8),
|
||||||
|
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
|
||||||
|
// -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
|
||||||
|
// -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
|
||||||
|
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
|
||||||
|
// Variables order must from small to big.
|
||||||
|
|
||||||
|
int nVisualVariables = 8;
|
||||||
|
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
|
||||||
|
|
||||||
|
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||||
|
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||||
|
// 1 -- limit for every variables, then get the min value.
|
||||||
|
// limitVector:
|
||||||
|
// 0 -- Each variable use the same limiter coefficient.
|
||||||
|
// 1 -- Each variable use the respective limiter coefficients.
|
||||||
|
// reconmeth:
|
||||||
|
// 0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
|
||||||
|
// 1 -- Q+, Q- use the min limiter coefficients of left and right cell.
|
||||||
|
|
||||||
|
int reconmeth = 1;
|
||||||
|
int limitVariables = 0;
|
||||||
|
int limitVector = 0;
|
|
@ -0,0 +1,30 @@
|
||||||
|
#########################################################################
|
||||||
|
# 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";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
//string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
|
@ -0,0 +1,21 @@
|
||||||
|
// pgridtype: The grid type.
|
||||||
|
// 0 -- unstruct grid.
|
||||||
|
// 1 -- struct grid.
|
||||||
|
// maxproc: The number of partition zones that want to be divided into,
|
||||||
|
// which is equal to the number of CPU processors you want.
|
||||||
|
// Usually, 50~100 thousands structured cells per CPU-Core is suggested.
|
||||||
|
// 30~70 thousands unstructured cells per CPU-Core is suggested.
|
||||||
|
// original_grid_file: Original grid file that want to be divided(PHengLEI type, *.fts).
|
||||||
|
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||||
|
|
||||||
|
int pgridtype = 1;
|
||||||
|
int maxproc = 4;
|
||||||
|
|
||||||
|
string original_grid_file = "./grid/rae2822_vis2d.fts";
|
||||||
|
string partition_grid_file = "./grid/rae2822_vis2d__4.fts";
|
||||||
|
|
||||||
|
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
|
||||||
|
// 1 -- single level.
|
||||||
|
// 2 -- 2 level.
|
||||||
|
// N -- N level,..., et al.
|
||||||
|
int numberOfMultigrid = 1;
|
Binary file not shown.
|
@ -0,0 +1,13 @@
|
||||||
|
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
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,196 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 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;
|
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// axisup: Type of Cartisien coordinates system, used in grid conversion.
|
||||||
|
// 1 -- Y upward. (default)
|
||||||
|
// 2 -- Z upward.
|
||||||
|
// from_gtype: Type of grid data type in grid conversion process.
|
||||||
|
// -1 -- MULTI_TYPE.
|
||||||
|
// 1 -- PHengLEI, *.fts.
|
||||||
|
// 2 -- CGNS, *.cgns.
|
||||||
|
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||||
|
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||||
|
// 5 -- Fluent, *.cas/*.msh.
|
||||||
|
// 6 -- Ustar, mgrid.in.
|
||||||
|
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||||
|
// 8 -- GMSH, *.msh.
|
||||||
|
// 9 -- Gridgen type of structured grid, *.dat/*.grd.
|
||||||
|
int gridtype = 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";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
|
@ -0,0 +1,24 @@
|
||||||
|
// pgridtype: The grid type.
|
||||||
|
// 0 -- unstruct grid.
|
||||||
|
// 1 -- struct grid.
|
||||||
|
// maxproc: The number of partition zones that want to be divided into,
|
||||||
|
// which is equal to the number of CPU processors you want.
|
||||||
|
// Usually, 50~100 thousands structured cells per CPU-Core is suggested.
|
||||||
|
// 30~70 thousands unstructured cells per CPU-Core is suggested.
|
||||||
|
// original_grid_file: Original grid file that want to be divided(PHengLEI type, *.fts).
|
||||||
|
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||||
|
// numberOfMultifile: The number of partition grid files that want to be dumped out.
|
||||||
|
|
||||||
|
int pgridtype = 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;
|
Binary file not shown.
|
@ -0,0 +1,161 @@
|
||||||
|
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
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,195 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 5000;
|
||||||
|
|
||||||
|
int intervalStepFlow = 500;
|
||||||
|
int intervalStepPlot = 500;
|
||||||
|
int intervalStepForce = 50;
|
||||||
|
int intervalStepRes = 10;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Inflow Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// refMachNumber: Mach number.
|
||||||
|
// attackd: Angle of attack.
|
||||||
|
// angleSlide: Angle of sideslip.
|
||||||
|
// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
|
||||||
|
// inflowParaType: The type of inflow parameters.
|
||||||
|
// 0 -- the nondimensional conditions.
|
||||||
|
// 1 -- the flight conditions.
|
||||||
|
// 2 -- the experiment conditions.
|
||||||
|
// 3 -- the subsonic boundary conditions.
|
||||||
|
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||||
|
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||||
|
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||||
|
// height: Fly height, unit of km.
|
||||||
|
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||||
|
// 1 dm = 0.1 m.
|
||||||
|
// 1 cm = 0.01 m.
|
||||||
|
// 1 mm = 0.001m.
|
||||||
|
// 1 inch = 0.0254m.
|
||||||
|
// 1 foot = 12 inches = 0.3048m.
|
||||||
|
// 1 yard = 3 feet = 0.9144m.
|
||||||
|
// 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 attackd = 0.00;
|
||||||
|
double angleSlide = 0.00;
|
||||||
|
|
||||||
|
double wallTemperature = -1.0;
|
||||||
|
|
||||||
|
int inflowParaType = 0;
|
||||||
|
double refReNumber = 1.0e5;
|
||||||
|
double refDimensionalTemperature = 293;
|
||||||
|
|
||||||
|
//int inflowParaType = 1;
|
||||||
|
//double height = 0.001;
|
||||||
|
|
||||||
|
//int inflowParaType = 2;
|
||||||
|
//double refDimensionalTemperature = 6051.024; // The total temperature, T*(1+(refGama-1)*M*M/2).
|
||||||
|
//double refDimensionalPressure = 4.299696E09; // The total pressure, p*(T0/T)^(refGama/(refGama-1)).
|
||||||
|
|
||||||
|
double gridScaleFactor = 1.0;
|
||||||
|
|
||||||
|
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||||
|
double forceReferenceLength = 1.0; // unit of meter.
|
||||||
|
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||||
|
double TorqueRefX = 0.0; // unit of meter.
|
||||||
|
double TorqueRefY = 0.0; // unit of meter.
|
||||||
|
double TorqueRefZ = 0.0; // unit of meter.
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Physical models #
|
||||||
|
#########################################################################
|
||||||
|
// viscousType : Viscous model.
|
||||||
|
// 0 -- Euler.
|
||||||
|
// 1 -- Lamilar.
|
||||||
|
// 3 -- 1eq turbulent.
|
||||||
|
// 4 -- 2eq turbulent.
|
||||||
|
// viscousName: Laminar or tubulent model.
|
||||||
|
// -- "1eq-sa", when viscousType = 3.
|
||||||
|
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||||
|
// DESType: Type of DES.
|
||||||
|
// 0 -- RANS.(default)
|
||||||
|
// 1 -- DES.
|
||||||
|
// 2 -- DDES.
|
||||||
|
// 3 -- IDDES.
|
||||||
|
|
||||||
|
int viscousType = 0;
|
||||||
|
string viscousName = "Euler";
|
||||||
|
|
||||||
|
//int viscousType = 1;
|
||||||
|
//string viscousName = "laminar";
|
||||||
|
|
||||||
|
//int viscousType = 3;
|
||||||
|
//string viscousName = "1eq-sa";
|
||||||
|
|
||||||
|
//int viscousType = 4;
|
||||||
|
//string viscousName = "2eq-kw-menter-sst";
|
||||||
|
|
||||||
|
int DESType = 0;
|
||||||
|
|
||||||
|
int roeEntropyFixMethod = 6;
|
||||||
|
double roeEntropyScale = 1.0;
|
||||||
|
#########################################################################
|
||||||
|
# Spatial Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
#*******************************************************************
|
||||||
|
# Struct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||||
|
// Using this when solve structered grid or hybrid.
|
||||||
|
// -- "vanleer", "steger", "ausmpw", "ausmpw+".
|
||||||
|
// str_limiter_name: Limiter of struct grid.
|
||||||
|
// -- "minmod", "3rd_minmod_smooth".
|
||||||
|
|
||||||
|
string inviscidSchemeName = "roe";
|
||||||
|
string str_limiter_name = "minvan";
|
||||||
|
|
||||||
|
#*******************************************************************
|
||||||
|
# UnStruct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||||
|
// Using this when solve Unstructered grid or hybrid.
|
||||||
|
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||||
|
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||||
|
// uns_limiter_name: Limiter of Unstruct grid.
|
||||||
|
// -- "vencat", "barth".
|
||||||
|
// -- "1st", meaning accuracy of first-order.
|
||||||
|
// -- "nolim", no limiter.
|
||||||
|
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||||
|
// The smaller the value, the more robust it is.
|
||||||
|
|
||||||
|
string uns_scheme_name = "vanleer";
|
||||||
|
string uns_limiter_name = "vencat";
|
||||||
|
double venkatCoeff = 0.5;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Temporal Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
// iunsteady: Steady or unsteady.
|
||||||
|
// 0 -- steady.
|
||||||
|
// 1 -- unsteay.
|
||||||
|
// CFLEnd: The CFL number, [0.1, 100].
|
||||||
|
// The bigger the value, the convergence faster but lower robustness.
|
||||||
|
// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
|
||||||
|
|
||||||
|
int iunsteady = 0;
|
||||||
|
|
||||||
|
double CFLStart = 0.1;
|
||||||
|
|
||||||
|
double CFLEnd = 20.0;
|
||||||
|
|
||||||
|
int CFLVaryStep = 100;
|
||||||
|
|
||||||
|
double ktmax = 1.0e10;
|
||||||
|
|
||||||
|
double LUSGSTolerance = 0.1;
|
||||||
|
#########################################################################
|
||||||
|
# File In or Out #
|
||||||
|
#########################################################################
|
||||||
|
// gridfile: The partitioned Grid file path, using relative path,
|
||||||
|
// which is relative to the working directory.
|
||||||
|
// IMPORTANT WARNING: The file index should be ignored,
|
||||||
|
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||||
|
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||||
|
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||||
|
|
||||||
|
string gridfile = "./grid/cylinder.fts";
|
||||||
|
int plotFieldType = 0;
|
||||||
|
|
||||||
|
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
|
||||||
|
// visualVariables: Variable types dumped, listed as following:
|
||||||
|
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
|
||||||
|
// -- viscosityLaminar(7), viscosityTurbulent(8),
|
||||||
|
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
|
||||||
|
// -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
|
||||||
|
// -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21).
|
||||||
|
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
|
||||||
|
// Variables order must from small to big.
|
||||||
|
|
||||||
|
int nVisualVariables = 8;
|
||||||
|
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
|
||||||
|
|
||||||
|
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||||
|
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||||
|
// 1 -- limit for every variables, then get the min value.
|
||||||
|
// limitVector:
|
||||||
|
// 0 -- Each variable use the same limiter coefficient.
|
||||||
|
// 1 -- Each variable use the respective limiter coefficients.
|
||||||
|
// reconmeth:
|
||||||
|
// 0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
|
||||||
|
// 1 -- Q+, Q- use the min limiter coefficients of left and right cell.
|
||||||
|
|
||||||
|
int reconmeth = 0;
|
||||||
|
int limitVariables = 0;
|
||||||
|
int limitVector = 1;
|
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// axisup: Type of Cartisien coordinates system, used in grid conversion.
|
||||||
|
// 1 -- Y upward. (default)
|
||||||
|
// 2 -- Z upward.
|
||||||
|
// from_gtype: Type of grid data type in grid conversion process.
|
||||||
|
// -1 -- MULTI_TYPE.
|
||||||
|
// 1 -- PHengLEI, *.fts.
|
||||||
|
// 2 -- CGNS, *.cgns.
|
||||||
|
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||||
|
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||||
|
// 5 -- Fluent, *.cas/*.msh.
|
||||||
|
// 6 -- Ustar, mgrid.in.
|
||||||
|
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||||
|
// 8 -- GMSH, *.msh.
|
||||||
|
// 9 -- Gridgen type of structured grid, *.dat/*.grd.
|
||||||
|
int gridtype = 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/cylinder.grd";
|
||||||
|
string out_gfile = "./grid/cylinder.fts";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
//string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
Binary file not shown.
|
@ -0,0 +1,9 @@
|
||||||
|
1
|
||||||
|
1
|
||||||
|
121 61
|
||||||
|
A
|
||||||
|
4
|
||||||
|
1 121 61 61 5
|
||||||
|
121 121 61 1 6
|
||||||
|
1 121 1 1 2
|
||||||
|
1 1 61 1 6
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,203 @@
|
||||||
|
#########################################################################
|
||||||
|
# General Control Parameter #
|
||||||
|
#########################################################################
|
||||||
|
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||||
|
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||||
|
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||||
|
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||||
|
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||||
|
|
||||||
|
int maxSimuStep = 10000;
|
||||||
|
|
||||||
|
int intervalStepFlow = 500;
|
||||||
|
int intervalStepPlot = 500;
|
||||||
|
int intervalStepForce = 50;
|
||||||
|
int intervalStepRes = 50;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# 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.73;
|
||||||
|
double attackd = 2.31;
|
||||||
|
double angleSlide = 0.00;
|
||||||
|
|
||||||
|
int inflowParaType = 0;
|
||||||
|
double refReNumber = 6.5e6;
|
||||||
|
double refDimensionalTemperature = 288.15;
|
||||||
|
|
||||||
|
//int inflowParaType = 1;
|
||||||
|
//double height = 0.001;
|
||||||
|
|
||||||
|
//int inflowParaType = 2;
|
||||||
|
//double refDimensionalTemperature = 6051.024; // The total temperature, T*(1+(refGama-1)*M*M/2).
|
||||||
|
//double refDimensionalPressure = 4.299696E09; // The total pressure, p*(T0/T)^(refGama/(refGama-1)).
|
||||||
|
|
||||||
|
double gridScaleFactor = 1.0;
|
||||||
|
|
||||||
|
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||||
|
double forceReferenceLength = 1.0; // unit of meter.
|
||||||
|
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||||
|
double TorqueRefX = 0.0; // unit of meter.
|
||||||
|
double TorqueRefY = 0.0; // unit of meter.
|
||||||
|
double TorqueRefZ = 0.0; // unit of meter.
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Physical models #
|
||||||
|
#########################################################################
|
||||||
|
// viscousType : Viscous model.
|
||||||
|
// 0 -- Euler.
|
||||||
|
// 1 -- Lamilar.
|
||||||
|
// 3 -- 1eq turbulent.
|
||||||
|
// 4 -- 2eq turbulent.
|
||||||
|
// viscousName: Laminar or tubulent model.
|
||||||
|
// -- "1eq-sa", when viscousType = 3.
|
||||||
|
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||||
|
// DESType: Type of DES.
|
||||||
|
// 0 -- RANS.(default)
|
||||||
|
// 1 -- DES.
|
||||||
|
// 2 -- DDES.
|
||||||
|
// 3 -- IDDES.
|
||||||
|
|
||||||
|
//int viscousType = 0;
|
||||||
|
//string viscousName = "Euler";
|
||||||
|
|
||||||
|
//int viscousType = 1;
|
||||||
|
//string viscousName = "laminar";
|
||||||
|
|
||||||
|
int viscousType = 3;
|
||||||
|
string viscousName = "1eq-sa";
|
||||||
|
|
||||||
|
//int viscousType = 4;
|
||||||
|
//string viscousName = "2eq-kw-menter-sst";
|
||||||
|
|
||||||
|
int DESType = 0;
|
||||||
|
|
||||||
|
int roeEntropyFixMethod = 6;
|
||||||
|
double roeEntropyScale = 1.0;
|
||||||
|
#########################################################################
|
||||||
|
# Spatial Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
#*******************************************************************
|
||||||
|
# Struct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||||
|
// Using this when solve structered grid or hybrid.
|
||||||
|
// -- "roe", "vanleer", "ausm+up", "ausmpw".
|
||||||
|
// str_limiter_name: Limiter of struct grid.
|
||||||
|
// -- "3rdsmooth", "smooth".
|
||||||
|
// -- "nolim", no limiter.
|
||||||
|
|
||||||
|
string inviscidSchemeName = "roe";
|
||||||
|
string str_limiter_name = "smooth";
|
||||||
|
double MUSCLCoefXk = 0.333333;
|
||||||
|
#*******************************************************************
|
||||||
|
# UnStruct Solver *
|
||||||
|
#*******************************************************************
|
||||||
|
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||||
|
// Using this when solve Unstructered grid or hybrid.
|
||||||
|
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||||
|
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||||
|
// uns_limiter_name: Limiter of Unstruct grid.
|
||||||
|
// -- "vencat", "barth".
|
||||||
|
// -- "1st", meaning accuracy of first-order.
|
||||||
|
// -- "nolim", no limiter.
|
||||||
|
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||||
|
// The smaller the value, the more robust it is.
|
||||||
|
|
||||||
|
string uns_scheme_name = "roe";
|
||||||
|
string uns_limiter_name = "vencat";
|
||||||
|
double venkatCoeff = 5.0;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Temporal Discretisation #
|
||||||
|
#########################################################################
|
||||||
|
// iunsteady: Steady or unsteady.
|
||||||
|
// 0 -- steady.
|
||||||
|
// 1 -- unsteay.
|
||||||
|
// CFLEnd: The CFL number, [0.1, 100].
|
||||||
|
// The bigger the value, the convergence faster but lower robustness.
|
||||||
|
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
|
||||||
|
// 1 -- is recommended for structured solver.
|
||||||
|
// 1-3 -- is recommended for unstructured solver.
|
||||||
|
|
||||||
|
int iunsteady = 0;
|
||||||
|
|
||||||
|
double CFLEnd = 10.0;
|
||||||
|
|
||||||
|
int nLUSGSSweeps = 1;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# Multi-Grid parameters #
|
||||||
|
#########################################################################
|
||||||
|
// nMGLevel: The number of Multi-Grid level.
|
||||||
|
// = 1 -- single-level.
|
||||||
|
// > 1 -- multi-level.
|
||||||
|
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||||
|
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||||
|
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||||
|
|
||||||
|
int nMGLevel = 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/rae2822_vis2d__4.fts";
|
||||||
|
int plotFieldType = 0;
|
||||||
|
|
||||||
|
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
|
||||||
|
// visualVariables: Variable types dumped, listed as following:
|
||||||
|
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
|
||||||
|
// -- viscosityLaminar(7), viscosityTurbulent(8),
|
||||||
|
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
|
||||||
|
// -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
|
||||||
|
// -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
|
||||||
|
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
|
||||||
|
// Variables order must from small to big.
|
||||||
|
|
||||||
|
int nVisualVariables = 8;
|
||||||
|
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
|
||||||
|
|
||||||
|
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||||
|
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||||
|
// 1 -- limit for every variables, then get the min value.
|
||||||
|
// limitVector:
|
||||||
|
// 0 -- Each variable use the same limiter coefficient.
|
||||||
|
// 1 -- Each variable use the respective limiter coefficients.
|
||||||
|
// reconmeth:
|
||||||
|
// 0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
|
||||||
|
// 1 -- Q+, Q- use the min limiter coefficients of left and right cell.
|
||||||
|
|
||||||
|
int reconmeth = 1;
|
||||||
|
int limitVariables = 0;
|
||||||
|
int limitVector = 0;
|
|
@ -0,0 +1,30 @@
|
||||||
|
#########################################################################
|
||||||
|
# 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";
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 2;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
//string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
|
@ -0,0 +1,21 @@
|
||||||
|
// pgridtype: The grid type.
|
||||||
|
// 0 -- unstruct grid.
|
||||||
|
// 1 -- struct grid.
|
||||||
|
// maxproc: The number of partition zones that want to be divided into,
|
||||||
|
// which is equal to the number of CPU processors you want.
|
||||||
|
// Usually, 50~100 thousands structured cells per CPU-Core is suggested.
|
||||||
|
// 30~70 thousands unstructured cells per CPU-Core is suggested.
|
||||||
|
// original_grid_file: Original grid file that want to be divided(PHengLEI type, *.fts).
|
||||||
|
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||||
|
|
||||||
|
int pgridtype = 1;
|
||||||
|
int maxproc = 4;
|
||||||
|
|
||||||
|
string original_grid_file = "./grid/rae2822_vis2d.fts";
|
||||||
|
string partition_grid_file = "./grid/rae2822_vis2d__4.fts";
|
||||||
|
|
||||||
|
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
|
||||||
|
// 1 -- single level.
|
||||||
|
// 2 -- 2 level.
|
||||||
|
// N -- N level,..., et al.
|
||||||
|
int numberOfMultigrid = 1;
|
Binary file not shown.
|
@ -0,0 +1,13 @@
|
||||||
|
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
|
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// axisup: Type of Cartisien coordinates system, used in grid conversion.
|
||||||
|
// 1 -- Y upward. (default)
|
||||||
|
// 2 -- Z upward.
|
||||||
|
// from_gtype: Type of grid data type in grid conversion process.
|
||||||
|
// -1 -- MULTI_TYPE.
|
||||||
|
// 1 -- PHengLEI, *.fts.
|
||||||
|
// 2 -- CGNS, *.cgns.
|
||||||
|
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||||
|
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||||
|
// 5 -- Fluent, *.cas/*.msh.
|
||||||
|
// 6 -- Ustar, mgrid.in.
|
||||||
|
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||||
|
// 8 -- GMSH, *.msh.
|
||||||
|
// 9 -- Gridgen type of structured grid, *.dat/*.grd.
|
||||||
|
int gridtype = 0;
|
||||||
|
int axisup = 1;
|
||||||
|
int from_gtype = 2;
|
||||||
|
|
||||||
|
#########################################################################
|
||||||
|
# File path #
|
||||||
|
#########################################################################
|
||||||
|
// from_gfile: path of original data file for unstructure grid convert from.
|
||||||
|
// out_gfile: path of target file for grid convert to, *.fts type of file usually.
|
||||||
|
string from_gfile = "./grid/M6_unstr.cgns";
|
||||||
|
string out_gfile = "./grid/M6_unstr.fts";
|
|
@ -0,0 +1,31 @@
|
||||||
|
#########################################################################
|
||||||
|
# Grid data type #
|
||||||
|
#########################################################################
|
||||||
|
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||||
|
// 0 -- Unstructured grid.
|
||||||
|
// 1 -- Structured grid.
|
||||||
|
// gridobj: Task type of grid treatment.
|
||||||
|
// 0 -- Grid generation of typical case, such as cylinder, flat plate, etc.
|
||||||
|
// 1 -- Grid conversion, from other grid data to PHenglEI, such as Fluent, CGNS.
|
||||||
|
// 2 -- Grid refinement.
|
||||||
|
// 3 -- Grid merging, merge two blocks into one block.
|
||||||
|
// 4 -- Grid deformation, achieve unstructured grid deformation.
|
||||||
|
// 5 -- Grid repairing, repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 6 -- Grid mirroring, mirror a symmetry grid to whole grid.
|
||||||
|
// 7 -- Grid type change, convert structured grid to unstructured grid.
|
||||||
|
// anisoRefine: If refine grid by anisoRefine type.
|
||||||
|
// 0 -- Not. (default)
|
||||||
|
// 1 -- Yes.
|
||||||
|
int gridtype = 0;
|
||||||
|
int gridobj = 2;
|
||||||
|
int anisoRefine = 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/M6_unstr__4.fts";
|
||||||
|
string out_gfile = "./grid/M6_unstr_refined__4.fts";
|
||||||
|
|
||||||
|
|
|
@ -0,0 +1,52 @@
|
||||||
|
string title = "PHengLEI Main Parameter Control File";
|
||||||
|
|
||||||
|
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||||
|
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
// ndim: Dimensional of the grid, 2 or 3.
|
||||||
|
// nparafile: the number of parameter files.
|
||||||
|
// nsimutask: simulation task type.
|
||||||
|
// 0 -- CFD Solver of NS or Turbulation.
|
||||||
|
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||||
|
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||||
|
// Grid reconstruction: such as grid adaptation.
|
||||||
|
// Grid merging: merge two blocks into one block.
|
||||||
|
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||||
|
// 2 -- Wall distance computation for turb-solver.
|
||||||
|
// 3 -- Grid partition.
|
||||||
|
int ndim = 3;
|
||||||
|
int nparafile = 1;
|
||||||
|
|
||||||
|
int nsimutask = 0;
|
||||||
|
//string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||||
|
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||||
|
|
||||||
|
int nsimutask = 1;
|
||||||
|
string parafilename = "./bin/grid_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 2;
|
||||||
|
//string parafilename = "./bin/cfd_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 3;
|
||||||
|
//string parafilename = "./bin/partition.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 1;
|
||||||
|
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 14;
|
||||||
|
//string parafilename = "./bin/integrative_solver.hypara";
|
||||||
|
|
||||||
|
//int nsimutask = 99;
|
||||||
|
//string parafilename = "./bin/post_processing.hypara";
|
||||||
|
|
||||||
|
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||||
|
int numberOfGridProcessor = 0;
|
||||||
|
// ATP read
|
||||||
|
//@string parafilename1 = ""
|
||||||
|
//@string parafilename2 = "";
|
|
@ -0,0 +1,21 @@
|
||||||
|
// pgridtype: The grid type.
|
||||||
|
// 0 -- unstruct grid.
|
||||||
|
// 1 -- struct grid.
|
||||||
|
// maxproc: The number of partition zones that want to be divided into,
|
||||||
|
// which is equal to the number of CPU processors you want.
|
||||||
|
// Usually, 50~100 thousands structured cells per CPU-Core is suggested.
|
||||||
|
// 30~70 thousands unstructured cells per CPU-Core is suggested.
|
||||||
|
// original_grid_file: Original grid file that want to be divided(PHengLEI type, *.fts).
|
||||||
|
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||||
|
|
||||||
|
int pgridtype = 0;
|
||||||
|
int maxproc = 4;
|
||||||
|
|
||||||
|
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.
|
||||||
|
// 2 -- 2 level.
|
||||||
|
// N -- N level,..., et al.
|
||||||
|
int numberOfMultigrid = 1;
|
Binary file not shown.
Binary file not shown.
29
README.md
29
README.md
|
@ -43,6 +43,12 @@ A 结构
|
||||||
| 三维结构激波干扰超声速层流平板 | A21_ThreeD_ShockwaveInteractionPlate_Ma2_Laminar_UDB_Struct_1CPU |
|
| 三维结构激波干扰超声速层流平板 | A21_ThreeD_ShockwaveInteractionPlate_Ma2_Laminar_UDB_Struct_1CPU |
|
||||||
| 三维结构弹体喷流(自定义边界) | A22_ThreeD_Jet_Ma3d33_SST_Struct_8CPU |
|
| 三维结构弹体喷流(自定义边界) | A22_ThreeD_Jet_Ma3d33_SST_Struct_8CPU |
|
||||||
| 二维结构S-K低速平板转捩模拟 | A23_TwoD_Plate_S-KSR_SST_Struct_4CPU |
|
| 二维结构S-K低速平板转捩模拟 | A23_TwoD_Plate_S-KSR_SST_Struct_4CPU |
|
||||||
|
| 二维结构30p30n翼型低速绕流(Matrix LUSGS) | A24_TwoD_30p30n_SA_MatrixLUSGS_Struct_4CPU |
|
||||||
|
| 二维结构低速湍流平板(Matrix LUSGS) | A25_TwoD_Plate_SST_LowMach_MatrixLUSGS_Struct_1CPU |
|
||||||
|
| 二维结构Rae2822翼型跨声速绕流(Matrix LUSGS) | A26_TwoD_Rae2822_SA_MatrixLUSGS_Struct_4CPU |
|
||||||
|
| 二维结构30p30n翼型低速绕流(熵修正) | A27_TwoD_30p30n_SA_Entropyfix6_Struct_4CPU |
|
||||||
|
| 二维结构圆柱超声速绕流(熵修正) | A28_TwoD_Cylinder_Inv_Ma3d0_Entropyfix6_Struct_1CPU |
|
||||||
|
| 二维结构Rae2822翼型跨声速绕流(熵修正) | A29_TwoD_Rae2822_SA_Entropyfix6_Struct_4CPU |
|
||||||
|
|
||||||
B 非结构
|
B 非结构
|
||||||
|
|
||||||
|
@ -81,6 +87,9 @@ D 高阶
|
||||||
| :------------------------- | :--------------------------------- |
|
| :------------------------- | :--------------------------------- |
|
||||||
| 二维结构高阶圆柱 | D01_TwoD_Cylinder_Laminar_HighOrder_Struct_1CPU |
|
| 二维结构高阶圆柱 | D01_TwoD_Cylinder_Laminar_HighOrder_Struct_1CPU |
|
||||||
| 三维结构高阶双椭球 | D02_ThreeD_DoubleEllipse_Laminar_HighOrder_Struct_48CPU |
|
| 三维结构高阶双椭球 | D02_ThreeD_DoubleEllipse_Laminar_HighOrder_Struct_48CPU |
|
||||||
|
| 二维WENN结构高阶30P30N翼型 | D03_TwoD_30P30N_SST_HighOrderWENN_Struct_6CPU |
|
||||||
|
| 三维WENN结构高阶M6机翼 | D04_ThreeD_ONERA_M6_SST_HighOrderWENN_Struct_8CPU |
|
||||||
|
| 三维WENN结构高阶双椭球 | D05_ThreeD_DoubleEllipse_Laminar_HighOrderWENN_Struct_48CPU |
|
||||||
|
|
||||||
E LES
|
E LES
|
||||||
|
|
||||||
|
@ -96,14 +105,15 @@ F 重叠网格、动网格
|
||||||
| :------------------------- | :--------------------------------- |
|
| :------------------------- | :--------------------------------- |
|
||||||
| 三维结构30p30n重叠网格装配 | F01_ThreeD_30p30n_OversetGrid_SA_Struct_4CPU |
|
| 三维结构30p30n重叠网格装配 | F01_ThreeD_30p30n_OversetGrid_SA_Struct_4CPU |
|
||||||
| 三维结构多弹体重叠网格装配 | F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU |
|
| 三维结构多弹体重叠网格装配 | F02_ThreeD_duodan_OversetGrid_SA_Struct_8CPU |
|
||||||
| 二维非结构双0012重叠网格装配 | F03_TwoD_Dual0012_OversetGrid_Invis_Unstruct_4CPU |
|
| 二维非结构双0012重叠网格装配 | F03_TwoD_Dual0012_OversetConfig_Unstruct_4CPU |
|
||||||
| 二维非结构双0012重叠网格计算 | F04_TwoD_Dual0012_OversetGrid_InvisCal_Unstruct_4CPU |
|
| 二维非结构双0012重叠网格计算 | F04_TwoD_Dual0012_OversetGrid_InvisCal_Unstruct_4CPU |
|
||||||
| 三维非结构AEDC重叠网格装配 | F05_ThreeD_AEDC_OversetGrid_SST_Unstruct_8CPU |
|
| 三维非结构AEDC吊舱分离重叠网格装配 | F05_ThreeD_AEDC_OversetConfig_Unstruct_8CPU |
|
||||||
| 三维非结构AEDC重叠网格计算 | F06_ThreeD_AEDC_OversetGrid_SSTCal_Unstruct_8CPU |
|
| 三维非结构AEDC吊舱分离非定常计算 | F06_ThreeD_AEDC_Separation_OversetGrid_Unsteady_InvisCal_Unstruct_8CPU |
|
||||||
| 二维结构NACA0012翼型俯仰振荡 | F07_TwoD_NACA0012_PitchingMovement_SA_Struct_4CPU |
|
| 二维结构NACA0012翼型俯仰振荡 | F07_TwoD_NACA0012_PitchingMovement_SA_Struct_4CPU |
|
||||||
| 三维结构Finner弹体俯仰振荡 | F08_ThreeD_Finner_Laminar_Ma2d5_Struct_16CPU |
|
| 三维结构Finner弹体俯仰振荡 | F08_ThreeD_Finner_Laminar_Ma2d5_Struct_16CPU |
|
||||||
| 鸟类扑翼网格变形(弹簧法) | F09_ThreeD_Brid_GridDeformationSPRING_Unstruct_1CPU |
|
| 鸟类扑翼网格变形(弹簧法) | F09_ThreeD_Brid_GridDeformationSPRING_Unstruct_1CPU |
|
||||||
| 鸟类扑翼网格变形(RBF方法) | F10_ThreeD_Brid_GridDeformationRBF_Unstruct_1CPU |
|
| 鸟类扑翼网格变形(RBF方法) | F10_ThreeD_Brid_GridDeformationRBF_Unstruct_1CPU |
|
||||||
|
| 非结构三维M6机翼并行加密 | F11_ThreeD_M6_GridParallelRefine_Unstruct_4CPU |
|
||||||
|
|
||||||
G 非平衡
|
G 非平衡
|
||||||
|
|
||||||
|
@ -113,6 +123,19 @@ G 非平衡
|
||||||
| 三维结构钝锥化学非平衡 | G02_ThreeD_Electre_Laminar_OneTemperMode_NCW_Struct |
|
| 三维结构钝锥化学非平衡 | G02_ThreeD_Electre_Laminar_OneTemperMode_NCW_Struct |
|
||||||
| 三维结构类返回舱化学非平衡 | G03_ |
|
| 三维结构类返回舱化学非平衡 | G03_ |
|
||||||
| 三维结构类航天飞机化学非平衡 | G04_ |
|
| 三维结构类航天飞机化学非平衡 | G04_ |
|
||||||
|
| 二维非结构高速圆柱两气体组分预混 | G05_TwoD_Cylinder_N2-O2_Premixing-Unstruct_1CPU |
|
||||||
|
| 三维非结构高速球头两气体组分预混 | G06_ThreeD_Sphere_N2-O2_Premixing-Unstruct_1CPU |
|
||||||
|
| 二维非结构管道喷流两气体组分预混 | G07_TwoD_Pipe_Air-CH4_Mixing-Unstruct_1CPU |
|
||||||
|
|
||||||
|
H 不可压
|
||||||
|
|
||||||
|
| 中文 | 英文 |
|
||||||
|
| :----------------------------------- | :-------------------------------------------------------- |
|
||||||
|
| 二维非结构不可压层流后台阶流动 | H01_PBSolver_TwoD_Backstep_Laminar_Re95_Unstruct_1CPU |
|
||||||
|
| 二维非结构不可压层流方腔自然对流 | H02_PBSolver_TwoD_Cavity_Boussinesq_Laminar_Unstruct_1CPU |
|
||||||
|
| 二维非结构不可压层流方腔流动 | H03_PBSolver_TwoD_Cavity_Laminar_Re1000_Unstruct_1CPU |
|
||||||
|
| 三维非结构不可压圆管湍流压降模拟 | H04_PBSolver_ThreeD_Pipe_Turb_Unstruct_1CPU |
|
||||||
|
| 三维非结构不可压AFF-1 半艇外流场模拟 | H05_PBSolver_ThreeD_Suboff_Laminar_Re225_1CPU |
|
||||||
|
|
||||||
X 用户分享算例
|
X 用户分享算例
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue