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