From bc500c88fc3aa4940a5f91f8d9ec341e0d8caff9 Mon Sep 17 00:00:00 2001 From: Mengmeng Zhang Date: Tue, 9 Jan 2024 18:54:54 +0100 Subject: [PATCH] updated the default.ainp.tmp file to the latest version, which allows the implicit scheme more efficient than the old multigrid scheme. Usage: NGRID=1 (multigrid not needed) STSTCONV=3 (turn on implicit) KSP_TOL, KSP_TOL_TURB=0.01 (convergence of the linear solver) --- ceasiompy/EdgeRun/files/default.ainp.tmp | 3943 +++++++++++++++------- 1 file changed, 2668 insertions(+), 1275 deletions(-) mode change 100755 => 100644 ceasiompy/EdgeRun/files/default.ainp.tmp diff --git a/ceasiompy/EdgeRun/files/default.ainp.tmp b/ceasiompy/EdgeRun/files/default.ainp.tmp old mode 100755 new mode 100644 index 979ac7f8e..30abac877 --- a/ceasiompy/EdgeRun/files/default.ainp.tmp +++ b/ceasiompy/EdgeRun/files/default.ainp.tmp @@ -1,11 +1,11 @@ -DEFAULT,N ,0,0,165 +DEFAULT,N ,0,0,216 * *************************************************************************** * * * Edge input parameter file * * This file uses the so called FFA-format * * * -MODIFICATIONS,L ,61,1,0 +MODIFICATIONS,L ,73,1,0 'Date Version Programmer Description ' '19991214 P.Eliasson Change def. VIS4,VIS2 ' '20010503 S.Wallin Introduced TURB_MOD_LIST ' @@ -39,7 +39,7 @@ MODIFICATIONS,L ,61,1,0 ' adapted to FBCWAL.F ' '20030930 3.2 S.-H. Peng Added Wilcoxs LRN model ' '20030930 3.2 P.Eliasson Variable ITHEDA removed ' -'20031110 3.2 M.Franke Added Rung LLR-model ' +'20031110 3.2 M.Franke Added Rung LLR-model (2D) ' '20040102 3.2 J.Smith PERT_CONST & IMOVGR inputs removed ' '20040210 3.2 S.-H. Peng Added S-A DES model ' '20040220 3.2 S.-H. Peng Added Smagorinsky LES model ' @@ -58,6 +58,12 @@ MODIFICATIONS,L ,61,1,0 '20040922 3.3 O.Amoignon Adjoint solver for inviscid flows ' '20041214 3.3 P.Eliasson Propeller file name added ' '20050929 3.3 S.-H. Peng Added hybrid HYB0 model ' +'20090630 4.1 St.Schmidt Added Spalart-Allmaras DDES/IDDES model ' +'20090630 4.1 St.Schmidt Added wall damping to Smagorinsky model ' +'20090630 4.1 St.Schmidt Added SST DES/DDES/IDDES models ' +'20090630 4.1 St.Schmidt Added Spalart-Allmaras Edwards model ' +'20090630 3.2 St.Schmidt Added Rung LLR-model (3D) ' +'20090630 4.1 St.Schmidt Added Rung SALSA RANS + DES/DDES/IDDES m' '20090909 4.1 P.Eliasson SA DES,LES model parameter changes ' '20090916 4.1 P.Eliasson Line-implicit time marching ' ' Weak boundary conditions ' @@ -66,17 +72,23 @@ MODIFICATIONS,L ,61,1,0 '20091108 4.1 L. Cavagna AERO-ELASTIC/THERMALS mostly uniformed ' '20091125 4.1 M. Tormalm Probe and inlet analysis in history file' '20100223 4.1 P.Eliasson Regional data supplied by file ' -'20100315 5.0 M.Tomac Added kw-sst + transition model ' +'20100630 4.1 St.Schmidt Added HD model ' +'20110911 5.2 St.Schmidt Added Rung RQEVM model ' +'20100315 5.2 M.Tomac Added kw-sst + transition model ' +'20160118 5.3 P.Eliasson Added option IPOSTMUL ' +'20170112 6.0 P.Eliasson Preprocessor variables updated ' +'20170905 6.0 P.Eliasson Low dispersion scheme introduced ' +'20180119 1.1 P.Eliasson M-Edge ' * *************************************************************************** * * * NAME,TYPE(I,R,L),NDIM(=1),NSIZE(=# ELEMENTS),NSUB(=0) * * * VERSION,S ,1,1,0 -'5.0.0 ' +'3.2.0 ' * TITLE,L ,1,1,0 -'Edge, an edge-based flow-solver from FFA. ( www.foi.se/edge ) ' +'M-Edge, an edge-based flow solver ' * *************************************************************************** * * @@ -84,13 +96,13 @@ TITLE,L ,1,1,0 * * *************************************************************************** * Initial solution option : -* 0 = Not available (created in code) +* 0 = Not available (created freestream in code) * 1 = Available (if full multigrid on coarsest level) * 2 = Restart of previous run * (solution available finest level incompatible with full MG) * INPRES,I ,1,1,0 -0 +0 * *************************************************************************** * * @@ -102,7 +114,7 @@ INPRES,I ,1,1,0 * 1 = N-S * INSEUL,I ,1,1,0 -__INSEUL__ +__INSEUL__ * *************************************************************************** * Gas model : @@ -184,6 +196,14 @@ IRTAB,R ,5,10,0 287.0 287.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 * *************************************************************************** +* Species discretization type : +* 0 = Central +* 1 = Upwind +* +IUPWMS,I ,1,1,0 +1 +* +*************************************************************************** * Entropy fix for upwind species equations (NSPEC>1) : * |SPEFIX| < 1. : Not valid (equivalent to SPEFIX=0) * |SPEFIX| > 1. : Entropy fix scaled isotropically with the spectral @@ -197,6 +217,28 @@ SPEFIX,R ,1,1,0 1.01 * *************************************************************************** +* Order of species equations upwind scheme : +* 1 = 1st order +* 2 = 2nd order +* +IORDMS,I ,1,1,0 +2 +* +*************************************************************************** +* 2nd order dissipation coefficient +* For species equations using central scheme only +* +VIS2MS,R ,1,1,0 +0.1 +* +*************************************************************************** +* 4th order dissipation coefficient +* For species equations using central scheme only +* +VIS4MS,R ,1,1,0 +0.1 +* +*************************************************************************** * * * --------------- Viscous data --------------- * * * @@ -212,7 +254,7 @@ ICOVIC,I ,1,1,0 * Viscosity at TFREE * RMU,R ,1,1,0 -1.7894E-5 +1.8478E-5 * *************************************************************************** * Prandtl number @@ -221,12 +263,13 @@ PRREF,R ,1,1,0 0.72 * *************************************************************************** -* Thin layer approximation or full N-S : -* 0 = Thin layer -* 1 = Full N-S +* Viscous type of operator : +* 1 = Original compact operator for normal derivatives +* 2 = Edge normal operator to account for skewness +* 3 = Coarse operator based on nodal gradients * -IFULNS,I ,1,1,0 -1 +IVISCOP,I ,1,1,0 +1 * *************************************************************************** * * @@ -248,19 +291,19 @@ __TFREE__ * Freestream velocity in x-direction * UFREE,R ,1,1,0 -__UFREE__ +__UFREE__ * *************************************************************************** * Freestream velocity in y-direction * VFREE,R ,1,1,0 --0.0 +__VFREE__ * *************************************************************************** * Freestream velocity in z-direction * WFREE,R ,1,1,0 -__WFREE__ +__WFREE__ * *************************************************************************** * Freestream turbulence intensity @@ -268,7 +311,7 @@ __WFREE__ * the values used in RANS computation) * TUFREE,R ,1,1,0 -0.0010 +0.001 * *************************************************************************** * Freestream viscosity ratio. @@ -301,19 +344,51 @@ SPMFFREE,R ,5,1,0 * (discrete adjoint implemented) * IUPWIN,I ,1,1,0 -0 +0 +* +*************************************************************************** +* Convection type of scheme +* 0 = *DEFAULT = conventional scheme according to IUPWIN +* 1 = Dispersion scheme, provide dispersion constant DISPCONST +* +CONVTYPE,I ,1,1,0 +0 +* +*************************************************************************** +* Dispersion constant >=0 if CONVTYPE=1 +* Default 1/3 which can increase accuracy & dispersion +* +DISPCONST,R ,1,1,0 +0.33333334 +* +*************************************************************************** +* Artificial dissipation type +* 0 = *DEFAULT = scalar dissipation +* 1 = Matrix dissipation; may require high VIS2 for good shock capturing +* 2 = Langer scalar dissipation; default values VIS2 = 8.0, VIS4 = 1/64 +* 3 = Langer matrix dissipation; default values VIS2 = 8.0, VIS4 = 1/64 +* +DISSTYPE,I ,1,1,0 +0 * *************************************************************************** * 2nd order dissipation coefficients (central schemes only) * VIS2,R ,1,1,0 -0.5 +0.5 * *************************************************************************** * 4th order dissipation coefficients (central schemes only) * VIS4,R ,1,1,0 -0.02 +0.02 +* +*************************************************************************** +* Damping coefficient of 4th order dissipation at different element types +* 1 = *DEFAULT = no damping. 0 = full damping +* +DAMPVIS4,R ,1,1,0 +1.0 * *************************************************************************** * Dissipation type (upwind schemes) : @@ -343,7 +418,7 @@ IORDER,I ,1,1,0 * 3 = SUPERBEE * LIMITE,I ,1,2,0 -1 1 +1 1 * *************************************************************************** * 1 = *DEFAULT = ROE averaging in upwind scheme @@ -361,7 +436,17 @@ IROEAV,I ,1,1,0 * radius (SR), = (ENTFIX-1.)*SR * ENTRFX,R ,1,2,0 -1.2 1.05 +1.2 1.05 +* +*************************************************************************** +* Method for gradient calculation +* 0 = *DEFAULT = Green-Gauss +* 1 = Weighted least squares +* 2 = LQ-factorization, Hessian second order. 27 weights(3d) stored/node +* Only applied on the finest grid, Green-Gauss on coarser +* +GRADTYPE,I ,1,1,0 +0 * *************************************************************************** * * @@ -371,7 +456,7 @@ ENTRFX,R ,1,2,0 * Number of grids : > 1 --> multigrid * NGRID,I ,1,1,0 -__NGRID__ +__NGRID__ * *************************************************************************** * Multigrid strategy : @@ -389,19 +474,20 @@ MGRSTR,I ,1,1,0 * Full multigrid : * 0 = No full multigrid * 1 = *DEFAULT = Full multigrid +* 2 = Grid sequencing * IFULMG,I ,1,1,0 -1 +1 * *************************************************************************** * Maximum number of full multigrid cycles (*DEFAULT = 500) * NFMGCY,I ,1,1,0 -500 +500 * *************************************************************************** * Number of sweeps (relaxation or Runge-Kutta) in -* fine-to-coarse (FC) and coarse-to-fine sweep of the MG cycle +* fine-to-coarse (FC) sweep of the MG cycle * NSWPFC,I ,1,10,0 1 1 1 1 1 1 1 1 1 1 @@ -433,7 +519,7 @@ MGSIMP,I ,1,1,0 * (central schemes only) * VIS0,R ,1,1,0 -0.15 +0.15 * *************************************************************************** * Reduction factor for CFL on coarser grids @@ -471,6 +557,177 @@ IRESTR,I ,1,1,0 * ---------------- Steady state time integrator ---------------- * * * *************************************************************************** +* Iterative approach towards steady state +* 1 = DEFAULT, Runge-Kutta explicit time marching with local time steps +* 2 = LU-SGS with global time steps +* 3 = Implicit Scheme: preconditioned Krylov method with local time steps +* +STSTCONV,I ,1,1,0 +3 +* +*************************************************************************** +* Number of LU-SGS sweeps/iteration >=1 +* +NSWPSGS,I ,1,1,0 +1 +* +*************************************************************************** +* Amount of LHS numerical dissipation for LU-SGS and GMRES +* LU-SGS, EPS2SGS>0: Fraction added to eigenvalues >= 0.15 recommended +* LU-SGS, EPS2SGS<0: Scalar LHS dissipation, <= -0.5 +* GMRES, EPS2SGS>0: Entropy fix for eigenvalues of first order upwind LHS +* GMRES, EPS2SGS<0: Scalar LHS dissipation, <= -0.3 recommended +* +EPS2SGS,R ,1,1,0 +0.3 +* +*************************************************************************** +* 0 LU-SGS, standard +* 1 LU-SGS combined with line-implicit +* +SGSLINES,I ,1,1,0 +1 +* +*************************************************************************** +* Node order for LU-SGS sweeps +* 1 original node ordering +* 2 neighbor to neighbor +* 3 implcit lines first + original ordering +* +SGSORDER,I ,1,1,0 +1 +* +*************************************************************************** +* Output LU-SGS convergence if NSWPSGS >=2 +* 0 Do not output linear convergene to screen +* 1 Output linear convergence to screen +* +SGSCONVOUT,I ,1,1,0 +0 +* +*************************************************************************** +* PETSc Krylov subspace solver +* 1 GMRES : Better convergence than BI-CGSTAB, higher memory usage +* 2 BI-CGSTAB : Slightly worse convergence than GMRES, lower memory usage +* +KSP_TYPE,I ,1,1,0 +1 +* +*************************************************************************** +* Convergence tolerance for the mean flow equations +* Reduce if non-linear iterations have convergence problems, e.g. 1.0E-2 +* A smaller value will increase the computational cost +* +KSP_TOL,R ,1,1,0 +0.01 +* +*************************************************************************** +* Convergence tolerance for the turbulence equations +* +KSP_TOL_TURB,R ,1,1,0 +0.01 +* +*************************************************************************** +* Maximum number of iterations of the Krylov subspace solver +* +KSP_MAXITER,I ,1,1,0 +100 +* +*************************************************************************** +* Minimum number of iterations of the mean flow +* +KSP_MINITER_M,I ,1,1,0 +1 +* +*************************************************************************** +* Minimum number of iterations of the turbulence +* +KSP_MINITER_T,I ,1,1,0 +1 +* +*************************************************************************** +* Minimum number of iterations of the turbulence omega equation +* +KSP_MINITER_TW,I ,1,1,0 +3 +* +*************************************************************************** +* To save memory, one can interrupt the GMRES iteration and make a restart +* Helps to reduce memory requirements, but can lead to worse convergence +* of the linear solver. Typically, a restart after around 30 iterations +* is used. With GMRES_RESTART=100, restart is basically not employed +* +GMRES_RESTART,I ,1,1,0 +100 +* +*************************************************************************** +* Level of fill-in for ILU(p) preconditioning, for the mean flow +* ILU(0) means: The preconditioning matrix has the same non-zero pattern +* as the original matrix +* ILU(p), p>0: Compared to the original matrix +* - new elements are allowed to be added to the preconditioning matrix +* - the preconditioning matrix has a more dense non-zero pattern +* More levels increase memory consumption. Dont go beyond 2 in 3d +* +ILU_LEVELS,I ,1,1,0 +0 +* +*************************************************************************** +* Level of fill-in for ILU(p) preconditioning, for the turbulence equation +* +ILU_LEVELS_TURB,I ,1,1,0 +0 +* +*************************************************************************** +* Iteration interval to update the KSP Jacobian and Preconditioner +* 1 : Update the Jacobian and the Preconditioner at every time step +* : Update the Jacobian and the Preconditioner at every time step +* Steady state: second value used when residuals reduced to JAC_PC_SWITCH +* Unsteady: only the second value used +* +JAC_PC_UPDATE,I ,1,2,0 +1 3 +* +*************************************************************************** +* Residual, when to switch the update interval JAC_PC_UPDATE +* of the Jacobian and Preconditioner Computation +* Used for steady state calculations only +* +JAC_PC_SWITCH,R ,1,1,0 +-4.0 +* +*************************************************************************** +* Additive Schwarz Overlap of domains for parallel ILU preconditioning +* Typical values: 1 or larger (1=PETSc default) +* +ASM_OVERLAP,I ,1,1,0 +1 +* +*************************************************************************** +* KSP_AUTO_DAMPING: Try to automatically dampen the Newton update +* (solution update of the linear Krylov solver) +* Damping is only applied for turbulent flow simulations +* +KSP_AUTO_DAMPING,I ,1,1,0 +0 +* +*************************************************************************** +* Update factor (damping) of the the Newton update for the mean flow +* (solution update of the linear Krylov solver) +* Only applied for turbulent flow simulations and KSP_AUTO_DAMPING=0 +* Range: 0 no damping +* +UPD_FACTOR_M,R ,1,1,0 +0.5 +* +*************************************************************************** +* Corresponding update factor (damping) for the turbulence equations +* (solution update of the linear Krylov solver) +* +UPD_FACTOR_T,R ,1,1,0 +0.5 +* +*************************************************************************** * Number of stages > 0 <=10 * NSTAGE,I ,1,1,0 @@ -523,6 +780,30 @@ IRKCOLI,R ,1,10,0 1.0 1.0 0.8 1.0 1.0 1.0 1.0 1.0 1.0 1.0 * *************************************************************************** +* Steady-state stabilization +* 0 = DEFAULT = No stabilization +* 1 = Stabilisation +* +ITSTAB,I ,1,1,0 +0 +* +*************************************************************************** +* Mean integration coefficient +* = + CSTM*(q-) +* DEFAULT = 0.01 +* +CSTM,R ,1,1,0 +0.01 +* +*************************************************************************** +* Feed-back stabilization coefficient +* dq/dt + (1-C)*R(q) = 0, where C=CSTS if R(q)*(q-)<0 else C=0 +* DEFAULT = 0.01 +* +CSTS,R ,1,1,0 +0.01 +* +*************************************************************************** * * * ---------------- Implicit residual smoothing ---------------- * * * @@ -556,7 +837,7 @@ NSMCYC,I ,1,1,0 * 2 = Explicit time accurate with global time step * ITIMAQ,I ,1,1,0 -0 +0 * *************************************************************************** * 0 = Compute global time step @@ -570,7 +851,7 @@ ISGLTS,I ,1,1,0 * Explicit (ITIMAQ = 2) global time step * DELTAT,R ,1,1,0 -0.005 +0.005 * *************************************************************************** * Convergence criteria, stop when MAX(LOG(RES)) < RESTAQ @@ -614,19 +895,27 @@ GAMN,R ,1,1,0 0.0 * *************************************************************************** -* 0 = Switch off time-averaging for steady computations -* 1 = Set for time-averaging for unsteady calculations (ITIMAQ/=0) -* This parameter can be used in update.ainp to intermediately -* start time-averaging during the running course without -* interruption by setting ITMEAN = 1 -* If INPRES=2 and T-averaging also restarts with a previous -* field, then one should set ITMEAN = 1 +* Extrapolation scheme for startguess for each global timestep +* 0 = *DEFAULT = No extrapolation +* 1 = Linear extrapolation +* 2 = Quadratic extrapolation +* +EXTPOLTA,I ,1,1,0 +0 +* +*************************************************************************** +* 0 = Switch off time-averaging +* 1 = Set for time-averaging for steady or unsteady calculations +* Ouputs primitive variables and mean of velocity squares +* plus pressure square. Intended for LES type calculations +* 2 = Time averaging but only output of primitive variables +* Restart will read in old mean flow file * ITMEAN,I ,1,1,0 0 * *************************************************************************** -* Time step at which time-averaging starts, function only if ITMEAN=1 +* Time step at which time-averaging starts, function only if ITMEAN>0 * NOTE: If INPRES=2 and ITMEAN=1, indicates that time-averaging * restarts from a previous run. In this case, ITERM < ITER. * @@ -634,30 +923,18 @@ ITERM,I ,1,1,0 2000 * *************************************************************************** -* Specification of points where time-series is recorded for unsteady -* computations. Specify a number of points NPTHIS each defined with -* coordinates and their local identities (boundary or interior points) -* Example: (p1,p2,p3,p4,p5) -* The data is specified in T_HIS_INFO +* Sample the mean solution to multiple (_it.bout) files every ITMULOU +* iterations * -NPTHIS,I ,1,1,0 +ITMULOU,I ,1,1,0 0 * -T_HIS_INFO,N ,0,0,2 -* -POINT_COOD,R ,5,3,0 -0.6858 0.0 -0.0254 -0.7874 -0.0254 0.0381 -1.2954 -0.01905 0.04445 -0.8128 -0.1016 -0.0254 -1.0414 0.0 0.0 +*************************************************************************** +* Specification of NPTHIS points where time-series is recorded for +* unsteady computations. Points input data file under CFIPOINTS * -POINT_INDEX,L ,5,2,0 -'point1 ' 'bound_point ' -'point2 ' 'bound_point ' -'point3 ' 'bound_point ' -'point4 ' 'bound_point ' -'point5 ' 'inter_point ' +NPTHIS,I ,1,1,0 +0 * *************************************************************************** * * @@ -668,50 +945,91 @@ POINT_INDEX,L ,5,2,0 * 1 = Rotation * IROT,I ,1,1,0 - 0 +0 * *************************************************************************** * Angular velocity and direction when rotation * The absolute value is the angular velocity * OMGROT,R ,1,3,0 -0.0 0.0 0.0 +0.0 0.0 0.0 * *************************************************************************** * Coordinates for one point on axis off rotation * OMGCOO,R ,1,3,0 -0.0 0.0 0.0 +0.0 0.0 0.0 * *************************************************************************** * * * ---------------- Maneuver option ---------------- * * * +MANEUVERS,N ,0,0,6 +* *************************************************************************** -* 0 = No maneuvers -* 1 = For specification of arbitrary maneuvers +* Maneuver solution options : +* 0 = No maneuver +* 1,2,3 = Prescribed motion through MOVDAT +* 202 = Coupled to external solver (full model) * -IMANEUV,I ,1,1,0 +ISOOPT,I ,1,1,0 0 * *************************************************************************** +* Coordinates for body reference point +* +CMGCOO,R ,1,3,0 +0.0 0.0 0.0 +* +*************************************************************************** * Input data format for the maneuver data given in MOVDAT * 0 = Time, v1, v2, v3, w1, w2, w3 * 1 = Time, v1, v2, v3, w1, w2, w3, W1, W2, W3 * 2 = Time, v1, v2, v3, w1, w2, w3, O11, O12, O13, O21.... O33 * -IMVDTF,I ,1,1,0 -1 -* MOVDAT,R ,2,10,0 0.0 0.0 0.0 0.0 0.0 0.0 125.664 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 125.664 0.0 0.0 125.664 * *************************************************************************** +* Controls where in the inner, ITIN, loop the maneuver functions +* are active : +* 0 = Never +* N>0 = Every N inner iterations (ignored for ISOOPT= 0, 1, 2, 3) +* +NITIN,I ,1,1,0 +1 +* +*************************************************************************** +* File CFIPRM (inp) is the update/comunication file for MANEUVERS +* - required to modify internal parameters and data +* +CFIPRM,L ,1,1,0 +'Edge_mv.ainp ' +* +*************************************************************************** +* Controls IO format for external coupled solution +* 0 = use same file +* 1 = append _in and _out to filename +* +IPRM,I ,1,1,0 +0 +* +*************************************************************************** * * * ---------------- On-line mesh-deformation ---------------- * * * -MSHD_CONSTANTS,N ,0,0,4 +MSHD_CONSTANTS,N ,0,0,11 +* +*************************************************************************** +* Mesh deformation method +* 1 = Laplacian smoothing, diagonal preconditioner +* 2 = Spring analogy, diagonal preconditioner +* 3 = Radial basis functions (RBF) +* 4 = Translation/rotation/scaling (TRS) +* +METHOD,I ,1,1,0 +1 * *************************************************************************** * Convergence tolerance : dimensionless @@ -732,12 +1050,46 @@ MAXIT,I ,1,1,0 1000 * *************************************************************************** -* Initialisation : -* 0 = Reference grid -* 1 = Previous timestep +* Exponent for Laplace edge length denominator * -INICO,I ,1,1,0 -0 +DEXP,R ,1,1,0 +2.0 +* +*************************************************************************** +* Exponent for spring analogy stiffness definition +* +BETA,R ,1,1,0 +1.0 +* +*************************************************************************** +* RBF function: 1(Gaussian), 2(IQ), 3(IMQ), 4(Wennland), 5(3rd poly.) +* +RBF_OPT,I ,1,1,0 +1 +* +*************************************************************************** +* RBF shape factor. Default 0.8 +* +RBF_SF_REL,R ,1,1,0 +0.8 +* +*************************************************************************** +* TRS inboard distance to start deformation decay. Estimated if <0 +* +TRS_DIST_IN,R ,1,1,0 +-1.0 +* +*************************************************************************** +* TRS outboard distance to end deformation decay. Estimated if <0 +* +TRS_DIST_OUT,R ,1,1,0 +-1.0 +* +*************************************************************************** +* TRS damping to damp out scaling/rotation. 0-not damped, 1-all damped +* +TRS_DAMP,R ,1,1,0 +0.0 * *************************************************************************** * * @@ -751,7 +1103,7 @@ INICO,I ,1,1,0 ISDRY,I ,1,1,0 0 * -AEROELASTICS,N ,0,0,13 +AEROELASTICS,N ,0,0,14 * *************************************************************************** * Aeroelastic solution options : @@ -760,19 +1112,20 @@ AEROELASTICS,N ,0,0,13 * 101 = Modal prescribed motion * 102 = Modal coupled : orthogonal modes central difference * 112 = Modal externally coupled +* 202 = Coupled to external solver (full model) * ISOOPT,I ,1,1,0 0 * *************************************************************************** * Mode type options : -* 0 = None, non-modal solutions PLACEHOLDER -* 1 = Grid perturbation -* 2 = Boundary surface -* 3 = Structural nodes PLACEHOLDER +* 0 = None, non-modal solutions PLACEHOLDER +* 1 = Volume grid perturbation +* 2 = Grid boundary surface - on-line mesh deformation +* 3 = Structural nodes PLACEHOLDER * MODOPT,I ,1,1,0 -0 +0 * *************************************************************************** * External force input on/off (coupled solutions only) : @@ -821,6 +1174,14 @@ CFIPRM,L ,1,1,0 'Edge_at.ainp ' * *************************************************************************** +* Controls IO format for external coupled solution +* 0 = use same file +* 1 = append _in and _out to filename +* +IPRM,I ,1,1,0 +0 +* +*************************************************************************** * File CFIMOD (amod) is the structural input file * - required only when the H coupling matrix is stored (in .bedg) * - contains structural data : nodes, label, coordinates and modal shapes @@ -857,13 +1218,13 @@ CFIMOS,L ,1,1,0 * Format: prefix.suffix / imin imax istep * CFIDIS,L ,1,1,0 -'bdis/Edge.bdis / 1001 1006 1 ' +'bdis/Edge.bdis / 1001 1006 1 ' * *************************************************************************** * * * ---------------- Aero-thermal solution options ---------------- * * * -AEROTHERMALS,N ,0,0,13 +AEROTHERMALS,N ,0,0,14 * *************************************************************************** * Thermal coupling option : @@ -901,6 +1262,14 @@ CFIPRM,L ,1,1,0 'Edge_at.ainp ' * *************************************************************************** +* Controls IO format for external coupled solution +* 0 = use same file +* 1 = append _in and _out to filename +* +IPRM,I ,1,1,0 +0 +* +*************************************************************************** * File CFIMOD (amod) is the structural thermal input file * - required only when the H coupling matrix is stored (in .bedg) * - contains structural data : nodes, label, coordinates and modal shapes @@ -974,90 +1343,205 @@ RELAXT,R ,1,1,0 * *************************************************************************** * * -* ------------- Trajectory solution options ---------------- * +* ---------------- Gusts options ---------------- * * * -TRAJECTORY_DATA,N ,0,0,7 +GUSTS,N ,0,0,6 * *************************************************************************** -* Trajectory solution options : -* 0 = No trajectory computations -* 1 = Quasi-steady trajectory computations -* 2 = Unsteady trajectory computations +* Gust option : +* 0 = No gust +* 1 = 1-cos gust * -ITRAJEC,I ,1,1,0 +ISOOPT,I ,1,1,0 0 * *************************************************************************** -* Total number of iterations for trajectory computations -* -ITMXTRA,I ,1,1,0 -1 -* -*************************************************************************** -* If ITRAJEC = 1, time step for trajectory computations +* Vector normal to gust plane +* cross-product of GUST_NORM x FREE_STREAM velocity vector +* determines direction of perturbation velocity * -DELTAT_QS,R ,1,1,0 -0.005 +GUST_NORM,R ,1,3,0 +0.0 1.0 0.0 * *************************************************************************** -* Gravitational acceleration vector in main frame coordinates +* Gust initial coordinate * -GRAV,R ,3,1,0 -0.0 -0.0 -9.80665 +ORIG,R ,1,3,0 +0.0 0.0 0.0 * *************************************************************************** -* Number of substeps for 6DoF computations +* Gust spatial length * -NSUBST,I ,1,1,0 -1 +LENGTH,R ,1,1,0 +1.0 * *************************************************************************** -* Parameter controlling the grid check ( > 1.) +* Gust perturbation velocity along PERT_DIR * -CHKMAX,R ,1,1,0 -10.0 +PERT_VEL,R ,1,1,0 +1.0 * *************************************************************************** -* The trajectory file ( .atra ) +* Initial time delay * -CFITRA,L ,1,1,0 -'Edge.atra ' +T_DELAY,R ,1,1,0 +0.0 * *************************************************************************** * * -* ---------------- Propeller options ---------------- * +* ---------------- MRF options ---------------- * * * +MRFS,N ,0,0,8 +* *************************************************************************** -* File CFIDISK contains the input to the propeller disk model -* The model is used as a boundary condition -* The data file is not in FFA-format, it is copied from Multnas +* MRF option : +* 0 = No moving reference frame +* 1 = Grid velocity from file +* 2 = Sinusoidal grid velocity +* 3 = Sinusoidal angular velocity * -CFIDISK,L ,1,1,0 -'prin.dat ' +IMRFOPT,I ,1,1,0 +0 * *************************************************************************** -* * -* --------------- Preconditioning options --------------- * -* * -*************************************************************************** -* 0 = No preconditioning -* 1 = Preconditioning for low Mach numbers used +* Normal to sinusodial grid veloctiy +* +SIN_NORM,R ,1,3,0 +0.0 1.0 0.0 +* +*************************************************************************** +* Frequency in hertz +* +SIN_FREQ,R ,1,1,0 +1.0 +* +*************************************************************************** +* Velocity amplitude of sinusoidal grid velocity +* +SIN_VEL,R ,1,1,0 +1.0 +* +*************************************************************************** +* Amplitude of sinusoidal motion in degrees +* +SIN_AMP,R ,1,1,0 +1.0 +* +*************************************************************************** +* Coordinates for one point on axis off rotation +* +OMGCOO_MRF,R ,1,3,0 +0.0 1.0 0.0 +* +*************************************************************************** +* Name of farfield BC to updated velocities for rotational motion +* +BNAME,L ,1,1,0 +'outer ' +* +*************************************************************************** +* Axis of rotation : +* 0 = x-axis +* 1 = y-axis +* 2 = z-axis +* +IAXIS,I ,1,1,0 +0 +* +*************************************************************************** +* * +* ------------- Trajectory solution options ---------------- * +* * +TRAJECTORY_DATA,N ,0,0,8 +* +*************************************************************************** +* Trajectory solution options : +* 0 = No trajectory computations +* 1 = Quasi-steady trajectory computations +* 2 = Unsteady trajectory computations +* +ITRAJEC,I ,1,1,0 +0 +* +*************************************************************************** +* Total number of iterations for trajectory computations +* +ITMXTRA,I ,1,1,0 +1 +* +*************************************************************************** +* If ITRAJEC = 1, time step for trajectory computations +* +DELTAT_QS,R ,1,1,0 +0.0 +* +*************************************************************************** +* Gravitational acceleration vector in main frame coordinates +* +GRAV,R ,3,1,0 +0.0 +0.0 +9.80665 +* +*************************************************************************** +* Number of substeps for 6DoF computations +* +NSUBST,I ,1,1,0 +1 +* +*************************************************************************** +* Increment for mesh modifications in the inner loop +* +NITIN,I ,1,1,0 +1 +* +*************************************************************************** +* Parameter controlling the grid check ( > 1.) +* +CHKMAX,R ,1,1,0 +10.0 +* +*************************************************************************** +* The trajectory file ( .atra ) +* +CFITRA,L ,1,1,0 +'Edge.atra ' +* +*************************************************************************** +* * +* ---------------- Propeller options ---------------- * +* * +*************************************************************************** +* File CFIDISK contains the input to the propeller disk model +* The model is used as a boundary condition +* The data file is not in FFA-format, it is copied from Multnas +* +CFIDISK,L ,1,1,0 +'prin.dat ' +* +*************************************************************************** +* * +* ----------- Low speed preconditioning options ----------- * +* * +*************************************************************************** +* 0 = No preconditioning +* 1 = Preconditioning for low Mach numbers used, Turkel/Choi&Merkle +* 2 = Preconditioning with art. diss. based on conservative variables +* Values fixed: RM0=1, RKBET1=1, RKBET2=1, ALPHAP=0, DELTA=0 * IPREPA,I ,1,1,0 -0 +0 * *************************************************************************** * RKBET2 preconditioning parameter, -* RKBET1 = 1, ALPHAP = 0 as default in Euranus +* RKBET1=1, ALPHAP=0, DELTA=1 as default in Euranus * RKBET2,R ,1,1,0 4.0 * *************************************************************************** * RM0 preconditioning parameter, -* RKBET1=1, ALPHAP=0 as default in Euranus +* RKBET1=1, ALPHAP=0, DELTA=1 as default in Euranus * RM0,R ,1,1,0 1.0 @@ -1077,13 +1561,13 @@ IPREVIS,I ,1,1,0 * Number of iterations and type of time step * ITMAX,I ,1,1,0 -1000 +__ITMAX__ * *************************************************************************** * Convergence criteria (orders of magnitude of reduction) * RESRED,R ,1,1,0 --6.0 +-7.0 * *************************************************************************** * Convergence criteria full multigrid (orders of magnitude @@ -1094,9 +1578,46 @@ RESFMG,R ,1,1,0 * *************************************************************************** * Inviscid CFL number +* <= 1.5 for explicit iterations (STSTCONV=1) +* Large for LU-SGS iterations (STSTCONV=2) +* Not used for implicit Krylov method (STSCONV=3) * CFL,R ,1,1,0 -__CFL__ +1.25 +* +*************************************************************************** +* Initial CFL number for implicit GMRES calculations (STSCONV=3) +* Applied at first iteration, then increased with factor CFL_FACTOR_KSP +* +CFL_MIN_KSP,R ,1,1,0 +1.5 +* +*************************************************************************** +* Factor to increase CFL for implicit Krylov calculations (STSCONV=3) +* +CFL_FACTOR_KSP,R ,1,1,0 +1.2 +* +*************************************************************************** +* Maximum CFL number for implicit Krylov calculations (STSCONV=3) +* First value applied for residuals >= CFL_MAX_SWITCH +* Second value applied for residuals < CFL_MAX_SWITCH +* +CFL_MAX_KSP,R ,1,2,0 +10000.0 10000.0 +* +*************************************************************************** +* Residual, when to switch the maximum allowed CFL number (STSCONV=3) +* +CFL_MAX_SWITCH,R ,1,1,0 +-4.0 +* +*************************************************************************** +* Maximum CFL number in largest cell for implicit calculations (STSCONV=3) +* Reduction -> 1 often required for supersonic calculations +* +CFL_MAX_MAXVOL,R ,1,1,0 +10000.0 * *************************************************************************** * Viscous CFL number : If > 0 --> Time step is minimum of inviscid @@ -1122,7 +1643,7 @@ CFLLI,R ,1,1,0 * Too high values will degrade rate of convergence * EPS2LI,R ,1,1,0 -0.2 +0.3 * *************************************************************************** * * @@ -1131,18 +1652,26 @@ EPS2LI,R ,1,1,0 *************************************************************************** * Post processing options. Additional fields written to * a separate solution file name in CFIPOST -* 0 = No post processing, default -* 1 = Cf, q, yplus for viscous calc. (INPRES=1) -* 2 = Mach number, Cp, energy, temperature and local time step -* 4 = Distance to the wall and laminar/turbulent field -* 8 = Residuals to all unknowns -* 16 = Reynolds stresses -* 32 = Gradients of the dependent variables -* 64 = Flow curvature vector +* 0 = No post processing, default +* 1 = Cf, q, yplus for viscous calc. (INPRES=1) +* 2 = Mach, Cp, energy, temperature, total states and local time step +* 4 = Distance to the wall and laminar/turbulent field +* 8 = Residuals to all unknowns +* 16 = Reynolds stresses +* 32 = Gradients of the dependent variables +* 64 = Flow curvature vector +* 128 = Q criterion and abs(rograd) * For several options sum up (1+4=5 for cf and wall distance etc.) * IPOST,I ,1,1,0 -15 +1 +* +*************************************************************************** +* Sample post solution files (IPOST>0) to multiple (_it.bout) +* files every IPMULOU iterations. +* +IPMULOU,I ,1,1,0 +0 * *************************************************************************** * EPS is used to avoid division by zero @@ -1155,7 +1684,7 @@ EPS,R ,1,1,0 * iterations. These files are used for restarts. See also INPRES. * IWRSOL,I ,1,1,0 -500 +100 * *************************************************************************** * Output solution to file, runtime parameter @@ -1177,6 +1706,13 @@ ISTOUT,I ,1,1,0 1 * *************************************************************************** +* Check the existence of new data in file update.ainp +* each BEDITFREQ iteration +* +BEDITFREQ,I ,1,1,0 +1 +* +*************************************************************************** * Number of the residual on standard output every ISTOUT iteration * 1 = *DEFAULT (density) * @@ -1185,8 +1721,8 @@ NRRES,I ,1,1,0 * *************************************************************************** * Region data file supplied (CFIREG) with regional data -* 0 = No data file supplied -* 1 = Data file supplied +* 0 = No data file supplied +* 1 = Data file supplied * IREGIONDATA,I ,1,1,0 0 @@ -1196,66 +1732,66 @@ IREGIONDATA,I ,1,1,0 * Direction for force 1 * IDCLP,R ,1,3,0 -__IDCLP__ +__IDCLP__ * *************************************************************************** * Direction for drag force * Direction for force 2 * IDCDP,R ,1,3,0 -__IDCDP__ +__IDCDP__ * *************************************************************************** * Direction for side force * Direction for force 3 * IDCCP,R ,1,3,0 -__IDCCP__ +__IDCCP__ * *************************************************************************** * Direction for pitch moment * Direction for moment 1 * IDCMP,R ,1,3,0 -__IDCMP__ +__IDCMP__ * *************************************************************************** * Direction for yaw moment * Direction for moment 2 * IDCNP,R ,1,3,0 -__IDCNP__ +__IDCNP__ * *************************************************************************** * Direction for roll moment * Direction for moment 3 * IDCRP,R ,1,3,0 -__IDCRP__ +__IDCRP__ * *************************************************************************** * Point for moment * IXMP,R ,1,3,0 -__IXMP__ +__IXMP__ * *************************************************************************** * Reference area to non-dimensionalize forces and moments * SREF,R ,1,1,0 -__SREF__ +__SREF__ * *************************************************************************** * Reference chord for non-dimensionalizing moment * CREF,R ,1,1,0 -__CREF__ +__CREF__ * *************************************************************************** * Reference width for non-dimensionalizing moment * BREF,R ,1,1,0 -__BREF__ +__BREF__ * *************************************************************************** * * @@ -1290,19 +1826,19 @@ VGKCONSTSCALE,R ,1,1,0 * The mesh file ( .bmsh ) * CFIMSH,L ,1,1,0 -'Edge.bmsh' +'../grid/Edge.bmsh ' * *************************************************************************** * Edge file name from the preprocessor ( .bedg ) * CFIEDG,L ,1,1,0 -'Edge.bedg' +'../grid/Edge.bedg ' * *************************************************************************** * The boundary condition file ( .aboc ) * CFIBOC,L ,1,1,0 -'Edge.aboc' +'../grid/Edge.aboc ' * *************************************************************************** * The residual file ( .bres ) @@ -1336,7 +1872,7 @@ CFIPOST,L ,1,1,0 'Post.bout ' * *************************************************************************** -* History recording file (PROBE_RECORD > 0 or AIP_RECORD=1) +* History recording file (IACHIS=1, IPROBEHIS=1, IAIPHIS=1 or ISAMMSH>=1) * CFIHIS,L ,1,1,0 'history.bhis ' @@ -1348,6 +1884,42 @@ CFIREG,L ,1,1,0 'Edge.areg ' * *************************************************************************** +* Boundary sampling input file (ISAMP = 1) +* +CFIBSI,L ,1,1,0 +'Edge_bsm.asmp ' +* +*************************************************************************** +* Aerodynamic coeff history input file (IACHIS = 1) +* +CFIACI,L ,1,1,0 +'Edge_ach.asmp ' +* +*************************************************************************** +* Probe coeff history input file (IPROBEHIS = 1) +* +CFIPRB,L ,1,1,0 +'Edge_probe.asmp ' +* +*************************************************************************** +* Inlet coeff history input file (IAIPHIS = 1) +* +CFIAIP,L ,1,1,0 +'Edge_aip.asmp ' +* +*************************************************************************** +* Points history input file (NPTHIS >= 1) +* +CFIPOINTS,L ,1,1,0 +'Edge_points.asmp ' +* +*************************************************************************** +* Embedded mesh file for sampling (ISAMMSH >= 1) +* +CFISAMMSH,L ,1,1,0 +'Edge_mshsamp.bmsh ' +* +*************************************************************************** * * * ---------------- Turbulence models ---------------- * * * @@ -1358,46 +1930,61 @@ CFIREG,L ,1,1,0 * 4 = LES models * ITURB,I ,1,1,0 -2 +2 * *************************************************************************** -* Choice of turbulence model : +* Choice of turbulence model : (*=DEFAULT , !=recommended) * --------------------------------------- -* Default RANS model: -* W&J EARSM + Hellsten k-omega -* Other recommended RANS models: +* One-equation models: * Spalart-Allmaras one-eq model -* Menter SST k-omega -* Menter SST k-omega + Transition -* W&J DRSM + Hellsten k-omega -* Other RANS models: +* Spalart-Allmaras Edwards +* Rung SALSA (!) +* Two-equation models: +* W&J EARSM + Hellsten k-omega (*) * std Wilcox k-omega -* Kok_TNT k-omega * Menter BSL k-omega +* Menter SST k-omega (!) +* Menter SST k-omega + Transition +* Kok_TNT k-omega * Wilcox LRN k-omega * PDH LRN k-omega -* W&J EARSM + std k-omega +* Rung LLR k-omega +* Rung RQEVM + std k-omega +* Algebraic Reynolds-Stress models: * GWJ Generalized EARSM + std k-omega +* W&J EARSM + std k-omega * W&J EARSM + Kok_TNT k-omega * W&J EARSM + Menter BSL k-omega 2009 * W&J CC-EARSM + std k-omega * W&J CC-EARSM + Menter BSL k-omega * W&J CC-EARSM + Hellsten k-omega -* W&J DRSM + std k-omega +* Reynolds-Stress Transport models: * SSG DRSM + std k-omega * SSG DRSM + Hellsten k-omega +* W&J DRSM + std k-omega +* W&J DRSM + Hellsten k-omega +* W&J DAMPED DRSM + Hellsten k-omega * --------------------------------------- * DES and hybrid RANS-LES models : * Spalart-Allmaras DES model +* Spalart-Allmaras DDES (!) +* Spalart-Allmaras IDDES +* Rung SALSA-DES +* Rung SALSA-DDES (!) +* Rung SALSA-IDDES +* Menter SST-DES k-omega +* Menter SST-DDES k-omega (!) +* Menter SST-IDDES k-omega * Peng hybrid HYB0 model * --------------------------------------- * LES models : * Smagorinsky SGS model +* Smagorinsky SGS with wall damping * Spalart-Allmaras SGS model * Yoshizawa k-eq SGS model * TURB_MOD_NAME,L ,1,1,0 -'W&J EARSM + Hellsten k-omega' +'Spalart-Allmaras one-eq model ' * *************************************************************************** * 2-eq. turb. discretization type : @@ -1405,21 +1992,21 @@ TURB_MOD_NAME,L ,1,1,0 * 1 = Upwind * IUPWKZ,I ,1,1,0 -1 +0 * *************************************************************************** * 2nd order dissipation coefficient -* For K-E turbulence model (central schemes only) +* For turbulence models using central scheme only * VIS2KZ,R ,1,1,0 -0.5 +1.0 * *************************************************************************** * 4th order dissipation coefficient -* For K-E turbulence model (central schemes only) +* For turbulence models using central scheme only * VIS4KZ,R ,1,1,0 -0.02 +0.05 * *************************************************************************** * Entropy fix for upwind 2-eq. turbulence : @@ -1432,7 +2019,7 @@ VIS4KZ,R ,1,1,0 * numbers to decrease the numerical dissipation. * TURFIX,R ,1,1,0 --1.1 +-1.1 * *************************************************************************** * Order of turbulent upwind scheme : @@ -1446,14 +2033,14 @@ IORDKZ,I ,1,1,0 * * * --------------- Turbulence model definitions --------------- * * * -TURB_MOD_LIST,N ,0,0,26 +TURB_MOD_LIST,N ,0,0,39 * *************************************************************************** * * -* --------------- Peng hybrid HYB0 model ---------- * +* --------------- Spalart-Allmaras one-eq model ---------- * * * -TURB_MOD,L ,1,1,12 -'Peng hybrid HYB0 model ' +TURB_MOD,L ,1,1,20 +'Spalart-Allmaras one-eq model ' * *************************************************************************** * * @@ -1467,51 +2054,67 @@ WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 -0 +1 * NAME_UNK,S ,1,1,0 -'zero_equation ' +'turb_work_vt ' * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -3 -* -NAME_EX,S ,1,3,0 -'turb_viscosity ' 'turb_kin_energy ' 'rans_viscosity ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * -CSGS,R ,1,1,0 -0.12 +SACB1,R ,1,1,0 +0.1355 * -CKSGS,R ,1,1,0 -0.0 +SACB2,R ,1,1,0 +0.622 * -SIGMA,R ,1,1,0 +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACT1,R ,1,1,0 1.0 * +SACT2,R ,1,1,0 +2.0 +* +SACT3,R ,1,1,0 +1.1 +* +SACT4,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* PRT,R ,1,1,0 -0.4 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 * *************************************************************************** * * @@ -1519,39 +2122,41 @@ PRT,R ,1,1,0 * * *************************************************************************** * * -* ----- Extra fields (turbulent viscosity) ----- * +* ----- Viscous terms ----- * * * -EXTRVA,S ,1,1,0 -'FXHYB0P ' +FVISF,S ,1,1,0 +'GVISSA ' * *************************************************************************** * * -* --------------- Spalart-Allmaras DES model ---------- * +* ----- Source terms ----- * * * -TURB_MOD,L ,1,1,20 -'Spalart-Allmaras DES model ' +GSOUKZ,S ,1,1,0 +'HSSAM ' * *************************************************************************** * * -* ----- Turbulence model program logics ----- * +* ----- Extra fields (turbulent viscosity) ----- * * * +EXTRVA,S ,1,1,0 +'FXSAM ' +* *************************************************************************** * * -* ----- If wall distance required ----- * +* --------------- Spalart-Allmaras Edwards ---------- * * * -WALL_DIST,N ,0,0,0 +TURB_MOD,L ,1,1,18 +'Spalart-Allmaras Edwards ' * *************************************************************************** * * -* ----- If MUT exists ----- * +* ----- Turbulence model program logics ----- * * * -MUT_AVAIL,N ,0,0,0 -* *************************************************************************** * * -* ----- turb_kin_energy exists ----- * +* ----- If wall distance required ----- * * * -K_AVAIL,N ,0,0,0 +WALL_DIST,N ,0,0,0 * *************************************************************************** * * @@ -1565,21 +2170,15 @@ NAME_UNK,S ,1,1,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_viscosity ' 'turb_kin_energy ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * -CDES,R ,1,1,0 -0.65 -* SACB1,R ,1,1,0 0.1355 * @@ -1595,6 +2194,12 @@ SACW2,R ,1,1,0 SACW3,R ,1,1,0 2.0 * +SACT1,R ,1,1,0 +1.0 +* +SACT2,R ,1,1,0 +2.0 +* SACV1,R ,1,1,0 7.1 * @@ -1626,88 +2231,22 @@ FVISF,S ,1,1,0 * * * ----- Source terms ----- * * * -GSOULDES,S ,1,1,0 -'HSDESSA ' +GSOUKZ,S ,1,1,0 +'HSSAE ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXDESSA ' -* -*************************************************************************** -* * -* --------------- Smagorinsky model ---------- * -* * -TURB_MOD,L ,1,1,11 -'Smagorinsky SGS model ' -* -*************************************************************************** -* * -* ----- Turbulence model program logics ----- * -* * -*************************************************************************** -* * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Number and names of turb eqs ----- * -* * -NTEQ,I ,1,1,0 -0 -* -NAME_UNK,S ,1,1,0 -'zero_equation ' -* -*************************************************************************** -* * -* ----- Number and names of additional turb fields ----- * -* * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_viscosity ' 'turb_kin_energy ' -* -*************************************************************************** -* * -* ----- Turbulence model constants ----- * -* * -CSMAG,R ,1,1,0 -0.2 -* -CKSMG,R ,1,1,0 -0.0 -* -SIGMA,R ,1,1,0 -1.0 -* -PRT,R ,1,1,0 -0.4 -* -*************************************************************************** -* * -* ----- Extra fields (turbulent viscosity and kin. energy) ----- * -* * -EXTRVA,S ,1,1,0 -'FXSGSSMG ' +'FXSAM ' * *************************************************************************** * * -* --------------- Spalart-Allmaras SGS model ---------- * +* --------------- Rung SALSA ---------- * * * -TURB_MOD,L ,1,1,20 -'Spalart-Allmaras SGS model ' +TURB_MOD,L ,1,1,16 +'Rung SALSA ' * *************************************************************************** * * @@ -1715,24 +2254,12 @@ TURB_MOD,L ,1,1,20 * * *************************************************************************** * * -* ----- If wall distance required ----- * +* ----- If wall distance required ----- * * * WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 @@ -1743,21 +2270,15 @@ NAME_UNK,S ,1,1,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_viscosity ' 'turb_kin_energy ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * -CDES,R ,1,1,0 -0.65 -* SACB1,R ,1,1,0 0.1355 * @@ -1804,22 +2325,22 @@ FVISF,S ,1,1,0 * * * ----- Source terms ----- * * * -GSOULDES,S ,1,1,0 -'HSLESSA ' +GSOUKZ,S ,1,1,0 +'HSSALSA ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSGSSA ' +'FXSAM ' * *************************************************************************** * * -* --------------- Yoshizawa SGS model ---------- * +* --------------- std Wilcox k-omega --------------- * * * -TURB_MOD,L ,1,1,14 -'Yoshizawa k-eq SGS model ' +TURB_MOD,L ,1,1,12 +'std Wilcox k-omega ' * *************************************************************************** * * @@ -1827,33 +2348,18 @@ TURB_MOD,L ,1,1,14 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 -1 +2 * -NAME_UNK,S ,1,1,0 -'turb_kin_energy ' +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -1861,17 +2367,23 @@ NAME_EX,S ,1,1,0 * * * ----- Turbulence model constants ----- * * * -CKMU,R ,1,1,0 -0.07 +ALPHA,R ,1,1,0 +0.556 * -CEPS,R ,1,1,0 -1.05 +BETA,R ,1,1,0 +0.075 * -SIGMA,R ,1,1,0 -1.0 +BETAS,R ,1,1,0 +0.09 +* +SIGMA,R ,1,2,0 +0.5 0.5 +* +WALL,R ,1,1,0 +10.0 * PRT,R ,1,1,0 -0.4 +0.9 * *************************************************************************** * Coefficient to multiply the limit of the turbulent production @@ -1886,31 +2398,24 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * -* ----- Viscous terms ----- * -* * -FVISF,S ,1,1,0 -'GVISLESYO ' -* -*************************************************************************** -* * * ----- Source terms ----- * * * -GSOULDES,S ,1,1,0 -'HSLESYO ' +GSOUKZ,S ,1,1,0 +'HSSTKW ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSGSYO ' +'FXSTKW ' * *************************************************************************** * * -* --------------- Spalart-Allmaras one-eq. model ---------- * +* --------------- Wilcox LRN k-omega --------------- * * * -TURB_MOD,L ,1,1,23 -'Spalart-Allmaras one-eq model ' +TURB_MOD,L ,1,1,12 +'Wilcox LRN k-omega ' * *************************************************************************** * * @@ -1918,79 +2423,40 @@ TURB_MOD,L ,1,1,23 * * *************************************************************************** * * -* ----- If wall distance required ----- * -* * -WALL_DIST,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 -1 +2 * -NAME_UNK,S ,1,1,0 -'turb_work_vt ' +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_viscosity ' 'turb_kin_energy ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * -SACB1,R ,1,1,0 -0.1355 -* -SACB2,R ,1,1,0 -0.622 +ALPHA,R ,1,1,0 +0.56 * -SACW1,R ,1,1,0 -3.24 +BETA,R ,1,1,0 +0.075 * -SACW2,R ,1,1,0 -0.3 +BETAS,R ,1,1,0 +0.09 * -SACW3,R ,1,1,0 -2.0 +SIGMA,R ,1,2,0 +0.5 0.5 * -SACT1,R ,1,1,0 +WALL,R ,1,1,0 1.0 * -SACT2,R ,1,1,0 -2.0 -* -SACT3,R ,1,1,0 -1.1 -* -SACT4,R ,1,1,0 -2.0 -* -SACV1,R ,1,1,0 -7.1 -* -SIGMA,R ,1,1,0 -1.5 -* PRT,R ,1,1,0 0.9 * @@ -2007,31 +2473,24 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * -* ----- Viscous terms ----- * -* * -FVISF,S ,1,1,0 -'GVISSA ' -* -*************************************************************************** -* * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSAM ' +'HSLRKWW ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSAM ' +'FXLRKWW ' * *************************************************************************** * * -* --------------- std Wilcox k-omega --------------- * +* --------------- Rung LLR k-omega --------------- * * * -TURB_MOD,L ,1,1,15 -'std Wilcox k-omega ' +TURB_MOD,L ,1,1,11 +'Rung LLR k-omega ' * *************************************************************************** * * @@ -2039,19 +2498,7 @@ TURB_MOD,L ,1,1,15 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Number and names of turb eqs ----- * +* ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 2 @@ -2061,11 +2508,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2073,14 +2517,11 @@ NAME_EX,S ,1,1,0 * * * ----- Turbulence model constants ----- * * * -ALPHA,R ,1,1,0 -0.556 -* BETA,R ,1,1,0 -0.075 +1.83 * BETAS,R ,1,1,0 -0.09 +0.27777 * SIGMA,R ,1,2,0 0.5 0.5 @@ -2089,7 +2530,7 @@ WALL,R ,1,1,0 10.0 * PRT,R ,1,1,0 -0.9 +1.0 * *************************************************************************** * Coefficient to multiply the limit of the turbulent production @@ -2107,21 +2548,21 @@ PROLIM,R ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSTKW ' +'HSSTLLR ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSTKW ' +'FXLLR ' * *************************************************************************** * * -* --------------- Wilcox LRN k-omega --------------- * +* --------------- Rung RQEVM + std k-omega --------------- * * * -TURB_MOD,L ,1,1,15 -'Wilcox LRN k-omega ' +TURB_MOD,L ,1,1,12 +'Rung RQEVM + std k-omega ' * *************************************************************************** * * @@ -2129,18 +2570,6 @@ TURB_MOD,L ,1,1,15 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 @@ -2151,11 +2580,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2164,7 +2590,7 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.56 +0.556 * BETA,R ,1,1,0 0.075 @@ -2176,10 +2602,10 @@ SIGMA,R ,1,2,0 0.5 0.5 * WALL,R ,1,1,0 -1.0 +10.0 * PRT,R ,1,1,0 -0.9 +1.0 * *************************************************************************** * Coefficient to multiply the limit of the turbulent production @@ -2197,21 +2623,21 @@ PROLIM,R ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSLRKWW ' +'HSSTRQEVM ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXLRKWW ' +'FXRQEVM ' * *************************************************************************** * * -* --------------- Rung LLR k-omega --------------- * +* --------------- PDH LRN k-omega --------------- * * * -TURB_MOD,L ,1,1,15 -'Rung LLR k-omega ' +TURB_MOD,L ,1,1,13 +'PDH LRN k-omega ' * *************************************************************************** * * @@ -2219,18 +2645,6 @@ TURB_MOD,L ,1,1,15 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 @@ -2241,11 +2655,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2254,22 +2665,25 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.43 +0.42 * BETA,R ,1,1,0 -0.83 +0.075 * BETAS,R ,1,1,0 -0.27667 +0.09 * SIGMA,R ,1,2,0 -0.5 0.5 +1.25 0.7407 +* +SIGMAD_CROSS,R ,1,1,0 +0.75 * WALL,R ,1,1,0 -10.0 +1.0 * PRT,R ,1,1,0 -1.0 +0.9 * *************************************************************************** * Coefficient to multiply the limit of the turbulent production @@ -2287,21 +2701,21 @@ PROLIM,R ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSTLLR ' +'HSLRKWP ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXLLR ' +'FXLRKWP ' * *************************************************************************** * * -* --------------- PDH LRN k-omega --------------- * +* --------------- Kok_TNT k-omega --------------- * * * -TURB_MOD,L ,1,1,16 -'PDH LRN k-omega ' +TURB_MOD,L ,1,1,13 +'Kok_TNT k-omega ' * *************************************************************************** * * @@ -2309,18 +2723,6 @@ TURB_MOD,L ,1,1,16 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 @@ -2331,11 +2733,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2344,7 +2743,7 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.42 +0.553 * BETA,R ,1,1,0 0.075 @@ -2353,13 +2752,13 @@ BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -1.25 0.7407 +0.66666 0.5 * SIGMAD_CROSS,R ,1,1,0 -0.75 +0.5 * WALL,R ,1,1,0 -1.0 +10.0 * PRT,R ,1,1,0 0.9 @@ -2380,21 +2779,21 @@ PROLIM,R ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSLRKWP ' +'HSSTCD ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXLRKWP ' +'FXSTKW ' * *************************************************************************** * * -* --------------- Kok_TNT k-omega --------------- * +* --------------- Menter BSL k-omega --------------- * * * -TURB_MOD,L ,1,1,16 -'Kok_TNT k-omega ' +TURB_MOD,L ,1,1,17 +'Menter BSL k-omega ' * *************************************************************************** * * @@ -2402,15 +2801,9 @@ TURB_MOD,L ,1,1,16 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * +* ----- If wall distance required ----- * * * -K_AVAIL,N ,0,0,0 +WALL_DIST,N ,0,0,0 * *************************************************************************** * * @@ -2424,11 +2817,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2437,19 +2827,25 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.553 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 * BETA,R ,1,1,0 0.075 * +BETA2,R ,1,1,0 +0.0828 +* BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -0.66666 0.5 +0.5 0.5 * -SIGMAD_CROSS,R ,1,1,0 -0.5 +SIGMA2,R ,1,2,0 +1.0 0.856 * WALL,R ,1,1,0 10.0 @@ -2470,10 +2866,17 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZBSL ' +* +*************************************************************************** +* * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSTCD ' +'HSKWBSL ' * *************************************************************************** * * @@ -2484,10 +2887,10 @@ EXTRVA,S ,1,1,0 * *************************************************************************** * * -* --------------- W&J EARSM + std k-omega --------------- * +* --------------- Menter SST k-omega --------------- * * * TURB_MOD,L ,1,1,18 -'W&J EARSM + std k-omega ' +'Menter SST k-omega ' * *************************************************************************** * * @@ -2495,21 +2898,9 @@ TURB_MOD,L ,1,1,18 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If Turb_extra_anis exists ----- * +* ----- If wall distance required ----- * * * -AEX_AVAIL,I ,1,0,0 +WALL_DIST,N ,0,0,0 * *************************************************************************** * * @@ -2523,11 +2914,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2536,16 +2924,28 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.556 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 * BETA,R ,1,1,0 0.075 * +BETA2,R ,1,1,0 +0.0828 +* BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -0.5 0.5 +0.85 0.5 +* +SIGMA2,R ,1,2,0 +1.0 0.856 +* +A1,R ,1,1,0 +0.31 * WALL,R ,1,1,0 10.0 @@ -2553,21 +2953,12 @@ WALL,R ,1,1,0 PRT,R ,1,1,0 0.9 * -C1,R ,1,1,0 -1.8 -* -*************************************************************************** -* Diffusion correction of EARSM model -* -CD,R ,1,1,0 -0.0 -* *************************************************************************** * Coefficient to multiply the limit of the turbulent production * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0E20 +1.0 * *************************************************************************** * * @@ -2575,24 +2966,31 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZBSL ' +* +*************************************************************************** +* * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSTKW ' +'HSKWSST ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXWJAR ' +'FXSST ' * *************************************************************************** * * -* ------------ GWJ Generalized EARSM + std k-omega ------------ * +* --------------- MENTER SST k-omega + Transition --------------- * * * -TURB_MOD,L ,1,1,19 -'GWJ Generalized EARSM + std k-omega ' +TURB_MOD,L ,1,1,40 +'Menter SST k-omega + Transition ' * *************************************************************************** * * @@ -2600,57 +2998,77 @@ TURB_MOD,L ,1,1,19 * * *************************************************************************** * * -* ----- If MUT exists ----- * +* ----- Transition model options ----- * * * -MUT_AVAIL,N ,0,0,0 +*************************************************************************** +* Governing correlation functions: +* 1 = *DEFAULT = Tomac-Pettersson +* +ICORRF,I ,1,1,0 +1 * *************************************************************************** * * -* ----- Turb_kin_energy exists ----- * +* ----- If wall distance required ----- * * * -K_AVAIL,N ,0,0,0 +WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If Turb_extra_anis exists ----- * +* ----- If MUT exists ----- * * * -AEX_AVAIL,I ,1,0,0 +MUT_AVAIL,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Turb_kin_energy exists ----- * +* * +K_AVAIL,N ,0,0,0 * *************************************************************************** * * * ----- Number and names of turb eqs ----- * * * NTEQ,I ,1,1,0 -2 +4 * -NAME_UNK,S ,1,2,0 -'turb_kin_energy ' 'turb_omega ' +NAME_UNK,S ,1,4,0 +'turb_kin_energy ' 'turb_omega ' 'intermittency ' 're_theta_t ' * *************************************************************************** * * * ----- Number and names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* -NAME_EX,S ,1,1,0 -'turb_viscosity ' +NAME_EX,S ,1,4,0 +'turb_viscosity ' 'laminar_turbulen' 'flt ' 'retc ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.556 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 * BETA,R ,1,1,0 0.075 * +BETA2,R ,1,1,0 +0.0828 +* BETAS,R ,1,1,0 0.09 * -SIGMA,R ,1,2,0 -0.5 0.5 +SIGMA,R ,1,4,0 +0.85 0.5 1.0 2.0 +* +SIGMA2,R ,1,2,0 +1.0 0.856 +* +A1,R ,1,1,0 +0.31 * WALL,R ,1,1,0 10.0 @@ -2658,24 +3076,73 @@ WALL,R ,1,1,0 PRT,R ,1,1,0 0.9 * -C1,R ,1,1,0 -1.8 +*************************************************************************** +* * +* ----- Transition model constants ----- * +* * +CA1,R ,1,1,0 +2.0 * -C,R ,1,1,0 -0.56 +CA2,R ,1,1,0 +0.06 * -*************************************************************************** -* Diffusion correction of EARSM model +CE1,R ,1,1,0 +1.0 * -CD,R ,1,1,0 +CE2,R ,1,1,0 +50.0 +* +CALPHA,R ,1,1,0 +0.5 +* +CTT,R ,1,1,0 +0.03 +* +S1,R ,1,1,0 +2.0 +* +FLTC1,R ,1,1,0 +0.3 +* +FLTC2,R ,1,1,0 +60.0 +* +FLTC3,R ,1,1,0 +53000.0 +* +FLTC4,R ,1,1,0 0.0 * +RETCC1,R ,1,1,0 +0.9 +* +RETCC2,R ,1,1,0 +0.1 +* +RETCC3,R ,1,1,0 +490.0 +* +RETCC4,R ,1,1,0 +390.0 +* +RETCC5,R ,1,1,0 +150.0 +* +LTLIM,R ,1,2,0 +-0.1 0.1 +* +TULIM,R ,1,1,0 +0.027 +* +RETTLIM,R ,1,2,0 +20.0 1450.0 +* *************************************************************************** * Coefficient to multiply the limit of the turbulent production * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0E20 +1.0 * *************************************************************************** * * @@ -2683,41 +3150,36 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * -* ----- Source terms ----- * +* ----- Viscous terms ----- * * * -GSOUKZ,S ,1,1,0 -'HSSTKW ' +FVISF,S ,1,1,0 +'GVISKBSLT ' * *************************************************************************** * * -* ----- Extra fields (turbulent viscosity) ----- * +* ----- Source terms ----- * * * -EXTRVA,S ,1,1,0 -'FXGWJA ' +GSOUKZ,S ,1,2,0 +'HSKWSST ' 'HSINTRETT ' * *************************************************************************** * * -* --------------- W&J EARSM + Kok_TNT k-omega --------------- * +* ----- Extra fields (turbulent viscosity) ----- * * * -TURB_MOD,L ,1,1,19 -'W&J EARSM + Kok_TNT k-omega ' +EXTRVA,S ,1,2,0 +'FXSST ' 'FXINTRETT ' * *************************************************************************** * * -* ----- Turbulence model program logics ----- * -* * -*************************************************************************** -* * -* ----- If MUT exists ----- * +* --------------- GWJ Generalized EARSM + std k-omega --------------- * * * -MUT_AVAIL,N ,0,0,0 +TURB_MOD,L ,1,1,16 +'GWJ Generalized EARSM + std k-omega ' * *************************************************************************** * * -* ----- Turb_kin_energy exists ----- * +* ----- Turbulence model program logics ----- * * * -K_AVAIL,N ,0,0,0 -* *************************************************************************** * * * ----- If Turb_extra_anis exists ----- * @@ -2736,11 +3198,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2749,7 +3208,7 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.525 +0.556 * BETA,R ,1,1,0 0.075 @@ -2758,10 +3217,7 @@ BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -0.66666 0.55 -* -SIGMAD_CROSS,R ,1,1,0 -0.3 +0.5 0.5 * WALL,R ,1,1,0 10.0 @@ -2772,6 +3228,9 @@ PRT,R ,1,1,0 C1,R ,1,1,0 1.8 * +C,R ,1,1,0 +0.56 +* *************************************************************************** * Diffusion correction of EARSM model * @@ -2794,21 +3253,21 @@ PROLIM,R ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSSTCD ' +'HSSTKW ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXWJAR ' +'FXGWJA ' * *************************************************************************** * * -* --------------- MENTER BSL k-omega --------------- * +* --------------- W&J EARSM + std k-omega --------------- * * * -TURB_MOD,L ,1,1,20 -'Menter BSL k-omega ' +TURB_MOD,L ,1,1,15 +'W&J EARSM + std k-omega ' * *************************************************************************** * * @@ -2816,21 +3275,9 @@ TURB_MOD,L ,1,1,20 * * *************************************************************************** * * -* ----- If wall distance required ----- * -* * -WALL_DIST,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * +* ----- If Turb_extra_anis exists ----- * * * -K_AVAIL,N ,0,0,0 +AEX_AVAIL,I ,1,0,0 * *************************************************************************** * * @@ -2844,11 +3291,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2857,38 +3301,38 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.55317 -* -ALPHA2,R ,1,1,0 -0.4403547 +0.556 * BETA,R ,1,1,0 0.075 * -BETA2,R ,1,1,0 -0.0828 -* BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 0.5 0.5 * -SIGMA2,R ,1,2,0 -1.0 0.856 -* WALL,R ,1,1,0 10.0 * PRT,R ,1,1,0 0.9 * +C1,R ,1,1,0 +1.8 +* +*************************************************************************** +* Diffusion correction of EARSM model +* +CD,R ,1,1,0 +0.0 +* *************************************************************************** * Coefficient to multiply the limit of the turbulent production * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0 +1.0E20 * *************************************************************************** * * @@ -2896,31 +3340,24 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * -* ----- Viscous terms ----- * -* * -FVISF,S ,1,1,0 -'GVISKBSL ' -* -*************************************************************************** -* * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSKWBSL ' +'HSSTKW ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSTKW ' +'FXWJAR ' * *************************************************************************** * * -* --------------- MENTER SST k-omega --------------- * +* --------------- W&J EARSM + Kok_TNT k-omega --------------- * * * -TURB_MOD,L ,1,1,21 -'Menter SST k-omega ' +TURB_MOD,L ,1,1,16 +'W&J EARSM + Kok_TNT k-omega ' * *************************************************************************** * * @@ -2928,21 +3365,9 @@ TURB_MOD,L ,1,1,21 * * *************************************************************************** * * -* ----- If wall distance required ----- * -* * -WALL_DIST,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * +* ----- If Turb_extra_anis exists ----- * * * -K_AVAIL,N ,0,0,0 +AEX_AVAIL,I ,1,0,0 * *************************************************************************** * * @@ -2956,11 +3381,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -2969,28 +3391,19 @@ NAME_EX,S ,1,1,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.55317 -* -ALPHA2,R ,1,1,0 -0.4403547 +0.525 * BETA,R ,1,1,0 0.075 * -BETA2,R ,1,1,0 -0.0828 -* BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -0.85 0.5 -* -SIGMA2,R ,1,2,0 -1.0 0.856 +0.66666 0.55 * -A1,R ,1,1,0 -0.31 +SIGMAD_CROSS,R ,1,1,0 +0.3 * WALL,R ,1,1,0 10.0 @@ -2998,12 +3411,21 @@ WALL,R ,1,1,0 PRT,R ,1,1,0 0.9 * +C1,R ,1,1,0 +1.8 +* +*************************************************************************** +* Diffusion correction of EARSM model +* +CD,R ,1,1,0 +0.0 +* *************************************************************************** * Coefficient to multiply the limit of the turbulent production * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0 +1.0E20 * *************************************************************************** * * @@ -3011,31 +3433,24 @@ PROLIM,R ,1,1,0 * * *************************************************************************** * * -* ----- Viscous terms ----- * -* * -FVISF,S ,1,1,0 -'GVISKBSL ' -* -*************************************************************************** -* * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSKWSST ' +'HSSTCD ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXSST ' +'FXWJAR ' * *************************************************************************** * * -* --------------- MENTER SST k-omega + Transition --------------- * +* --------------- W&J EARSM + Menter BSL k-omega 2009 ---------- * * * -TURB_MOD,L ,1,1,41 -'Menter SST k-omega + Transition ' +TURB_MOD,L ,1,1,23 +'W&J EARSM + Menter BSL k-omega 2009 ' * *************************************************************************** * * @@ -3043,82 +3458,61 @@ TURB_MOD,L ,1,1,41 * * *************************************************************************** * * -* ----- Transition model options ----- * -* * -*************************************************************************** -* Governing correlation functions: -* 0 = *DEFAULT = Tomac-Pettersson -* 1 = Menter-Langtry -* -ICORRF,I ,1,1,0 -0 -* -*************************************************************************** -* * * ----- If wall distance required ----- * * * WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * +* ----- If Turb_extra_anis exists ----- * * * -MUT_AVAIL,N ,0,0,0 +AEX_AVAIL,I ,1,0,0 * *************************************************************************** * * -* ----- Turb_kin_energy exists ----- * +* ----- Number and names of turb eqs ----- * * * -K_AVAIL,N ,0,0,0 +NTEQ,I ,1,1,0 +2 * -*************************************************************************** -* * -* ----- Number and names of turb eqs ----- * -* * -NTEQ,I ,1,1,0 -4 -* -NAME_UNK,S ,1,4,0 -'turb_kin_energy ' 'turb_omega ' 'intermittency ' 're_tetha_t ' +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* NAME_EX,S ,1,2,0 -'turb_viscosity ' 'laminar_turbulen' +'turb_viscosity ' 'mut_simp_diff ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.55317 +0.556 * ALPHA2,R ,1,1,0 -0.4403547 +0.44035 * BETA,R ,1,1,0 0.075 * +BETA1,R ,1,1,0 +0.075 +* BETA2,R ,1,1,0 0.0828 * BETAS,R ,1,1,0 0.09 * -SIGMA,R ,1,4,0 -0.85 0.5 1.0 2.0 +SIGMA,R ,1,2,0 +0.5 0.5 * SIGMA2,R ,1,2,0 1.0 0.856 * -A1,R ,1,1,0 -0.31 -* WALL,R ,1,1,0 10.0 * @@ -3126,66 +3520,11 @@ PRT,R ,1,1,0 0.9 * *************************************************************************** -* * -* ----- Transition model constants ----- * -* * -CA1,R ,1,1,0 -2.0 -* -CA2,R ,1,1,0 -0.06 -* -CE1,R ,1,1,0 -1.0 -* -CE2,R ,1,1,0 -50.0 -* -CALPHA,R ,1,1,0 -0.5 -* -CTT,R ,1,1,0 -0.03 -* -S1,R ,1,1,0 -2.0 -* -FLTC1,R ,1,1,0 -0.3188 -* -FLTC2,R ,1,1,0 -100.0 +* Diffusion correction of EARSM model * -FLTC3,R ,1,1,0 +CD,R ,1,1,0 0.0 * -FLTC4,R ,1,1,0 -73000.0 -* -RETCC1,R ,1,1,0 -0.92345 -* -RETCC2,R ,1,1,0 -187.0 -* -RETCC3,R ,1,1,0 -0.615 -* -RETCC4,R ,1,1,0 -30000.0 -* -RETCC5,R ,1,1,0 -0.0080 -* -LTLIM,R ,1,2,0 --0.1 0.1 -* -TULIM,R ,1,1,0 -0.027 -* -RETTLIM,R ,1,1,0 -20.0 -* *************************************************************************** * Coefficient to multiply the limit of the turbulent production * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) @@ -3193,6 +3532,15 @@ RETTLIM,R ,1,1,0 PROLIM,R ,1,1,0 1.0 * +A1,R ,1,1,0 +1.245 +* +A3,R ,1,1,0 +1.8 +* +A4,R ,1,1,0 +2.25 +* *************************************************************************** * * * ----- Turbulence model subroutines ----- * @@ -3202,28 +3550,28 @@ PROLIM,R ,1,1,0 * ----- Viscous terms ----- * * * FVISF,S ,1,1,0 -'GVISKBSLT ' +'GVISKZBSL ' * *************************************************************************** * * * ----- Source terms ----- * * * -GSOUKZ,S ,1,2,0 -'HSKWSSTT ' 'HSINTRETT ' +GSOUKZ,S ,1,1,0 +'HSKWBSL ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * -EXTRVA,S ,1,2,0 -'FXSST ' 'FXINTRETT ' +EXTRVA,S ,1,1,0 +'FXWJAR ' * *************************************************************************** * * -* --------------- W&J EARSM + MENTER BSL k-omega 2009 ---------- * +* --------------- W&J EARSM + Hellsten k-omega --------------- * * * -TURB_MOD,L ,1,1,26 -'W&J EARSM + Menter BSL k-omega 2009 ' +TURB_MOD,L ,1,1,24 +'W&J EARSM + Hellsten k-omega ' * *************************************************************************** * * @@ -3237,18 +3585,6 @@ WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- If Turb_extra_anis exists ----- * * * AEX_AVAIL,I ,1,0,0 @@ -3265,29 +3601,26 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_viscosity ' 'mut_simp_diff ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.556 +0.518 * ALPHA2,R ,1,1,0 -0.44035 +0.44 * BETA,R ,1,1,0 -0.075 +0.0747 * BETA1,R ,1,1,0 -0.075 +0.0747 * BETA2,R ,1,1,0 0.0828 @@ -3296,22 +3629,31 @@ BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,2,0 -0.5 0.5 +1.1 0.53 * SIGMA2,R ,1,2,0 -1.0 0.856 +1.1 1.0 * -WALL,R ,1,1,0 -10.0 +SIGMAD1,R ,1,1,0 +1.0 * -PRT,R ,1,1,0 -0.9 +SIGMAD2,R ,1,1,0 +0.4 +* +CMIX,R ,1,1,0 +1.5 * *************************************************************************** * Diffusion correction of EARSM model * CD,R ,1,1,0 -0.0 +2.2 +* +WALL,R ,1,1,0 +10.0 +* +PRT,R ,1,1,0 +0.9 * *************************************************************************** * Coefficient to multiply the limit of the turbulent production @@ -3320,15 +3662,9 @@ CD,R ,1,1,0 PROLIM,R ,1,1,0 1.0 * -A1,R ,1,1,0 -1.245 -* -A3,R ,1,1,0 +C1,R ,1,1,0 1.8 * -A4,R ,1,1,0 -2.25 -* *************************************************************************** * * * ----- Turbulence model subroutines ----- * @@ -3338,14 +3674,14 @@ A4,R ,1,1,0 * ----- Viscous terms ----- * * * FVISF,S ,1,1,0 -'GVISKBSL ' +'GVISKZHELLST ' * *************************************************************************** * * * ----- Source terms ----- * * * GSOUKZ,S ,1,1,0 -'HSKWBSL ' +'HSSHELLST ' * *************************************************************************** * * @@ -3358,7 +3694,7 @@ EXTRVA,S ,1,1,0 * * * --------------- W&J CC-EARSM + std k-omega --------------- * * * -TURB_MOD,L ,1,1,20 +TURB_MOD,L ,1,1,17 'W&J CC-EARSM + std k-omega ' * *************************************************************************** @@ -3367,18 +3703,6 @@ TURB_MOD,L ,1,1,20 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- If Turb_extra_anis exists ----- * * * AEX_AVAIL,I ,1,0,0 @@ -3401,11 +3725,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -3472,7 +3793,7 @@ EXTRVA,S ,1,1,0 * * * --------------- W&J CC-EARSM + Menter BSL k-omega --------------- * * * -TURB_MOD,L ,1,1,26 +TURB_MOD,L ,1,1,23 'W&J CC-EARSM + Menter BSL k-omega ' * *************************************************************************** @@ -3487,18 +3808,6 @@ WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- If Turb_extra_anis exists ----- * * * AEX_AVAIL,I ,1,0,0 @@ -3521,11 +3830,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -3591,7 +3897,7 @@ A0,R ,1,1,0 * ----- Viscous terms ----- * * * FVISF,S ,1,1,0 -'GVISKBSL ' +'GVISKZBSL ' * *************************************************************************** * * @@ -3609,10 +3915,10 @@ EXTRVA,S ,1,1,0 * *************************************************************************** * * -* --------------- W&J EARSM + Hellsten k-omega --------------- * +* --------------- W&J CC-EARSM + Hellsten k-omega --------------- * * * -TURB_MOD,L ,1,1,27 -'W&J EARSM + Hellsten k-omega ' +TURB_MOD,L ,1,1,26 +'W&J CC-EARSM + Hellsten k-omega ' * *************************************************************************** * * @@ -3626,21 +3932,15 @@ WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * +* ----- If Turb_extra_anis exists ----- * * * -K_AVAIL,N ,0,0,0 +AEX_AVAIL,I ,1,0,0 * *************************************************************************** * * -* ----- If Turb_extra_anis exists ----- * +* ----- If curvature correction ----- * * * -AEX_AVAIL,I ,1,0,0 +CURV_CORR,I ,1,0,0 * *************************************************************************** * * @@ -3654,11 +3954,8 @@ NAME_UNK,S ,1,2,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -1 -* NAME_EX,S ,1,1,0 'turb_viscosity ' * @@ -3721,6 +4018,9 @@ PROLIM,R ,1,1,0 C1,R ,1,1,0 1.8 * +A0,R ,1,1,0 +-0.72 +* *************************************************************************** * * * ----- Turbulence model subroutines ----- * @@ -3748,10 +4048,10 @@ EXTRVA,S ,1,1,0 * *************************************************************************** * * -* --------------- W&J-CC EARSM + Hellsten k-omega --------------- * +* --------------- W&J DRSM + std k-omega --------------- * * * -TURB_MOD,L ,1,1,29 -'W&J CC-EARSM + Hellsten k-omega ' +TURB_MOD,L ,1,1,20 +'W&J DRSM + std k-omega ' * *************************************************************************** * * @@ -3759,96 +4059,42 @@ TURB_MOD,L ,1,1,29 * * *************************************************************************** * * -* ----- If wall distance required ----- * -* * -WALL_DIST,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If MUT exists ----- * +* ----- If Reynolds_stress exists ----- * * * -MUT_AVAIL,N ,0,0,0 +RS_AVAIL,N ,0,0,0 * *************************************************************************** * * -* ----- Turb_kin_energy exists ----- * +* ----- Number and names of turb eqs ----- * * * -K_AVAIL,N ,0,0,0 +NTEQ,I ,1,2,0 +5 7 * -*************************************************************************** -* * -* ----- If Turb_extra_anis exists ----- * -* * -AEX_AVAIL,I ,1,0,0 +NAME_UNK,S ,1,12,0 +'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'turb_omega ' 'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'RS_13 ' 'RS_23 ' 'turb_omega ' * *************************************************************************** * * -* ----- If curvature correction ----- * +* ----- Names of additional turb fields ----- * * * -CURV_CORR,I ,1,0,0 +NAME_EX,S ,1,2,0 +'turb_kin_energy ' 'turb_viscosity ' * *************************************************************************** * * -* ----- Number and names of turb eqs ----- * +* ----- Turbulence model constants ----- * * * -NTEQ,I ,1,1,0 -2 +ALPHA,R ,1,1,0 +0.556 * -NAME_UNK,S ,1,2,0 -'turb_kin_energy ' 'turb_omega ' -* -*************************************************************************** -* * -* ----- Number and names of additional turb fields ----- * -* * -NTEX,I ,1,1,0 -1 -* -NAME_EX,S ,1,1,0 -'turb_viscosity ' -* -*************************************************************************** -* * -* ----- Turbulence model constants ----- * -* * -ALPHA,R ,1,1,0 -0.518 -* -ALPHA2,R ,1,1,0 -0.44 -* -BETA,R ,1,1,0 -0.0747 -* -BETA1,R ,1,1,0 -0.0747 -* -BETA2,R ,1,1,0 -0.0828 +BETA,R ,1,1,0 +0.075 * BETAS,R ,1,1,0 0.09 * -SIGMA,R ,1,2,0 -1.1 0.53 -* -SIGMA2,R ,1,2,0 -1.1 1.0 -* -SIGMAD1,R ,1,1,0 -1.0 -* -SIGMAD2,R ,1,1,0 -0.4 -* -CMIX,R ,1,1,0 -1.5 -* -*************************************************************************** -* Diffusion correction of EARSM model -* -CD,R ,1,1,0 -2.2 +SIGMA,R ,1,12,0 +0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 * WALL,R ,1,1,0 10.0 @@ -3861,45 +4107,57 @@ PRT,R ,1,1,0 * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0 -* -C1,R ,1,1,0 -1.8 -* -A0,R ,1,1,0 --0.72 +1.0E20 * *************************************************************************** * * -* ----- Turbulence model subroutines ----- * +* ----- Pressure-strain model constants ----- * * * +PSC10,R ,1,1,0 +4.6 +* +PSC11,R ,1,1,0 +1.24 +* +PSC2,R ,1,1,0 +0.47 +* +PSC2S,R ,1,1,0 +0.0 +* +PSC3,R ,1,1,0 +2.0 +* +PSC4,R ,1,1,0 +0.56 +* +PSC5,R ,1,1,0 +0.0 +* *************************************************************************** * * -* ----- Viscous terms ----- * +* ----- Turbulence model subroutines ----- * * * -FVISF,S ,1,1,0 -'GVISKZHELLST ' -* *************************************************************************** * * * ----- Source terms ----- * * * -GSOUKZ,S ,1,1,0 -'HSSHELLST ' +GSOUKZ,S ,1,2,0 +'HSSTKW ' 'HSSRST ' * *************************************************************************** * * * ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXWJAR ' +'FXXRST ' * *************************************************************************** * * -* --------------- W&J DRSM + std k-omega --------------- * +* --------------- W&J DRSM + Hellsten k-omega --------------- * * * -TURB_MOD,L ,1,1,23 -'W&J DRSM + std k-omega ' +TURB_MOD,L ,1,1,28 +'W&J DRSM + Hellsten k-omega ' * *************************************************************************** * * @@ -3907,15 +4165,9 @@ TURB_MOD,L ,1,1,23 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * +* ----- If wall distance required ----- * * * -K_AVAIL,N ,0,0,0 +WALL_DIST,N ,0,0,0 * *************************************************************************** * * @@ -3935,11 +4187,8 @@ NAME_UNK,S ,1,12,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* NAME_EX,S ,1,2,0 'turb_kin_energy ' 'turb_viscosity ' * @@ -3948,16 +4197,34 @@ NAME_EX,S ,1,2,0 * ----- Turbulence model constants ----- * * * ALPHA,R ,1,1,0 -0.556 +0.518 +* +ALPHA2,R ,1,1,0 +0.44 * BETA,R ,1,1,0 -0.075 +0.0747 +* +BETA2,R ,1,1,0 +0.0828 * BETAS,R ,1,1,0 0.09 * SIGMA,R ,1,12,0 -0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 +1.1 1.1 1.1 1.1 0.53 1.1 1.1 1.1 1.1 1.1 1.1 0.53 +* +SIGMA2,R ,1,12,0 +1.1 1.1 1.1 1.1 1.0 1.1 1.1 1.1 1.1 1.1 1.1 1.0 +* +SIGMAD1,R ,1,1,0 +1.0 +* +SIGMAD2,R ,1,1,0 +0.4 +* +CMIX,R ,1,1,0 +1.5 * WALL,R ,1,1,0 10.0 @@ -4003,10 +4270,17 @@ PSC5,R ,1,1,0 * * *************************************************************************** * * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZHELLST ' +* +*************************************************************************** +* * * ----- Source terms ----- * * * GSOUKZ,S ,1,2,0 -'HSSTKW ' 'HSSRST ' +'HSSHELLST ' 'HSSRST ' * *************************************************************************** * * @@ -4017,10 +4291,10 @@ EXTRVA,S ,1,1,0 * *************************************************************************** * * -* --------------- W&J DRSM + Hellsten k-omega --------------- * +* --------------- W&J DAMPED DRSM + Hellsten k-omega --------------- * * * TURB_MOD,L ,1,1,31 -'W&J DRSM + Hellsten k-omega ' +'W&J DAMPED DRSM + Hellsten k-omega ' * *************************************************************************** * * @@ -4034,18 +4308,6 @@ WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- Turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- If Reynolds_stress exists ----- * * * RS_AVAIL,N ,0,0,0 @@ -4062,11 +4324,8 @@ NAME_UNK,S ,1,12,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* NAME_EX,S ,1,2,0 'turb_kin_energy ' 'turb_viscosity ' * @@ -4142,6 +4401,15 @@ PSC4,R ,1,1,0 PSC5,R ,1,1,0 0.0 * +PSDC10,R ,1,1,0 +0.57 +* +PSDC11,R ,1,1,0 +0.9 +* +PSAPLUS,R ,1,1,0 +26.0 +* *************************************************************************** * * * ----- Turbulence model subroutines ----- * @@ -4158,7 +4426,7 @@ FVISF,S ,1,1,0 * ----- Source terms ----- * * * GSOUKZ,S ,1,2,0 -'HSSHELLST ' 'HSSRST ' +'HSSHELLST ' 'DAMPED_HSSRST ' * *************************************************************************** * * @@ -4171,7 +4439,7 @@ EXTRVA,S ,1,1,0 * * * --------------- SSG DRSM + std k-omega --------------- * * * -TURB_MOD,L ,1,1,23 +TURB_MOD,L ,1,1,20 'SSG DRSM + std k-omega ' * *************************************************************************** @@ -4180,18 +4448,6 @@ TURB_MOD,L ,1,1,23 * * *************************************************************************** * * -* ----- If MUT exists ----- * -* * -MUT_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * -* ----- If turb_kin_energy exists ----- * -* * -K_AVAIL,N ,0,0,0 -* -*************************************************************************** -* * * ----- If Reynolds_stress exists ----- * * * RS_AVAIL,N ,0,0,0 @@ -4208,11 +4464,8 @@ NAME_UNK,S ,1,12,0 * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* NAME_EX,S ,1,2,0 'turb_kin_energy ' 'turb_viscosity ' * @@ -4276,114 +4529,1290 @@ PSC5,R ,1,1,0 * *************************************************************************** * * -* ----- Turbulence model subroutines ----- * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOUKZ,S ,1,2,0 +'HSSTKW ' 'HSSRST ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity and K) ----- * +* * +EXTRVA,S ,1,1,0 +'FXXRST ' +* +*************************************************************************** +* * +* --------------- SSG DRSM + Hellsten k-omega --------------- * +* * +TURB_MOD,L ,1,1,28 +'SSG DRSM + Hellsten k-omega ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- If Reynolds_stress exists ----- * +* * +RS_AVAIL,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs (2D and 3D) ----- * +* * +NTEQ,I ,1,2,0 +5 7 +* +NAME_UNK,S ,1,12,0 +'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'turb_omega ' 'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'RS_13 ' 'RS_23 ' 'turb_omega ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,2,0 +'turb_kin_energy ' 'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +ALPHA,R ,1,1,0 +0.518 +* +ALPHA2,R ,1,1,0 +0.44 +* +BETA,R ,1,1,0 +0.0747 +* +BETA2,R ,1,1,0 +0.0828 +* +BETAS,R ,1,1,0 +0.09 +* +SIGMA,R ,1,12,0 +1.1 1.1 1.1 1.1 0.53 1.1 1.1 1.1 1.1 1.1 1.1 0.53 +* +SIGMA2,R ,1,12,0 +1.1 1.1 1.1 1.1 1.0 1.1 1.1 1.1 1.1 1.1 1.1 1.0 +* +SIGMAD1,R ,1,1,0 +1.0 +* +SIGMAD2,R ,1,1,0 +0.4 +* +CMIX,R ,1,1,0 +1.5 +* +WALL,R ,1,1,0 +10.0 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0E20 +* +*************************************************************************** +* * +* ----- Pressure-strain model constants ----- * +* * +PSC10,R ,1,1,0 +3.4 +* +PSC11,R ,1,1,0 +1.8 +* +PSC2,R ,1,1,0 +0.8 +* +PSC2S,R ,1,1,0 +1.3 +* +PSC3,R ,1,1,0 +1.25 +* +PSC4,R ,1,1,0 +0.4 +* +PSC5,R ,1,1,0 +4.2 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZHELLST ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOUKZ,S ,1,2,0 +'HSSHELLST ' 'HSSRST ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity and K) ----- * +* * +EXTRVA,S ,1,1,0 +'FXXRST ' +* +*************************************************************************** +* * +* --------------- Peng hybrid HYB0 model ---------- * +* * +TURB_MOD,L ,1,1,9 +'Peng hybrid HYB0 model ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +0 +* +NAME_UNK,S ,1,1,0 +'zero_equation ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,2,0 +'turb_viscosity ' 'rans_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CSGS,R ,1,1,0 +0.12 +* +CKSGS,R ,1,1,0 +0.0 +* +SIGMA,R ,1,1,0 +1.0 +* +PRT,R ,1,1,0 +0.4 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXHYB0P ' +* +*************************************************************************** +* * +* --------------- Spalart-Allmaras DES model ---------- * +* * +TURB_MOD,L ,1,1,17 +'Spalart-Allmaras DES model ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDESSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Spalart-Allmaras DDES ---------- * +* * +TURB_MOD,L ,1,1,17 +'Spalart-Allmaras DDES ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDDESSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Spalart-Allmaras IDDES ---------- * +* * +TURB_MOD,L ,1,1,17 +'Spalart-Allmaras IDDES ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSIDDESSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Menter SST-DES k-omega --------------- * +* * +TURB_MOD,L ,1,1,20 +'Menter SST-DES k-omega ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +2 +* +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDESKW,R ,1,1,0 +0.78 +* +CDESKE,R ,1,1,0 +0.61 +* +ALPHA,R ,1,1,0 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 +* +BETA,R ,1,1,0 +0.075 +* +BETA2,R ,1,1,0 +0.0828 +* +BETAS,R ,1,1,0 +0.09 +* +SIGMA,R ,1,2,0 +0.85 0.5 +* +SIGMA2,R ,1,2,0 +1.0 0.856 +* +A1,R ,1,1,0 +0.31 +* +WALL,R ,1,1,0 +10.0 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZBSL ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDESSST ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSST ' +* +*************************************************************************** +* * +* --------------- Menter SST-DDES k-omega --------------- * +* * +TURB_MOD,L ,1,1,20 +'Menter SST-DDES k-omega ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +2 +* +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDESKW,R ,1,1,0 +0.78 +* +CDESKE,R ,1,1,0 +0.61 +* +ALPHA,R ,1,1,0 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 +* +BETA,R ,1,1,0 +0.075 +* +BETA2,R ,1,1,0 +0.0828 +* +BETAS,R ,1,1,0 +0.09 +* +SIGMA,R ,1,2,0 +0.85 0.5 +* +SIGMA2,R ,1,2,0 +1.0 0.856 +* +A1,R ,1,1,0 +0.31 +* +WALL,R ,1,1,0 +10.0 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZBSL ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDDESSST ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSST ' +* +*************************************************************************** +* * +* --------------- Menter SST-IDDES k-omega --------------- * +* * +TURB_MOD,L ,1,1,20 +'Menter SST-IDDES k-omega ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +2 +* +NAME_UNK,S ,1,2,0 +'turb_kin_energy ' 'turb_omega ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDESKW,R ,1,1,0 +0.78 +* +CDESKE,R ,1,1,0 +0.61 +* +ALPHA,R ,1,1,0 +0.55317 +* +ALPHA2,R ,1,1,0 +0.4403547 +* +BETA,R ,1,1,0 +0.075 +* +BETA2,R ,1,1,0 +0.0828 +* +BETAS,R ,1,1,0 +0.09 +* +SIGMA,R ,1,2,0 +0.85 0.5 +* +SIGMA2,R ,1,2,0 +1.0 0.856 +* +A1,R ,1,1,0 +0.31 +* +WALL,R ,1,1,0 +10.0 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISKZBSL ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSIDDESSST ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSST ' +* +*************************************************************************** +* * +* --------------- Rung SALSA-DES ---------- * +* * +TURB_MOD,L ,1,1,17 +'Rung SALSA-DES ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDESSALSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Rung SALSA-DDES ---------- * +* * +TURB_MOD,L ,1,1,17 +'Rung SALSA-DDES ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSDDESSALSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Rung SALSA-IDDES ---------- * +* * +TURB_MOD,L ,1,1,17 +'Rung SALSA-IDDES ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 +* +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 +* +NAME_UNK,S ,1,1,0 +'turb_work_vt ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CDES,R ,1,1,0 +0.65 +* +SACB1,R ,1,1,0 +0.1355 +* +SACB2,R ,1,1,0 +0.622 +* +SACW1,R ,1,1,0 +3.24 +* +SACW2,R ,1,1,0 +0.3 +* +SACW3,R ,1,1,0 +2.0 +* +SACV1,R ,1,1,0 +7.1 +* +SIGMA,R ,1,1,0 +1.5 +* +PRT,R ,1,1,0 +0.9 +* +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 +* +*************************************************************************** +* * +* ----- Turbulence model subroutines ----- * +* * +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' +* +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSIDDESSALSA ' +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSAM ' +* +*************************************************************************** +* * +* --------------- Smagorinsky SGS model ---------- * +* * +TURB_MOD,L ,1,1,8 +'Smagorinsky SGS model ' +* +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +0 +* +NAME_UNK,S ,1,1,0 +'zero_equation ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CSMAG,R ,1,1,0 +0.2 +* +CKSMG,R ,1,1,0 +0.0 +* +SIGMA,R ,1,1,0 +1.0 +* +PRT,R ,1,1,0 +0.4 +* +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity and kin. energy) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSGSSMG ' +* +*************************************************************************** +* * +* --------------- Smagorinsky SGS with wall damping ---------- * * * +TURB_MOD,L ,1,1,9 +'Smagorinsky SGS with wall damping ' +* *************************************************************************** * * -* ----- Source terms ----- * +* ----- Turbulence model program logics ----- * * * -GSOUKZ,S ,1,2,0 -'HSSTKW ' 'HSSRST ' -* *************************************************************************** * * -* ----- Extra fields (turbulent viscosity and K) ----- * +* ----- If wall distance required ----- * * * -EXTRVA,S ,1,1,0 -'FXXRST ' +WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* --------------- SSG DRSM + Hellsten k-omega --------------- * +* ----- Number and names of turb eqs ----- * * * -TURB_MOD,L ,1,1,31 -'SSG DRSM + Hellsten k-omega ' +NTEQ,I ,1,1,0 +0 +* +NAME_UNK,S ,1,1,0 +'zero_equation ' * *************************************************************************** * * -* ----- Turbulence model program logics ----- * +* ----- Names of additional turb fields ----- * * * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* *************************************************************************** * * -* ----- If wall distance required ----- * +* ----- Turbulence model constants ----- * * * -WALL_DIST,N ,0,0,0 +CSMAG,R ,1,1,0 +0.2 +* +CKSMG,R ,1,1,0 +0.0 +* +SIGMA,R ,1,1,0 +1.0 +* +PRT,R ,1,1,0 +0.4 * *************************************************************************** * * -* ----- If MUT exists ----- * +* ----- Extra fields (turbulent viscosity and kin. energy) ----- * * * -MUT_AVAIL,N ,0,0,0 +EXTRVA,S ,1,1,0 +'FXSGSSMG2 ' * *************************************************************************** * * -* ----- If turb_kin_energy exists ----- * +* --------------- Spalart-Allmaras SGS model ---------- * * * -K_AVAIL,N ,0,0,0 +TURB_MOD,L ,1,1,17 +'Spalart-Allmaras SGS model ' * *************************************************************************** * * -* ----- If Reynolds_stress exists ----- * +* ----- Turbulence model program logics ----- * * * -RS_AVAIL,N ,0,0,0 +*************************************************************************** +* * +* ----- If wall distance required ----- * +* * +WALL_DIST,N ,0,0,0 * *************************************************************************** * * -* ----- Number and names of turb eqs (2D and 3D) ----- * +* ----- Number and names of turb eqs ----- * * * -NTEQ,I ,1,2,0 -5 7 +NTEQ,I ,1,1,0 +1 * -NAME_UNK,S ,1,12,0 -'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'turb_omega ' 'RS_11 ' 'RS_22 ' 'RS_33 ' 'RS_12 ' 'RS_13 ' 'RS_23 ' 'turb_omega ' +NAME_UNK,S ,1,1,0 +'turb_work_vt ' * *************************************************************************** * * -* ----- Number and names of additional turb fields ----- * +* ----- Names of additional turb fields ----- * * * -NTEX,I ,1,1,0 -2 -* -NAME_EX,S ,1,2,0 -'turb_kin_energy ' 'turb_viscosity ' +NAME_EX,S ,1,1,0 +'turb_viscosity ' * *************************************************************************** * * * ----- Turbulence model constants ----- * * * -ALPHA,R ,1,1,0 -0.518 -* -ALPHA2,R ,1,1,0 -0.44 -* -BETA,R ,1,1,0 -0.0747 +CDES,R ,1,1,0 +0.65 * -BETA2,R ,1,1,0 -0.0828 +SACB1,R ,1,1,0 +0.1355 * -BETAS,R ,1,1,0 -0.09 +SACB2,R ,1,1,0 +0.622 * -SIGMA,R ,1,12,0 -1.1 1.1 1.1 1.1 0.53 1.1 1.1 1.1 1.1 1.1 1.1 0.53 +SACW1,R ,1,1,0 +3.24 * -SIGMA2,R ,1,12,0 -1.1 1.1 1.1 1.1 1.0 1.1 1.1 1.1 1.1 1.1 1.1 1.0 +SACW2,R ,1,1,0 +0.3 * -SIGMAD1,R ,1,1,0 -1.0 +SACW3,R ,1,1,0 +2.0 * -SIGMAD2,R ,1,1,0 -0.4 +SACV1,R ,1,1,0 +7.1 * -CMIX,R ,1,1,0 +SIGMA,R ,1,1,0 1.5 * -WALL,R ,1,1,0 -10.0 -* PRT,R ,1,1,0 0.9 * @@ -4392,32 +5821,83 @@ PRT,R ,1,1,0 * PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) * PROLIM,R ,1,1,0 -1.0E20 +1.0 * *************************************************************************** * * -* ----- Pressure-strain model constants ----- * +* ----- Turbulence model subroutines ----- * * * -PSC10,R ,1,1,0 -3.4 +*************************************************************************** +* * +* ----- Viscous terms ----- * +* * +FVISF,S ,1,1,0 +'GVISSA ' * -PSC11,R ,1,1,0 -1.8 +*************************************************************************** +* * +* ----- Source terms ----- * +* * +GSOULDES,S ,1,1,0 +'HSLESSA ' * -PSC2,R ,1,1,0 -0.8 +*************************************************************************** +* * +* ----- Extra fields (turbulent viscosity) ----- * +* * +EXTRVA,S ,1,1,0 +'FXSGSSA ' * -PSC2S,R ,1,1,0 -1.3 +*************************************************************************** +* * +* --------------- Yoshizawa SGS ---------- * +* * +TURB_MOD,L ,1,1,11 +'Yoshizawa k-eq SGS model ' * -PSC3,R ,1,1,0 -1.25 +*************************************************************************** +* * +* ----- Turbulence model program logics ----- * +* * +*************************************************************************** +* * +* ----- Number and names of turb eqs ----- * +* * +NTEQ,I ,1,1,0 +1 * -PSC4,R ,1,1,0 +NAME_UNK,S ,1,1,0 +'turb_kin_energy ' +* +*************************************************************************** +* * +* ----- Names of additional turb fields ----- * +* * +NAME_EX,S ,1,1,0 +'turb_viscosity ' +* +*************************************************************************** +* * +* ----- Turbulence model constants ----- * +* * +CKMU,R ,1,1,0 +0.07 +* +CEPS,R ,1,1,0 +1.05 +* +SIGMA,R ,1,1,0 +1.0 +* +PRT,R ,1,1,0 0.4 * -PSC5,R ,1,1,0 -4.2 +*************************************************************************** +* Coefficient to multiply the limit of the turbulent production +* PLIM = PROLIM*RHO*RK*SQRT(PROD/MUTC) +* +PROLIM,R ,1,1,0 +1.0 * *************************************************************************** * * @@ -4428,21 +5908,21 @@ PSC5,R ,1,1,0 * ----- Viscous terms ----- * * * FVISF,S ,1,1,0 -'GVISKZHELLST ' +'GVISLESYO ' * *************************************************************************** * * * ----- Source terms ----- * * * -GSOUKZ,S ,1,2,0 -'HSSHELLST ' 'HSSRST ' +GSOULDES,S ,1,1,0 +'HSLESYO ' * *************************************************************************** * * -* ----- Extra fields (turbulent viscosity and K) ----- * +* ----- Extra fields (turbulent viscosity) ----- * * * EXTRVA,S ,1,1,0 -'FXXRST ' +'FXSGSYO ' * *************************************************************************** * * @@ -4461,43 +5941,45 @@ NLEVEL,I ,1,1,0 5 * *************************************************************************** -* Removement of singletons (YES = 1, NO = 0) +* Coarsening ratio in the prismatic layer ( CRATIOB > 1.0) * -RMSING,I ,1,1,0 -1 +CRATIOB,R ,1,1,0 +4.0 * *************************************************************************** -* Number of prismatic layers +* Reduction of target number of agglomerated cells close to boundaries, +* >=1 if active * -NPLAYER,I ,1,1,0 -40 +NAGGLTBC,I ,1,1,0 +0 * *************************************************************************** -* Coarsening ratio in the prismatic layer ( CRATIOB > 1.0) +* Number of fine layers (nodes) where NAGGLTBC apply, >=1 if active * -CRATIOB,R ,1,1,0 -4.0 +NLAYERTBC,I ,1,1,0 +0 * *************************************************************************** * Number of partitions. NPART > 1 specifies a parallel calculation. * Used by both flow solver and preprocessor * NPART,I ,1,1,0 -__NPART__ +__NPART__ * *************************************************************************** -* Partitioning option -* 1 = Original option with algorithm based on agglomeration -* 2 = Partitioning based on METIS +* Node weight in implicit lines. Default 2000, can be reduced to 1000 +* if LINEIM=0. No effect for Euler meshes or if STR_LINES_OPT>1 * -PTYPE,I ,1,1,0 -2 +METIS_L_V_WEIGHT,I ,1,1,0 +2000 * *************************************************************************** -* Flag for printing the dual grid (YES = 1, NO = 0) +* 1 = Account for stretched lines in partitioning and agglomeration +* 2 = Account for stretched lines in agglomeration only +* 3 = Do not account for stretched lines at all * -PDUAL,I ,1,1,0 -0 +STR_LINES_OPT,I ,1,1,0 +1 * *************************************************************************** * Type of edge coloring : @@ -4540,8 +6022,8 @@ IWALLDIS,I ,1,1,0 * * *************************************************************************** * Adaption parameters taken from file: refinement_data.aref -* 0 = Data in this file used -* 1 = Data in refinement_data.aref used +* 0 = Data in this file used +* 1 = Data in refinement_data.aref used * EXTADAPT,I ,1,1,0 0 @@ -4724,134 +6206,45 @@ CFIBLAD,L ,1,1,0 * specified point in the volume or on a boundary * * * *************************************************************************** -* Number of active probes if >0 +* Probe data history on/off (IPROBEHIS=1/0) +* Probe input data file under CFIPRB * -PROBE_RECORD,I ,1,1,3 +IPROBEHIS,I ,1,1,0 0 * *************************************************************************** -* Coordinates of probes -* -PROBE_COORD,R ,8,3,0 -0.8939345 0.0105114 0.0759227 -0.8939345 -0.0299341 0.0678775 -0.8939345 -0.0528447 0.0335895 -0.8939345 -0.0447995 -0.006856 -0.8939345 -0.0105114 -0.0297666 -0.8939345 0.0299341 -0.0217215 -0.8939345 0.0528447 0.0125665 -0.8939345 0.0447995 0.0530121 -* +* * +* --------------- Sampling surface options ------------- * +* * *************************************************************************** -* Name and location of probes: -* interior = inside the volume -* boundary = on boundary surface +* Aerodynamic coeff history on/off (IACHIS=1/0) +* Aerodynamic coeff history input file name supplied to CFIACI * -PROBE_NAME,L ,8,2,0 -'PD1 ' 'interior ' -'PD2 ' 'interior ' -'PD3 ' 'interior ' -'PD4 ' 'interior ' -'PD5 ' 'interior ' -'PD6 ' 'interior ' -'PD7 ' 'interior ' -'PD8 ' 'interior ' +IACHIS,I ,1,1,0 +0 * *************************************************************************** -* Write probe info to a mesh file: 1=yes,0=no +* Sampling boundary data on/off (ISAMP=1/0) +* Sampling boundary data input file name supplied to CFIBSI * -PMESH,I ,1,1,1 -1 +ISAMP,I ,1,1,0 +0 * *************************************************************************** -* Name of mesh file +* Sampling on supplied mesh (ISAMMSH=1) +* Sampling mesh file name supplied to CFISAMMSH * -CFIPROBE,L ,1,1,0 -'Probe.bmsh ' +ISAMMSH,I ,1,1,0 +0 * *************************************************************************** * * * --------------- Inlet analysis options --------------- * * * *************************************************************************** -* 0 = No inlet analysis -* 1 = Calculate pressure recovery and distortion at specified AIP location -* (AIP=Aerodynamic Interface Plane must be circular) -* -* Note: Extra data is added to the bedg-file in the preprocessor step. -* If any setting is changed, make sure to run the preprocessor again. +* Inlet analysis option on/off (IAIPHIS=1/0) +* Inlet analysis input data file under CFIAIP * -AIP_RECORD,I ,1,1,7 +IAIPHIS,I ,1,1,0 0 * -*************************************************************************** -* Write Aip mesh to file: 1=yes,0=no -* -AIPMESH,I ,1,1,1 -1 -* -*************************************************************************** -* Name of mesh file -* -CFIAIP,L ,1,1,0 -'Aip.bmsh ' -* -*************************************************************************** -* Mid-point of AIP and probe rake -* -AIP_POINT,R ,1,3,0 -4.7217 0.0 2.1 -* -*************************************************************************** -* Radius of AIP and probe rake -* -AIP_RADIUS,R ,1,1,0 -0.33 -* -*************************************************************************** -* AIP and probe rake plane normal direction -* -AIP_NORM,R ,1,3,0 -1.0 0.0 0.0 -* -*************************************************************************** -* Start distortion analysis from this direction and rotate -* clockwise looking from the AIP_NORM vector. This will determine -* the angle of maximum distortion. -* -AIP_ROT,R ,1,3,0 -0.0 0.0 1.0 -* -*************************************************************************** -* Pos 1: Number of radial position of AIP rake -* Pos 2: Number of arms of AIP rake -* Pos 3: Full circle=1, half circle=0 (if using a symmetry plane at AIP) -* -AIP_GRID,I ,1,3,0 -5 5 0 -* -*************************************************************************** -* Evaluate freestream total pressure for AIP reference -* 1 = At specified point in volume (PT_POINT) -* 2 = At specified boundary name (PT_BC, mean values evaluated) -* 3 = Use Free stream values from this input file (PFREE,TFREE...) -* -PT_PROBE,I ,1,1,2 -1 -* -*************************************************************************** -* Coordinates of points for free stream references -* -PT_POINT,R ,1,3,0 -3.0717 0.15 2.4 -* -*************************************************************************** -* Boundary condition name for free stream reference -* -PT_BC,L ,1,1,0 -'BC1_on_INLET ' -* -*************************************************************************** -* * -* ---------------- The end of regular Edge parameters ---------------- * -* *