新功能算例添加，部分已有功能参数及文档修改
This commit is contained in:
hechao
parent
17098ab5f3
commit
e75f62e76c
|
|
@ -0,0 +1,205 @@
|
|||
#########################################################################
|
||||
# General Control Parameter #
|
||||
#########################################################################
|
||||
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||
|
||||
int maxSimuStep = 10000;
|
||||
|
||||
int intervalStepFlow = 1000;
|
||||
int intervalStepPlot = 1000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 0.8395;
|
||||
double attackd = 3.06;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 1.801e7;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//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 = 0.53703; // unit of meter.
|
||||
double forceReferenceArea = 0.5204; // unit of meter^2.
|
||||
double TorqueRefX = 0.0; // unit of meter.
|
||||
double TorqueRefY = 0.0; // unit of meter.
|
||||
double TorqueRefZ = 0.0; // unit of meter.
|
||||
|
||||
#########################################################################
|
||||
# Physical models #
|
||||
#########################################################################
|
||||
// viscousType : Viscous model.
|
||||
// 0 -- Euler.
|
||||
// 1 -- Lamilar.
|
||||
// 3 -- 1eq turbulent.
|
||||
// 4 -- 2eq turbulent.
|
||||
// viscousName: Laminar or tubulent model.
|
||||
// -- "1eq-sa", when viscousType = 3.
|
||||
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||
// DESType: Type of DES.
|
||||
// 0 -- RANS.(default)
|
||||
// 1 -- DES.
|
||||
// 2 -- DDES.
|
||||
// 3 -- IDDES.
|
||||
|
||||
//int viscousType = 0;
|
||||
//string viscousName = "Euler";
|
||||
|
||||
//int viscousType = 1;
|
||||
//string viscousName = "laminar";
|
||||
|
||||
int viscousType = 3;
|
||||
string viscousName = "1eq-sa";
|
||||
|
||||
//int viscousType = 4;
|
||||
//string viscousName = "2eq-kw-menter-sst";
|
||||
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 2;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||
// Using this when solve structered grid or hybrid.
|
||||
// -- "roe", "vanleer", "ausm+up", "ausmpw".
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string inviscidSchemeName = "roe";
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||
// Using this when solve Unstructered grid or hybrid.
|
||||
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat", "barth".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_scheme_name = "roe";
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 5.0;
|
||||
|
||||
#########################################################################
|
||||
# Temporal Discretisation #
|
||||
#########################################################################
|
||||
// iunsteady: Steady or unsteady.
|
||||
// 0 -- steady.
|
||||
// 1 -- unsteay.
|
||||
// CFLEnd: The CFL number, [0.1, 100].
|
||||
// The bigger the value, the convergence faster but lower robustness.
|
||||
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
|
||||
// 1 -- is recommended for structured solver.
|
||||
// 1-3 -- is recommended for unstructured solver.
|
||||
|
||||
int iunsteady = 0;
|
||||
|
||||
double CFLEnd = 10.0;
|
||||
|
||||
int nLUSGSSweeps = 1;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 2;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/m6_str__4.fts";
|
||||
int plotFieldType = 0;
|
||||
|
||||
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
|
||||
// visualVariables: Variable types dumped, listed as following:
|
||||
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
|
||||
// -- viscosityLaminar(7), viscosityTurbulent(8),
|
||||
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
|
||||
// -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
|
||||
// -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
|
||||
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
|
||||
// Variables order must from small to big.
|
||||
|
||||
int nVisualVariables = 8;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
|
||||
|
||||
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||
// 1 -- limit for every variables, then get the min value.
|
||||
// limitVector:
|
||||
// 0 -- Each variable use the same limiter coefficient.
|
||||
// 1 -- Each variable use the respective limiter coefficients.
|
||||
// reconmeth:
|
||||
// 0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
|
||||
// 1 -- Q+, Q- use the min limiter coefficients of left and right cell.
|
||||
|
||||
int reconmeth = 1;
|
||||
int limitVariables = 0;
|
||||
int limitVector = 0;
|
||||
|
||||
double MUSCLCoefXk = 0.333333;
|
||||
|
|
@ -0,0 +1,35 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 1;
|
||||
int nAxisRotateTimes = 1;
|
||||
int axisRotateOrder[] = [1];
|
||||
double axisRotateAngles[]= [90.0];
|
||||
int from_gtype = 2;
|
||||
|
||||
#########################################################################
|
||||
# File path #
|
||||
#########################################################################
|
||||
// from_gfile: path of original data file for unstructure grid convert from.
|
||||
// out_gfile: path of target file for grid convert to, *.fts type of file usually.
|
||||
string from_gfile = "./grid/m6_str.cgns";
|
||||
string out_gfile = "./grid/m6_str.fts";
|
||||
|
|
@ -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 = 1;
|
||||
int maxproc = 4;
|
||||
|
||||
string original_grid_file = "./grid/m6_str.fts";
|
||||
string partition_grid_file = "./grid/m6_str__4.fts";
|
||||
|
||||
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
|
||||
// 1 -- single level.
|
||||
// 2 -- 2 level.
|
||||
// N -- N level,..., et al.
|
||||
int numberOfMultigrid = 2;
|
||||
|
|
@ -44,8 +44,8 @@ double attackd = 3.06;
|
|||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 1.171e7;
|
||||
double refDimensionalTemperature = 288;
|
||||
double refReNumber = 1.801e7;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//int inflowParaType = 1;
|
||||
//double height = 0.001;
|
||||
|
|
@ -57,8 +57,8 @@ double refDimensionalTemperature = 288;
|
|||
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 forceReferenceLength = 0.53703; // unit of meter.
|
||||
double forceReferenceArea = 0.5204; // 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.
|
||||
|
|
|
|||
|
|
@ -0,0 +1,56 @@
|
|||
# nBoundaryConditions: Number of global boundary conditions.
|
||||
# bcName : Boundary condition name.
|
||||
# bcType(in PHengLEI): Boundary condition type.
|
||||
|
||||
# How to set boundary condition, for example:
|
||||
# string bcName = "Wall";
|
||||
# {
|
||||
# int bcType = 2;
|
||||
# int viscousType = 1;
|
||||
# double wallTemperature = -1.0;
|
||||
# double uWall = 0.0;
|
||||
# double vWall = 0.0;
|
||||
# double wWall = 0.0;
|
||||
# }
|
||||
# string bcName = "Inflow";
|
||||
# {
|
||||
# int bcType = 5;
|
||||
# int inflowParaType = 0;
|
||||
# double refMachNumber = 0.73;
|
||||
# double attackd = 2.79;
|
||||
# double angleSlide = 0.0;
|
||||
# double refReNumber = 6.5e6;
|
||||
# double refDimensionalTemperature = 288.15;
|
||||
# }
|
||||
|
||||
# For more information, see examples/bin/boundary_condition.hypara file!!!
|
||||
|
||||
int nBoundaryConditions = 3;
|
||||
string bcName = "wall";
|
||||
{
|
||||
string bodyName = "body";
|
||||
int bcType = 2;
|
||||
}
|
||||
string bcName = "farfield";
|
||||
{
|
||||
int bcType = 4;
|
||||
}
|
||||
string bcName = "outflow";
|
||||
{
|
||||
int bcType = 6;
|
||||
}
|
||||
|
||||
# 'bcType' is defined as following:
|
||||
# -2: WAKE
|
||||
# -1: INTERFACE
|
||||
# 0 : NO_BOUNDARY_CONDITION
|
||||
# 1 : EXTRAPOLATION
|
||||
# 2 : SOLID_SURFACE
|
||||
# 3 : SYMMETRY
|
||||
# 4 : FARFIELD
|
||||
# 5 : INFLOW
|
||||
# 6 : OUTFLOW
|
||||
# 52: PRESSURE_INLET
|
||||
# 62: PRESSURE_OUTLET
|
||||
# 61: OUTFLOW_CONFINED
|
||||
# 7 : POLE
|
||||
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,193 @@
|
|||
#########################################################################
|
||||
# General Control Parameter #
|
||||
#########################################################################
|
||||
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||
|
||||
int maxSimuStep = 20000;
|
||||
|
||||
int intervalStepFlow = 2000;
|
||||
int intervalStepPlot = 2000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 2.85;
|
||||
double attackd = 0.00;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
double wallTemperature = 276;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 7.391e7;
|
||||
double refDimensionalTemperature = 100.3;
|
||||
|
||||
//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 wallFunctionType = 1;
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 2;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*************************************************
|
||||
******************
|
||||
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||
// Using this when solve structered grid or hybrid.
|
||||
// -- "vanleer", "steger", "ausmpw".
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "minmod", "3rd_minmod_smooth".
|
||||
|
||||
string inviscidSchemeName = "vanleer";
|
||||
string str_limiter_name = "minmod";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||
// Using this when solve Unstructered grid or hybrid.
|
||||
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat", "barth".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_scheme_name = "vanleer";
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 0.5;
|
||||
|
||||
#########################################################################
|
||||
# Temporal Discretisation #
|
||||
#########################################################################
|
||||
// iunsteady: Steady or unsteady.
|
||||
// 0 -- steady.
|
||||
// 1 -- unsteay.
|
||||
// CFLEnd: The CFL number, [0.1, 100].
|
||||
// The bigger the value, the convergence faster but lower robustness.
|
||||
// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
|
||||
|
||||
int iunsteady = 0;
|
||||
|
||||
double CFLEnd = 30.0;
|
||||
|
||||
#########################################################################
|
||||
# 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/str_1__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];
|
||||
|
||||
int nVisualWallVariables = 6;
|
||||
int visualWallVariables[] = [0, 1, 2, 3, 4, 5];
|
||||
|
||||
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||
// 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,35 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 1;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[]= [0.0, 0.0, 0.0];
|
||||
int from_gtype = 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/str_1.cgns";
|
||||
string out_gfile = "./grid/str_1.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/str_1.fts";
|
||||
string partition_grid_file = "./grid/str_1__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;
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
|
|
@ -0,0 +1,50 @@
|
|||
# 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 = "blade";
|
||||
{
|
||||
string bodyName = "body";
|
||||
int bcType = 2;
|
||||
}
|
||||
string bcName = "inlet";
|
||||
{
|
||||
int bcType = 52;
|
||||
double totalPressure = 107882.3;
|
||||
double totalTemperature = 293;
|
||||
int directionMethod = 0;
|
||||
double direction_inlet[] = 1, 0, 0;
|
||||
}
|
||||
string bcName = "outlet";
|
||||
{
|
||||
int bcType = 62;
|
||||
double staticPressure = 101325.0;
|
||||
}
|
||||
|
||||
# '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,202 @@
|
|||
#########################################################################
|
||||
# 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 = 100;
|
||||
int intervalStepPlot = 100;
|
||||
int intervalStepForce = 20;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
|
||||
// 0 -- no precondition process. (default, mach > 0.3)
|
||||
// 1 -- carry out precondition process. (mach number <= 0.3)
|
||||
int ifLowSpeedPrecon = 0;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 0.2;
|
||||
double attackd = 0.00;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 5000000.0;
|
||||
double refDimensionalTemperature = 300.00;
|
||||
//int inflowParaType = 1;
|
||||
//double height = 0.001;
|
||||
|
||||
//int inflowParaType = 2;
|
||||
//double refDimensionalTemperature = 302.4; // The total temperature, T*(1+(refGama-1)*M*M/2).
|
||||
//double refDimensionalPressure = 118309.784; // The total pressure, p*(T0/T)^(refGama/(refGama-1)).
|
||||
|
||||
double gridScaleFactor = 1.0;
|
||||
|
||||
double 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.
|
||||
|
||||
int periodicType = 1;
|
||||
double translationLength[] = [0.0,0.1,0.0];
|
||||
double rotationAngle = 0.0;
|
||||
|
||||
#########################################################################
|
||||
# Physical models #
|
||||
#########################################################################
|
||||
// viscousType : Viscous model.
|
||||
// 0 -- Euler.
|
||||
// 1 -- Lamilar.
|
||||
// 3 -- 1eq turbulent.
|
||||
// 4 -- 2eq turbulent.
|
||||
// viscousName: Laminar or tubulent model.
|
||||
// -- "1eq-sa", when viscousType = 3.
|
||||
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||
// DESType: Type of DES.
|
||||
// 0 -- RANS.(default)
|
||||
// 1 -- DES.
|
||||
// 2 -- DDES.
|
||||
// 3 -- IDDES.
|
||||
|
||||
//int viscousType = 0;
|
||||
//string viscousName = "Euler";
|
||||
|
||||
//int viscousType = 1;
|
||||
//string viscousName = "laminar";
|
||||
|
||||
int viscousType = 3;
|
||||
string viscousName = "1eq-sa";
|
||||
|
||||
//int viscousType = 4;
|
||||
//string viscousName = "2eq-kw-menter-sst";
|
||||
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 3;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_limiter_name = "nolim";
|
||||
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 = 5;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 1;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/Durham-2d__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 = 10;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15,39];
|
||||
|
||||
// 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,39 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 0;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[] = [0.0, 0.0, 0.0];
|
||||
int from_gtype = 2;
|
||||
|
||||
int periodicType = 1;
|
||||
double translationLength[] = [0.0,0.1,0.0];
|
||||
double rotationAngle = 0.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/Durham-2d.cgns";
|
||||
string out_gfile = "./grid/Durham-2d.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).
|
||||
|
||||
int pgridtype = 0;
|
||||
int maxproc = 4;
|
||||
|
||||
string original_grid_file = "./grid/Durham-2d.fts";
|
||||
string partition_grid_file = "./grid/Durham-2d__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;
|
||||
int periodicType = 1;
|
||||
double translationLength[] = [0.0,0.1,0.0];
|
||||
double rotationAngle = 0.0;
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,205 @@
|
|||
#########################################################################
|
||||
# General Control Parameter #
|
||||
#########################################################################
|
||||
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||
|
||||
int maxSimuStep = 10000;
|
||||
|
||||
int intervalStepFlow = 1000;
|
||||
int intervalStepPlot = 1000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 0.8395;
|
||||
double attackd = 3.06;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 1.801e7;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//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 = 0.53703; // unit of meter.
|
||||
double forceReferenceArea = 0.5204; // unit of meter^2.
|
||||
double TorqueRefX = 0.0; // unit of meter.
|
||||
double TorqueRefY = 0.0; // unit of meter.
|
||||
double TorqueRefZ = 0.0; // unit of meter.
|
||||
|
||||
#########################################################################
|
||||
# Physical models #
|
||||
#########################################################################
|
||||
// viscousType : Viscous model.
|
||||
// 0 -- Euler.
|
||||
// 1 -- Lamilar.
|
||||
// 3 -- 1eq turbulent.
|
||||
// 4 -- 2eq turbulent.
|
||||
// viscousName: Laminar or tubulent model.
|
||||
// -- "1eq-sa", when viscousType = 3.
|
||||
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||
// DESType: Type of DES.
|
||||
// 0 -- RANS.(default)
|
||||
// 1 -- DES.
|
||||
// 2 -- DDES.
|
||||
// 3 -- IDDES.
|
||||
|
||||
//int viscousType = 0;
|
||||
//string viscousName = "Euler";
|
||||
|
||||
//int viscousType = 1;
|
||||
//string viscousName = "laminar";
|
||||
|
||||
int viscousType = 3;
|
||||
string viscousName = "1eq-sa";
|
||||
|
||||
//int viscousType = 4;
|
||||
//string viscousName = "2eq-kw-menter-sst";
|
||||
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 2;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||
// Using this when solve structered grid or hybrid.
|
||||
// -- "roe", "vanleer", "ausm+up", "ausmpw".
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string inviscidSchemeName = "roe";
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||
// Using this when solve Unstructered grid or hybrid.
|
||||
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat", "barth".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_scheme_name = "roe";
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 5.0;
|
||||
|
||||
#########################################################################
|
||||
# Temporal Discretisation #
|
||||
#########################################################################
|
||||
// iunsteady: Steady or unsteady.
|
||||
// 0 -- steady.
|
||||
// 1 -- unsteay.
|
||||
// CFLEnd: The CFL number, [0.1, 100].
|
||||
// The bigger the value, the convergence faster but lower robustness.
|
||||
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
|
||||
// 1 -- is recommended for structured solver.
|
||||
// 1-3 -- is recommended for unstructured solver.
|
||||
|
||||
int iunsteady = 0;
|
||||
|
||||
double CFLEnd = 10.0;
|
||||
|
||||
int nLUSGSSweeps = 1;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 2;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/m6_unstr__4.fts";
|
||||
int plotFieldType = 0;
|
||||
|
||||
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||
// nVisualVariables: Number of variables want to be dumped for tecplot visualization.
|
||||
// visualVariables: Variable types dumped, listed as following:
|
||||
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
|
||||
// -- viscosityLaminar(7), viscosityTurbulent(8),
|
||||
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
|
||||
// -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
|
||||
// -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), iblank(81).
|
||||
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
|
||||
// Variables order must from small to big.
|
||||
|
||||
int nVisualVariables = 8;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
|
||||
|
||||
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||
// 1 -- limit for every variables, then get the min value.
|
||||
// limitVector:
|
||||
// 0 -- Each variable use the same limiter coefficient.
|
||||
// 1 -- Each variable use the respective limiter coefficients.
|
||||
// reconmeth:
|
||||
// 0 -- When reconstruct face value, Q+, Q- use respective limiter coefficients.
|
||||
// 1 -- Q+, Q- use the min limiter coefficients of left and right cell.
|
||||
|
||||
int reconmeth = 1;
|
||||
int limitVariables = 0;
|
||||
int limitVector = 0;
|
||||
|
||||
double MUSCLCoefXk = 0.333333;
|
||||
|
|
@ -0,0 +1,35 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 0;
|
||||
int nAxisRotateTimes = 1;
|
||||
int axisRotateOrder[] = [1];
|
||||
double axisRotateAngles[]= [90.0];
|
||||
int from_gtype = 2;
|
||||
|
||||
#########################################################################
|
||||
# File path #
|
||||
#########################################################################
|
||||
// from_gfile: path of original data file for unstructure grid convert from.
|
||||
// out_gfile: path of target file for grid convert to, *.fts type of file usually.
|
||||
string from_gfile = "./grid/m6_unstr.cgns";
|
||||
string out_gfile = "./grid/m6_unstr.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;
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
|
|
@ -0,0 +1,61 @@
|
|||
# nBoundaryConditions: Number of global boundary conditions.
|
||||
# bcName : Boundary condition name.
|
||||
# bcType(in PHengLEI): Boundary condition type.
|
||||
|
||||
# How to set boundary condition, for example:
|
||||
# string bcName = "Wall";
|
||||
# {
|
||||
# int bcType = 2;
|
||||
# int viscousType = 1;
|
||||
# double wallTemperature = -1.0;
|
||||
# double uWall = 0.0;
|
||||
# double vWall = 0.0;
|
||||
# double wWall = 0.0;
|
||||
# }
|
||||
# string bcName = "Inflow";
|
||||
# {
|
||||
# int bcType = 5;
|
||||
# int inflowParaType = 0;
|
||||
# double refMachNumber = 0.73;
|
||||
# double attackd = 2.79;
|
||||
# double angleSlide = 0.0;
|
||||
# double refReNumber = 6.5e6;
|
||||
# double refDimensionalTemperature = 288.15;
|
||||
# }
|
||||
|
||||
# For more information, see examples/bin/boundary_condition.hypara file!!!
|
||||
|
||||
int nBoundaryConditions = 3;
|
||||
string bcName = "wall";
|
||||
{
|
||||
string bodyName = "body";
|
||||
int bcType = 2;
|
||||
}
|
||||
string bcName = "inlet";
|
||||
{
|
||||
int bcType = 52;
|
||||
double totalPressure = 101625;
|
||||
double totalTemperature = 300;
|
||||
double direction_inlet[] = 1, 0, 0;
|
||||
int directionMethod = 1;
|
||||
}
|
||||
string bcName = "outlet";
|
||||
{
|
||||
int bcType = 62;
|
||||
double staticPressure = 101325;
|
||||
}
|
||||
|
||||
# 'bcType' is defined as following:
|
||||
# -2: WAKE
|
||||
# -1: INTERFACE
|
||||
# 0 : NO_BOUNDARY_CONDITION
|
||||
# 1 : EXTRAPOLATION
|
||||
# 2 : SOLID_SURFACE
|
||||
# 3 : SYMMETRY
|
||||
# 4 : FARFIELD
|
||||
# 5 : INFLOW
|
||||
# 6 : OUTFLOW
|
||||
# 52: PRESSURE_INLET
|
||||
# 62: PRESSURE_OUTLET
|
||||
# 61: OUTFLOW_CONFINED
|
||||
# 7 : POLE
|
||||
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 = 5000;
|
||||
|
||||
int intervalStepFlow = 1000;
|
||||
int intervalStepPlot = 1000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
|
||||
// 0 -- no precondition process. (default, mach > 0.3)
|
||||
// 1 -- carry out precondition process. (mach number <= 0.3)
|
||||
int ifLowSpeedPrecon = 1;
|
||||
|
||||
#########################################################################
|
||||
# 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 attackd = 0;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 5;
|
||||
double refDimensionalVelocity = 10;
|
||||
double refDimensionalPressure = 101325;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//double refMachNumber = 0.1;
|
||||
//double attackd = 0.00;
|
||||
//double angleSlide = 0.00;
|
||||
//
|
||||
//int inflowParaType = 0;
|
||||
//double refReNumber = 6.0e5;
|
||||
//double refDimensionalTemperature = 288.15;
|
||||
|
||||
double gridScaleFactor = 1.0;
|
||||
|
||||
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||
double forceReferenceLength = 1.0; // unit of meter.
|
||||
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||
double TorqueRefX = 0.0; // unit of meter.
|
||||
double TorqueRefY = 0.0; // unit of meter.
|
||||
double TorqueRefZ = 0.0; // unit of meter.
|
||||
|
||||
#########################################################################
|
||||
# Physical models #
|
||||
#########################################################################
|
||||
// viscousType : Viscous model.
|
||||
// 0 -- Euler.
|
||||
// 1 -- Laminar.
|
||||
// 3 -- 1eq turbulent.
|
||||
// 4 -- 2eq turbulent.
|
||||
// viscousName: Laminar or tubulent model.
|
||||
// -- "1eq-sa", when viscousType = 3.
|
||||
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||
// DESType: Type of DES.
|
||||
// 0 -- RANS.(default)
|
||||
// 1 -- DES.
|
||||
// 2 -- DDES.
|
||||
// 3 -- IDDES.
|
||||
|
||||
//int viscousType = 0;
|
||||
//string viscousName = "Euler";
|
||||
|
||||
int viscousType = 1;
|
||||
string viscousName = "laminar";
|
||||
|
||||
//int viscousType = 3;
|
||||
//string viscousName = "1eq-sa";
|
||||
|
||||
//int viscousType = 4;
|
||||
//string viscousName = "2eq-kw-menter-sst";
|
||||
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 3;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 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 = 100.0;
|
||||
|
||||
int nLUSGSSweeps = 1;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 1;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/filter1.fts";
|
||||
int plotFieldType = 1;
|
||||
|
||||
// ----------------- 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 = 14;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 38, 39, 40, 53, 54, 55, 63];
|
||||
|
||||
// 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;
|
||||
|
||||
int isPorousZone = 1;
|
||||
double viscousResistanceCoeff[] = [4.2e6 0.0 0.0];
|
||||
double inertialResistanceCoeff[] = [6.349 0.0 0.0];
|
||||
|
|
@ -0,0 +1,41 @@
|
|||
#########################################################################
|
||||
# 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.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order. The size of array "axisRotateOrder" is equal to nAxisRotateTimes.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
|
||||
int gridtype = 0;
|
||||
int axisup = 1;
|
||||
int from_gtype = 2;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[] = [0.0, 0.0, 0.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/filter1.cgns";
|
||||
string out_gfile = "./grid/filter1.fts";
|
||||
|
|
@ -0,0 +1,62 @@
|
|||
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/incompressible.hypara";
|
||||
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_para.hypara";
|
||||
|
||||
//int nsimutask = 2;
|
||||
//string parafilename = "./bin/cfd_para.hypara";
|
||||
|
||||
//int nsimutask = 3;
|
||||
//string parafilename = "./bin/partition.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||
|
||||
//int nsimutask = 5;
|
||||
//string parafilename = "./bin/overset_grid_view.hypara";
|
||||
|
||||
//int nsimutask = 13;
|
||||
//string parafilename = "./bin/lbm_para.hypara";
|
||||
|
||||
//int nsimutask = 14;
|
||||
//string parafilename = "./bin/integrative_solver.hypara";
|
||||
|
||||
//int nsimutask = 15;
|
||||
//string parafilename = "./bin/gridcheck_solver.hypara";
|
||||
|
||||
//int nsimutask = 99;
|
||||
//string parafilename = "./bin/post_processing.hypara";
|
||||
|
||||
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||
int numberOfGridProcessor = 0;
|
||||
// ATP read
|
||||
//@string parafilename1 = ""
|
||||
//@string parafilename2 = "";
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
Binary file not shown.
|
|
@ -0,0 +1,61 @@
|
|||
# nBoundaryConditions: Number of global boundary conditions.
|
||||
# bcName : Boundary condition name.
|
||||
# bcType(in PHengLEI): Boundary condition type.
|
||||
|
||||
# How to set boundary condition, for example:
|
||||
# string bcName = "Wall";
|
||||
# {
|
||||
# int bcType = 2;
|
||||
# int viscousType = 1;
|
||||
# double wallTemperature = -1.0;
|
||||
# double uWall = 0.0;
|
||||
# double vWall = 0.0;
|
||||
# double wWall = 0.0;
|
||||
# }
|
||||
# string bcName = "Inflow";
|
||||
# {
|
||||
# int bcType = 5;
|
||||
# int inflowParaType = 0;
|
||||
# double refMachNumber = 0.73;
|
||||
# double attackd = 2.79;
|
||||
# double angleSlide = 0.0;
|
||||
# double refReNumber = 6.5e6;
|
||||
# double refDimensionalTemperature = 288.15;
|
||||
# }
|
||||
|
||||
# For more information, see examples/bin/boundary_condition.hypara file!!!
|
||||
|
||||
int nBoundaryConditions = 3;
|
||||
string bcName = "wall";
|
||||
{
|
||||
string bodyName = "body";
|
||||
int bcType = 2;
|
||||
}
|
||||
string bcName = "inlet";
|
||||
{
|
||||
int bcType = 52;
|
||||
double totalPressure = 102095;
|
||||
double totalTemperature = 300;
|
||||
double direction_inlet[] = 1, 0, 0;
|
||||
int directionMethod = 1;
|
||||
}
|
||||
string bcName = "outlet";
|
||||
{
|
||||
int bcType = 62;
|
||||
double staticPressure = 101325;
|
||||
}
|
||||
|
||||
# 'bcType' is defined as following:
|
||||
# -2: WAKE
|
||||
# -1: INTERFACE
|
||||
# 0 : NO_BOUNDARY_CONDITION
|
||||
# 1 : EXTRAPOLATION
|
||||
# 2 : SOLID_SURFACE
|
||||
# 3 : SYMMETRY
|
||||
# 4 : FARFIELD
|
||||
# 5 : INFLOW
|
||||
# 6 : OUTFLOW
|
||||
# 52: PRESSURE_INLET
|
||||
# 62: PRESSURE_OUTLET
|
||||
# 61: OUTFLOW_CONFINED
|
||||
# 7 : POLE
|
||||
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 = 20000;
|
||||
|
||||
int intervalStepFlow = 1000;
|
||||
int intervalStepPlot = 1000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
|
||||
// 0 -- no precondition process. (default, mach > 0.3)
|
||||
// 1 -- carry out precondition process. (mach number <= 0.3)
|
||||
int ifLowSpeedPrecon = 0;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceReferenceLength, forceReferenceLengthSpanWise, forceReferenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double attackd = 0;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 5;
|
||||
double refDimensionalVelocity = 10;
|
||||
double refDimensionalPressure = 101325;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//double refMachNumber = 0.1;
|
||||
//double attackd = 0.00;
|
||||
//double angleSlide = 0.00;
|
||||
//
|
||||
//int inflowParaType = 0;
|
||||
//double refReNumber = 6.0e5;
|
||||
//double refDimensionalTemperature = 288.15;
|
||||
|
||||
double gridScaleFactor = 1.0;
|
||||
|
||||
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||
double forceReferenceLength = 1.0; // unit of meter.
|
||||
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||
double TorqueRefX = 0.0; // unit of meter.
|
||||
double TorqueRefY = 0.0; // unit of meter.
|
||||
double TorqueRefZ = 0.0; // unit of meter.
|
||||
|
||||
#########################################################################
|
||||
# Physical models #
|
||||
#########################################################################
|
||||
// viscousType : Viscous model.
|
||||
// 0 -- Euler.
|
||||
// 1 -- Laminar.
|
||||
// 3 -- 1eq turbulent.
|
||||
// 4 -- 2eq turbulent.
|
||||
// viscousName: Laminar or tubulent model.
|
||||
// -- "1eq-sa", when viscousType = 3.
|
||||
// -- "2eq-kw-menter-sst", when viscousType = 4.
|
||||
// DESType: Type of DES.
|
||||
// 0 -- RANS.(default)
|
||||
// 1 -- DES.
|
||||
// 2 -- DDES.
|
||||
// 3 -- IDDES.
|
||||
|
||||
//int viscousType = 0;
|
||||
//string viscousName = "Euler";
|
||||
|
||||
int viscousType = 1;
|
||||
string viscousName = "laminar";
|
||||
|
||||
//int viscousType = 3;
|
||||
//string viscousName = "1eq-sa";
|
||||
|
||||
//int viscousType = 4;
|
||||
//string viscousName = "2eq-kw-menter-sst";
|
||||
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 3;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 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 = 100.0;
|
||||
|
||||
int nLUSGSSweeps = 1;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 1;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/filter3_12W.fts";
|
||||
int plotFieldType = 1;
|
||||
|
||||
// ----------------- 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 = 14;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 38, 39, 40, 53, 54, 55, 63];
|
||||
|
||||
// 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;
|
||||
|
||||
int isPorousZone = 1;
|
||||
double viscousResistanceCoeff[] = [3.846e6 0.0 0.0];
|
||||
double inertialResistanceCoeff[] = [20.414 0.0 0.0];
|
||||
|
|
@ -0,0 +1,41 @@
|
|||
#########################################################################
|
||||
# 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.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order. The size of array "axisRotateOrder" is equal to nAxisRotateTimes.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
|
||||
int gridtype = 0;
|
||||
int axisup = 1;
|
||||
int from_gtype = 2;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[] = [0.0, 0.0, 0.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/filter3_12W.cgns";
|
||||
string out_gfile = "./grid/filter3_12W.fts";
|
||||
|
|
@ -0,0 +1,62 @@
|
|||
string title = "PHengLEI Main Parameter Control File";
|
||||
|
||||
// IMPORTANT NOTICE: DON NOT MODIFY THE FOWLLOWING LINE.
|
||||
string defaultParaFile = "./bin/cfd_para.hypara";
|
||||
|
||||
// ndim: Dimensional of the grid, 2 or 3.
|
||||
// nparafile: the number of parameter files.
|
||||
// nsimutask: simulation task type.
|
||||
// 0 -- CFD Solver of NS or Turbulation.
|
||||
// 1 -- Grid generation: for special typical cases, such as cylinder, flat plate, etc.
|
||||
// Grid conversion: from other format to PHengLEI format (.fts).
|
||||
// Grid reconstruction: such as grid adaptation.
|
||||
// Grid merging: merge two blocks into one block.
|
||||
// Grid repairing: repair the original grid in order to remove the negative volume cells.
|
||||
// 2 -- Wall distance computation for turb-solver.
|
||||
// 3 -- Grid partition.
|
||||
int ndim = 3;
|
||||
int nparafile = 1;
|
||||
|
||||
int nsimutask = 0;
|
||||
string parafilename = "./bin/cfd_para_subsonic.hypara";
|
||||
//string parafilename = "./bin/cfd_para_transonic.hypara";
|
||||
//string parafilename = "./bin/cfd_para_supersonic.hypara";
|
||||
//string parafilename = "./bin/cfd_para_hypersonic.hypara";
|
||||
//string parafilename = "./bin/incompressible.hypara";
|
||||
//string parafilename = "./bin/cfd_para_incompressible.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_para.hypara";
|
||||
|
||||
//int nsimutask = 2;
|
||||
//string parafilename = "./bin/cfd_para.hypara";
|
||||
|
||||
//int nsimutask = 3;
|
||||
//string parafilename = "./bin/partition.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_deform_para.hypara";
|
||||
|
||||
//int nsimutask = 1;
|
||||
//string parafilename = "./bin/grid_refine_para.hypara";
|
||||
|
||||
//int nsimutask = 5;
|
||||
//string parafilename = "./bin/overset_grid_view.hypara";
|
||||
|
||||
//int nsimutask = 13;
|
||||
//string parafilename = "./bin/lbm_para.hypara";
|
||||
|
||||
//int nsimutask = 14;
|
||||
//string parafilename = "./bin/integrative_solver.hypara";
|
||||
|
||||
//int nsimutask = 15;
|
||||
//string parafilename = "./bin/gridcheck_solver.hypara";
|
||||
|
||||
//int nsimutask = 99;
|
||||
//string parafilename = "./bin/post_processing.hypara";
|
||||
|
||||
// ---------------- Advanced Parameters, DO NOT care it ----------------
|
||||
int numberOfGridProcessor = 0;
|
||||
// ATP read
|
||||
//@string parafilename1 = ""
|
||||
//@string parafilename2 = "";
|
||||
|
|
@ -0,0 +1,13 @@
|
|||
# nVolumeConditions: number of global volume conditions.
|
||||
# vcName : Volume Condition Name.
|
||||
# vcType(in PHengLEI): Volume Condition Type.
|
||||
|
||||
int nVolumeConditions = 2;
|
||||
string vcName = "porous";
|
||||
{
|
||||
int vcType = 0;
|
||||
}
|
||||
string vcName = "airways";
|
||||
{
|
||||
int vcType = 1;
|
||||
}
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
Binary file not shown.
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,193 @@
|
|||
#########################################################################
|
||||
# General Control Parameter #
|
||||
#########################################################################
|
||||
// maxSimuStep: The max simulation step, don't care simulation is restart or not.
|
||||
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
|
||||
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
|
||||
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
|
||||
// intervalStepRes: The step intervals for residual 'res.dat' saved.
|
||||
|
||||
int maxSimuStep = 20000;
|
||||
|
||||
int intervalStepFlow = 2000;
|
||||
int intervalStepPlot = 2000;
|
||||
int intervalStepForce = 100;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// wallTemperature: Temprature of the solid wall, minus value is for adiabatic boundary condition.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 2.85;
|
||||
double attackd = 0.00;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
double wallTemperature = 276;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 7.391e7;
|
||||
double refDimensionalTemperature = 100.3;
|
||||
|
||||
//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 wallFunctionType = 1;
|
||||
int DESType = 0;
|
||||
|
||||
int roeEntropyFixMethod = 2;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*************************************************
|
||||
******************
|
||||
// inviscidSchemeName: Spatial discretisation scheme of struct grid.
|
||||
// Using this when solve structered grid or hybrid.
|
||||
// -- "vanleer", "steger", "ausmpw".
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "minmod", "3rd_minmod_smooth".
|
||||
|
||||
string inviscidSchemeName = "vanleer";
|
||||
string str_limiter_name = "minmod";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_scheme_name: Spatial discretisation scheme of Unstruct grid.
|
||||
// Using this when solve Unstructered grid or hybrid.
|
||||
// -- "vanleer", "roe", "steger", "kfvs", "lax_f", "hlle".
|
||||
// -- "ausm+", "ausmdv", "ausm+w", "ausmpw", "ausmpwplus".
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat", "barth".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_scheme_name = "vanleer";
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 0.5;
|
||||
|
||||
#########################################################################
|
||||
# Temporal Discretisation #
|
||||
#########################################################################
|
||||
// iunsteady: Steady or unsteady.
|
||||
// 0 -- steady.
|
||||
// 1 -- unsteay.
|
||||
// CFLEnd: The CFL number, [0.1, 100].
|
||||
// The bigger the value, the convergence faster but lower robustness.
|
||||
// ktmax: The lower the value, the more robustness, 1.0e5 - 1.0e10.
|
||||
|
||||
int iunsteady = 0;
|
||||
|
||||
double CFLEnd = 30.0;
|
||||
|
||||
#########################################################################
|
||||
# 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/unstr_1__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];
|
||||
|
||||
int nVisualWallVariables = 6;
|
||||
int visualWallVariables[] = [0, 1, 2, 3, 4, 5];
|
||||
|
||||
// limitVariables: Limit model (It is useful only if limitVector is 0).
|
||||
// 0 -- limit only for pressure and denstiny, then get the min value.
|
||||
// 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,35 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 0;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[]= [0.0, 0.0, 0.0];
|
||||
int from_gtype = 2;
|
||||
|
||||
#########################################################################
|
||||
# File path #
|
||||
#########################################################################
|
||||
// from_gfile: path of original data file for unstructure grid convert from.
|
||||
// out_gfile: path of target file for grid convert to, *.fts type of file usually.
|
||||
string from_gfile = "./grid/unstr_1.cgns";
|
||||
string out_gfile = "./grid/unstr_1.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 = 0;
|
||||
int maxproc = 4;
|
||||
|
||||
string original_grid_file = "./grid/unstr_1.fts";
|
||||
string partition_grid_file = "./grid/unstr_1__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;
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
|
|
@ -0,0 +1,69 @@
|
|||
# nBoundaryConditions: Number of global boundary conditions.
|
||||
# bcName : Boundary condition name.
|
||||
# bcType(in PHengLEI): Boundary condition type.
|
||||
|
||||
# How to set boundary condition, for example:
|
||||
# string bcName = "Wall";
|
||||
# {
|
||||
# int bcType = 2;
|
||||
# int viscousType = 1;
|
||||
# double wallTemperature = -1.0;
|
||||
# double uWall = 0.0;
|
||||
# double vWall = 0.0;
|
||||
# double wWall = 0.0;
|
||||
# }
|
||||
# string bcName = "Inflow";
|
||||
# {
|
||||
# int bcType = 5;
|
||||
# int inflowParaType = 0;
|
||||
# double refMachNumber = 0.73;
|
||||
# double attackd = 2.79;
|
||||
# double angleSlide = 0.0;
|
||||
# double refReNumber = 6.5e6;
|
||||
# double refDimensionalTemperature = 288.15;
|
||||
# }
|
||||
|
||||
# For more information, see examples/bin/boundary_condition.hypara file!!!
|
||||
|
||||
int nBoundaryConditions = 5;
|
||||
string bcName = "Wall";
|
||||
{
|
||||
string bodyName = "body";
|
||||
int bcType = 2;
|
||||
}
|
||||
string bcName = "Symmetry";
|
||||
{
|
||||
int bcType = 3;
|
||||
}
|
||||
string bcName = "Farfield";
|
||||
{
|
||||
int bcType = 4;
|
||||
}
|
||||
string bcName = "PressureInlet";
|
||||
{
|
||||
int bcType = 52;
|
||||
int directionMethod = 1;
|
||||
double totalPressure = 118309.784;
|
||||
double totalTemperature = 302.4;
|
||||
double direction_inlet[] = 1, 0, 0;
|
||||
}
|
||||
string bcName = "PressureOutlet";
|
||||
{
|
||||
int bcType = 62;
|
||||
double staticPressure = 115056.0;
|
||||
}
|
||||
|
||||
# 'bcType' is defined as following:
|
||||
# -2: WAKE
|
||||
# -1: INTERFACE
|
||||
# 0 : NO_BOUNDARY_CONDITION
|
||||
# 1 : EXTRAPOLATION
|
||||
# 2 : SOLID_SURFACE
|
||||
# 3 : SYMMETRY
|
||||
# 4 : FARFIELD
|
||||
# 5 : INFLOW
|
||||
# 6 : OUTFLOW
|
||||
# 52: PRESSURE_INLET
|
||||
# 62: PRESSURE_OUTLET
|
||||
# 61: OUTFLOW_CONFINED
|
||||
# 7 : POLE
|
||||
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 = 50000;
|
||||
|
||||
int intervalStepFlow = 2000;
|
||||
int intervalStepPlot = 2000;
|
||||
int intervalStepForce = 200;
|
||||
int intervalStepRes = 10;
|
||||
|
||||
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
|
||||
// 0 -- no precondition process. (default, mach > 0.3)
|
||||
// 1 -- carry out precondition process. (mach number <= 0.3)
|
||||
int ifLowSpeedPrecon = 1;
|
||||
|
||||
#########################################################################
|
||||
# Inflow Parameter #
|
||||
#########################################################################
|
||||
// refMachNumber: Mach number.
|
||||
// attackd: Angle of attack.
|
||||
// angleSlide: Angle of sideslip.
|
||||
// inflowParaType: The type of inflow parameters.
|
||||
// 0 -- the nondimensional conditions.
|
||||
// 1 -- the flight conditions.
|
||||
// 2 -- the experiment conditions.
|
||||
// 3 -- the subsonic boundary conditions.
|
||||
// refReNumber: Reynolds number, which is based unit length, unit of 1/m.
|
||||
// refDimensionalTemperature: Dimensional reference temperature, or the total temperature only for the experiment condition.
|
||||
// refDimensionalPressure: Dimensional reference pressure, or the total pressure only for the experiment condition.
|
||||
// height: Fly height, unit of km.
|
||||
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter.Common dimensions like:
|
||||
// 1 dm = 0.1 m.
|
||||
// 1 cm = 0.01 m.
|
||||
// 1 mm = 0.001m.
|
||||
// 1 inch = 0.0254m.
|
||||
// 1 foot = 12 inches = 0.3048m.
|
||||
// 1 yard = 3 feet = 0.9144m.
|
||||
// forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
|
||||
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
|
||||
|
||||
double refMachNumber = 0.2;
|
||||
double attackd = 0.00;
|
||||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 5e6;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//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 = 3;
|
||||
double roeEntropyScale = 1.0;
|
||||
#########################################################################
|
||||
# Spatial Discretisation #
|
||||
#########################################################################
|
||||
#*******************************************************************
|
||||
# Struct Solver *
|
||||
#*******************************************************************
|
||||
// str_limiter_name: Limiter of struct grid.
|
||||
// -- "3rdsmooth", "smooth".
|
||||
// -- "nolim", no limiter.
|
||||
|
||||
string str_limiter_name = "smooth";
|
||||
|
||||
#*******************************************************************
|
||||
# UnStruct Solver *
|
||||
#*******************************************************************
|
||||
// uns_limiter_name: Limiter of Unstruct grid.
|
||||
// -- "vencat".
|
||||
// -- "1st", meaning accuracy of first-order.
|
||||
// -- "nolim", no limiter.
|
||||
// venkatCoeff: Coefficient of vencat limiter, when uns_limiter_name = 'vencat'.
|
||||
// The smaller the value, the more robust it is.
|
||||
|
||||
string uns_limiter_name = "vencat";
|
||||
double venkatCoeff = 50.0;
|
||||
|
||||
#########################################################################
|
||||
# Temporal Discretisation #
|
||||
#########################################################################
|
||||
// iunsteady: Steady or unsteady.
|
||||
// 0 -- steady.
|
||||
// 1 -- unsteay.
|
||||
// CFLEnd: The CFL number, [0.1, 100].
|
||||
// The bigger the value, the convergence faster but lower robustness.
|
||||
// nLUSGSSweeps: Number of Sub-iteration of LU-SGS.
|
||||
// 1 -- is recommended for structured solver.
|
||||
// 1-3 -- is recommended for unstructured solver.
|
||||
|
||||
int iunsteady = 0;
|
||||
|
||||
double CFLEnd = 100.0;
|
||||
|
||||
int nLUSGSSweeps = 4;
|
||||
|
||||
#########################################################################
|
||||
# Multi-Grid parameters #
|
||||
#########################################################################
|
||||
// nMGLevel: The number of Multi-Grid level.
|
||||
// = 1 -- single-level.
|
||||
// > 1 -- multi-level.
|
||||
// flowInitStep: Flow initialization step, 0 - 500 is suggested.
|
||||
// Multi-Grid : Number of steps computing on coarse grid, during flow initialization.
|
||||
// Single-Grid: Number of steps computing using first-order with vanleer, during flow initialization.
|
||||
|
||||
int nMGLevel = 1;
|
||||
int flowInitStep = 100;
|
||||
|
||||
#########################################################################
|
||||
# File In or Out #
|
||||
#########################################################################
|
||||
// gridfile: The partitioned Grid file path, using relative path,
|
||||
// which is relative to the working directory.
|
||||
// IMPORTANT WARNING: The file index should be ignored,
|
||||
// e.g. if the partitioned grid is rae2822_hybrid2d__4_0.fts,
|
||||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/2d_plate_545x385_unstr__4.fts";
|
||||
int plotFieldType = 1;
|
||||
|
||||
// ----------------- 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 = 10;
|
||||
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 7, 8, 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,35 @@
|
|||
#########################################################################
|
||||
# Grid data type #
|
||||
#########################################################################
|
||||
// gridtype: Grid type for generation, conversion, reconstruction, merging.
|
||||
// 0 -- Unstructured grid.
|
||||
// 1 -- Structured grid.
|
||||
// nAxisRotateTimes: number of axis rotating times, zero (default) meaning without rotating.
|
||||
// axisRotateOrder : axis rotating order.
|
||||
// 1 -- X-axis.
|
||||
// 2 -- Y-axis.
|
||||
// 3 -- Z-axis.
|
||||
// axisRotateAngles: axis rotating angles (degree), which are corresponding to the axis rotating order.
|
||||
// from_gtype: Type of grid data type in grid conversion process.
|
||||
// -1 -- MULTI_TYPE.
|
||||
// 1 -- PHengLEI, *.fts.
|
||||
// 2 -- CGNS, *.cgns.
|
||||
// 3 -- Plot3D type of structured grid, *.dat/*.grd.
|
||||
// 4 -- Fieldview type of unstructured grid, *.dat/*.inp.
|
||||
// 5 -- Fluent, *.cas/*.msh.
|
||||
// 6 -- Ustar, mgrid.in.
|
||||
// 7 -- Hybrid, include both of unstructured and structured grid, *.fts.
|
||||
// 8 -- GMSH, *.msh.
|
||||
int gridtype = 0;
|
||||
int nAxisRotateTimes = 0;
|
||||
int axisRotateOrder[] = [1, 2, 3];
|
||||
double axisRotateAngles[] = [0.0, 0.0, 0.0];
|
||||
int from_gtype = 2;
|
||||
|
||||
#########################################################################
|
||||
# File path #
|
||||
#########################################################################
|
||||
// from_gfile: path of original data file for unstructure grid convert from.
|
||||
// out_gfile: path of target file for grid convert to, *.fts type of file usually.
|
||||
string from_gfile = "./grid/2d_plate_545x385_unstr.cgns";
|
||||
string out_gfile = "./grid/2d_plate_545x385_unstr.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 = 0;
|
||||
int maxproc = 4;
|
||||
|
||||
string original_grid_file = "./grid/2d_plate_545x385_unstr.fts";
|
||||
string partition_grid_file = "./grid/2d_plate_545x385_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;
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
Binary file not shown.
|
|
@ -44,8 +44,8 @@ double attackd = 3.06;
|
|||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 1.814e7;
|
||||
double refDimensionalTemperature = 270;
|
||||
double refReNumber = 1.801e7;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//int inflowParaType = 1;
|
||||
//double height = 0.001;
|
||||
|
|
@ -58,9 +58,9 @@ double wallTemperature = -1.0;
|
|||
|
||||
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 forceReferenceLengthSpanWise = 1.1963; // unit of meter.
|
||||
double forceReferenceLength = 0.64607; // unit of meter.
|
||||
double forceReferenceArea = 0.7532; // 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.
|
||||
|
|
|
|||
Binary file not shown.
|
|
@ -44,7 +44,7 @@ double attackd = -0.304;
|
|||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 2.12464589e7;
|
||||
double refReNumber = 2.12e7;
|
||||
double refDimensionalTemperature = 274.1;
|
||||
|
||||
//int inflowParaType = 1;
|
||||
|
|
@ -58,9 +58,9 @@ double wallTemperature = -1.0;
|
|||
|
||||
double gridScaleFactor = 0.001;
|
||||
|
||||
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||
double forceReferenceLength = 1.0; // unit of meter.
|
||||
double forceReferenceArea = 1.0; // unit of meter^2.
|
||||
double forceReferenceLengthSpanWise = 0.58565; // unit of meter.
|
||||
double forceReferenceLength = 0.1412; // unit of meter.
|
||||
double forceReferenceArea = 0.0727; // 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.
|
||||
|
|
|
|||
|
|
@ -44,8 +44,8 @@ double attackd = 3.06;
|
|||
double angleSlide = 0.00;
|
||||
|
||||
int inflowParaType = 0;
|
||||
double refReNumber = 1.171e7;
|
||||
double refDimensionalTemperature = 288;
|
||||
double refReNumber = 1.801e7;
|
||||
double refDimensionalTemperature = 300;
|
||||
|
||||
//int inflowParaType = 1;
|
||||
//double height = 0.001;
|
||||
|
|
@ -56,9 +56,9 @@ double refDimensionalTemperature = 288;
|
|||
|
||||
double gridScaleFactor = 1.0;
|
||||
|
||||
double forceRefenenceLengthSpanWise = 1.0; // unit of meter.
|
||||
double forceRefenenceLength = 1.0; // unit of meter.
|
||||
double forceRefenenceArea = 1.0; // unit of meter^2.
|
||||
double forceReferenceLengthSpanWise = 1.0; // unit of meter.
|
||||
double forceReferenceLength = 0.53703; // unit of meter.
|
||||
double forceReferenceArea = 0.52041; // 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.
|
||||
|
|
|
|||
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
|
@ -265,7 +265,7 @@ int flowInitStep = 0;
|
|||
// Please use 'rae2822_hybrid2d__4.fts' here!
|
||||
// plotFieldType: If dump out the whole field results to tecplot or not, 0 / 1.
|
||||
|
||||
string gridfile = "./grid/naca0012_str_26m_yup__400.fts";
|
||||
string gridfile = "./grid/naca0012_str_26m_yup__200.fts";
|
||||
int plotFieldType = 1;
|
||||
|
||||
// ----------------- Advanced Parameters, DO NOT care it ----------------
|
||||
|
|
@ -9,10 +9,10 @@
|
|||
// partition_grid_file: Target partition grid file(PHengLEI type, *.fts).
|
||||
|
||||
int pgridtype = 1;
|
||||
int maxproc = 400;
|
||||
int maxproc = 200;
|
||||
|
||||
string original_grid_file = "./grid/naca0012_str_26m_yup.fts";
|
||||
string partition_grid_file = "./grid/naca0012_str_26m_yup__400.fts";
|
||||
string partition_grid_file = "./grid/naca0012_str_26m_yup__200.fts";
|
||||
|
||||
// numberOfMultigrid: Number of multi-grid levels, ONLY used for structured grid.
|
||||
// 1 -- single level.
|
||||
|
|
@ -0,0 +1,5 @@
|
|||
红山开源风雷算例库原始网格获取百度网盘链接：
|
||||
链接：http://pan.baidu.com/s/1aZ9cdkp6CkT9il4fEpnTcA
|
||||
提取码：w47m
|
||||
|
||||
注：plot3D格式网格需同时下载.grd和.inp文件
|
||||
Binary file not shown.
Binary file not shown.
Binary file not shown.
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,34 @@
|
|||
0 ！MRT set 1, BGK set 0
|
||||
0 solver16600a.plt ！CONTI, resume :1 , new 0 //follows input datafile name
|
||||
0 ！LES yes: 1 , no: 0
|
||||
1 mesh_2d.dat ！GEO read complex geo? yes:1, no:0 //follows input datafile name
|
||||
1 ! MB multiblock LBM? yes:1, no:0
|
||||
0 2.e-3 !CONV using converge criterion? yes:1, no:0
|
||||
600 ！NX
|
||||
400 ! NY
|
||||
1 ! NZ
|
||||
240 ！ NX2
|
||||
160 ！ NY2
|
||||
1 ! NZ2
|
||||
200 ! LowX
|
||||
160 ! LowY
|
||||
0 ! LowZ
|
||||
100 ! framerate
|
||||
100000 ！total steps to terminate
|
||||
1.0 ! density
|
||||
0.54 ！ tau
|
||||
(0.0 0.0, 0.0000000) ！ Volumetric force vector (gx, gy, gz)
|
||||
D2Q9 ! velocity set
|
||||
Poiseuille ! boundary condition
|
||||
! note: in the following, boundary types: periodic,velocity, pressure
|
||||
! are allowed. if "velocity" is applied, please strictly follow the format "velocity (%lf, %lf, %lf)".
|
||||
velocity (0.1, 0.00, 0.0) ！xmin face
|
||||
pressure 1.0 ！xmax face
|
||||
velocity (0.1, 0.00, 0.0) ！ymin face
|
||||
velocity (0.1, 0.00, 0.0) ！ymax face,
|
||||
nonslip ！zmin face
|
||||
nonslip ！zmax face
|
||||
-------------------！ initial velocity field
|
||||
(0.1, 0.00, 0.0)
|
||||
1 ! Multiplt geo number
|
||||
0 ! use IBPeskin or not
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
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 = 17;
|
||||
string parafilename = "./bin/cfd_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,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
1 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
1 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.3 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
4.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.004 10000 ! kbend、kstrech
|
||||
0.0 0.4 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
2.0 4.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,321 @@
|
|||
RectPlat
|
||||
157 157 1
|
||||
END
|
||||
157 X Y Z
|
||||
1 0.5000000000 0.0000000000 0.0000000000
|
||||
2 0.4996000533 0.0199946671 0.0000000000
|
||||
3 0.4984008532 0.0399573470 0.0000000000
|
||||
4 0.4964043179 0.0598561036 0.0000000000
|
||||
5 0.4936136417 0.0796591033 0.0000000000
|
||||
6 0.4900332889 0.0993346654 0.0000000000
|
||||
7 0.4856689874 0.1188513132 0.0000000000
|
||||
8 0.4805277192 0.1381778243 0.0000000000
|
||||
9 0.4746177090 0.1572832803 0.0000000000
|
||||
10 0.4679484118 0.1761371166 0.0000000000
|
||||
11 0.4605304970 0.1947091712 0.0000000000
|
||||
12 0.4523758316 0.2129697325 0.0000000000
|
||||
13 0.4434974614 0.2308895878 0.0000000000
|
||||
14 0.4339095898 0.2484400689 0.0000000000
|
||||
15 0.4236275555 0.2655930990 0.0000000000
|
||||
16 0.4126678075 0.2823212367 0.0000000000
|
||||
17 0.4010478789 0.2985977207 0.0000000000
|
||||
18 0.3887863594 0.3143965120 0.0000000000
|
||||
19 0.3759028646 0.3296923360 0.0000000000
|
||||
20 0.3624180054 0.3444607226 0.0000000000
|
||||
21 0.3483533547 0.3586780454 0.0000000000
|
||||
22 0.3337314129 0.3723215600 0.0000000000
|
||||
23 0.3185755721 0.3853694394 0.0000000000
|
||||
24 0.3029100783 0.3978008100 0.0000000000
|
||||
25 0.2867599930 0.4095957842 0.0000000000
|
||||
26 0.2701511529 0.4207354924 0.0000000000
|
||||
27 0.2531101286 0.4312021136 0.0000000000
|
||||
28 0.2356641821 0.4409789034 0.0000000000
|
||||
29 0.2178412231 0.4500502211 0.0000000000
|
||||
30 0.1996697647 0.4584015544 0.0000000000
|
||||
31 0.1811788772 0.4660195430 0.0000000000
|
||||
32 0.1623981422 0.4728919997 0.0000000000
|
||||
33 0.1433576048 0.4790079301 0.0000000000
|
||||
34 0.1240877258 0.4843575501 0.0000000000
|
||||
35 0.1046193329 0.4889323012 0.0000000000
|
||||
36 0.0849835715 0.4927248650 0.0000000000
|
||||
37 0.0652118544 0.4957291741 0.0000000000
|
||||
38 0.0453358122 0.4979404223 0.0000000000
|
||||
39 0.0253872425 0.4993550720 0.0000000000
|
||||
40 0.0053980585 0.4999708601 0.0000000000
|
||||
41 -0.0145997612 0.4997868015 0.0000000000
|
||||
42 -0.0345742243 0.4988031907 0.0000000000
|
||||
43 -0.0544933761 0.4970216011 0.0000000000
|
||||
44 -0.0743253501 0.4944448830 0.0000000000
|
||||
45 -0.0940384194 0.4910771586 0.0000000000
|
||||
46 -0.1136010473 0.4869238154 0.0000000000
|
||||
47 -0.1329819378 0.4819914981 0.0000000000
|
||||
48 -0.1521500855 0.4762880971 0.0000000000
|
||||
49 -0.1710748256 0.4698227368 0.0000000000
|
||||
50 -0.1897258824 0.4626057604 0.0000000000
|
||||
51 -0.2080734183 0.4546487134 0.0000000000
|
||||
52 -0.2260880811 0.4459643255 0.0000000000
|
||||
53 -0.2437410512 0.4365664898 0.0000000000
|
||||
54 -0.2610040876 0.4264702408 0.0000000000
|
||||
55 -0.2778495731 0.4156917304 0.0000000000
|
||||
56 -0.2942505586 0.4042482019 0.0000000000
|
||||
57 -0.3101808060 0.3921579625 0.0000000000
|
||||
58 -0.3256148303 0.3794403541 0.0000000000
|
||||
59 -0.3405279403 0.3661157220 0.0000000000
|
||||
60 -0.3548962782 0.3522053829 0.0000000000
|
||||
61 -0.3686968578 0.3377315903 0.0000000000
|
||||
62 -0.3819076010 0.3227174992 0.0000000000
|
||||
63 -0.3945073736 0.3071871289 0.0000000000
|
||||
64 -0.4064760185 0.2911653248 0.0000000000
|
||||
65 -0.4177943886 0.2746777182 0.0000000000
|
||||
66 -0.4284443767 0.2577506859 0.0000000000
|
||||
67 -0.4384089452 0.2404113075 0.0000000000
|
||||
68 -0.4476721529 0.2226873223 0.0000000000
|
||||
69 -0.4562191807 0.2046070848 0.0000000000
|
||||
70 -0.4640363551 0.1861995197 0.0000000000
|
||||
71 -0.4711111703 0.1674940751 0.0000000000
|
||||
72 -0.4774323082 0.1485206757 0.0000000000
|
||||
73 -0.4829896562 0.1293096748 0.0000000000
|
||||
74 -0.4877743238 0.1098918061 0.0000000000
|
||||
75 -0.4917786565 0.0902981339 0.0000000000
|
||||
76 -0.4949962483 0.0705600040 0.0000000000
|
||||
77 -0.4974219517 0.0507089932 0.0000000000
|
||||
78 -0.4990518860 0.0307768587 0.0000000000
|
||||
79 -0.4998834439 0.0107954879 0.0000000000
|
||||
80 -0.4999152948 -0.0092031535 0.0000000000
|
||||
81 -0.4991473879 -0.0291870717 0.0000000000
|
||||
82 -0.4975809517 -0.0491242969 0.0000000000
|
||||
83 -0.4952184920 -0.0689829336 0.0000000000
|
||||
84 -0.4920637885 -0.0887312124 0.0000000000
|
||||
85 -0.4881218878 -0.1083375402 0.0000000000
|
||||
86 -0.4833990963 -0.1277705510 0.0000000000
|
||||
87 -0.4779029693 -0.1469991562 0.0000000000
|
||||
88 -0.4716422995 -0.1659925941 0.0000000000
|
||||
89 -0.4646271027 -0.1847204793 0.0000000000
|
||||
90 -0.4568686016 -0.2031528511 0.0000000000
|
||||
91 -0.4483792082 -0.2212602216 0.0000000000
|
||||
92 -0.4391725037 -0.2390136231 0.0000000000
|
||||
93 -0.4292632169 -0.2563846537 0.0000000000
|
||||
94 -0.4186672005 -0.2733455235 0.0000000000
|
||||
95 -0.4074014059 -0.2898690989 0.0000000000
|
||||
96 -0.3954838560 -0.3059289455 0.0000000000
|
||||
97 -0.3829336162 -0.3214993710 0.0000000000
|
||||
98 -0.3697707644 -0.3365554662 0.0000000000
|
||||
99 -0.3560163582 -0.3510731444 0.0000000000
|
||||
100 -0.3416924018 -0.3650291804 0.0000000000
|
||||
101 -0.3268218104 -0.3784012477 0.0000000000
|
||||
102 -0.3114283739 -0.3911679536 0.0000000000
|
||||
103 -0.2955367184 -0.4033088743 0.0000000000
|
||||
104 -0.2791722672 -0.4148045868 0.0000000000
|
||||
105 -0.2623611999 -0.4256367005 0.0000000000
|
||||
106 -0.2451304107 -0.4357878862 0.0000000000
|
||||
107 -0.2275074651 -0.4452419043 0.0000000000
|
||||
108 -0.2095205561 -0.4539836303 0.0000000000
|
||||
109 -0.1911984589 -0.4619990794 0.0000000000
|
||||
110 -0.1725704849 -0.4692754284 0.0000000000
|
||||
111 -0.1536664350 -0.4758010369 0.0000000000
|
||||
112 -0.1345165516 -0.4815654653 0.0000000000
|
||||
113 -0.1151514703 -0.4865594916 0.0000000000
|
||||
114 -0.0956021713 -0.4907751265 0.0000000000
|
||||
115 -0.0758999292 -0.4942056260 0.0000000000
|
||||
116 -0.0560762635 -0.4968455018 0.0000000000
|
||||
117 -0.0361628876 -0.4986905308 0.0000000000
|
||||
118 -0.0161916589 -0.4997377614 0.0000000000
|
||||
119 0.0038054731 -0.4999855182 0.0000000000
|
||||
120 0.0237965171 -0.4994334047 0.0000000000
|
||||
121 0.0437494917 -0.4980823044 0.0000000000
|
||||
122 0.0636324765 -0.4959343787 0.0000000000
|
||||
123 0.0834136629 -0.4929930637 0.0000000000
|
||||
124 0.1030614053 -0.4892630650 0.0000000000
|
||||
125 0.1225442713 -0.4847503497 0.0000000000
|
||||
126 0.1418310927 -0.4794621373 0.0000000000
|
||||
127 0.1608910146 -0.4734068878 0.0000000000
|
||||
128 0.1796935452 -0.4665942882 0.0000000000
|
||||
129 0.1982086045 -0.4590352373 0.0000000000
|
||||
130 0.2164065722 -0.4507418280 0.0000000000
|
||||
131 0.2342583357 -0.4417273279 0.0000000000
|
||||
132 0.2517353357 -0.4320061582 0.0000000000
|
||||
133 0.2688096129 -0.4215938709 0.0000000000
|
||||
134 0.2854538521 -0.4105071234 0.0000000000
|
||||
135 0.3016414260 -0.3987636520 0.0000000000
|
||||
136 0.3173464380 -0.3863822438 0.0000000000
|
||||
137 0.3325437633 -0.3733827064 0.0000000000
|
||||
138 0.3472090896 -0.3597858364 0.0000000000
|
||||
139 0.3613189554 -0.3456133858 0.0000000000
|
||||
140 0.3748507880 -0.3308880275 0.0000000000
|
||||
141 0.3877829393 -0.3156333189 0.0000000000
|
||||
142 0.4000947205 -0.2998736644 0.0000000000
|
||||
143 0.4117664356 -0.2836342760 0.0000000000
|
||||
144 0.4227794122 -0.2669411332 0.0000000000
|
||||
145 0.4331160320 -0.2498209416 0.0000000000
|
||||
146 0.4427597585 -0.2323010897 0.0000000000
|
||||
147 0.4516951638 -0.2144096057 0.0000000000
|
||||
148 0.4599079532 -0.1961751120 0.0000000000
|
||||
149 0.4673849880 -0.1776267800 0.0000000000
|
||||
150 0.4741143065 -0.1587942830 0.0000000000
|
||||
151 0.4800851433 -0.1397077491 0.0000000000
|
||||
152 0.4852879463 -0.1203977126 0.0000000000
|
||||
153 0.4897143921 -0.1008950654 0.0000000000
|
||||
154 0.4933573993 -0.0812310076 0.0000000000
|
||||
155 0.4962111399 -0.0614369976 0.0000000000
|
||||
156 0.4982710485 -0.0415447014 0.0000000000
|
||||
157 0.4995338298 -0.0215859426 0.0000000000
|
||||
END
|
||||
157 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
51 51 52 52
|
||||
52 52 53 53
|
||||
53 53 54 54
|
||||
54 54 55 55
|
||||
55 55 56 56
|
||||
56 56 57 57
|
||||
57 57 58 58
|
||||
58 58 59 59
|
||||
59 59 60 60
|
||||
60 60 61 61
|
||||
61 61 62 62
|
||||
62 62 63 63
|
||||
63 63 64 64
|
||||
64 64 65 65
|
||||
65 65 66 66
|
||||
66 66 67 67
|
||||
67 67 68 68
|
||||
68 68 69 69
|
||||
69 69 70 70
|
||||
70 70 71 71
|
||||
71 71 72 72
|
||||
72 72 73 73
|
||||
73 73 74 74
|
||||
74 74 75 75
|
||||
75 75 76 76
|
||||
76 76 77 77
|
||||
77 77 78 78
|
||||
78 78 79 79
|
||||
79 79 80 80
|
||||
80 80 81 81
|
||||
81 81 82 82
|
||||
82 82 83 83
|
||||
83 83 84 84
|
||||
84 84 85 85
|
||||
85 85 86 86
|
||||
86 86 87 87
|
||||
87 87 88 88
|
||||
88 88 89 89
|
||||
89 89 90 90
|
||||
90 90 91 91
|
||||
91 91 92 92
|
||||
92 92 93 93
|
||||
93 93 94 94
|
||||
94 94 95 95
|
||||
95 95 96 96
|
||||
96 96 97 97
|
||||
97 97 98 98
|
||||
98 98 99 99
|
||||
99 99 100 100
|
||||
100 100 101 101
|
||||
101 101 102 102
|
||||
102 102 103 103
|
||||
103 103 104 104
|
||||
104 104 105 105
|
||||
105 105 106 106
|
||||
106 106 107 107
|
||||
107 107 108 108
|
||||
108 108 109 109
|
||||
109 109 110 110
|
||||
110 110 111 111
|
||||
111 111 112 112
|
||||
112 112 113 113
|
||||
113 113 114 114
|
||||
114 114 115 115
|
||||
115 115 116 116
|
||||
116 116 117 117
|
||||
117 117 118 118
|
||||
118 118 119 119
|
||||
119 119 120 120
|
||||
120 120 121 121
|
||||
121 121 122 122
|
||||
122 122 123 123
|
||||
123 123 124 124
|
||||
124 124 125 125
|
||||
125 125 126 126
|
||||
126 126 127 127
|
||||
127 127 128 128
|
||||
128 128 129 129
|
||||
129 129 130 130
|
||||
130 130 131 131
|
||||
131 131 132 132
|
||||
132 132 133 133
|
||||
133 133 134 134
|
||||
134 134 135 135
|
||||
135 135 136 136
|
||||
136 136 137 137
|
||||
137 137 138 138
|
||||
138 138 139 139
|
||||
139 139 140 140
|
||||
140 140 141 141
|
||||
141 141 142 142
|
||||
142 142 143 143
|
||||
143 143 144 144
|
||||
144 144 145 145
|
||||
145 145 146 146
|
||||
146 146 147 147
|
||||
147 147 148 148
|
||||
148 148 149 149
|
||||
149 149 150 150
|
||||
150 150 151 151
|
||||
151 151 152 152
|
||||
152 152 153 153
|
||||
153 153 154 154
|
||||
154 154 155 155
|
||||
155 155 156 156
|
||||
156 156 157 157
|
||||
157 157 1 1
|
||||
END
|
||||
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,34 @@
|
|||
0 ！MRT set 1, BGK set 0
|
||||
0 solver16600a.plt ！CONTI, resume :1 , new 0 //follows input datafile name
|
||||
0 ！LES yes: 1 , no: 0
|
||||
0 mesh_2d.dat ！GEO read complex geo? yes:1, no:0 //follows input datafile name
|
||||
0 ! MB multiblock LBM? yes:1, no:0
|
||||
0 2.e-3 !CONV using converge criterion? yes:1, no:0
|
||||
800 ！NX
|
||||
600 ! NY
|
||||
1 ! NZ
|
||||
240 ！ NX2
|
||||
160 ！ NY2
|
||||
1 ! NZ2
|
||||
200 ! LowX
|
||||
160 ! LowY
|
||||
0 ! LowZ
|
||||
150 ! framerate
|
||||
100000 ！total steps to terminate
|
||||
1.0 ! density
|
||||
0.54 ！ tau
|
||||
(0.0 0.0, 0.0000000) ！ Volumetric force vector (gx, gy, gz)
|
||||
D2Q9 ! velocity set
|
||||
Poiseuille ! boundary condition
|
||||
! note: in the following, boundary types: periodic,velocity, pressure
|
||||
! are allowed. if "velocity" is applied, please strictly follow the format "velocity (%lf, %lf, %lf)".
|
||||
velocity (0.1, 0.00, 0.0) ！xmin face
|
||||
pressure 1.0 ！xmax face
|
||||
velocity (0.1, 0.00, 0.0) ！ymin face
|
||||
velocity (0.1, 0.00, 0.0) ！ymax face,
|
||||
nonslip ！zmin face
|
||||
nonslip ！zmax face
|
||||
-------------------！ initial velocity field
|
||||
(0.1, 0.00, 0.0)
|
||||
2 ! Multiplt geo number
|
||||
1 ! use IBPeskin or not
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
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 = 17;
|
||||
string parafilename = "./bin/cfd_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,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
1 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
1 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.3 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
4.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.004 10000 ! kbend、kstrech
|
||||
0.0 0.4 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
2.0 4.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
2 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
2 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.3 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
4.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.004 10000 ! kbend、kstrech
|
||||
0.0 0.0 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
2.0 2.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,321 @@
|
|||
RectPlat
|
||||
157 157 1
|
||||
END
|
||||
157 X Y Z
|
||||
1 0.5000000000 0.0000000000 0.0000000000
|
||||
2 0.4996000533 0.0199946671 0.0000000000
|
||||
3 0.4984008532 0.0399573470 0.0000000000
|
||||
4 0.4964043179 0.0598561036 0.0000000000
|
||||
5 0.4936136417 0.0796591033 0.0000000000
|
||||
6 0.4900332889 0.0993346654 0.0000000000
|
||||
7 0.4856689874 0.1188513132 0.0000000000
|
||||
8 0.4805277192 0.1381778243 0.0000000000
|
||||
9 0.4746177090 0.1572832803 0.0000000000
|
||||
10 0.4679484118 0.1761371166 0.0000000000
|
||||
11 0.4605304970 0.1947091712 0.0000000000
|
||||
12 0.4523758316 0.2129697325 0.0000000000
|
||||
13 0.4434974614 0.2308895878 0.0000000000
|
||||
14 0.4339095898 0.2484400689 0.0000000000
|
||||
15 0.4236275555 0.2655930990 0.0000000000
|
||||
16 0.4126678075 0.2823212367 0.0000000000
|
||||
17 0.4010478789 0.2985977207 0.0000000000
|
||||
18 0.3887863594 0.3143965120 0.0000000000
|
||||
19 0.3759028646 0.3296923360 0.0000000000
|
||||
20 0.3624180054 0.3444607226 0.0000000000
|
||||
21 0.3483533547 0.3586780454 0.0000000000
|
||||
22 0.3337314129 0.3723215600 0.0000000000
|
||||
23 0.3185755721 0.3853694394 0.0000000000
|
||||
24 0.3029100783 0.3978008100 0.0000000000
|
||||
25 0.2867599930 0.4095957842 0.0000000000
|
||||
26 0.2701511529 0.4207354924 0.0000000000
|
||||
27 0.2531101286 0.4312021136 0.0000000000
|
||||
28 0.2356641821 0.4409789034 0.0000000000
|
||||
29 0.2178412231 0.4500502211 0.0000000000
|
||||
30 0.1996697647 0.4584015544 0.0000000000
|
||||
31 0.1811788772 0.4660195430 0.0000000000
|
||||
32 0.1623981422 0.4728919997 0.0000000000
|
||||
33 0.1433576048 0.4790079301 0.0000000000
|
||||
34 0.1240877258 0.4843575501 0.0000000000
|
||||
35 0.1046193329 0.4889323012 0.0000000000
|
||||
36 0.0849835715 0.4927248650 0.0000000000
|
||||
37 0.0652118544 0.4957291741 0.0000000000
|
||||
38 0.0453358122 0.4979404223 0.0000000000
|
||||
39 0.0253872425 0.4993550720 0.0000000000
|
||||
40 0.0053980585 0.4999708601 0.0000000000
|
||||
41 -0.0145997612 0.4997868015 0.0000000000
|
||||
42 -0.0345742243 0.4988031907 0.0000000000
|
||||
43 -0.0544933761 0.4970216011 0.0000000000
|
||||
44 -0.0743253501 0.4944448830 0.0000000000
|
||||
45 -0.0940384194 0.4910771586 0.0000000000
|
||||
46 -0.1136010473 0.4869238154 0.0000000000
|
||||
47 -0.1329819378 0.4819914981 0.0000000000
|
||||
48 -0.1521500855 0.4762880971 0.0000000000
|
||||
49 -0.1710748256 0.4698227368 0.0000000000
|
||||
50 -0.1897258824 0.4626057604 0.0000000000
|
||||
51 -0.2080734183 0.4546487134 0.0000000000
|
||||
52 -0.2260880811 0.4459643255 0.0000000000
|
||||
53 -0.2437410512 0.4365664898 0.0000000000
|
||||
54 -0.2610040876 0.4264702408 0.0000000000
|
||||
55 -0.2778495731 0.4156917304 0.0000000000
|
||||
56 -0.2942505586 0.4042482019 0.0000000000
|
||||
57 -0.3101808060 0.3921579625 0.0000000000
|
||||
58 -0.3256148303 0.3794403541 0.0000000000
|
||||
59 -0.3405279403 0.3661157220 0.0000000000
|
||||
60 -0.3548962782 0.3522053829 0.0000000000
|
||||
61 -0.3686968578 0.3377315903 0.0000000000
|
||||
62 -0.3819076010 0.3227174992 0.0000000000
|
||||
63 -0.3945073736 0.3071871289 0.0000000000
|
||||
64 -0.4064760185 0.2911653248 0.0000000000
|
||||
65 -0.4177943886 0.2746777182 0.0000000000
|
||||
66 -0.4284443767 0.2577506859 0.0000000000
|
||||
67 -0.4384089452 0.2404113075 0.0000000000
|
||||
68 -0.4476721529 0.2226873223 0.0000000000
|
||||
69 -0.4562191807 0.2046070848 0.0000000000
|
||||
70 -0.4640363551 0.1861995197 0.0000000000
|
||||
71 -0.4711111703 0.1674940751 0.0000000000
|
||||
72 -0.4774323082 0.1485206757 0.0000000000
|
||||
73 -0.4829896562 0.1293096748 0.0000000000
|
||||
74 -0.4877743238 0.1098918061 0.0000000000
|
||||
75 -0.4917786565 0.0902981339 0.0000000000
|
||||
76 -0.4949962483 0.0705600040 0.0000000000
|
||||
77 -0.4974219517 0.0507089932 0.0000000000
|
||||
78 -0.4990518860 0.0307768587 0.0000000000
|
||||
79 -0.4998834439 0.0107954879 0.0000000000
|
||||
80 -0.4999152948 -0.0092031535 0.0000000000
|
||||
81 -0.4991473879 -0.0291870717 0.0000000000
|
||||
82 -0.4975809517 -0.0491242969 0.0000000000
|
||||
83 -0.4952184920 -0.0689829336 0.0000000000
|
||||
84 -0.4920637885 -0.0887312124 0.0000000000
|
||||
85 -0.4881218878 -0.1083375402 0.0000000000
|
||||
86 -0.4833990963 -0.1277705510 0.0000000000
|
||||
87 -0.4779029693 -0.1469991562 0.0000000000
|
||||
88 -0.4716422995 -0.1659925941 0.0000000000
|
||||
89 -0.4646271027 -0.1847204793 0.0000000000
|
||||
90 -0.4568686016 -0.2031528511 0.0000000000
|
||||
91 -0.4483792082 -0.2212602216 0.0000000000
|
||||
92 -0.4391725037 -0.2390136231 0.0000000000
|
||||
93 -0.4292632169 -0.2563846537 0.0000000000
|
||||
94 -0.4186672005 -0.2733455235 0.0000000000
|
||||
95 -0.4074014059 -0.2898690989 0.0000000000
|
||||
96 -0.3954838560 -0.3059289455 0.0000000000
|
||||
97 -0.3829336162 -0.3214993710 0.0000000000
|
||||
98 -0.3697707644 -0.3365554662 0.0000000000
|
||||
99 -0.3560163582 -0.3510731444 0.0000000000
|
||||
100 -0.3416924018 -0.3650291804 0.0000000000
|
||||
101 -0.3268218104 -0.3784012477 0.0000000000
|
||||
102 -0.3114283739 -0.3911679536 0.0000000000
|
||||
103 -0.2955367184 -0.4033088743 0.0000000000
|
||||
104 -0.2791722672 -0.4148045868 0.0000000000
|
||||
105 -0.2623611999 -0.4256367005 0.0000000000
|
||||
106 -0.2451304107 -0.4357878862 0.0000000000
|
||||
107 -0.2275074651 -0.4452419043 0.0000000000
|
||||
108 -0.2095205561 -0.4539836303 0.0000000000
|
||||
109 -0.1911984589 -0.4619990794 0.0000000000
|
||||
110 -0.1725704849 -0.4692754284 0.0000000000
|
||||
111 -0.1536664350 -0.4758010369 0.0000000000
|
||||
112 -0.1345165516 -0.4815654653 0.0000000000
|
||||
113 -0.1151514703 -0.4865594916 0.0000000000
|
||||
114 -0.0956021713 -0.4907751265 0.0000000000
|
||||
115 -0.0758999292 -0.4942056260 0.0000000000
|
||||
116 -0.0560762635 -0.4968455018 0.0000000000
|
||||
117 -0.0361628876 -0.4986905308 0.0000000000
|
||||
118 -0.0161916589 -0.4997377614 0.0000000000
|
||||
119 0.0038054731 -0.4999855182 0.0000000000
|
||||
120 0.0237965171 -0.4994334047 0.0000000000
|
||||
121 0.0437494917 -0.4980823044 0.0000000000
|
||||
122 0.0636324765 -0.4959343787 0.0000000000
|
||||
123 0.0834136629 -0.4929930637 0.0000000000
|
||||
124 0.1030614053 -0.4892630650 0.0000000000
|
||||
125 0.1225442713 -0.4847503497 0.0000000000
|
||||
126 0.1418310927 -0.4794621373 0.0000000000
|
||||
127 0.1608910146 -0.4734068878 0.0000000000
|
||||
128 0.1796935452 -0.4665942882 0.0000000000
|
||||
129 0.1982086045 -0.4590352373 0.0000000000
|
||||
130 0.2164065722 -0.4507418280 0.0000000000
|
||||
131 0.2342583357 -0.4417273279 0.0000000000
|
||||
132 0.2517353357 -0.4320061582 0.0000000000
|
||||
133 0.2688096129 -0.4215938709 0.0000000000
|
||||
134 0.2854538521 -0.4105071234 0.0000000000
|
||||
135 0.3016414260 -0.3987636520 0.0000000000
|
||||
136 0.3173464380 -0.3863822438 0.0000000000
|
||||
137 0.3325437633 -0.3733827064 0.0000000000
|
||||
138 0.3472090896 -0.3597858364 0.0000000000
|
||||
139 0.3613189554 -0.3456133858 0.0000000000
|
||||
140 0.3748507880 -0.3308880275 0.0000000000
|
||||
141 0.3877829393 -0.3156333189 0.0000000000
|
||||
142 0.4000947205 -0.2998736644 0.0000000000
|
||||
143 0.4117664356 -0.2836342760 0.0000000000
|
||||
144 0.4227794122 -0.2669411332 0.0000000000
|
||||
145 0.4331160320 -0.2498209416 0.0000000000
|
||||
146 0.4427597585 -0.2323010897 0.0000000000
|
||||
147 0.4516951638 -0.2144096057 0.0000000000
|
||||
148 0.4599079532 -0.1961751120 0.0000000000
|
||||
149 0.4673849880 -0.1776267800 0.0000000000
|
||||
150 0.4741143065 -0.1587942830 0.0000000000
|
||||
151 0.4800851433 -0.1397077491 0.0000000000
|
||||
152 0.4852879463 -0.1203977126 0.0000000000
|
||||
153 0.4897143921 -0.1008950654 0.0000000000
|
||||
154 0.4933573993 -0.0812310076 0.0000000000
|
||||
155 0.4962111399 -0.0614369976 0.0000000000
|
||||
156 0.4982710485 -0.0415447014 0.0000000000
|
||||
157 0.4995338298 -0.0215859426 0.0000000000
|
||||
END
|
||||
157 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
51 51 52 52
|
||||
52 52 53 53
|
||||
53 53 54 54
|
||||
54 54 55 55
|
||||
55 55 56 56
|
||||
56 56 57 57
|
||||
57 57 58 58
|
||||
58 58 59 59
|
||||
59 59 60 60
|
||||
60 60 61 61
|
||||
61 61 62 62
|
||||
62 62 63 63
|
||||
63 63 64 64
|
||||
64 64 65 65
|
||||
65 65 66 66
|
||||
66 66 67 67
|
||||
67 67 68 68
|
||||
68 68 69 69
|
||||
69 69 70 70
|
||||
70 70 71 71
|
||||
71 71 72 72
|
||||
72 72 73 73
|
||||
73 73 74 74
|
||||
74 74 75 75
|
||||
75 75 76 76
|
||||
76 76 77 77
|
||||
77 77 78 78
|
||||
78 78 79 79
|
||||
79 79 80 80
|
||||
80 80 81 81
|
||||
81 81 82 82
|
||||
82 82 83 83
|
||||
83 83 84 84
|
||||
84 84 85 85
|
||||
85 85 86 86
|
||||
86 86 87 87
|
||||
87 87 88 88
|
||||
88 88 89 89
|
||||
89 89 90 90
|
||||
90 90 91 91
|
||||
91 91 92 92
|
||||
92 92 93 93
|
||||
93 93 94 94
|
||||
94 94 95 95
|
||||
95 95 96 96
|
||||
96 96 97 97
|
||||
97 97 98 98
|
||||
98 98 99 99
|
||||
99 99 100 100
|
||||
100 100 101 101
|
||||
101 101 102 102
|
||||
102 102 103 103
|
||||
103 103 104 104
|
||||
104 104 105 105
|
||||
105 105 106 106
|
||||
106 106 107 107
|
||||
107 107 108 108
|
||||
108 108 109 109
|
||||
109 109 110 110
|
||||
110 110 111 111
|
||||
111 111 112 112
|
||||
112 112 113 113
|
||||
113 113 114 114
|
||||
114 114 115 115
|
||||
115 115 116 116
|
||||
116 116 117 117
|
||||
117 117 118 118
|
||||
118 118 119 119
|
||||
119 119 120 120
|
||||
120 120 121 121
|
||||
121 121 122 122
|
||||
122 122 123 123
|
||||
123 123 124 124
|
||||
124 124 125 125
|
||||
125 125 126 126
|
||||
126 126 127 127
|
||||
127 127 128 128
|
||||
128 128 129 129
|
||||
129 129 130 130
|
||||
130 130 131 131
|
||||
131 131 132 132
|
||||
132 132 133 133
|
||||
133 133 134 134
|
||||
134 134 135 135
|
||||
135 135 136 136
|
||||
136 136 137 137
|
||||
137 137 138 138
|
||||
138 138 139 139
|
||||
139 139 140 140
|
||||
140 140 141 141
|
||||
141 141 142 142
|
||||
142 142 143 143
|
||||
143 143 144 144
|
||||
144 144 145 145
|
||||
145 145 146 146
|
||||
146 146 147 147
|
||||
147 147 148 148
|
||||
148 148 149 149
|
||||
149 149 150 150
|
||||
150 150 151 151
|
||||
151 151 152 152
|
||||
152 152 153 153
|
||||
153 153 154 154
|
||||
154 154 155 155
|
||||
155 155 156 156
|
||||
156 156 157 157
|
||||
157 157 1 1
|
||||
END
|
||||
|
|
@ -0,0 +1,108 @@
|
|||
Filament
|
||||
51 50 1
|
||||
END
|
||||
51 X Y Z
|
||||
1 0.5400000000 0.0000000000 0.0000000000
|
||||
2 0.5600000000 0.0000000000 0.0000000000
|
||||
3 0.5800000000 0.0000000000 0.0000000000
|
||||
4 0.6000000000 0.0000000000 0.0000000000
|
||||
5 0.6200000000 0.0000000000 0.0000000000
|
||||
6 0.6400000000 0.0000000000 0.0000000000
|
||||
7 0.6600000000 0.0000000000 0.0000000000
|
||||
8 0.6800000000 0.0000000000 0.0000000000
|
||||
9 0.7000000000 0.0000000000 0.0000000000
|
||||
10 0.7200000000 0.0000000000 0.0000000000
|
||||
11 0.7400000000 0.0000000000 0.0000000000
|
||||
12 0.7600000000 0.0000000000 0.0000000000
|
||||
13 0.7800000000 0.0000000000 0.0000000000
|
||||
14 0.8000000000 0.0000000000 0.0000000000
|
||||
15 0.8200000000 0.0000000000 0.0000000000
|
||||
16 0.8400000000 0.0000000000 0.0000000000
|
||||
17 0.8600000000 0.0000000000 0.0000000000
|
||||
18 0.8800000000 0.0000000000 0.0000000000
|
||||
19 0.9000000000 0.0000000000 0.0000000000
|
||||
20 0.9200000000 0.0000000000 0.0000000000
|
||||
21 0.9400000000 0.0000000000 0.0000000000
|
||||
22 0.9600000000 0.0000000000 0.0000000000
|
||||
23 0.9800000000 0.0000000000 0.0000000000
|
||||
24 1.0000000000 0.0000000000 0.0000000000
|
||||
25 1.0200000000 0.0000000000 0.0000000000
|
||||
26 1.0400000000 0.0000000000 0.0000000000
|
||||
27 1.0600000000 0.0000000000 0.0000000000
|
||||
28 1.0800000000 0.0000000000 0.0000000000
|
||||
29 1.1000000000 0.0000000000 0.0000000000
|
||||
30 1.1200000000 0.0000000000 0.0000000000
|
||||
31 1.1400000000 0.0000000000 0.0000000000
|
||||
32 1.1600000000 0.0000000000 0.0000000000
|
||||
33 1.1800000000 0.0000000000 0.0000000000
|
||||
34 1.2000000000 0.0000000000 0.0000000000
|
||||
35 1.2200000000 0.0000000000 0.0000000000
|
||||
36 1.2400000000 0.0000000000 0.0000000000
|
||||
37 1.2600000000 0.0000000000 0.0000000000
|
||||
38 1.2800000000 0.0000000000 0.0000000000
|
||||
39 1.3000000000 0.0000000000 0.0000000000
|
||||
40 1.3200000000 0.0000000000 0.0000000000
|
||||
41 1.3400000000 0.0000000000 0.0000000000
|
||||
42 1.3600000000 0.0000000000 0.0000000000
|
||||
43 1.3800000000 0.0000000000 0.0000000000
|
||||
44 1.4000000000 0.0000000000 0.0000000000
|
||||
45 1.4200000000 0.0000000000 0.0000000000
|
||||
46 1.4400000000 0.0000000000 0.0000000000
|
||||
47 1.4600000000 0.0000000000 0.0000000000
|
||||
48 1.4800000000 0.0000000000 0.0000000000
|
||||
49 1.5000000000 0.0000000000 0.0000000000
|
||||
50 1.5200000000 0.0000000000 0.0000000000
|
||||
51 1.5400000000 0.0000000000 0.0000000000
|
||||
END
|
||||
50 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
END
|
||||
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,34 @@
|
|||
0 ！MRT set 1, BGK set 0
|
||||
0 solver16600a.plt ！CONTI, resume :1 , new 0 //follows input datafile name
|
||||
0 ！LES yes: 1 , no: 0
|
||||
0 mesh_2d.dat ！GEO read complex geo? yes:1, no:0 //follows input datafile name
|
||||
0 ! MB multiblock LBM? yes:1, no:0
|
||||
0 2.e-3 !CONV using converge criterion? yes:1, no:0
|
||||
800 ！NX
|
||||
600 ! NY
|
||||
1 ! NZ
|
||||
240 ！ NX2
|
||||
160 ！ NY2
|
||||
1 ! NZ2
|
||||
200 ! LowX
|
||||
160 ! LowY
|
||||
0 ! LowZ
|
||||
250 ! framerate
|
||||
100000 ！total steps to terminate
|
||||
1.0 ! density
|
||||
0.54 ！ tau
|
||||
(0.0 0.0, 0.0000000) ！ Volumetric force vector (gx, gy, gz)
|
||||
D2Q9 ! velocity set
|
||||
Poiseuille ! boundary condition
|
||||
! note: in the following, boundary types: periodic,velocity, pressure
|
||||
! are allowed. if "velocity" is applied, please strictly follow the format "velocity (%lf, %lf, %lf)".
|
||||
velocity (0.1, 0.00, 0.0) ！xmin face
|
||||
pressure 1.0 ！xmax face
|
||||
velocity (0.1, 0.00, 0.0) ！ymin face
|
||||
velocity (0.1, 0.00, 0.0) ！ymax face,
|
||||
nonslip ！zmin face
|
||||
nonslip ！zmax face
|
||||
-------------------！ initial velocity field
|
||||
(0.1, 0.00, 0.0)
|
||||
1 ! Multiplt geo number
|
||||
1 ! use IBPeskin or not
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
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 = 17;
|
||||
string parafilename = "./bin/cfd_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,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
2 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
2 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.3 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
4.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.01 10000 ! kbend、kstrech
|
||||
0.4 0.0 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
5.0 4.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,108 @@
|
|||
Filament
|
||||
51 50 1
|
||||
END
|
||||
51 X Y Z
|
||||
1 0.0000000000 0.0000000000 0.0000000000
|
||||
2 0.0200000000 0.0000000000 0.0000000000
|
||||
3 0.0400000000 0.0000000000 0.0000000000
|
||||
4 0.0600000000 0.0000000000 0.0000000000
|
||||
5 0.0800000000 0.0000000000 0.0000000000
|
||||
6 0.1000000000 0.0000000000 0.0000000000
|
||||
7 0.1200000000 0.0000000000 0.0000000000
|
||||
8 0.1400000000 0.0000000000 0.0000000000
|
||||
9 0.1600000000 0.0000000000 0.0000000000
|
||||
10 0.1800000000 0.0000000000 0.0000000000
|
||||
11 0.2000000000 0.0000000000 0.0000000000
|
||||
12 0.2200000000 0.0000000000 0.0000000000
|
||||
13 0.2400000000 0.0000000000 0.0000000000
|
||||
14 0.2600000000 0.0000000000 0.0000000000
|
||||
15 0.2800000000 0.0000000000 0.0000000000
|
||||
16 0.3000000000 0.0000000000 0.0000000000
|
||||
17 0.3200000000 0.0000000000 0.0000000000
|
||||
18 0.3400000000 0.0000000000 0.0000000000
|
||||
19 0.3600000000 0.0000000000 0.0000000000
|
||||
20 0.3800000000 0.0000000000 0.0000000000
|
||||
21 0.4000000000 0.0000000000 0.0000000000
|
||||
22 0.4200000000 0.0000000000 0.0000000000
|
||||
23 0.4400000000 0.0000000000 0.0000000000
|
||||
24 0.4600000000 0.0000000000 0.0000000000
|
||||
25 0.4800000000 0.0000000000 0.0000000000
|
||||
26 0.5000000000 0.0000000000 0.0000000000
|
||||
27 0.5200000000 0.0000000000 0.0000000000
|
||||
28 0.5400000000 0.0000000000 0.0000000000
|
||||
29 0.5600000000 0.0000000000 0.0000000000
|
||||
30 0.5800000000 0.0000000000 0.0000000000
|
||||
31 0.6000000000 0.0000000000 0.0000000000
|
||||
32 0.6200000000 0.0000000000 0.0000000000
|
||||
33 0.6400000000 0.0000000000 0.0000000000
|
||||
34 0.6600000000 0.0000000000 0.0000000000
|
||||
35 0.6800000000 0.0000000000 0.0000000000
|
||||
36 0.7000000000 0.0000000000 0.0000000000
|
||||
37 0.7200000000 0.0000000000 0.0000000000
|
||||
38 0.7400000000 0.0000000000 0.0000000000
|
||||
39 0.7600000000 0.0000000000 0.0000000000
|
||||
40 0.7800000000 0.0000000000 0.0000000000
|
||||
41 0.8000000000 0.0000000000 0.0000000000
|
||||
42 0.8200000000 0.0000000000 0.0000000000
|
||||
43 0.8400000000 0.0000000000 0.0000000000
|
||||
44 0.8600000000 0.0000000000 0.0000000000
|
||||
45 0.8800000000 0.0000000000 0.0000000000
|
||||
46 0.9000000000 0.0000000000 0.0000000000
|
||||
47 0.9200000000 0.0000000000 0.0000000000
|
||||
48 0.9400000000 0.0000000000 0.0000000000
|
||||
49 0.9600000000 0.0000000000 0.0000000000
|
||||
50 0.9800000000 0.0000000000 0.0000000000
|
||||
51 1.0000000000 0.0000000000 0.0000000000
|
||||
END
|
||||
50 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
END
|
||||
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,34 @@
|
|||
0 ！MRT set 1, BGK set 0
|
||||
0 solver16600a.plt ！CONTI, resume :1 , new 0 //follows input datafile name
|
||||
0 ！LES yes: 1 , no: 0
|
||||
0 mesh_2d.dat ！GEO read complex geo? yes:1, no:0 //follows input datafile name
|
||||
0 ! MB multiblock LBM? yes:1, no:0
|
||||
0 2.e-3 !CONV using converge criterion? yes:1, no:0
|
||||
1600 ！NX
|
||||
800 ! NY
|
||||
1 ! NZ
|
||||
240 ！ NX2
|
||||
160 ！ NY2
|
||||
1 ! NZ2
|
||||
200 ! LowX
|
||||
160 ! LowY
|
||||
0 ! LowZ
|
||||
250 ! framerate
|
||||
100000 ！total steps to terminate
|
||||
1.0 ! density
|
||||
0.54 ！ tau
|
||||
(0.0 0.0, 0.0000000) ！ Volumetric force vector (gx, gy, gz)
|
||||
D2Q9 ! velocity set
|
||||
Poiseuille ! boundary condition
|
||||
! note: in the following, boundary types: periodic,velocity, pressure
|
||||
! are allowed. if "velocity" is applied, please strictly follow the format "velocity (%lf, %lf, %lf)".
|
||||
velocity (0.1, 0.00, 0.0) ！xmin face
|
||||
pressure 1.0 ！xmax face
|
||||
velocity (0.1, 0.00, 0.0) ！ymin face
|
||||
velocity (0.1, 0.00, 0.0) ！ymax face,
|
||||
nonslip ！zmin face
|
||||
nonslip ！zmax face
|
||||
-------------------！ initial velocity field
|
||||
(0.1, 0.00, 0.0)
|
||||
1 ! Multiplt geo number
|
||||
1 ! use IBPeskin or not
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
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 = 17;
|
||||
string parafilename = "./bin/cfd_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,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
2 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
2 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.3 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
4.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.0005 10000 ! kbend、kstrech
|
||||
0.0 0.4 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
2.0 5.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,108 @@
|
|||
Filament
|
||||
51 50 1
|
||||
END
|
||||
51 X Y Z
|
||||
1 0.0000000000 0.0000000000 0.0000000000
|
||||
2 0.0200000000 0.0000000000 0.0000000000
|
||||
3 0.0400000000 0.0000000000 0.0000000000
|
||||
4 0.0600000000 0.0000000000 0.0000000000
|
||||
5 0.0800000000 0.0000000000 0.0000000000
|
||||
6 0.1000000000 0.0000000000 0.0000000000
|
||||
7 0.1200000000 0.0000000000 0.0000000000
|
||||
8 0.1400000000 0.0000000000 0.0000000000
|
||||
9 0.1600000000 0.0000000000 0.0000000000
|
||||
10 0.1800000000 0.0000000000 0.0000000000
|
||||
11 0.2000000000 0.0000000000 0.0000000000
|
||||
12 0.2200000000 0.0000000000 0.0000000000
|
||||
13 0.2400000000 0.0000000000 0.0000000000
|
||||
14 0.2600000000 0.0000000000 0.0000000000
|
||||
15 0.2800000000 0.0000000000 0.0000000000
|
||||
16 0.3000000000 0.0000000000 0.0000000000
|
||||
17 0.3200000000 0.0000000000 0.0000000000
|
||||
18 0.3400000000 0.0000000000 0.0000000000
|
||||
19 0.3600000000 0.0000000000 0.0000000000
|
||||
20 0.3800000000 0.0000000000 0.0000000000
|
||||
21 0.4000000000 0.0000000000 0.0000000000
|
||||
22 0.4200000000 0.0000000000 0.0000000000
|
||||
23 0.4400000000 0.0000000000 0.0000000000
|
||||
24 0.4600000000 0.0000000000 0.0000000000
|
||||
25 0.4800000000 0.0000000000 0.0000000000
|
||||
26 0.5000000000 0.0000000000 0.0000000000
|
||||
27 0.5200000000 0.0000000000 0.0000000000
|
||||
28 0.5400000000 0.0000000000 0.0000000000
|
||||
29 0.5600000000 0.0000000000 0.0000000000
|
||||
30 0.5800000000 0.0000000000 0.0000000000
|
||||
31 0.6000000000 0.0000000000 0.0000000000
|
||||
32 0.6200000000 0.0000000000 0.0000000000
|
||||
33 0.6400000000 0.0000000000 0.0000000000
|
||||
34 0.6600000000 0.0000000000 0.0000000000
|
||||
35 0.6800000000 0.0000000000 0.0000000000
|
||||
36 0.7000000000 0.0000000000 0.0000000000
|
||||
37 0.7200000000 0.0000000000 0.0000000000
|
||||
38 0.7400000000 0.0000000000 0.0000000000
|
||||
39 0.7600000000 0.0000000000 0.0000000000
|
||||
40 0.7800000000 0.0000000000 0.0000000000
|
||||
41 0.8000000000 0.0000000000 0.0000000000
|
||||
42 0.8200000000 0.0000000000 0.0000000000
|
||||
43 0.8400000000 0.0000000000 0.0000000000
|
||||
44 0.8600000000 0.0000000000 0.0000000000
|
||||
45 0.8800000000 0.0000000000 0.0000000000
|
||||
46 0.9000000000 0.0000000000 0.0000000000
|
||||
47 0.9200000000 0.0000000000 0.0000000000
|
||||
48 0.9400000000 0.0000000000 0.0000000000
|
||||
49 0.9600000000 0.0000000000 0.0000000000
|
||||
50 0.9800000000 0.0000000000 0.0000000000
|
||||
51 1.0000000000 0.0000000000 0.0000000000
|
||||
END
|
||||
50 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
END
|
||||
File diff suppressed because it is too large
Load Diff
|
|
@ -0,0 +1,34 @@
|
|||
0 ！MRT set 1, BGK set 0
|
||||
0 solver16600a.plt ！CONTI, resume :1 , new 0 //follows input datafile name
|
||||
0 ！LES yes: 1 , no: 0
|
||||
0 mesh_2d.dat ！GEO read complex geo? yes:1, no:0 //follows input datafile name
|
||||
0 ! MB multiblock LBM? yes:1, no:0
|
||||
0 2.e-3 !CONV using converge criterion? yes:1, no:0
|
||||
1600 ！NX
|
||||
800 ! NY
|
||||
1 ! NZ
|
||||
240 ！ NX2
|
||||
160 ！ NY2
|
||||
1 ! NZ2
|
||||
200 ! LowX
|
||||
160 ! LowY
|
||||
0 ! LowZ
|
||||
250 ! framerate
|
||||
80000 ！total steps to terminate
|
||||
1.0 ! density
|
||||
0.54 ！ tau
|
||||
(0.0 0.0, 0.0000000) ！ Volumetric force vector (gx, gy, gz)
|
||||
D2Q9 ! velocity set
|
||||
Poiseuille ! boundary condition
|
||||
! note: in the following, boundary types: periodic,velocity, pressure
|
||||
! are allowed. if "velocity" is applied, please strictly follow the format "velocity (%lf, %lf, %lf)".
|
||||
velocity (0.1, 0.00, 0.0) ！xmin face
|
||||
pressure 1.0 ！xmax face
|
||||
velocity (0.1, 0.00, 0.0) ！ymin face
|
||||
velocity (0.1, 0.00, 0.0) ！ymax face,
|
||||
nonslip ！zmin face
|
||||
nonslip ！zmax face
|
||||
-------------------！ initial velocity field
|
||||
(0.1, 0.00, 0.0)
|
||||
1 ! Multiplt geo number
|
||||
1 ! use IBPeskin or not
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
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 = 17;
|
||||
string parafilename = "./bin/cfd_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,25 @@
|
|||
Title = "IBFSI input file"
|
||||
!------------basic paramters------------
|
||||
0.02 ! ds: ds = dt
|
||||
(0.2 0.5) ! Position of the zero point of solid structure
|
||||
20 ! OpenMP using cores
|
||||
100 ! Re: Reynolds number, used to calculate tau in LBMSolver
|
||||
6 ! Supp: The distribution function affects the region
|
||||
1 ! iMotion: 1 prescribed rigid motion、2 flexilble fillament model、3 free rigid motion；
|
||||
2 ! subMotionType: Detailed requirements with different iMotion
|
||||
(0.0 0.0 0.0) ! gravity: have dimension
|
||||
!------------when iMotion is equal to 1 and subMotionType1 is equal to 2------------
|
||||
1 ! Is rotate move or not
|
||||
0.0 0.0 ! Initial Angle and initial phase of the Angle
|
||||
0.4 ! Amplitude of rotation: the amplitude of the trailing edge
|
||||
5.0 ! Period of pitching: a multiple of the dimensionless period
|
||||
!------------when iMotion is equal to 2------------
|
||||
0.004 10000 ! kbend、kstrech
|
||||
0.0 0.4 ! Amplitude of heaving(displacement)、pitching(angle)
|
||||
2.0 4.0 ! Period of heaving、pitching: a multiple of the dimensionless period(Cannot equal zero)
|
||||
(0.0 0.0) ! Leading edge fixed point accelerated speed
|
||||
!------------when iMotion is equal to 3------------
|
||||
10 ! Mass Ratio
|
||||
0.2 0.2 ! Dimensionless natural frequency: y、theta
|
||||
0.0 0.0 ! Damp ratio: y、theta
|
||||
(0.0 0.0) ! Centroid position: (x, y)
|
||||
|
|
@ -0,0 +1,108 @@
|
|||
Filament
|
||||
51 50 1
|
||||
END
|
||||
51 X Y Z
|
||||
1 0.0000000000 0.0000000000 0.0000000000
|
||||
2 0.0200000000 0.0000000000 0.0000000000
|
||||
3 0.0400000000 0.0000000000 0.0000000000
|
||||
4 0.0600000000 0.0000000000 0.0000000000
|
||||
5 0.0800000000 0.0000000000 0.0000000000
|
||||
6 0.1000000000 0.0000000000 0.0000000000
|
||||
7 0.1200000000 0.0000000000 0.0000000000
|
||||
8 0.1400000000 0.0000000000 0.0000000000
|
||||
9 0.1600000000 0.0000000000 0.0000000000
|
||||
10 0.1800000000 0.0000000000 0.0000000000
|
||||
11 0.2000000000 0.0000000000 0.0000000000
|
||||
12 0.2200000000 0.0000000000 0.0000000000
|
||||
13 0.2400000000 0.0000000000 0.0000000000
|
||||
14 0.2600000000 0.0000000000 0.0000000000
|
||||
15 0.2800000000 0.0000000000 0.0000000000
|
||||
16 0.3000000000 0.0000000000 0.0000000000
|
||||
17 0.3200000000 0.0000000000 0.0000000000
|
||||
18 0.3400000000 0.0000000000 0.0000000000
|
||||
19 0.3600000000 0.0000000000 0.0000000000
|
||||
20 0.3800000000 0.0000000000 0.0000000000
|
||||
21 0.4000000000 0.0000000000 0.0000000000
|
||||
22 0.4200000000 0.0000000000 0.0000000000
|
||||
23 0.4400000000 0.0000000000 0.0000000000
|
||||
24 0.4600000000 0.0000000000 0.0000000000
|
||||
25 0.4800000000 0.0000000000 0.0000000000
|
||||
26 0.5000000000 0.0000000000 0.0000000000
|
||||
27 0.5200000000 0.0000000000 0.0000000000
|
||||
28 0.5400000000 0.0000000000 0.0000000000
|
||||
29 0.5600000000 0.0000000000 0.0000000000
|
||||
30 0.5800000000 0.0000000000 0.0000000000
|
||||
31 0.6000000000 0.0000000000 0.0000000000
|
||||
32 0.6200000000 0.0000000000 0.0000000000
|
||||
33 0.6400000000 0.0000000000 0.0000000000
|
||||
34 0.6600000000 0.0000000000 0.0000000000
|
||||
35 0.6800000000 0.0000000000 0.0000000000
|
||||
36 0.7000000000 0.0000000000 0.0000000000
|
||||
37 0.7200000000 0.0000000000 0.0000000000
|
||||
38 0.7400000000 0.0000000000 0.0000000000
|
||||
39 0.7600000000 0.0000000000 0.0000000000
|
||||
40 0.7800000000 0.0000000000 0.0000000000
|
||||
41 0.8000000000 0.0000000000 0.0000000000
|
||||
42 0.8200000000 0.0000000000 0.0000000000
|
||||
43 0.8400000000 0.0000000000 0.0000000000
|
||||
44 0.8600000000 0.0000000000 0.0000000000
|
||||
45 0.8800000000 0.0000000000 0.0000000000
|
||||
46 0.9000000000 0.0000000000 0.0000000000
|
||||
47 0.9200000000 0.0000000000 0.0000000000
|
||||
48 0.9400000000 0.0000000000 0.0000000000
|
||||
49 0.9600000000 0.0000000000 0.0000000000
|
||||
50 0.9800000000 0.0000000000 0.0000000000
|
||||
51 1.0000000000 0.0000000000 0.0000000000
|
||||
END
|
||||
50 I J K
|
||||
1 1 2 2
|
||||
2 2 3 3
|
||||
3 3 4 4
|
||||
4 4 5 5
|
||||
5 5 6 6
|
||||
6 6 7 7
|
||||
7 7 8 8
|
||||
8 8 9 9
|
||||
9 9 10 10
|
||||
10 10 11 11
|
||||
11 11 12 12
|
||||
12 12 13 13
|
||||
13 13 14 14
|
||||
14 14 15 15
|
||||
15 15 16 16
|
||||
16 16 17 17
|
||||
17 17 18 18
|
||||
18 18 19 19
|
||||
19 19 20 20
|
||||
20 20 21 21
|
||||
21 21 22 22
|
||||
22 22 23 23
|
||||
23 23 24 24
|
||||
24 24 25 25
|
||||
25 25 26 26
|
||||
26 26 27 27
|
||||
27 27 28 28
|
||||
28 28 29 29
|
||||
29 29 30 30
|
||||
30 30 31 31
|
||||
31 31 32 32
|
||||
32 32 33 33
|
||||
33 33 34 34
|
||||
34 34 35 35
|
||||
35 35 36 36
|
||||
36 36 37 37
|
||||
37 37 38 38
|
||||
38 38 39 39
|
||||
39 39 40 40
|
||||
40 40 41 41
|
||||
41 41 42 42
|
||||
42 42 43 43
|
||||
43 43 44 44
|
||||
44 44 45 45
|
||||
45 45 46 46
|
||||
46 46 47 47
|
||||
47 47 48 48
|
||||
48 48 49 49
|
||||
49 49 50 50
|
||||
50 50 51 51
|
||||
END
|
||||
File diff suppressed because it is too large
Load Diff
Some files were not shown because too many files have changed in this diff Show More
Loading…
Reference in New Issue