Optional Parameters

Optional parameters can be specified in a list at the end of the input file.
 
This is done by specifying the name of the parameter, and its value on the same line.
 
The value must be in floating point decimal form, and must line up with the header: VALUE
 
The listed parameters and their values are output at the start of the run, providing the user with a check of the values read in.
 
The following is a list of all of the available optional parameters sorted by function with the typical values used for each. Click on the parameter names for complete descriptions. An alphabetized list is also included below for convenience.
Hover over a parameter for a short definition or CLICK on it to be taken to it's listing below.

OPTIONAL PARAMETERS BY FUNCTION:

DES: Detached Eddy Simulation

DISCRETIZATION TYPE (CONVECTIVE)
IFLUX_TYPEIFLUX_TYPE: Specifies the type of flux or dissipation for the spatial discretization scheme.
(click for more)
..............:(=0,1,2)
IGRAD_TYPEIGRAD_TYPE: Defines the values used in gradient reconstruction with IFLUX_TYPE=1,2 (IGRAD_TYPE is inactive for IFLUX_TYPE=0)
(click for more)
..............:(=1,2)
IFREEZE_GRADIFREEZE_GRAD: Used to enable freezing of gradient calculation within the stages of the Runge-Kutta multi-stage scheme. This speeds up execution, but may be less stable. Gradients are recomputed at each time-step in the first stage of the Rung-Kutta scheme, so the final converged result will be identical regardless of the value of IFREEZE_GRAD
(click for more)
............:(=0,1)
ILIM_TYPEILIM_TYPE: Select Type of Limiter for Upwind Reconstruction Schemes (IFLUX_TYPE=1,2...Not active for IFLUX_TYPE=0)
(click for more)
...............:(=-1,0,1,2)
FK_LIMITFK_LIMIT: Parameter in Venkataskrishnan TVB limiter. Low values tend to weaker enforcement of monotonicity, high values tend toward original monotone limiter.
(click for more)
................:(0.0 < 100.0)
IFREEZE_LIMIFREEZE_LIM: Specifies if and how to freeze limiters (for IFLUX_TYPE=1 or 2 only, inactive for IFLUX_TYPE=0) Convergence to steady-state can be disrupted by limiters which switch back and forth at each iteration. Freezing these limiters may enable or enhance otherwise poor convergence, but the final steady-state result may depend on the manner in which the limiters have been frozen and the convergence history (i.e. a restart solution may converge to slightly different results). In general these differences should be acceptably small.
(click for more)
.............:(=1,2,3)
NFREEZE_LIMNFREEZE_LIM: Controls limiter freezing. Freeze limiters (IFREEZE_LIM=1) or begin freezing limiters (IFREEZE_LIM=2,3) after NFREEZE_LIM cycles. Only active for IFREEZE_LIM > 0
(click for more)
.............:(-1 &.< NCYC)
BIHFACTRBIHFACTR: is a factor which scales the linearization of the biharmonic dissipation terms in the Jacobian.
(click for more)
................:(=20.0)
BIH_BNDY_FACTRBIH_BNDY_FACTR: Specifies the solid wall boundary treatment of artificial dissipation (Active for IFLUX_TYPE=0 only). Second differences at the wall are scaled by BIH_BNDY_FACTR prior to calculating the full biharmonic (4th difference) dissipation terms.
(click for more)
..........:(0.0 < 1.0)
IFINE_BNDY_DISSIPIFINE_BNDY_DISSIP: Controls modification of boundary dissipation for schemes other than biharmonic dissipation (IFLUX_TYPE=0) IFINE_BNDY_DISSIP = -1 : Zero out dissipation for tangential boundary condition points (see innerbc_0.f)
(click for more)
.......:(=-1,0,1,2)

DISCRETIZATION TYPE (VISCOUS TERMS)
DISCRETIZATION TYPE (TURBULENCE MODEL)

ENTROPY FIXES

TIME STEP LIMITING FOR ROBUSTNESS
IMESSAGE_LEVELIMESSAGE_LEVEL: Select Level of output messages when time-step limiting or other limiting occurs.
(click for more)
..........:(=0,1,2)
ISAFE_LEVELISAFE_LEVEL: Controls level of checking for unphysical states in time stepping procedure.
(click for more)
.............:(=0,1,2,3,4)
ISAFE_MGISAFE_MG: ISAFE_MG ... (=0,1,2)
(click for more)
................:(=0,1,2)
FACTR_MGFACTR_MG: Can be used to enhance robustness of multigrid by decreasing size of transfered residuals to coarse grids. And in turn increasing corrections back to fine grid. Overall convergence history should be identical, but non-linear instabilities on coarse grids may be avoided.
(click for more)
................:(1.0 < 10.0)
IGEOM_LIMITIGEOM_LIMIT: Option to use (1) or Discard (0) Grid Based Limiting
(click for more)
.............:(=0,1)
FGEOM_FRACFGEOM_FRAC: Sets Maximum fraction of limit grid points. When limiting based on grid cell geometry is established (prior to flow calculations), if more than FGEOM_FRAC * nnode grid points are found to be limited (nnode = total number of grid points), execution is halted.
(click for more)
..............:(0.0 < 1.0)
PMINPMIN: Minimum pressure value before limiting occurs.
....................:(0.0 < 1.0)
PMIN_MGPMIN_MG: Minimum pressure value for omitting multigrid updates
.................:(0.0 < 1.0)
RHOMINRHOMIN: Minimum density value before limiting occurs.
..................:(0.0 < 1.0)
RHOMIN_MGRHOMIN_MG: Minimum densuty value for omitting multigrid updates.
...............:(0.0 < 1.0)
 
LOW MACH NUMBER PRECONDITIONING
 
INPUT/OUTPUT CONTROL PARAMETERS
 
RESTART FILE AND CHECKPOINT FILE PARAMETERS
 
IBL SPECIFIC PARAMETERS (ACTIVE ONLY WHEN RUNNING INTERACTIVE BOUNDARY LAYER)
 

OPTIONAL PARAMETERS BY NAME:

Click on a parameter to toggle a full description and additional information.

BETA_MIN ..................: (0.0 < 1.0)
BETA_MIN: relates exclusively to low Mach number preconditioning.

This is an optional parameter and is only active for IPC_LOW_MACH = -1 In this case, BETA_MIN may take on values between 0 and 1. Generally (i.e. IPC_LOW_MACH = +1) low Mach number preconditioning uses the value BETA_MIN = min(1,3 x Mach**2) by default. Lower values of BETA_MIN may be less stable while providing more low Mach number preconditioning effect. Higher numbers produce the opposite, with the effect of low-Mach number preconditioning vanishes for BETA_MIN=1.

BIH_BNDY_FACTR ............: (0.0 < 1.0)
BIH_BNDY_FACTR: specifies the solid wall boundary treatment of artificial dissipation (Active for IFLUX_TYPE=0 only). Second differences at the wall are scaled by BIH_BNDY_FACTR prior to calculating the full biharmonic (4th difference) dissipation terms. The two extreme values correspond to:

BIH_BNDY_FACTR = 0: Omit second differences on wall in construction of biharmonic dissipation (Better skin friction, possibly lower robustness)

BIH_BNDY_FACTR = 1: Do not modify second differences at wall.

BIHFACTR ..................: (=20.0)
BIHFACTR: is a factor which scales the linearization of the biharmonic dissipation terms in the Jacobian.

Generally, a 1st order Jacobian is used for the point or line implicit solution algorithm.

When using a biharmonic (artificial) dissipation scheme, there is no simple correspondance between the linearization of these terms and a 1st order (2nd difference dissipation) Jacobian. Therefore, we use the nearest neighbor entries from the biharmonic construction, and scale these by the factor BIHFACTR.
BIHFACTR = 20.0 : Standard Value. Higher value produces more diagonally dominant matrix (slower and more robust convergence) while lower values produce the opposite.

C_DES .....................: (0 < C_DES <= 1) default 0.65
C_DES: parameter for DES length scale, set to 0.65 according to DES97

CP_FACTOR .................: (0.0 < 1.0)
CP_FACTOR: For non-converging boundary layer code, scale Cps by factor CP_FACTOR and rerun IBL station.

Not active on Nash IBL code.

Recommended value: 0.5

EFACTRC1 ..................: (0.0 < 1.0)
EFACTRC1: Value of entropy fix for COARSE GRID (in multigrid) Discretization Limit value of minimum eigenvalue by EFACTRC1 of maximum eigenvalue.

EFACTRC1 = 0.0 : No limiting
EFACTRC1 = 1.0 : Scalar Dissipation Scheme (slower, more robust) Affects only convergence, not final solution (fine grid discretization)

EFACTRC2 ..................: (0.0 < 1.0)
EFACTRC2: Value of entropy fix for COARSE GRID (in multigrid) Jacobian Limit value of minimum eigenvalue by EFACTRC2 of maximum eigenvalue.

EFACTRC2 = 0.0 : No limiting
EFACTRC2 = 1.0 : Scalar Dissipation Scheme (slower convergence, more robust) Affects only convergence, not final solution (fine grid discretization)

EFACTRF1 ..................: (0.0 < 1.0)
EFACTRF1: Value of entropy fix for FINE GRID Discretization Limit value of minimum eigenvalue by EFACTRF1 of maximum eigenvalue.

EFACTRF1 = 0.0 : No limiting
EFACTRF1 = 1.0 : Scalar Dissipation Scheme (slower convergence, more diffusion, more robust)
Affects convergence, AND final solution (fine grid discretization)

EFACTRF2 ..................: (0.0 < 1.0)
EFACTRF2: Value of entropy fix for FINE GRID Jacobian Limit value of minimum eigenvalue by EFACTRF2 of maximum eigenvalue.

EFACTRF2 = 0.0 : No limiting
EFACTRF2 = 1.0 : Scalar Dissipation Scheme (slower convergence, more robust) Affects only convergence, not final solution (fine grid Jacobian only)

FACTR_MG ..................: (1.0 < 10.0)
FACTR_MG: Can be used to enhance robustness of multigrid by decreasing size of transfered residuals to coarse grids. And in turn increasing corrections back to fine grid. Overall convergence history should be identical, but non-linear instabilities on coarse grids may be avoided.

NOTE: FACTR_MG=1 recommended. Other values seldom used.
FACTR_MG > 1: Rescale for robustness.
FACTR_MG = 1: Baseline Value, no rescaling of multigrid terms

FD_DES.....................: (< 0,0,1) default 1
FD_DES: parameter for selecting DES or DDES
FD_DES = 1 will run DES97 (default)
FD_DES = -1 will run DDES2006
FD_DES = 0 will run the RANS

FGEOM_FRAC ................: (0.0 < 1.0)
FGEOM_FRAC: Sets Maximum fraction of limit grid points. When limiting based on grid cell geometry is established (prior to flow calculations), if more than FGEOM_FRAC * nnode grid points are found to be limited (nnode = total number of grid points), execution is halted.

FGEOM_FRAC > 0 : Use values between 0 and 1.0
FGEOM_FRAC > 1.0 : Limit is Inactive

FK_LIMIT ..................: (0.0 < 100.0)
FK_LIMIT: Parameter in Venkataskrishnan TVB limiter. Low values tend to weaker enforcement of monotonicity, high values tend toward original monotone limiter.

FK_LIMIT = 0.0: Removes effect of limiter, acts as unlimited discretization
LIMIT_inf -> Monotone : Limiter FK
Typical Values:
5.0 < FK_LIMIT < 100.0
Only active for : IFLUX_TYPE>0, and ILIM_TYPE = 1

FNSFACTR ..................: (=0.0,1.0)
FNSFACTR: Controls Navier-Stokes Discretization.

FNSFACTR = 0.0 : Use thin Layer Navier-Stokes using a Pseudo- Laplacian Operator Based on Mesh Edges This is the prefered option. The result is actually thin layer in all three directions, or full Navier-Stokes under the assumption of incompressibility and constant viscosity. This is also the most robust option.

FNSFACTR = 1.0 : Use Full Navier-Stokes terms using the above Pseudo-Laplacian complemented by second differences computed as gradients of gradients for corrections to Laplacian and inclusion of cross terms. (involves stencils of neighbors of neighbors) This may be less robust as exact Jacobians of these terms are not available.

FNSFACTT ..................: (=0.0,1.0)
FNSFACTT: Controls Turbulence model dissipation discretization.

FNSFACTT = 0.0 : Use a Pseudo- Laplacian Operator Based on Mesh Edges for Dissipation terms. This is the prefered option.

FNSFACTT = 1.0 : Correct the above Pseudo- Laplacian Operator Based on Mesh Edges for Dissipation terms using second derivatives computed as gradients of gradients (involves stencils of neighbors of neighbors) This may be less robust as exact Jacobians of these terms are not available.

Click on a parameter to toggle a full description and additional information.

FNSGRAD ...................: (=1.0,2.0)
FNSGRAD: Gradient construction options for Navier-Stokes Terms (only active when FNSFACTR = 1.0)

FNSGRAD = 1.0 : Use Green-Gauss Construction for Navier-Stokes Terms when FNSFACTR = 1.0.
FNSGRAD = 2.0 : Use Least-Squares Cosntruction for Navier-Stokes Terms when FNSFACTR = 1.0

FNSGRAT ...................: (=1.0,2.0)
FNSGRAT: Gradient construction options for Turbulence Model Dissipative Terms (only active when FNSFACTT = 1.0)

FNSGRAT = 1.0 : Use Green-Gauss Construction for Turbulence Model Dissipative Terms when FNSFACTT = 1.0

FNSGRAT = 2.0 : Use Least-Squares Cosntruction for Turbulence Model Dissipative Terms when FNSFACTT = 1.0

ICHK_PARALLEL .............: (=1,2)
ICHK_PARALLEL: Enable/Disable Parallel output for Checkpoint Files Checkpoint files are written at regular intervals. Parallel output is implemented as all processors simultaneously writing to a common file system. On some systems, this causes the Network File System (NFS) to overload. To avoid this, ICHK_PARALLEL=0 instructs each processor to write to the common file system in sequence, one at a time, while the other processors wait their turn. This is slower but more robust for many NFS systems.

ICHK_PARALLEL = 0 : Sequential output to common file system using I/O from one processor at a time.
ICHK_PARALLEL = 1 : Parallel output to common file system using I/O from all processors simultaneously.

IFINE_BNDY_DISSIP .........: (=-1,0,1,2)
IFINE_BNDY_DISSIP: Controls modification of boundary dissipation for schemes other than biharmonic dissipation (IFLUX_TYPE=0) IFINE_BNDY_DISSIP = -1 : Zero out dissipation for tangential boundary condition points (see innerbc_0.f)

IFINE_BNDY_DISSIP = 0 : No change, as expected on all grid levels
IFINE_BNDY_DISSIP = 1 : Use coarse grid dissipation values at ghost edge in rfluxb_0.f rfluxb_1.f
IFINE_BNDY_DISSIP = 2 : Use LIM_GEOM values at all boundary edges:
eg. values used for LIM_GEOM points (values set in set_lim_values.f)
BUT this OPTION requires nsmoo_limit(LIM_GEOM) = 0.

IFLUX_TYPE ................: (=0,1,2)
IFLUX_TYPE: Specifies the type of flux or dissipation for the spatial discretization scheme. Currently 3 options are implemented:

IFLUX_TYPE = 0 : Original biharmonic matrix dissipation

IFLUX_TYPE = 1 : Roe Rieman Flux Difference Splitting

IFLUX_TYPE = 2 : Van-Leer Flux Vector Splitting

For Navier-Stokes flows, IFLUX_TYPE = 0 is recommended, as there are still issues of accuracy with the other schemes.

For Euler (inviscid) flows, IFLUX_TYPE = 1 is relatively robust for supersonic flows especially when using a limiter (see ILIM_TYPE).

IFREEZE_LIM ...............: (=1,2,3)
IFREEZE_LIM: Specifies if and how to freeze limiters (for IFLUX_TYPE=1 or 2 only, inactive for IFLUX_TYPE=0) Convergence to steady-state can be disrupted by limiters which switch back and forth at each iteration. Freezing these limiters may enable or enhance otherwise poor convergence, but the final steady-state result may depend on the manner in which the limiters have been frozen and the convergence history (i.e. a restart solution may converge to slightly different results). In general these differences should be acceptably small.

IFREEZE_LIM = 0 : Never Freeze limiters
IFREEZE_LIM = 1 : Freeze all limiters at current values after NFREEZE_LIM iterations

IFREEZE_LIM = 2 : After NFREEZE_LIM iterations, begin a moving average of limiter values

IFREEZE_LIM = 3 : After NFREEZE_LIM iterations, use minimum limiter value of previous and next iteration (produces minimum value over all remaining iterations)

IFREEZE_GRAD ..............: (=0,1)
IFREEZE_GRAD: Used to enable freezing of gradient calculation within the stages of the Runge-Kutta multi-stage scheme. This speeds up execution, but may be less stable. Gradients are recomputed at each time-step in the first stage of the Rung-Kutta scheme, so the final converged result will be identical regardless of the value of IFREEZE_GRAD

IFREEZE_GRAD = 0 : Do not freeze any gradient calculations.

IFREEZE_GRAD = 1 : Freeze at all Runge-Kutta stages after 1st stage (but recompute at start of each new time step).

IGRAD_TYPE ................: (=1,2)
IGRAD_TYPE: Defines the values used in gradient reconstruction with IFLUX_TYPE=1,2 (IGRAD_TYPE is inactive for IFLUX_TYPE=0)

IGRAD_TYPE = 1 : Gradients are computed using Primitive Variables

IGRAD_TYPE = 2 : Gradients are computed using Conserved Variables

IGEOM_LIMIT ...............: (=0,1)
IGEOM_LIMIT: Option to use (1) or Discard (0) Grid Based Limiting

IGEOM_LIMIT = 0 : Discard

IGEOM_LIMIT = 1 : Enforce

(Recommend =0 Discard for this version)

ILIM_TYPE .................: (= -1,0,1,2)
ILIM_TYPE: Select Type of Limiter for Upwind Reconstruction Schemes (IFLUX_TYPE=1,2...Not active for IFLUX_TYPE=0)
ILIM_TYPE < 0 : First Order Discretization: Set All Gradients = 0.0

ILIM_TYPE = 0 : No Limiting, allow non-monotone solutions (least diffusive)

ILIM_TYPE = 1 : Venkatakrishnan Smooth (TVB) Limiter (less diffusive) --> Also select FK_LIMIT value

ILIM_TYPE = 2 : Barth Monotone Limiter (most diffusive)

ILINE_SOLVE ...............: (=0,1)
ILINE_SOLVE: Select / Omit Line Solver. Line solver increases convergence for Navier-Stokes flows. Robustness problems have been encountered for IFLUX_TYPE=1,2. For IFLUX_TYPE= 0, line solver should always be used.

ILINE_SOLVE = 0 : Omit Line Solver. Use Point Jacobi only.

ILINE_SOLVE = 1 : Use Line Solver in Boundary Layer Regions.

IMESSAGE_LEVEL ............: (=0,1,2)
IMESSAGE_LEVEL: Select Level of output messages when time-step limiting or other limiting occurs.

IMESSAGE_LEVEL=0 : No output messages

IMESSAGE_LEVEL=1 : Moderate output Messages

IMESSAGE_LEVEL=2 : Full Output Message (recommended)

IN_PARALLEL ...............: (=1,2)
IN_PARALLEL: Enable/Disable Parallel Input for Reading Grid Files. Parallel input is implemented as all processors simultaneously reading from a common file system. On some systems, this causes the (Network File System) NFS to overload. To avoid this, IN_PARALLEL=0 instructs each processor to read from the common file system in sequence, one at a time, while the other processors wait their turn. This is slower but more robust for many NFS systems.

IN_PARALLEL = 0 : Sequential input from common file system using I/O from one processor at a time.

IN_PARALLEL = 1 : Parallel input from common file system using I/O from all processors simultaneously.

Click on a parameter to toggle a full description and additional information.

INRES_PARALLEL ............: (=1,2)
INRES_PARALLEL: Enable/Disable Parallel Input for Reading Restart Files. Parallel input is implemented as all processors simultaneously reading from a common file system. On some systems, this causes the (Network File System) NFS to overload. To avoid this, INRES_PARALLEL=0 instructs each processor to read from the common file system in sequence, one at a time, while the other processors wait their turn. This is slower but more robust for many NFS systems.

INRES_PARALLEL = 0 : Sequential input from common file system using I/O from one processor at a time.

INRES_PARALLEL = 1 : Parallel input from common file system using I/O from all processors simultaneously.

IOUT_PARALLEL .............: (=1,2)
IOUT_PARALLEL: Enable/Disable Parallel Output for Writing out Restart Files. Parallel output is implemented as all processors simultaneously writing to a common file system. On some systems, this causes the (Network File System) NFS to overload. To avoid this, IOUT_PARALLEL=0 instructs each processor to write to the common file system in sequence, one at a time, while the other processors wait their turn. This is slower but more robust for many NFS systems.

IOUT_PARALLEL = 0 : Sequential output to common file system using I/O from one processor at a time.

IOUT_PARALLEL = 1 : Parallel output to common file system using I/O from all processors simultaneously.

IPC_LOW_MACH ..............: (=-1,0,1)
IPC_LOW_MACH: Select or Omit Low Mach Number Preconditioning

IPC_LOW_MACH = 0 : No Preconditioning, use regular base scheme

IPC_LOW_MACH = 1 : Standard Low-Mach Preconditioning (use preset value BETA_MIN=3Mach**2)

IPC_LOW_MACH = -1 : Custome Low-Mach Preconditioning: Also Set value for : BETA_MIN

IPC_RAMP ..................: (0 < NCYC)
IPC_RAMP: Ability to Ramp in the Low_Mach number preconditioning for increased robustness.

IPC_RAMP = 0 : No ramping, apply full preconditioning from first iteration.

IPC_RAMP = 100 : Ramp in low-Mach number preconditioning over first 100 cycles and apply full preconditioning thereafter.

IRESTART_AUX_TYPE .........: (=0,1,2,3)
IRESTART_AUX_TYPE: Determines content of restart.aux.out auxiliary restart file (or directory with partitioned data)

IRESTART_AUX_TYPE = 0 : No auxiliary restart file
IRESTART_AUX_TYPE = 1 : Cf,Y+,on surface, eddy viscosity elsewhere

IRESTART_AUX_TYPE = 2 : Oil Flow surface velocities(1st 3 entries), Cf, Y+ on surface

IRESTART_AUX_TYPE = 3 : Transition Mask, last iteration number for limited time step at given grid points

IRESTART_AUX_TYPE = 5 : Pressure Coefficient and Skin Friction vector (Cp,CF,CFx,CFy,CFz)

ISAFE_LEVEL ...............: (=0,1,2,3,4)
ISAFE_LEVEL: Controls level of checking for unphysical states in time stepping procedure.

ISAFE_LEVEL = 0 : No Checks for Negative Density/Pressures (slightly faster execution)

ISAFE_LEVEL = 1 : Checks for Negative Density/Pressures (reduces time step accordingly)

ISAFE_LEVEL = 2 : After above fixes, if negative values still occur omit update at these points

ISAFE_LEVEL = 3 : Perturb Energy Values for STUCK pts (where dt -> 0.0)

ISAFE_LEVEL = 4 : More Robust itbc(but O(h) flux through walls)

(requires IFINE_BNDY_DISSIP .ne. -1)

Recommended Value: ISAFE_LEVEL = 2

ISAFE_MG ..................: (=0,1,2)
ISAFE_MG:

ISAFE_MG = 0 : Regular collc and addw Multigrid transfer routines (i.e. no change to base algorithm)

ISAFE_MG = 1 : Omit transfer or residuals (collc) to coarse grids at (time-step) limited points

ISAFE_MG = 2 : Omit transfer of residuals (collc) AND prolongation back to fine grid (addw) at (time-step) limited points

ITE_IBL ...................: (= -2,-1,0,1,2,3)
ITE_IBL: Trailing edge treatment for IBL code

ITE_IBL = 0 : No special Trailing Edge Treatment

ITE_IBL = 1 : Extrapolate Velocities at x>xclip_ibl (%chord)

ITE_IBL = -1 : Smooth Velocities at x>xclip_ibl (%chord)

ITE_IBL = 2 : Extrapolate Velocities at last NTE_IBL Points

ITE_IBL = -2 : Smooth Velocities at last NTE_IBL Points

ITE_IBL = 3 : IBL not computed on last NTE_IBL Points and blowing velocities extrapolated

NCYC_CHECKPT ..............: (0 < NCYC)
NCYC_CHECKPT: Write out checkpoint file aftyer NCYC_CHECKPT cycles. Alternate checkpt files are written out: checkpt.1 every odd multiple of NCYC_CHECKPT and checkpt.2 every even multiple of NCYC_CHECKPT. Rather than overwriting the latest checkpt file, this approach avoids possible loss of latest checkpt file if failure occurs during checkpt write.

NCYC_CHECKPT < 0 or = 0 : Omit Checkpoint files

NFREEZE_LIM ...............: (-1 < NCYC)
NFREEZE_LIM: Controls limiter freezing. Freeze limiters (IFREEZE_LIM=1) or begin freezing limiters (IFREEZE_LIM=2,3) after NFREEZE_LIM cycles. Only active for IFREEZE_LIM > 0

Click on a parameter to toggle a full description and additional information.

NTE_IBL ...................: (0 < NPTS)
NTE_IBL: Number of pts when ITE_IBL = 2,3

PMIN ......................: (0.0 < 1.0)
PMIN: Minimum pressure value before limiting occurs.

PMIN_BCWALL_IBL ...........: (0.0 < 1.0)
PMIN_BCWALL_IBL: Limit on Minimum Wall Pressure at blowing points

PMIN_MG ...................: (0.0 < 1.0)
PMIN_MG: Minimum pressure value for omitting multigrid updates

RHOMIN ....................: (0.0 < 1.0)
RHOMIN: Minimum density value before limiting occurs.

RHOMIN_MG .................: (0.0 < 1.0)
RHOMIN_MG: Minimum densuty value for omitting multigrid updates.

SIGMA_DES..................: (0 < SIGMA_DES <= 1) default 1
SIGMA DES: parameter for the scaling of the artificial damping.

  • 0 < SIGMA_DES < 1 is for the fixed scaling value.
  • SIGMA_DES = 1 will run the hybrid scheme(equation 6). However, the parameters in equation (6), (7), (8) and (9) need to further investigation to ensue the robustness of the scheme. This option is not recommended right now.

VNMAX_IBL .................: (0.0 < 1.0)
VNMAX_IBL: Maximum Limit on Blowing Velocities

XCLIP_IBL .................: (0.0 < 1.0)
XCLIP_IBL: Fraction of chord after which extrapolation is done when ITE_IBL = 1 (i.e. XCLIP_IBL = 0.99)