parameter_example.yml 8.48 KB
Newer Older
1
# Define the system of units to use internally.
2
InternalUnitSystem:
3
4
5
6
7
  UnitMass_in_cgs:     1   # Grams
  UnitLength_in_cgs:   1   # Centimeters
  UnitVelocity_in_cgs: 1   # Centimeters per second
  UnitCurrent_in_cgs:  1   # Amperes
  UnitTemp_in_cgs:     1   # Kelvin
Matthieu Schaller's avatar
Matthieu Schaller committed
8

9
# Parameters for the task scheduling
Matthieu Schaller's avatar
Matthieu Schaller committed
10
Scheduler:
11
12
13
14
15
16
  nr_queues:              0         # (Optional) The number of task queues to use. Use 0  to let the system decide.
  cell_max_size:          8000000   # (Optional) Maximal number of interactions per task if we force the split (this is the default value).
  cell_sub_size_pair:     256000000 # (Optional) Maximal number of interactions per sub-pair task  (this is the default value).
  cell_sub_size_self:     32000     # (Optional) Maximal number of interactions per sub-self task  (this is the default value).
  cell_split_size:        400       # (Optional) Maximal number of particles per cell (this is the default value).
  max_top_level_cells:    12        # (Optional) Maximal number of top-level cells in any dimension. The number of top-level cells will be the cube of this (this is the default value).
17
  tasks_per_cell:         0         # (Optional) The average number of tasks per cell. If not large enough the simulation will fail (means guess...).
Matthieu Schaller's avatar
Matthieu Schaller committed
18

19
# Parameters governing the time integration (Set dt_min and dt_max to the same value for a fixed time-step run.)
Matthieu Schaller's avatar
Matthieu Schaller committed
20
TimeIntegration:
21
22
23
24
25
  time_begin: 0.    # The starting time of the simulation (in internal units).
  time_end:   1.    # The end time of the simulation (in internal units).
  dt_min:     1e-6  # The minimal time-step size of the simulation (in internal units).
  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).

26
27
# Parameters governing the snapshots
Snapshots:
28
29
30
  basename:   output      # Common part of the name of output files
  time_first: 0.          # Time of the first output (in internal units)
  delta_time: 0.01        # Time difference between consecutive outputs (in internal units)
31
  compression: 0          # (Optional) Set the level of compression of the HDF5 datasets [0-9]. 0 does no compression.
32
33
34
35
36
  UnitMass_in_cgs:     1  # (Optional) Unit system for the outputs (Grams)
  UnitLength_in_cgs:   1  # (Optional) Unit system for the outputs (Centimeters)
  UnitVelocity_in_cgs: 1  # (Optional) Unit system for the outputs (Centimeters per second)
  UnitCurrent_in_cgs:  1  # (Optional) Unit system for the outputs (Amperes)
  UnitTemp_in_cgs:     1  # (Optional) Unit system for the outputs (Kelvin)
37

38
39
# Parameters governing the conserved quantities statistics
Statistics:
Matthieu Schaller's avatar
Matthieu Schaller committed
40
41
  delta_time:          1e-2      # Time between statistics output
  energy_file_name:    energy    # (Optional) File name for energy output
42
  timestep_file_name:  timesteps # (Optional) File name for timing information output. Note: No underscores "_" allowed in file name
43

44
# Parameters for the hydrodynamics scheme
Matthieu Schaller's avatar
Matthieu Schaller committed
45
SPH:
46
47
  resolution_eta:        1.2348   # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
48
  h_tolerance:           1e-4     # (Optional) Relative accuracy of the Netwon-Raphson scheme for the smoothing lengths.
49
  h_max:                 10.      # (Optional) Maximal allowed smoothing length in internal units. Defaults to FLT_MAX if unspecified.
50
51
  max_volume_change:     1.4      # (Optional) Maximal allowed change of kernel volume over one time-step.
  max_ghost_iterations:  30       # (Optional) Maximal number of iterations allowed to converge towards the smoothing length.
Matthieu Schaller's avatar
Matthieu Schaller committed
52

53
54
# Parameters for the self-gravity scheme
Gravity:
55
56
57
58
59
60
  eta:          0.025    # Constant dimensionless multiplier for time integration.
  theta:        0.7      # Opening angle (Multipole acceptance criterion)
  epsilon:      0.1      # Softening length (in internal units).
  a_smooth:     1.25     # (Optional) Smoothing scale in top-level cell sizes to smooth the long-range forces over (this is the default value).
  r_cut_max:    4.5      # (Optional) Cut-off in number of top-level cells beyond which no FMM forces are computed (this is the default value).
  r_cut_min:    0.1      # (Optional) Cut-off in number of top-level cells below which no truncation of FMM forces are performed (this is the default value).
61

62
# Parameters related to the initial conditions
Matthieu Schaller's avatar
Matthieu Schaller committed
63
InitialConditions:
64
  file_name:  SedovBlast/sedov.hdf5 # The file to read
65
  cleanup_h:   0                    # (Optional) Clean the values of h that are read in. Set to 1 to activate.
66
67
  h_scaling:  1.                    # (Optional) A scaling factor to apply to all smoothing lengths in the ICs.
  shift_x:    0.                    # (Optional) A shift to apply to all particles read from the ICs (in internal units).
Matthieu Schaller's avatar
Matthieu Schaller committed
68
69
  shift_y:    0.
  shift_z:    0.
70
  replicate:  2                     # (Optional) Replicate all particles along each axis a given number of times. Default 1.
Matthieu Schaller's avatar
Matthieu Schaller committed
71

72
# Parameters governing domain decomposition
Matthieu Schaller's avatar
Matthieu Schaller committed
73
DomainDecomposition:
74
75
76
77
78
79
80
81
82
83
84
  initial_type:     simple_metis # (Optional) The initial decomposition strategy: "grid",
                                 #            "simple_metis", "weighted_metis", or "vectorized".
  initial_grid_x:   10      # (Optional) Grid size if the "grid" strategy is chosen.
  initial_grid_y:   10      # ""
  initial_grid_z:   10      # ""
  repartition_type: task_weights # (Optional) The re-decomposition strategy: "none", "task_weights", "particle_weights",
                                 #            "edge_task_weights" or "hybrid_weights".
  trigger:          0.05    # (Optional) Fractional (<1) CPU time difference between MPI ranks required to trigger a
                            # new decomposition, or number of steps (>1) between decompositions
  minfrac:          0.9     # (Optional) Fractional of all particles that should be updated in previous step when
                            # using CPU time trigger
85

86
# Parameters related to external potentials --------------------------------------------
87

88
# Point mass external potentials
89
PointMassPotential:
90
  position_x:      50.      # location of external point mass (internal units)
91
  position_y:      50.
92
  position_z:      50.
93
94
  mass:            1e10     # mass of external point mass (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
95

96
# Isothermal potential parameters
97
IsothermalPotential:
98
  position_x:      100.     # Location of centre of isothermal potential with respect to centre of the box (internal units)
99
  position_y:      100.
100
  position_z:      100.
101
102
  vrot:            200.     # Rotation speed of isothermal potential (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
103
  epsilon:         0.1      # Softening size (internal units)
104

105
# Disk-patch potential parameters
106
107
108
DiscPatchPotential:
  surface_density: 10.      # Surface density of the disc (internal units)
  scale_height:    100.     # Scale height of the disc (internal units)
109
110
111
  z_disc:          400.     # Position of the disc along the z-axis (internal units)
  z_trunc:         300.     # (Optional) Distance from the disc along z-axis above which the potential gets truncated.
  z_max:           380.     # (Optional) Distance from the disc along z-axis above which the potential is set to 0.
112
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
113
  growth_time:     5.       # (Optional) Time for the disc to grow to its final size (multiple of the dynamical time)
114

115
116
117
118
119
120
# Sine Wave potential
SineWavePotential:
  amplitude:        10.     # Amplitude of the sine wave (internal units)
  timestep_limit:   1.      # Time-step dimensionless pre-factor.
  growth_time:      0.      # (Optional) Time for the potential to grow to its final size.
 
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
# Parameters related to cooling function  ----------------------------------------------

# Constant du/dt cooling function
ConstCooling:
  cooling_rate: 1.          # Cooling rate (du/dt) (internal units)
  min_energy:   1.          # Minimal internal energy per unit mass (internal units)
  cooling_tstep_mult: 1.    # Dimensionless pre-factor for the time-step condition

# Constant lambda cooling function
LambdaCooling:
  lambda:                      2.0   # Cooling rate (in cgs units)
  minimum_temperature:         1.0e4 # Minimal temperature (Kelvin)
  mean_molecular_weight:       0.59  # Mean molecular weight
  hydrogen_mass_abundance:     0.75  # Hydrogen mass abundance (dimensionless)
  cooling_tstep_mult:          1.0   # Dimensionless pre-factor for the time-step condition