1# PHASTA Version 1.5 Input File 2# 3# Basic format is 4# 5# Key Phrase : Acceptable Value (integer, double, logical, or phrase 6# list of integers, list of doubles ) 7# 8 9# To find the Key Phrases, first look in input.config. There you will 10# find the default values for everything that is allowed to have a 11# default. Some things are not allowed to default and must be entered 12# in this file. In most cases, the acceptable inputs also appear in 13# the input.config file. To add acceptable inputs you should only have 14# to modify: 15 16# input_asci.cc :which matches the string and translates it to the 17# parameter change. If it is a new parameter, you must 18# modify: 19# 20# common.h (to carry new parameter through code) 21# common_c.h (to carry the parameter from C to Fortran) 22# 23# 24# In case it is not clear by now, # allows you to comment either from the 25# beginning of the line or to the right as shown below. 26# 27#SOLUTION CONTROL <--- These are for your organizational clarity (not required) 28#{ <--- 29 Equation of State: Incompressible # sets ipress=-1 matflag(1,n) 30 Number of Timesteps: 20 #replaces nsteps(1) (ntseq wired =1) 31 Time Step Size: 0.1 # Delt(1) 32# Turbulence Model: RANS # No-Model iturb=0, RANS =-1 LES=1 33#} 34 35Print Error Indicators: True 36#MATERIAL PROPERTIES 37#{ 38 Viscosity: 0.01 # fills datmat (2 values REQUIRED if iLset=1) 39 Density: 1.0 # ditto 40# Scalar Diffusivity: 0.1 # fills scdiff(1:nsclrS) 41#} 42 43 44#LINEAR SOLVER 45#{ 46 Number of Solves per Left-hand-side Formation: 2 #nupdat/LHSupd(1) 47#} 48 49#DISCRETIZATION CONTROL 50#{ 51 Basis Function Order: 1 # ipord 52 Quadrature Rule on Interior: 2 #int(1) 53 Quadrature Rule on Boundary: 2 #intb(1) 54 Include Viscous Correction in Stabilization: True # if p=1 idiff=1 55 # if p=2 idiff=2 56 Lumped Mass Fraction on Left-hand-side: 1.0 # flmpl 57 Lumped Mass Fraction on Right-hand-side: 1.0 # flmpr 58#} 59 Surface ID for Integrated Mass: 2 60 Number of Force Surfaces: 1 61 Surface ID's for Force Calculation: 1 62 63TURBULENCE MODELING PARAMETERS 64{ # lines below are only read if ||| is true 65 Dynamic Model Type : Standard # adds zero to iturb LES 66 Filter Integration Rule: 1 #ifrule adds ifrule-1 to iturb LES 67# Turbulence Wall Model Type: None #itwmod=0 RANSorLES 68# Turbulence Wall Model Type: Slip Velocity #itwmod=1 RANSorLES 69 Turbulence Wall Model Type: Effective Viscosity #itwmod=2 RANSorLES 70 Velocity Averaging Steps : 500. # wtavei= 1/this RANSorLES 71 Dynamic Model Averaging Steps : 500. # dtavei= 1/this LES 72# negative values to the two previous entries make their value ISTEP in code 73# Anil...leave as any negative value 74 Filter Width Ratio : 3. # fwr1 LES 75 } 76# 77# 78#This last one is brand new. It allows you to construct your step 79#from elementary operations. It works under the premise that a step is 80#constructed from from a series of solves and updates. The table goes like 81#this: 82# solve flow = 0; update flow = 1 83# solve scalar 1 = 10; update scalar 1 = 11 84# solve scalar 2 = 20; update scalar 2 = 21 85# solve scalar 3 = 30; update scalar 3 = 31 86# : : 87# solve scalar n = n*10; update scalar n = n*10+1 88# solve heat = (n+1)*10; update T = (n+1)*10+1 89# 90 91# Below we have an example of solving the flow with two iterations 92# (solve, update, solve,update) what would have been achieve before by 93# setting niter=2 94 95# 96#STEP SEQUENCE 97#{ 98# Step Construction : 0 1 0 1 # this is the standard two iteration 99# Step Construction : 0 1 0 1 0 1 0 1 0 1 100 Step Construction : 0 1 # this is the standard one iteration 101# Step Construction : 0 1 10 11 102# Step Construction : 0 1 10 11 0 1 10 11 0 1 10 11 20 21 20 21 20 21 103# This one is one Tony might like where the solver would solve the flow with 104# the first scalar (3 times with an update immediately after each solve) 105# followed by 3 successive solves of the second scalar (with an update after 106# each solve) 107# 108 109 110# NOTE: An update consists of adding the delta from the appropriate 111# solve to the appropriate part of the Y vector, followed by 112# reapplication of the boundary conditions. You have to ask for an 113# update. It is not implicit that it will follow each solve. This is to 114# enable you to solve successive pieces BEFORE updating the solution as 115# is sometimes convenient/necessary. Here is an example of that 116 117# Step Construction : 0 1 10 20 11 21 0 1 10 20 11 21 118# here we solve the flow, update the flow, solve scalar 1, solve scalar 2 (BEFORE UPDATING SCALAR 1) then update both scalars, then repeat the process 119 120 121#} 122