Changes between Version 7 and Version 8 of PulsedJets


Ignore:
Timestamp:
09/06/13 14:15:36 (11 years ago)
Author:
ehansen
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • PulsedJets

    v7 v8  
    1 Fiducial:
     1== Fiducial Parameters ==
     2[[CollapsibleStart(global.data)]]
     3{{{
     4nDim     = 2                            ! number of dimensions for this problem (1-3)
     5GmX      = 42,420,1                     ! Base grid resolution [x,y,z]
     6MaxLevel = 7                            ! Maximum level for this simulation (0 is fixed grid)
     7
     8GxBounds = 0d0,0d0,0d0,6d0,60d0,0d0     ! Problem boundaries in computational units,format:
     9                                        ! (xlower, ylower, zlower, xupper, yupper, zupper)
     10                                        ! For 2D problems, set zlower and zupper to 0.d0.
     11
     12Gmthbc   = 5,4,1,1,1,1                  ! Sets the physical boundary conditions at the edge of the problem domain
     13                                        ! format:  (x1, y1, z1, x2, y2, z2)
     14                                        ! 1-Extrapolated,
     15                                        ! 2-Periodic,
     16                                        ! 3-ReflectingWall (Field lines do not penetrate)
     17                                        ! 4-Reflect_BParallel (Field lines held normal
     18                                        ! 5-Reflect_Cylindrical (Like reflect wall, but also changes sign of phi
     19                                        !                        components of velocity and magnetic fields)
     20                                        ! [1,1,1,1,1,1]
     21
     22!=============================================================================================
     23! Parameters Related to Time
     24!=============================================================================================
     25start_time        = 0d0                 ! The initial time in computational units [0d0]. Do NOT change this on restarts.
     26final_time        = .2866183232d0                       ! The final time in computational units.
     27final_frame       = 10                  ! The final frame [10]
     28lRestart          = F                   ! restart flag [F]
     29lPostProcess      = F                   ! Post Process (preempts lrestart) [F]
     30restart_frame     = 47                  ! The frame from which restarts should begin.  Ignored unless lRestart or lPostProcess = .true.
     31initial_maxspeed  = 1000.d0             ! A lower bound on the initial maximum wave speed in computational units.
     32cfl_vars          = .8, .4, .1          ! Maximum allowed CFL before restarting step, desired CFL to run at, relaxation parameter.  Smaller numbers give more gradual relaxation [1, .3, .5]
     33
     34
     35!=============================================================================================
     36! Parameters Related to Performance/Algorithms
     37!=============================================================================================
     38lSkipProfile            = T                             ! Flag for skipping profiling [T] - should set this to false for good performance on production runs.
     39lUseOriginalNewSubGrids = F                             ! T for Old Subgrid Algo, F for New Subgrid Algo. [F]
     40MinimumGridPoints       = 2                             ! The minimum number of cells a grid can have along one dimension [2] (ignored unless lUseOriginalNewSubGrids=T).
     41iThreaded               = -1                                    ! -1 is non-threaded, 0 is pseudo-threaded, 1 is threaded [0].  If you are using threaded you must compile with pThreadFlag=2
     42LevelBalance            = 0d0,0d0                       ! Try to balance load on every level? [0d0,0d0]
     43
     44/
     45
     46!============================================================================================
     47!Level-specific variables (All LevelData lists should be initialized from 0 to MaxLevel).
     48!============================================================================================
     49
     50&LevelData
     51qTolerance              = .10,.25,.40,1d30,1d30,1d30,1d30,1d30          ! field variable tolerances on each level
     52DesiredFillRatios       = .70,.70,.70,.70,.70,.70,.70,.70       ! Determines how much of the new subgrids must be filled with refined cells (also ignored unless lUseOriginalNewSubGrids=T)
     53buffers                 = 4, 4, 4, 4, 4, 4, 4, 4                ! ErrFlag buffering
     54/
     55
     56!============================================================================================
     57!Level-specific variables (All LevelData lists should be initialized from 0 to MaxLevel).
     58!============================================================================================
     59
     60&PoissonData
     61solver                  = 1             ! Which solver to use for linear system solve? StructPCG=1, StructGMRes=2 [1]
     62tolerance               = 1e-6          ! Solver convergence tolerance [1e-6]
     63printlevel              = 0             ! Sets print level for hypre solver.  Increasing this will make solvers more verbose. [0]
     64MaxIters                = 1000          ! Maximum number of iterations to attempt convergence.  [1000]
     65hverbosity              = 0             ! Anything > 0 will cause linear systems to be outputted as matrix, rhs, lhs files. [0]
     66mthbc                   = 2 2 2 2 2 2   ! What type of boundary conditions to use for gravitational potential
     67                                        ! format:  (x1, y1, z1, x2, y2, z2)
     68                                        ! 0-Zero Slope,
     69                                        ! 1-Reflecting,
     70                                        ! 2-Periodic,
     71                                        ! 3-Multipole Expansion,
     72                                        ! 4-User Specified (User is reponsible for setting Phi at physical boundaries)
     73}}}
     74[[CollapsibleEnd]]
     75[[CollapsibleStart(physics.data)]]
     76{{{
     77iEOS            = 1            ! Equation of State [0]-ideal gas, 1-MultiSpecies gas, 3-Van der Waals gas, 4-Isothermal
     78gamma           = 1.666667d0   ! Adiabatic index used for ideal gas and Van der Waals gas [1.66666666667]
     79IsoTemp         = 10d0         ! Isothermal temperature for Isothermal EOS
     80CriticalDens    = 30.693677103 ! maximum density in Van der Waals EOS
     81AttractionParam = 1.38d6       ! attraction pram in Van der Waals EOS
     82/
     83&AbundanceData
     84HeFrac = 0.085d0
     85xH     = 0.01d0
     86xHe    = 0d0
     87xHeII  = 0d0
     88ZFrac  = 0.001165d0
     89Zmu    = 16.48d0
     90Xmu    = 0d0
     91xZ     = 4.13127d-4
     92/
     93
     94&PhysicsData
     95
     96!================================================================================
     97! Field based refinement control
     98!================================================================================
     99InterpOpts = 0,0,0,0,0,0,0,0,0,0,0,0,0  ! [0]-constant, 1-minmod, 2-superbee, 3-vanLeer, 4-mc, 5-Parabolic (not conservative), 6-Linear (not limited).  It is recommended to use constant interpolation for any fields in q that correspond to aux fields (iBx, iBy, iBz, etc...)
     100refineVariableFactor = 1d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0 ! weight factors for each field used in triggering refinement.  Combined with qtolerance for that level. [0d0]
     101
     102!================================================================================
     103! Global source term switches
     104!================================================================================
     105lSelfGravity   = .false.    ! Turns on self-gravity if true [.false.]
     106UniformGravity = 0d0        ! Gravitational acceleration in the -y direction [0d0]
     107OmegaRot       = 0d0        ! Shifts to a rotating reference frame (enables coriolis and centripetal force terms)
     108iCylindrical   = 2          ! [0]-No cylindrical geometry, 1-Cylindrical with no angular momentum, 2-Cylindrical with angular momentum
     109iCooling       = 4          ! Cooling module [NoCool=0, AnalyticCool=1, DMCool=2, IICool=3, ZCool=4]
     110SrcPrecision   = 1d-3       ! Relative precision desired for source terms.  [1d-3]
     111lTrackHydrogen = .true.
     112lTrackHelium   = .true.
     113ne_min=1d-3
     114FloorTemp=1d3
     115!================================================================================
     116! Explicit solver switches
     117!================================================================================
     118lResistive     = .false.     ! Turns on resistivity [.false.]
     119ResType        = 2           ! Resistivity type. 1-constant, 2-user defined distribution, 3-Spitzer. [1]
     120resistivity    = 1.2d-1      ! resistivity value [0d0]
     121lConductive    = .false.     ! Turns on mhd heat conduction [.false.]
     122ConType        = 1           ! Conduction type. 1-constant, 2-Spitzer. [1]
     123conductivity   = 8d-2        ! conductivity value [0d0]
     124lThermalFluxLimiter= .false. ! Turns on saturation heat flux. (only works for ConType==2). [.false.]
     125lIsoDiff       = .false.     ! Turns on isotropic heat conduction. Allows isotropicity to be in the range (0,1). [.false.]
     126lViscous       = .false.     ! Turns on viscosity (2D only) [.false.]
     127VisType        = 1           ! viscosity type. 1-classic, 2-Braginskii [1]
     128bviscosity     = 0d0         ! viscosity value [0d0]
     129
     130!================================================================================
     131! Density Protection Options
     132!================================================================================
     133lRestartOnDensityProtections = .false.  ! Do density protections trigger restarts?  [.false.]
     134iDensityProtect  = 2                    ! 0-Set to MinDensity, 1-Set to minimum nearby density, [2]-Average nearby densities
     135iMomentumProtect = 2                    ! 0-Conserve momentum, 1-Set to zero, [2]-Average nearby velocities
     136MinDensity       = 1d-4                ! Minimum computational density before protection is triggered [1d-10]
     137
     138!================================================================================
     139! Pressure Protection Options
     140!================================================================================
     141lRestartOnPressureProtections = .false. ! Do pressure protections trigger restarts? [.false.]
     142iPressureProtect = 2                    ! 0-Set to MinTemp, 1-Set to minimim nearby pressure, [2]-Set to average nearby pressure, 3-Set to minimum nearby temperature, 4-Set to average nearby temperature
     143MinTemp          = 1d-10                ! [1d-10] minimum allowed temperature for the system in Kelvin before protection is triggered
     144
     145!================================================================================
     146! Other Protection Options
     147!================================================================================
     148lRestartOnHypreMaxIters = T             ! Restart steps if hypre does not converge
     149
     150!================================================================================
     151! Description of various scaling parameters
     152! Define one of each of the following: [nScale/rScale], [TempScale/pScale], and set the other to 0d0.
     153!================================================================================
     154nScale          =       1d2,              ! number density scale parameter [particles/cc]
     155rScale          =       0,                ! density scale [g/cc]
     156TempScale       =       1d4,                ! temperature scale parameter [Kelvin]
     157pScale          =       0d0,              ! pressure scale [dynes/cm^2]
     158lScale          =       5e15,       ! length scale parameter [cm] (AU=1.49598e13, pc=3.08568025e18, R_sun=6.955e10
     159
     160!================================================================================
     161! MHD related section
     162!================================================================================
     163lMHD             =  .false.      ! Magnetic Fields present? [.false.]
     164lCheckDivergence =  .false.      ! Turn on divergence checking [.false.]
     165}}}
     166[[CollapsibleEnd]]
     167[[CollpasibleStart(solver.data)]]
     168{{{
     169        iScheme = 1     ! Sweep = [1]
     170        iSolver = 0     ! [0]-Default HLL Solver, 1-iExactRS, 2-HLLC, 3-HLL, 4-HLLC_ISO, 6-HLLD, 8-HLLD_ISO
     171/
     172&SweepData
     173  InterpOrder = 2          ! Interpolation order 1-Piecewise constant, 2-Piecewise linear, [3]-Piecewise Parabolic
     174  lChar_Limiters = .true.  ! Limit using characteristic variables instead of primitive? [.true.]
     175  lCautious=.true.         ! Protect quantities throughout hyperbolic step [.true.]
     176  ViscCD = 0               ! [0] = NO_VISCOSITY, 1 = ROE_VISCOSITY, 2 = H_VISCOSITY
     177  iTracer = 0              ! [0] = NONLAGRANGIAN, 1 = LAGRANGIAN
     178  lApplyDiffusion=.true.   ! Apply additional diffusion at converging flows [.true.]
     179  DIFF_ALPHA = .1          ! Diffusion coefficient [.1]
     180  lApplyLOF = .false.      ! Apply local oscillation filter (reduces carbuncles) [.false.]
     181  LOF_ALPHA = .075         ! Strength of local oscillation filter. [.075]
     182  lCTU=.true.              ! Use corner transport upwind?  [.true.]
     183  lLimiter=.true.          ! Use additional slope limiters [.true.]
     184}}}
     185[[CollapsibleEnd]]
     186[[CollapsibleStart(problem.data)]]
     187{{{
     188namb = 1d2           ! ambient density [cm^-3]
     189tamb = 1d4           ! ambient temperature [K]
     190njet = 5d2           ! jet density [cm^-3]
     191vjet = 200d5          ! initial jet velocity [cm/s]
     192vrange = 50d5         ! (+/-) range of jet velocities [cm/s]
     193Rjet = 1d0           ! jet radius [cu]
     194Rm = 0.6d0          ! radius at which maximum B field is reached
     195Betam = 1d0        ! plasma beta sets the maximum B field
     196Bgeom = 1           ! jet magnetic field geometry (1 = toroidal, 2 = helical)
     197vtype = 1            ! type of veloicty pulse (0 = constant, 1 = sinusoidal, 2 = random)
     198freq = 10d0             ! frequency of velocity pulses in units of 1/final_time
     199                     ! so if sinusoidal, this is how many sine waves the velocity will go through
     200                     ! if random, this is how many times the the velocity will change to a new value
     201theta = 0d0            ! angle of precession of the injection velocity [deg]
     202precfreq = 30d0        ! frequency of the precession in units of 1/final_time
     203seed = 1             ! number used to seed the random number generator (allows runs to be repeatable
     204                     ! and ensures that each processor uses the same random velocity)
     205simtype = 0          ! 0 = single jet, 1 = colliding jets (fields rotating in same direction)
     206                     ! 2 = colliding jets (fields rotating in opposite direction)
     207N = 1               ! number of sub-sampling iterations
     208CoolingRes = 10      ! number of cells per cooling length (using a velocity of 50 km/s)
     209
     210}}}
     211[[CollapsibleEnd]]
     212
    2213Resolution: 42X420 + 7 AMR[[BR]]
    3214