Changes between Version 11 and Version 12 of PulsedJets


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

—

Legend:

Unmodified
Added
Removed
Modified

  • PulsedJets

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