blob: d2caef31356184932dc647d8ad5a9ab64d4a6a85 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
##
ActiveThorns = "time MoL cartgrid3d carpetioascii ioutil boundary hydrobase grhydro grhydro_initdata nanchecker eos_omni CoordBase LoopControl CarpetLib Carpet CarpetReduce spacemask ADMCoupling ADMMacros ADMBASE STATICCONFORMAL TMUNUBASE SymBase localreduce CarpetIOBasic initbase "
ActiveThorns = "CarpetIOHDF5"
driver::ghost_size=2 # for TVD
grhydro::grhydro_stencil=2
Time::dtfac = 0.25 # 0.5 see p 409 (0.25 from David Neilsen, Eric WHirschmann, and R Steven Millward arXiv:gr-qc/0512147v1)
# rk2 dies quite quickly
methodoflines::ODE_Method = "rk3"
methodoflines::MoL_Intermediate_Steps=3
#methodoflines::MoL_Num_Scratch_Levels=1
ADMBase::initial_shift = "zero"
HydroBase::initial_hydro="rotor"
HydroBase::timelevels=2
HydroBase::evolution_method="grhydro"
HydroBase::Bvec_evolution_method="grhydro"
# rotor problem from DelZanna
# needs to be before EOS_Omni
grhydro_initdata::rotor_use_smoothing = no
SpaceMask::use_mask="yes"
EOS_Omni::gl_gamma=5./3. # relativistic rotor problem of DelZanna & Buccianti & Londrille AaA_400_397_2003
# try to mimic the core fluid
EOS_Omni::poly_k = "grhydro_initdata::rotor_pressin/grhydro_initdata::rotor_rhoin**EOS_Omni::gl_gamma"
EOS_Omni::poly_gamma="EOS_Omni::gl_gamma"
grhydro::recon_method="tvd"
grhydro::riemann_solver="HLLE"
grhydro::tvd_limiter="vanleerMC"
grhydro::bound="flat"
grhydro::grhydro_eos_type = "General"
grhydro::grhydro_eos_table="Ideal_Fluid"
grhydro::Grhydro_MaxNumConstrainedVars=13
grhydro::GRHydro_MaxNumEvolvedVars=10
grhydro::GRHydro_MaxNumSandRVars=0 # ADMBase is not evolved
grhydro::rho_abs_min = 1e-14 # need some reasonable value here, even though no density ever gets this low. This is about the same as GRHydro's default value from rho_central and rho_rel_min
grhydro::clean_divergence = "yes"
grhydro::transport_constraints = "no"
grhydro::decouple_normal_Bfield = "yes"
grhydro::kap_dc = 5.0
grhydro::track_divB = "yes"
grid::type = "coordbase"
grid::domain = "full"
CoordBase::domainsize = "minmax"
CoordBase::xmin = 0.
CoordBase::ymin = 0.
CoordBase::zmin = -0.0125
CoordBase::xmax = 1.0
CoordBase::ymax = 1.0
CoordBase::zmax = 0.0125
CoordBase::boundary_size_x_lower = carpet::ghost_size
CoordBase::boundary_size_y_lower = carpet::ghost_size
CoordBase::boundary_size_z_lower = carpet::ghost_size
CoordBase::boundary_size_x_upper = carpet::ghost_size
CoordBase::boundary_size_y_upper = carpet::ghost_size
CoordBase::boundary_size_z_upper = carpet::ghost_size
CoordBase::spacing = "numcells"
CoordBase::ncells_x = 400
CoordBase::ncells_y = 400
CoordBase::ncells_z = 10
carpet::domain_from_coordbase = "yes"
Cactus::terminate = "time"
Cactus::cctk_final_time = 0.4 # see DelZanna p 410
IO::out_dir = $parfile
CarpetIOBasic::outInfo_vars="
hydrobase::rho
hydrobase::press
hydrobase::eps
hydrobase::vel[0] hydrobase::vel[1]
hydrobase::w_lorentz
GRHydro::bcons[0] GRHydro::bcons[1]
GRHydro::divB GRHydro::psidc
"
CarpetIOBasic::outInfo_every=1
CarpetIOASCII::out1D_every = 10
CarpetIOASCII::out1D_z = "no"
CarpetIOASCII::out1D_d = "no"
CarpetIOASCII::out1D_xline_y = 0.5
CarpetIOASCII::out1D_yline_x = 0.5
CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec grhydro::Bcons grhydro::divB grhydro::psidc"
CarpetIOHDF5::one_file_per_group = "yes"
CarpetIOHDF5::open_one_input_file_at_a_time = "yes"
CarpetIOHDF5::out2D_every = 1
CarpetIOHDF5::out2D_xy = "yes"
CarpetIOHDF5::out2D_xz = "no"
CarpetIOHDF5::out2D_yz = "no"
CarpetIOHDF5::out2D_xyplane_z = 0.0
CarpetIOHDF5::out2D_vars = "
GRHydro::Bcons
GrHydro::dens
GrHydro::divB
GRHydro::scon
GRHydro::tau
HydroBase::Bvec
HydroBase::eps
HydroBase::press
HydroBase::rho
HydroBase::vel
HydroBase::w_lorentz
GRHydro::psidc
"
|
