1# PHASTA Version 1.6 Input File 2# Here you will 3# find the default values for everything that is allowed to have a 4# default. Some things are not allowed to default and must be entered 5# in this file. In most cases, the acceptable inputs also appear in 6# the input.config file. To add acceptable inputs you should only have 7# to modify: 8# input_asci.cc :which matches the string and translates it to the 9# parameter change. If it is a new parameter, you must 10# modify: 11# 12# common.h (to carry new parameter through code) 13# common_c.h (to carry the parameter from C to Fortran) 14# 15## 16 17FLOW CONTROL PARAMETERS 18{ 19#inlet boundary conditions 20 Take BC from IC at Inlet : 0 # iI2Binlet, give the surf ID which want to use taking BC from IC 21 Inlet Bulk x Velocity: NODEFAULT # inletVelX 22#outlet pressure 23 Set Outlet Pressure: 0 # isetOutPres, give the surf ID which want to be set 24 Uniform Outlet Pressure: NODEFAULT # outPres1 25#initial condition 26 Specify Initial Conditions: 0 #isetInitial 27 Initial X Velocity: NODEFAULT #xvel_ini 28 Initial Y Velocity: NODEFAULT #yvel_ini 29 Initial Z Velocity: NODEFAULT #zvel_ini 30 Initial Temp: NODEFAULT #temp_ini 31 Initial Pressure: NODEFAULT #pres_ini 32 Initial Scalar 1: NODEFAULT #evis_ini 33} 34 35SOLUTION CONTROL 36{ 37 38 Equation of State: NODEFAULT 39# Equation of State: Compressible # sets ipress=0 matflag(1,n) 40# Equation of State: Incompressible # sets ipress=-1 matflag(1,n) 41 Viscous Control: Viscous #replaces navier (0 no, 1 yes) 42 Turbulence Model: No-Model # No-Model, RANS-SA, 43 # RANS, RANS-KE, LES, DES97, DDES 44 Number of Timesteps: NODEFAULT #replaces nsteps(1) (ntseq wired =1) 45 Time Step Size: NODEFAULT # Delt(1) 46 Solve Flow : True # solveflow 47 Solve Heat : False # solveheat 48 Solve Scalars : 0 # nsclrS 49 Solve Level Set : 0 # iLSet 50 # total number of scalars must be <=4 51 # impl=10*logical(solflow)+nsclr*100 + 52 # 1*logical(solveheat) 53} Control 54 55MATERIAL PROPERTIES 56{ 57 Shear Law: Constant Viscosity # ishear=0 => matflag(2,n) 58 Bulk Viscosity Law: Constant Bulk Viscosity # ibulk=0 => matflag(3,n) 59 Conductivity Law: Constant Conductivity # icond=0 => matflag(4,n) 60 Viscosity: NODEFAULT # fills datmat (2 values REQUIRED if iLset=1) 61 Density: 1.0 # ditto 62 Thermal Conductivity: 0.2 # ditto 63 Prandtl Number: 0.72 # Compressible code sets diffusivity with this 64 Scalar Diffusivity: 0.2 # fills scdiff(1:nsclrS) 65 Body Force Option: None # ibody=0 => matflag(5,n) 66# Body Force Option: Vector # ibody=1 => matflag(5,n) 67# Body Force Option: Boussinesq # ibody=2 => matflag(5,n) 68# Body Force Option: User e3source.f # ibody=3 => matflag(5,n) 69# Body Force Option: Cooling Analytic # ibody=4 => matflag(5,n) 70# Body Force Option: Cooling Initial Condition # ibody=5 => matflag(5,n) 71 Body Force: 0.0 0.0 0.0 # (datmat(i,5,n),i=1,nsd) 72 Body Force Pressure Gradient: 0.0 0.0 0.0 # (datmat(i,7,n),i=1,nsd) 73 Zero Mean Pressure : False # True causes pressure avg set to zero 74 # ONLY appropriate if no pressure BCs. 75 Rotating Frame of Reference: False 76 Rotating Frame of Reference Rotation Rate: 0. 0. 0. 77 Surface Tension Option: No #isurf=0 78# Surface Tension Option: Yes #isurf=1 79 Maximum Value of Sponge Parameter: 1.0 80 Inflow Cooling Sponge Ends at z: 12.0 81 Outflow Cooling Sponge Begins at z: 24.0 82 Radial Cooling Sponge Begins at r: 7.5 83 Sponge Growth Coefficient Outflow : .1 84 Sponge Growth Coefficient Inflow : .5 85 Sponge for Continuity Equation : False 86 Sponge for x Momentum Equation : False 87 Sponge for y Momentum Equation : False 88 Sponge for z Momentum Equation : False 89 Sponge for Energy Equation : False 90 91# damp vortices near outlet by ramping up viscosity 92 Ramp Up Viscosity Near Outlet: 0 #irampViscOutlet 93 Stretch X Coordinate Near Outlet: 0 #istretchOutlet 94} 95 96OUTPUT CONTROL 97{ 98 Number of Timesteps between Restarts: 500 #replaces nout/ntout 99 Number of SyncIO Files: 1 # nsynciofiles 100 Verbosity Level: 3 #replaces necho 101 Print Statistics: False #False=> ioform=1, True=> ioform=2 102 Print Wall Fluxes: False #No current action but it will come later 103 Print Residual at End of Step: False # T lstres=1 F lstres=0 104 Print Error Indicators: False # F ierrcalc = 0, T ierrcalc = 1 105 Print FieldView: False # outpar.iofieldv 106 Print ybar: False # F ioybar = 0, T ioybar = 1 107 Print vorticity: False # F ivort = 0, T ivort = 1 108 Number of Steps in a Cycle: 0 # nstepsincycle 109 Number of Phases in a Cycle: 0 # nphasesincycle 110 Number of Initial Cycles to Skip in Phase Average: 0 # ncycles_startphaseavg 111 Print Velocity Hessian: False # F ihessian = 0, T ihessian = 1 112 Number of Error Smoothing Iterations: 3 # ierrsmooth 113 Surface ID for Integrated Mass: 1 # isrfIM 114 Number of Force Surfaces: 0 # nsrfCM 115 Surface ID's for Force Calculation: NODEFAULT # nreadlist(j),j=1,nsrfCM 116# this must be processed as in input.f or passed to input.f for processing 117# not read if nsrfCM=0 118 119 Ranks per core: 1 # for varts only 120 Cores per node: 1 # for varts only 121# By default, with these 2 settings, 122# probe 1 will be treated by rank (numpe-1) 123# probe 2 will be treated by rank (numpe-1)-1 124# probe 3 will be treated by rank (numpe-1)-2 125# On BGQ, with 4 and 16 respectively, 126# probe 1 will be treated by rank (numpe-1) 127# probe 2 will be treated by rank (numpe-1)-1*(4*16) on another node 128# probe 3 will be treated by rank (numpe-1)-2*(4*16) on another node 129 130 Data Block Format : binary #iotype, options 'binary','ascii' 131} 132 133LINEAR SOLVER 134{ 135# Solver Type: ACUSIM # iprjFlag=0 ipresPrjFlag=0 136 Solver Type: ACUSIM with P Projection # iprjFlag=0 ipresPrjFlag=1 137# Solver Type: ACUSIM with Velocity Projection # iprjFlag=1 ipresPrjFlag=0 138# Solver Type: ACUSIM with Full Projection # iprjFlag=1 ipresPrjFlag=1 139# The above 4 are for incompressible flow. 140 141# The next two are for compresible flow. 142# Solver Type: GMRES EBE 143# Solver Type: GMRES Matrix Free 144 145 Number of GMRES Sweeps per Solve: 1 # replaces nGMRES 146 Number of Krylov Vectors per GMRES Sweep: 50 # replaces Kspace 147 Number of Solves per Left-hand-side Formation: 1 #nupdat/LHSupd(1) 148 Tolerance on Momentum Equations: 0.1 # epstol(1) 149 Number of Solves of Temperature per Left-hand-side Formation: 1 150 Temperature Solver Tolerance: 0.001 151 Number of Solves of Scalar 1 per Left-hand-side Formation: 1 152 Number of Solves of Scalar 2 per Left-hand-side Formation: 1 153 Number of Solves of Scalar 3 per Left-hand-side Formation: 1 154 Number of Solves of Scalar 4 per Left-hand-side Formation: 1 155 Scalar 1 Solver Tolerance: 0.001 156 Scalar 2 Solver Tolerance: 0.001 157 Scalar 3 Solver Tolerance: 0.001 158 Scalar 4 Solver Tolerance: 0.001 159 Tolerance on ACUSIM Pressure Projection: 0.1 # prestol 160 Minimum Number of Iterations per Nonlinear Iteration: 10 # minIters 161 Maximum Number of Iterations per Nonlinear Iteration: 200 # maxIters 162 Velocity Delta Ratio :0. #utol deltol(1,1) Stop factor for steady solve 163 Pressure Delta Ratio :0. #ptol deltol(1,2) Stop factor for steady solve 164 Number of Velocity Projection Vectors: 10 #nPrjs 165 Number of Pressure Projection Vectors: 10 #nPresPrjs 166 ACUSIM Verbosity Level : 0 #iverbose 167} 168 169DISCRETIZATION CONTROL 170{ 171 Basis Function Order: 1 # ipord 172 Time Integration Rule: First Order # 1st Order sets rinf(1) -1 173# Time Integration Rule: Second Order # Second Order sets rinf next 174# Time Integration Rho Infinity: 0.5 # rinf(1) Only used for 2nd order 175 Predictor at Start of Step : Same Velocity # ipred=1 (more options later) 176 Weak Form: SUPG # alternate is Galerkin only for compressible 177 Flow Advection Form: Convective # iconvflow=2 178# Flow Advection Form: Conservative # iconvflow=1 179 Scalar Advection Form: Convective # iconvsclr=2 180# Scalar Advection Form: Conservative # iconvsclr=1 181# Use Conservative Scalar Convection Velocity: True 182 Use Conservative Scalar Convection Velocity: False 183 Tau Matrix: Diagonal-Shakib #itau=0 184# Tau Matrix: Diagonal-Franca #itau=1 185# Tau Matrix: Diagonal-Jansen(dev) #itau=2 186# Tau Matrix: Diagonal-Compressible #itau=3 187# Tau Matrix: Matrix-Mallet #itau=10 188 Tau Time Constant: 1. #dtsfct 189 Tau C Scale Factor: 1.0 # taucfct best value depends # on Tau Matrix chosen 190 Remove Time Term from Stabilization: 1 #iremoveStabTimeTerm 191 Discontinuity Capturing: Off # Sets IDC to 0 for now 192# Discontinuity Capturing: "DC-mallet" #Sets IDC to 1 193 Scalar Discontinuity Capturing: 0 0 #Sets idcsclr to [0 0], no DC 194 #on any scalar 195# idcsclr(1)--> Type of DC (see flow), idcsclr(2)---> on which scalar DC acting 196# Scalar Discontinuity Capturing: 1 1 #Sets DC=1 on first scalar 197# Scalar Discontinuity Capturing: 1 2 #Sets DC=1 on second scalar 198 Include Viscous Correction in Stabilization: True # if p=1 idiff=1 199 # if p=2 idiff=2 200 Lumped Mass Fraction on Left-hand-side: 0. # flmpl 201 Lumped Mass Fraction on Right-hand-side: 0. # flmpr 202 Dump CFL: False #iCFLworst=0 203 Quadrature Rule on Interior: 2 #int(1) 204 Quadrature Rule on Boundary: 2 #intb(1) 205 Number of Elements Per Block: 64 #ibksiz 206 Entropy Form of Pressure Constraint on Weight Space: 0 # 1 turns it on 207} 208 209SOLUTION SCALING PARAMETERS 210{ 211 Density: 1. #ro This is used in sponge 212 Velocity: 1. #vel This affects tau currently 213 Pressure: 1. #pres This is used in sponge 214 Temperature: 1. #temper This scales diagonal energy tau (see e3tau.f) 215 Entropy: 1. #entrop 216} 217 218TURBULENCE MODELING PARAMETERS 219{ # lines below are only read if ||| is true 220 Turn Off Source Terms for Scalars: False 221 Decay Multiplier for Scalars : 1.0 (this number multiplies scalars each step) 222 Number of Homogenous Directions : 0 223 Dynamic Model Type : Standard # adds zero to iles LES 224# Dynamic Model Type : Bardina # adds 10 to iles LES 225# Dynamic Model Type : Projection # adds 20 to iles LES 226 Filter Integration Rule: 1 #ifrule adds ifrule-1 to iles LES 227 Double Filter : False #Filter applied twice to create wide filter? 228 Model Statistics: False #Collect statistics on model? 229 Model/SUPG Dissipation : False # Get diss. due to model and SUPG stresses? 230 231 DES Edge Length: NODEFAULT # shorter edges are pure les while twice longer edges are pure RANS, finally in between is a blend 232 DES SA Minimum Edge Length: 0.0 # limit edge length for DES based on SA model 233 234# 235# The following are models under development (at various stages of validation) 236# 237 Dynamic Sub-Model Type: None 238# Dynamic Sub-Model Type: DFWR # dmc with dynamic filter width ratio 239# Dynamic Sub-Model Type: SUPG # dmc w/ SUPG interation 240# Dynamic Sub-Model Type: ConsistentProj # Consistent projection filter dmc 241 Projection Filter Type: Linear # Projection onto lins. as the filter 242# Projection Filter Type: Quadratic # Projection onto quads as the filter. 243 Lumping Factor for Filter : 0.5 # Lumping factor for projection filter 244 Turbulence Wall Model Type: None #itwmod=0 RANSorLES 245# Turbulence Wall Model Type: Slip Velocity #itwmod=1 RANSorLES 246# Turbulence Wall Model Type: Effective Viscosity #itwmod=2 RANSorLES 247 Velocity Averaging Steps : 500. # wtavei= 1/this RANSorLES 248 Dynamic Model Averaging Steps : 500. # dtavei= 1/this LES 249# negative values to the two previous entries make their value ISTEP in code 250# Anil...leave as any negative value 251 Filter Width Ratio : 3. # fwr1 LES 252 Target Viscosity For Step NSTEP: 0 # uses linear ramp between start 253 # with a zero here the code disables 254 # that feature and holds value 255 Limit u1 : 0. 0. 0. # switch min max change switch from zero to activate 256 Limit u2 : 0 0 0 257 Limit u3 : 0 0 0 258 Limit Pressure : 0 0 0 259 Limit Temperature : 0 0 0 260 Limit Scalar 1 : 0 0 0 261 Limit Scalar 2 : 0 0 0 262 Limit Scalar 3 : 0 0 0 263 Limit Scalar 4 : 0 0 0 264 265 Number of Father Nodes: 0 # value of SONFATH 266 267 Ramp Inflow : False 268 Mdot Ramp Inflow Start and Stop : 0.0 0.0 269 Mdot Ramp Lower FC Start and Stop : 0.0 0.0 270 Mdot Ramp Upper FC Start and Stop : 0.0 0.0 271 272} 273 274SPEBC MODELING PARAMETERS 275{ 276 SPEBC Model Active:-1 # irscale=-1 means not active, see genscale.f 277# Next lines required only if above not non-negative 278 Interpolate Pressure: 0 # intpres 279 Distance between Planes: NODEFAULT # plandist 280 Theta Angle of Arc: 0.0 # thetag = 0.0 for cartisian case 281 Distance for Velocity Averaging: NODEFAULT # ds 282# Tolerance for cylindrical case to account for discratisation of 283# the curved surface 284 SPEBC Cylindrical Tolerance: 0.01 # tolerence 285 Radius of recycle plane: 0.5 #radcyl only needed in cylindrical case 286 Inlet Boundary Layer Thickness: NODEFAULT # rbltin 287# Velocity Averaging Steps needed for SPEBC 288 Vertical Velocity Scale Factor: NODEFAULT # rvscal 289} 290CARDIOVASCULAR MODELING PARAMETERS 291{ 292 Time Varying Boundary Conditions From File: False # F itvn=0 T itvn=1 293 BCT Time Scale Factor : 1.0 294 Number of Coupled Surfaces: 0 # icardio 295 Pressure Coupling: None # Explicit, Implicit, P-Implicit 296 # none ipvsq=0, expl ipvsq=1, 297 # impl ipvsq=2, P-Imp ipvsq=3 298 Number of Resistance Surfaces: 0 # numResistSrfs 299 List of Resistance Surfaces: NODEFAULT # nsrflistResist(j), j=0,MAXSURF 300 Resistance Values : NODEFAULT # ValueListResist(j),j=1,icardio 301 Number of Impedance Surfaces: 0 # numImpSrfs 302 List of Impedance Surfaces: NODEFAULT # nsrflistImp(j), j=0,MAXSURF 303 Impedance From File: False #False impfile=0, True impfile=1 304 Number of RCR Surfaces: 0 # numRCRSrfs 305 List of RCR Surfaces: NODEFAULT # nsrflistRCR(j), j=0,MAXSURF 306 RCR Values From File: False #False ircrfile=0, True ircrfile=1 307 Deformable Wall: False #False ideformwall=0, True ideformwall=1 308 Density of Vessel Wall: NODEFAULT # rhovw 309 Thickness of Vessel Wall: NODEFAULT # thicknessvw 310 Young Mod of Vessel Wall: NODEFAULT # evw 311 Poisson Ratio of Vessel Wall: 0.5 # rnuvw 312 Shear Constant of Vessel Wall: NODEFAULT # rshearconstantvw 313 Wall Mass Matrix for LHS: True # iwallmassfactor=1 314# Wall Mass Matrix for LHS: False # iwallmassfactor=0 315 Wall Stiffness Matrix for LHS: True # iwallstiffactor=1 316# Wall Stiffness Matrix for LHS: False # iwallstiffactor=0 317 Viscous Flux Flag: True # iviscflux=1 318# Viscous Flux Flag: False # iviscflux=0 319 } 320 321 322LEVEL SET MODELING PARAMETERS 323{ 324#(Only read if Solve Level Set is true=> iLS.ne.0) 325 Number of Elements Across Interface: 3 #epsilon_ls 326 Number of Elements Across Interface for Redistancing: 3 #epsilon_lsd 327 Apply Volume Constraint: False #ivconstraint=0 328# Apply Volume Constraint: True #ivconstraint=1 329 Pseudo Time step for Redistancing: 0.001 #dtset 330 Explicit Solve for Redistance Field: 0 #iExpLSSclr2 (0-implicit, 1-explicit) 331 Explicit Solve for Scalar 1 Field: 0 #iExpLSSclr1 (0-implicit, 1-explicit) 332 } 333 334DISABLED FEATURES 335{ 336 iALE : 0 337 icoord: 0 338 irs : 2 339 iexec : 1 340 ntseq : 1 341 imap : 0 342# ivart : 2 # the old practice of using ivart to 343# # set advective form is depricated 344 CFLfl : 1 #CFLfl(1) 345 CFLsl : 10 #CFLsl(1) 346 iepstm: 10 347} 348 349 350 351STEP SEQUENCE 352{ 353 354 Step Construction : 0 1 355} 356 357#AMG PARAMETERS < Note: you will need to compile the code with AMG=1 as an exported environment variable to link in PARPACK and have the following options available to you> 358{ 359# 360 Employ AMG: False # True-1/False-0 (irun_amg = 1, see solfar.f) 361 Run AMG As CG-preconditioner: 0 # 0: no-run. 362 # 1: always preconditioned by AMG. 363 # 2: first try plain CG, if hit plateau, restart CG with AMG. 364 # 3: first try plain CG, if hit maximum iteration, restart CG with AMG. 365 Strong Criterion Eps: 0.25 # "strong criterion" Stuben constant (strong_eps) 366 AMG Convergence Eps: 1.0E-7 # AMG convergence eps 367 AMG Verbosity: 10 # set to 1 for output/0 for no 368 AMG Neg_Sten: 1 # set to 1 for only negs, 0 for pos and neg coeff's 369 AMG Nlevel: 10 # number of levels for tAMG cycle, 2-V cycle (MAX IS 10) 370 AMG Coarsest Solver : 1 # 0: smoother 1: smoother to solve , 2:direct 371 AMG Relaxation Omega: -1 # Omega Relaxation Factor GS / Jac, <0 for GS, >0 for 1/omega of Jacobi 372 AMG Freeze Setup: 10000 # how many solves per setup of AMG 373 AMG Interpolation Type: Direct # Standard or Direct 374 AMG Truncation Set: 0.5 # Interpolation/Prolongation truncate if less than value times diagonal 375 AMG GGB nev: 10 # number of eigenvalues in GGB projection (-1 means no GGB) 376 AMG GGB ncv: 30 # number of memory slots allocated for parpack setup 377# AMG Smoother Type: Gauss-Seidel # Specify smoother 378# AMG Smoother Type: ChebyShev # Specify smoother 379 AMG Smoother Type: MLS # Specify smoother 380 AMG Chebyshev Eigenvalue ratio: 0.1 # Eigen ratio for Chebyshev smoothing: smoothing for ev from beta-ratio*beta 381 AMG MLS Degree: 2 # polynomial smoothing degree {1,2,3,4}, for MLS and Chebyshev only 382 AMG PPE Scale: 2 # 0: original, 1: by diag(C), 2: by diag(PPE) 383 AMG Run Reduced Serial: 0 # <DEPRICATED> 0: no run, N: reduced serial if run serial; output information if parallel. 384} 385 386