This commit is contained in:
hechao 2021-03-26 15:53:14 +08:00
parent 7f8e0990db
commit 7774cb2a9a
16 changed files with 731 additions and 600 deletions

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// viscousType : Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,18 +215,20 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
// 0 -- no precondition process. (default, mach > 0.3) // 0 -- no precondition process. (default, mach > 0.3)
// 1 -- carry out precondition process. (mach number <= 0.3) // 1 -- carry out precondition process. (mach number <= 0.3)
int maxSimuStep = 3000; int maxSimuStep = 20000;
int intervalStepFlow = 100; int intervalStepFlow = 1000;
int intervalStepPlot = 100; int intervalStepPlot = 1000;
int intervalStepForce = 1; int intervalStepSample = 1000;
int intervalStepRes = 1; int intervalStepForce = 100;
int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
// compressible: // compressible:
@ -247,18 +259,24 @@ int compressible = 1;
// 2 -- dump out wall Q & the typical position Q of ball. // 2 -- dump out wall Q & the typical position Q of ball.
// 3 -- dump out wall Q & the typical position Q of cone. // 3 -- dump out wall Q & the typical position Q of cone.
// 4 -- dump out wall Q & the typical position Q of double sphere. // 4 -- dump out wall Q & the typical position Q of double sphere.
// gridScaleFactor: The customizable unit of the grid, default value is 1.0 for meter. // 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. // forceRefenenceLength, forceRefenenceLengthSpanWise, forceRefenenceArea: Reference length, SpanWise length and area, independent of grid unit.
// TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit. // TorqueRefX, TorqueRefY, TorqueRefZ: Reference point, independent of grid unit.
// radiationCoef: the radiation coefficient on wall, it is used to compute the radiation heat flux on wall when the boundary // radiationCoef: The radiation coefficient on wall, it is used to compute the radiation heat flux on wall when the boundary
// condition is radiation equilibrium temperature, and 0.8 is the default value. // condition is radiation equilibrium temperature, and 0.8 is the default value.
double refMachNumber = 0.20; double refMachNumber = 0.73;
double attackd = 0.00; double attackd = 2.79;
double angleSlide = 0.00; double angleSlide = 0.00;
int inflowParaType = 0; int inflowParaType = 0;
double refReNumber = 3900; double refReNumber = 6.5e6;
double refDimensionalTemperature = 288.15; double refDimensionalTemperature = 288.15;
double freestream_vibration_temperature = 10000.00; double freestream_vibration_temperature = 10000.00;
@ -290,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The total pressure of subsonic inlet.
// subsonicInletTotalTemperature: The total temperature of subsonic inlet.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The total pressure of subsonic outlet.
// subsonicOutletTotalTemperature: The total temperature of subsonic outlet.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -329,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -369,7 +369,7 @@ string str_limiter_name = "vanalbada";
// 4 -- Ustar limiter model, without grid size unitary. // 4 -- Ustar limiter model, without grid size unitary.
// 7 -- default used. // 7 -- default used.
// venkatCoeff: Cofficient of vencat, when using vencat limter. // venkatCoeff: Cofficient of vencat, when using vencat limter.
// limitVariables: Limit model. // limitVariables: Limit model (It is useful only if limitVector is 0).
// 0 -- limit only for pressure and denstiny, then get the min value. // 0 -- limit only for pressure and denstiny, then get the min value.
// 1 -- limit for every variables, then get the min value. // 1 -- limit for every variables, then get the min value.
// limitVector: // limitVector:
@ -386,16 +386,16 @@ string str_limiter_name = "vanalbada";
// roeEntropyScale: Entropy fix (correction) coefficient scale, default is 1.0. // roeEntropyScale: Entropy fix (correction) coefficient scale, default is 1.0.
// It is used to scale the default Roe entropy fix coefficients. // It is used to scale the default Roe entropy fix coefficients.
//int iviscous = 0; //int viscousType = 0;
//string viscousName = "Euler"; //string viscousName = "Euler";
int viscousType = 1; //int viscousType = 1;
string viscousName = "laminar"; //string viscousName = "laminar";
//int iviscous = 3; int viscousType = 3;
//string viscousName = "1eq-sa"; string viscousName = "1eq-sa";
//int iviscous = 4; //int viscousType = 4;
//string viscousName = "2eq-kw-menter-sst"; //string viscousName = "2eq-kw-menter-sst";
int DESType = 0; int DESType = 0;
@ -466,14 +466,14 @@ double roeEntropyScale = 1.0;
// RKStage: The number of Runge-Kutta step. // RKStage: The number of Runge-Kutta step.
// lamda: Cofficient of Runge-Kutta step. // lamda: Cofficient of Runge-Kutta step.
int iunsteady = 1; int iunsteady = 0;
double physicalTimeStep = 0.005; double physicalTimeStep = 0.01;
int ifStartFromSteadyResults = 0; int ifStartFromSteadyResults = 0;
int ifStaticsFlowField = 0; int ifStaticsFlowField = 0;
int startStatisticStep = 0; int startStatisticStep = 800000;
int min_sub_iter = 20; int min_sub_iter = 50;
int max_sub_iter = 20; int max_sub_iter = 50;
double tol_sub_iter = 0.01; double tol_sub_iter = 0.01;
int tscheme = 4; int tscheme = 4;
@ -502,10 +502,10 @@ int maxale = 10;
double dtau = 0.001; double dtau = 0.001;
double dtau_max = 1E-01; double dtau_max = 1E-01;
int wallFunctionType = 1; int wallFunctionType = 0;
int RKStage = 2; int RKStage = 2;
double [] = 0.5, 1.0; double lamda[] = 0.5, 1.0;
//int RKStage = 1; //int RKStage = 1;
//double lamda[] = 1.0; //double lamda[] = 1.0;
@ -535,7 +535,7 @@ double [] = 0.5, 1.0;
// visualfile: The visualization file path of flowfield , write data for every default (intervalStepPlot) steps. // visualfile: The visualization file path of flowfield , write data for every default (intervalStepPlot) steps.
// wall_aircoefile: The file path to save flowfield variables of wall, write data for every default steps. // wall_aircoefile: The file path to save flowfield variables of wall, write data for every default steps.
int numberOfGridGroups = 1; int numberOfGridGroups = 1;
string gridfile = "./grid/cylinder_Re3900_270w_split120.fts"; string gridfile = "./grid/rae2822_hybrid2d__4.fts";
int walldistMethod = 1; int walldistMethod = 1;
@ -549,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -561,13 +561,13 @@ string componentforcefile = "results/component_aircoef.dat";
string overset_gridfile = "iblank.ovs"; string overset_gridfile = "iblank.ovs";
int isPlotVolumeField = 1; int isPlotVolumeField = 0;
// visualfileType: The file type of visualfile. // visualfileType: The file type of visualfile.
// 0 -- Tecplot binary. // 0 -- Tecplot binary.
// 1 -- Tecplot ASCII. // 1 -- Tecplot ASCII.
int visualfileType = 0; int visualfileType = 1;
// visualSlice: The slice of tecflow. // visualSlice: The slice of tecflow.
// 0 -- Do not save slice data. // 0 -- Do not save slice data.
@ -588,17 +588,35 @@ double lowerPlotFieldBox[] = [0.0 0.0 0.0];
double upperPlotFieldBox[] = [1.0 1.0 1.0]; double upperPlotFieldBox[] = [1.0 1.0 1.0];
// nVisualVariables: Number of variables want to be dumped for tecplot visualization. // nVisualVariables: Number of variables want to be dumped for tecplot visualization.
// visualVariables : variable types dumped, listed as following: // visualVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6) // -- density(0), u(1), v(2), w(3), pressure(4), temperature(5), mach(6),
// -- viscosityLaminar(7), viscosityTurbulent(8), // -- viscosityLaminar(7), viscosityTurbulent(8),
// -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12), strain_rate(13), Q_criteria(14) // -- vorticity_x(9), vorticity_y(10), vorticity_z(11), vorticityMagnitude(12),
// -- Cp(15), timeStep(16), volume(17) // -- strain_rate(13), Q_criteria(14), Cp(15), timeStep(16), volume(17),
// -- modeledTKE(18),modeleddissipationrate(19), SSTF1(20), SSTF2(21) ,vibration temperature(33), electron temperature(34) // -- modeledTKE(18), modeleddissipationrate(19), SSTF1(20), SSTF2(21), vibration temperature(33), electron temperature(34),
// Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!. // -- gradientUx(41), gradientUy(42), gradientVx(43), gradientVy(44).
// Arriables order must from small to larger. // Important Warning: Array size of visualVariables MUST be equal to nVisualVariables!!!
int nVisualVariables = 10; // Variables order must from small to big.
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 11, 12, 15];; int nVisualVariables = 8;
int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// dumpStandardModel: Dump many standard model data.
// 1 -- Turbulent flat plate.
int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -620,51 +638,11 @@ double muoo = 1.0e-1;
double kwoo = 1.0; double kwoo = 1.0;
# maximum eddy viscosity (myt/my) max. # maximum eddy viscosity (myt/my) max.
double eddyViscosityLimit = 1.0e5; double eddyViscosityLimit = 1.0e10;
double sdilim = 1.0e20; double sdilim = 1.0e20;
double coef_kvist = 1.0; double coef_kvist = 1.0;
int monitor_vistmax = 0; int monitor_vistmax = 0;
//-----------------------------------------------------------------------
# LES Parameter #
//-----------------------------------------------------------------------
// iLES : Create LESSolver or not.
// >= 1 - Create LESSolver;
// < 1 - not;
// sgsmodel : subgrid scale model
// = "smagorinsky" ;
// = "dsm" ;
// = "wale" ;
// deltaFunctionType : = 1 - MAX(deltai, deltaj, deltak)
// = 2 - pow(deltai * deltaj *deltak, 1/3)
// = 3 - Devloped by Scotti
// wallDampingFunctionType : = 0 - no wall function
// = 1 - van Driest
// = 2 - developed by Dr. Deng Xiaobing
// = 3 - developed by Piomelli
// turbViscousCutType : turbulent viscosity cut type
// = 0 - mu_total = mut + mul
// = 1 - mu_total = max(mut-mul,0)+ mul
// = 2 - mu_total = max(mut ,0)+ mul
// smagConstant : constant of smagorinsky model.
// waleConstant : constant of wale model
// filterDirection[3] : filter variables in i, j, k direction or not;
// averageDirection[3] : average variables in i, j, k direction or not;
// isotropicConstant : constant of isotropic part of SGS stress
int iLES = 1;
string sgsmodel = "smagorinsky";
int deltaFunctionType = 2;
int wallDampingFunctionType = 1;
int turbViscousCutType = 2;
double smagConstant = 0.135;
double isotropicConstant = 0.0;
double waleConstant = 0.6;
int filterDirection[] = [1, 1, 0];
int averageDirection[] = [1, 1, 0];
double testFilterScale = 2.0;
int averageWidth = 1;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Other Parameter # # Other Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -805,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.
@ -822,7 +803,7 @@ string structhighordergradient = "conservation";
// 0 -- the effect of the limiter is cancelled, means the first-order interpolations. // 0 -- the effect of the limiter is cancelled, means the first-order interpolations.
// allReduceStep: Iteration intervals for MPI AllReduce operation, default is 1. // allReduceStep: Iteration intervals for MPI AllReduce operation, default is 1.
string outtimesc = "impbd2"; string outtimesc = "impbd2";
double MUSCLCoefXk = 0.8; double MUSCLCoefXk = -1;
double MUSCLCoefXb = 1.0; double MUSCLCoefXb = 1.0;
int allReduceStep = 1; int allReduceStep = 1;

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// viscousType : Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// viscousType : Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// viscousType : Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// viscousType : Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.

View File

@ -146,6 +146,16 @@ int numberOfReferenceCP = 10;
double influenceRadius = 20; double influenceRadius = 20;
int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0; int symmetryPlane = 3; // 1: plane of x=0; 2: plane of y=0; 3: plane of z=0;
// ----------------- Periodic Parameters --------------------------------
// Notice:Rotational periodicity only support rotation along the X axis!
// periodicType: Which periodic boundary is used.
// 0 -- without Periodic Boundary.
// 1 -- Translational periodicity.
// 2 -- Rotational periodicity.
int periodicType = 0;
double translationLength[] = [0.0,0.0,0.0];
double rotationAngle = 0.0;
######################################################################### #########################################################################
# Default parameters for Partition # # Default parameters for Partition #
######################################################################### #########################################################################
@ -205,6 +215,7 @@ int numberOfMultigrid = 1;
// maxSimuStep: The max simulation step, don't care simulation is restart or not. // maxSimuStep: The max simulation step, don't care simulation is restart or not.
// intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved. // intervalStepFlow: The step intervals for flow variables file 'flow.dat' saved.
// intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved. // intervalStepPlot: The step intervals for tecplot visual file 'tecflow.dat' saved.
// intervalStepSample: The step intervals for monitored probes variables file 'sample.dat' saved.
// intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved. // intervalStepForce: The step intervals for aerodynamics coefficients file 'aircoef.dat' saved.
// intervalStepRes: The step intervals for residual file 'res.dat' saved. // intervalStepRes: The step intervals for residual file 'res.dat' saved.
// ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow. // ifLowSpeedPrecon: Precondition process to accelerate convergence for low speed flow.
@ -215,6 +226,7 @@ int maxSimuStep = 20000;
int intervalStepFlow = 1000; int intervalStepFlow = 1000;
int intervalStepPlot = 1000; int intervalStepPlot = 1000;
int intervalStepSample = 1000;
int intervalStepForce = 100; int intervalStepForce = 100;
int intervalStepRes = 10; int intervalStepRes = 10;
int ifLowSpeedPrecon = 0; int ifLowSpeedPrecon = 0;
@ -296,24 +308,6 @@ double TorqueRefX = 0.0; // unit of meter.
double TorqueRefY = 0.0; // unit of meter. double TorqueRefY = 0.0; // unit of meter.
double TorqueRefZ = 0.0; // unit of meter. double TorqueRefZ = 0.0; // unit of meter.
// directionMethod: The method of determining direction.
// 1 -- using direciton.
// 2 -- using face normal.
// subsonicInletTotalPressure: The dimensional total pressure of subsonic inlet if height > 0, else non-dimensional.
// subsonicInletTotalTemperature: The dimensional total temperature of subsonic inlet if height > 0, else non-dimensional.
// direction_inlet: The direction of inlet.
// subsonicOutletTotalPressure: The dimensional total pressure of subsonic outlet if height > 0, else non-dimensional.
// subsonicOutletTotalTemperature: The dimensional total temperature of subsonic outlet if height > 0, else non-dimensional.
// direction_outlet: The direction of outlet.
int directionMethod = 2;
double subsonicInletTotalPressure = 1.2e6;
double subsonicInletTotalTemperature = 1300;
double direction_inlet[] = 1, 0, 0;
double subsonicOutletTotalPressure = 17.8571428;
double subsonicOutletTotalTemperature = 1.0;
double direction_outlet[] = 1, 0, 0;
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Spatial Discretisation # # Spatial Discretisation #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -335,7 +329,7 @@ string str_limiter_name = "vanalbada";
#******************************************************************* #*******************************************************************
# UnStruct Solver or Common * # UnStruct Solver or Common *
#******************************************************************* #*******************************************************************
// iviscous: Viscous model. // viscousType: Viscous model.
// 0 -- Euler. // 0 -- Euler.
// 1 -- Lamilar. // 1 -- Lamilar.
// 2 -- Algebraic. // 2 -- Algebraic.
@ -555,7 +549,7 @@ string turbfile = "results/turb.dat";
string visualfile = "results/tecflow.plt"; string visualfile = "results/tecflow.plt";
string Qwall_file = "results/Qwall.dat"; string Qwall_file = "results/Qwall.dat";
string wall_aircoefile = "results/wall_aircoef.dat"; string wall_aircoefile = "results/wall_aircoef.dat";
string probesflowfile = "results/sample.dat";
string surfacefile = ""; string surfacefile = "";
string wall_varfile = ""; string wall_varfile = "";
@ -610,6 +604,19 @@ int visualVariables[] = [0, 1, 2, 3, 4, 5, 6, 15];
// 1 -- Turbulent flat plate. // 1 -- Turbulent flat plate.
int dumpStandardModel = 0; int dumpStandardModel = 0;
// ifSetProbesToMonitor: Set probes location to Monitor.
// 0 -- Do not monitor.
// 1 -- To monitor and save data to files.
// probesDefineFile: probes location information file.
// nProbeVariables: Number of variables want to be dumped for probes monitered.
// probeVariables : Variable types dumped, listed as following:
// -- density(0), u(1), v(2), w(3), pressure(4), temperature(5).
// Important Warning: Array size of probeVariables MUST be equal to nProbeVariables!!!
// probeVariables order must from small to big.
int ifSetProbesToMonitor = 0;
string probesDefineFile = "bin/probes_location.hypara";
int nProbeVariables = 6;
int probeVariables[] = [0, 1, 2, 3, 4, 5];
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
# Turbulence Parameter # # Turbulence Parameter #
//----------------------------------------------------------------------- //-----------------------------------------------------------------------
@ -776,11 +783,14 @@ string zoneInverseFileName = "./grid/zoneInverseMapping.inp";
// structhighordergradient: // structhighordergradient:
// -- "conservation", "chain_rule". // -- "conservation", "chain_rule".
int ifvfd = 0; int ifvfd = 0;
string str_highorder_solver = "WCNS";
int SolverStructOrder = 0; int SolverStructOrder = 0;
double str_highorder_interpolation_epsilon = 1.0e-6; double str_highorder_interpolation_epsilon = 1.0e-6;
string str_highorder_interpolation_type = "test"; string str_highorder_interpolation_type = "test";
string str_highorder_flux_name = "steger"; string str_highorder_flux_name = "steger";
string structhighordergradient = "conservation"; string structhighordergradient = "conservation";
double coefofstrflux = 0.5;
double limitcoefofinterface = 0.0;
// ----------------- Advanced choices ----------------------------------- // ----------------- Advanced choices -----------------------------------
// outtimesc: Time stepping scheme for the outer loop. // outtimesc: Time stepping scheme for the outer loop.