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 Verbosity Level: 3 #replaces necho 100 Print Statistics: False #False=> ioform=1, True=> ioform=2 101 Print Wall Fluxes: False #No current action but it will come later 102 Print Residual at End of Step: False # T lstres=1 F lstres=0 103 Print Error Indicators: False # F ierrcalc = 0, T ierrcalc = 1 104 Print FieldView: False # outpar.iofieldv 105 Print ybar: False # F ioybar = 0, T ioybar = 1 106 Print vorticity: False # F ivort = 0, T ivort = 1 107 Number of Steps in a Cycle: 0 # nstepsincycle 108 Number of Phases in a Cycle: 0 # nphasesincycle 109 Number of Initial Cycles to Skip in Phase Average: 0 # ncycles_startphaseavg 110 Print Velocity Hessian: False # F ihessian = 0, T ihessian = 1 111 Number of Error Smoothing Iterations: 3 # ierrsmooth 112 Surface ID for Integrated Mass: 1 # isrfIM 113 Number of Force Surfaces: 0 # nsrfCM 114 Surface ID's for Force Calculation: NODEFAULT # nreadlist(j),j=1,nsrfCM 115# this must be processed as in input.f or passed to input.f for processing 116# not read if nsrfCM=0 117 118 Ranks per core: 1 # for varts only 119 Cores per node: 1 # for varts only 120# By default, with these 2 settings, 121# probe 1 will be treated by rank (numpe-1) 122# probe 2 will be treated by rank (numpe-1)-1 123# probe 3 will be treated by rank (numpe-1)-2 124# On BGQ, with 4 and 16 respectively, 125# probe 1 will be treated by rank (numpe-1) 126# probe 2 will be treated by rank (numpe-1)-1*(4*16) on another node 127# probe 3 will be treated by rank (numpe-1)-2*(4*16) on another node 128 129 Data Block Format : binary #iotype, options 'binary','ascii' 130} 131 132LINEAR SOLVER 133{ 134# Solver Type: ACUSIM # iprjFlag=0 ipresPrjFlag=0 135 Solver Type: ACUSIM with P Projection # iprjFlag=0 ipresPrjFlag=1 136# Solver Type: ACUSIM with Velocity Projection # iprjFlag=1 ipresPrjFlag=0 137# Solver Type: ACUSIM with Full Projection # iprjFlag=1 ipresPrjFlag=1 138# The above 4 are for incompressible flow. 139 140# The next two are for compresible flow. 141# Solver Type: GMRES EBE 142# Solver Type: GMRES Matrix Free 143 144 Number of GMRES Sweeps per Solve: 1 # replaces nGMRES 145 Number of Krylov Vectors per GMRES Sweep: 50 # replaces Kspace 146 Number of Solves per Left-hand-side Formation: 1 #nupdat/LHSupd(1) 147 Tolerance on Momentum Equations: 0.1 # epstol(1) 148 Number of Solves of Temperature per Left-hand-side Formation: 1 149 Temperature Solver Tolerance: 0.001 150 Number of Solves of Scalar 1 per Left-hand-side Formation: 1 151 Number of Solves of Scalar 2 per Left-hand-side Formation: 1 152 Number of Solves of Scalar 3 per Left-hand-side Formation: 1 153 Number of Solves of Scalar 4 per Left-hand-side Formation: 1 154 Scalar 1 Solver Tolerance: 0.001 155 Scalar 2 Solver Tolerance: 0.001 156 Scalar 3 Solver Tolerance: 0.001 157 Scalar 4 Solver Tolerance: 0.001 158 Tolerance on ACUSIM Pressure Projection: 0.1 # prestol 159 Minimum Number of Iterations per Nonlinear Iteration: 10 # minIters 160 Maximum Number of Iterations per Nonlinear Iteration: 200 # maxIters 161 Velocity Delta Ratio :0. #utol deltol(1,1) Stop factor for steady solve 162 Pressure Delta Ratio :0. #ptol deltol(1,2) Stop factor for steady solve 163 Number of Velocity Projection Vectors: 10 #nPrjs 164 Number of Pressure Projection Vectors: 10 #nPresPrjs 165 ACUSIM Verbosity Level : 0 #iverbose 166} 167 168DISCRETIZATION CONTROL 169{ 170 Basis Function Order: 1 # ipord 171 Time Integration Rule: First Order # 1st Order sets rinf(1) -1 172# Time Integration Rule: Second Order # Second Order sets rinf next 173# Time Integration Rho Infinity: 0.5 # rinf(1) Only used for 2nd order 174 Predictor at Start of Step : Same Velocity # ipred=1 (more options later) 175 Weak Form: SUPG # alternate is Galerkin only for compressible 176 Flow Advection Form: Convective # iconvflow=2 177# Flow Advection Form: Conservative # iconvflow=1 178 Scalar Advection Form: Convective # iconvsclr=2 179# Scalar Advection Form: Conservative # iconvsclr=1 180# Use Conservative Scalar Convection Velocity: True 181 Use Conservative Scalar Convection Velocity: False 182 Tau Matrix: Diagonal-Shakib #itau=0 183# Tau Matrix: Diagonal-Franca #itau=1 184# Tau Matrix: Diagonal-Jansen(dev) #itau=2 185# Tau Matrix: Diagonal-Compressible #itau=3 186# Tau Matrix: Matrix-Mallet #itau=10 187 Tau Time Constant: 1. #dtsfct 188 Tau C Scale Factor: 1.0 # taucfct best value depends # on Tau Matrix chosen 189 Remove Time Term from Stabilization: 1 #iremoveStabTimeTerm 190 Discontinuity Capturing: Off # Sets IDC to 0 for now 191# Discontinuity Capturing: "DC-mallet" #Sets IDC to 1 192 Scalar Discontinuity Capturing: 0 0 #Sets idcsclr to [0 0], no DC 193 #on any scalar 194# idcsclr(1)--> Type of DC (see flow), idcsclr(2)---> on which scalar DC acting 195# Scalar Discontinuity Capturing: 1 1 #Sets DC=1 on first scalar 196# Scalar Discontinuity Capturing: 1 2 #Sets DC=1 on second scalar 197 Include Viscous Correction in Stabilization: True # if p=1 idiff=1 198 # if p=2 idiff=2 199 Lumped Mass Fraction on Left-hand-side: 0. # flmpl 200 Lumped Mass Fraction on Right-hand-side: 0. # flmpr 201 Dump CFL: False #iCFLworst=0 202 Quadrature Rule on Interior: 2 #int(1) 203 Quadrature Rule on Boundary: 2 #intb(1) 204 Number of Elements Per Block: 64 #ibksiz 205 Entropy Form of Pressure Constraint on Weight Space: 0 # 1 turns it on 206} 207 208SOLUTION SCALING PARAMETERS 209{ 210 Density: 1. #ro This is used in sponge 211 Velocity: 1. #vel This affects tau currently 212 Pressure: 1. #pres This is used in sponge 213 Temperature: 1. #temper This scales diagonal energy tau (see e3tau.f) 214 Entropy: 1. #entrop 215} 216 217TURBULENCE MODELING PARAMETERS 218{ # lines below are only read if ||| is true 219 Turn Off Source Terms for Scalars: False 220 Decay Multiplier for Scalars : 1.0 (this number multiplies scalars each step) 221 Number of Homogenous Directions : 0 222 Dynamic Model Type : Standard # adds zero to iles LES 223# Dynamic Model Type : Bardina # adds 10 to iles LES 224# Dynamic Model Type : Projection # adds 20 to iles LES 225 Filter Integration Rule: 1 #ifrule adds ifrule-1 to iles LES 226 Double Filter : False #Filter applied twice to create wide filter? 227 Model Statistics: False #Collect statistics on model? 228 Model/SUPG Dissipation : False # Get diss. due to model and SUPG stresses? 229 230 DES Edge Length: NODEFAULT # shorter edges are pure les while twice longer edges are pure RANS, finally in between is a blend 231 DES SA Minimum Edge Length: 0.0 # limit edge length for DES based on SA model 232 233# 234# The following are models under development (at various stages of validation) 235# 236 Dynamic Sub-Model Type: None 237# Dynamic Sub-Model Type: DFWR # dmc with dynamic filter width ratio 238# Dynamic Sub-Model Type: SUPG # dmc w/ SUPG interation 239# Dynamic Sub-Model Type: ConsistentProj # Consistent projection filter dmc 240 Projection Filter Type: Linear # Projection onto lins. as the filter 241# Projection Filter Type: Quadratic # Projection onto quads as the filter. 242 Lumping Factor for Filter : 0.5 # Lumping factor for projection filter 243 Turbulence Wall Model Type: None #itwmod=0 RANSorLES 244# Turbulence Wall Model Type: Slip Velocity #itwmod=1 RANSorLES 245# Turbulence Wall Model Type: Effective Viscosity #itwmod=2 RANSorLES 246 Velocity Averaging Steps : 500. # wtavei= 1/this RANSorLES 247 Dynamic Model Averaging Steps : 500. # dtavei= 1/this LES 248# negative values to the two previous entries make their value ISTEP in code 249# Anil...leave as any negative value 250 Filter Width Ratio : 3. # fwr1 LES 251 Target Viscosity For Step NSTEP: 0 # uses linear ramp between start 252 # with a zero here the code disables 253 # that feature and holds value 254 Limit u1 : 0. 0. 0. # switch min max change switch from zero to activate 255 Limit u2 : 0 0 0 256 Limit u3 : 0 0 0 257 Limit Pressure : 0 0 0 258 Limit Temperature : 0 0 0 259 Limit Scalar 1 : 0 0 0 260 Limit Scalar 2 : 0 0 0 261 Limit Scalar 3 : 0 0 0 262 Limit Scalar 4 : 0 0 0 263 264 Number of Father Nodes: 0 # value of SONFATH 265 266 Ramp Inflow : False 267 Mdot Ramp Inflow Start and Stop : 0.0 0.0 268 Mdot Ramp Lower FC Start and Stop : 0.0 0.0 269 Mdot Ramp Upper FC Start and Stop : 0.0 0.0 270 271} 272 273SPEBC MODELING PARAMETERS 274{ 275 SPEBC Model Active:-1 # irscale=-1 means not active, see genscale.f 276# Next lines required only if above not non-negative 277 Interpolate Pressure: 0 # intpres 278 Distance between Planes: NODEFAULT # plandist 279 Theta Angle of Arc: 0.0 # thetag = 0.0 for cartisian case 280 Distance for Velocity Averaging: NODEFAULT # ds 281# Tolerance for cylindrical case to account for discratisation of 282# the curved surface 283 SPEBC Cylindrical Tolerance: 0.01 # tolerence 284 Radius of recycle plane: 0.5 #radcyl only needed in cylindrical case 285 Inlet Boundary Layer Thickness: NODEFAULT # rbltin 286# Velocity Averaging Steps needed for SPEBC 287 Vertical Velocity Scale Factor: NODEFAULT # rvscal 288} 289CARDIOVASCULAR MODELING PARAMETERS 290{ 291 Time Varying Boundary Conditions From File: False # F itvn=0 T itvn=1 292 BCT Time Scale Factor : 1.0 293 Number of Coupled Surfaces: 0 # icardio 294 Pressure Coupling: None # Explicit, Implicit, P-Implicit 295 # none ipvsq=0, expl ipvsq=1, 296 # impl ipvsq=2, P-Imp ipvsq=3 297 Number of Resistance Surfaces: 0 # numResistSrfs 298 List of Resistance Surfaces: NODEFAULT # nsrflistResist(j), j=0,MAXSURF 299 Resistance Values : NODEFAULT # ValueListResist(j),j=1,icardio 300 Number of Impedance Surfaces: 0 # numImpSrfs 301 List of Impedance Surfaces: NODEFAULT # nsrflistImp(j), j=0,MAXSURF 302 Impedance From File: False #False impfile=0, True impfile=1 303 Number of RCR Surfaces: 0 # numRCRSrfs 304 List of RCR Surfaces: NODEFAULT # nsrflistRCR(j), j=0,MAXSURF 305 RCR Values From File: False #False ircrfile=0, True ircrfile=1 306 Deformable Wall: False #False ideformwall=0, True ideformwall=1 307 Density of Vessel Wall: NODEFAULT # rhovw 308 Thickness of Vessel Wall: NODEFAULT # thicknessvw 309 Young Mod of Vessel Wall: NODEFAULT # evw 310 Poisson Ratio of Vessel Wall: 0.5 # rnuvw 311 Shear Constant of Vessel Wall: NODEFAULT # rshearconstantvw 312 Wall Mass Matrix for LHS: True # iwallmassfactor=1 313# Wall Mass Matrix for LHS: False # iwallmassfactor=0 314 Wall Stiffness Matrix for LHS: True # iwallstiffactor=1 315# Wall Stiffness Matrix for LHS: False # iwallstiffactor=0 316 Viscous Flux Flag: True # iviscflux=1 317# Viscous Flux Flag: False # iviscflux=0 318 } 319 320 321LEVEL SET MODELING PARAMETERS 322{ 323#(Only read if Solve Level Set is true=> iLS.ne.0) 324 Number of Elements Across Interface: 3 #epsilon_ls 325 Number of Elements Across Interface for Redistancing: 3 #epsilon_lsd 326 Apply Volume Constraint: False #ivconstraint=0 327# Apply Volume Constraint: True #ivconstraint=1 328 Pseudo Time step for Redistancing: 0.001 #dtset 329 Explicit Solve for Redistance Field: 0 #iExpLSSclr2 (0-implicit, 1-explicit) 330 Explicit Solve for Scalar 1 Field: 0 #iExpLSSclr1 (0-implicit, 1-explicit) 331 } 332 333DISABLED FEATURES 334{ 335 iALE : 0 336 icoord: 0 337 irs : 2 338 iexec : 1 339 ntseq : 1 340 imap : 0 341# ivart : 2 # the old practice of using ivart to 342# # set advective form is depricated 343 CFLfl : 1 #CFLfl(1) 344 CFLsl : 10 #CFLsl(1) 345 iepstm: 10 346} 347 348 349 350STEP SEQUENCE 351{ 352 353 Step Construction : 0 1 354} 355 356#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> 357{ 358# 359 Employ AMG: False # True-1/False-0 (irun_amg = 1, see solfar.f) 360 Run AMG As CG-preconditioner: 0 # 0: no-run. 361 # 1: always preconditioned by AMG. 362 # 2: first try plain CG, if hit plateau, restart CG with AMG. 363 # 3: first try plain CG, if hit maximum iteration, restart CG with AMG. 364 Strong Criterion Eps: 0.25 # "strong criterion" Stuben constant (strong_eps) 365 AMG Convergence Eps: 1.0E-7 # AMG convergence eps 366 AMG Verbosity: 10 # set to 1 for output/0 for no 367 AMG Neg_Sten: 1 # set to 1 for only negs, 0 for pos and neg coeff's 368 AMG Nlevel: 10 # number of levels for tAMG cycle, 2-V cycle (MAX IS 10) 369 AMG Coarsest Solver : 1 # 0: smoother 1: smoother to solve , 2:direct 370 AMG Relaxation Omega: -1 # Omega Relaxation Factor GS / Jac, <0 for GS, >0 for 1/omega of Jacobi 371 AMG Freeze Setup: 10000 # how many solves per setup of AMG 372 AMG Interpolation Type: Direct # Standard or Direct 373 AMG Truncation Set: 0.5 # Interpolation/Prolongation truncate if less than value times diagonal 374 AMG GGB nev: 10 # number of eigenvalues in GGB projection (-1 means no GGB) 375 AMG GGB ncv: 30 # number of memory slots allocated for parpack setup 376# AMG Smoother Type: Gauss-Seidel # Specify smoother 377# AMG Smoother Type: ChebyShev # Specify smoother 378 AMG Smoother Type: MLS # Specify smoother 379 AMG Chebyshev Eigenvalue ratio: 0.1 # Eigen ratio for Chebyshev smoothing: smoothing for ev from beta-ratio*beta 380 AMG MLS Degree: 2 # polynomial smoothing degree {1,2,3,4}, for MLS and Chebyshev only 381 AMG PPE Scale: 2 # 0: original, 1: by diag(C), 2: by diag(PPE) 382 AMG Run Reduced Serial: 0 # <DEPRICATED> 0: no run, N: reduced serial if run serial; output information if parallel. 383} 384 385