diff --git a/examples/grhd_tov/avmda_spherical1d.toml b/examples/grhd_tov/avmda_spherical1d.toml index d9bfc250d6bee3f09115d76d259c1e8341b74e43..805f9fbe5cb9fae56edd522eca9f5223bd78c2de 100644 --- a/examples/grhd_tov/avmda_spherical1d.toml +++ b/examples/grhd_tov/avmda_spherical1d.toml @@ -10,13 +10,13 @@ id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" atm_factor = 1e-8 formulation = "spherical1d" av_drag = 0.99 -c2p_freeze_atm_reset = false +# c2p_freeze_atm_reset = false # training_wheels = true [Mesh] range = [ 0.0, 20.0 ] n = 5 -k = 24 +k = 16 basis = "lgl" periodic = false diff --git a/examples/grhd_tov/avmda_spherical1d_k17.toml b/examples/grhd_tov/avmda_spherical1d_k17.toml new file mode 100644 index 0000000000000000000000000000000000000000..3ab403f831562e5957a02ec2a22e6921441332b5 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k17.toml @@ -0,0 +1,85 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 17 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k17_pt.toml b/examples/grhd_tov/avmda_spherical1d_k17_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..8bcc37473615f2b03526b5ce0ca228e10cc1e1af --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k17_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 17 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k18.toml b/examples/grhd_tov/avmda_spherical1d_k18.toml new file mode 100644 index 0000000000000000000000000000000000000000..5aa66c8a61e04b3e3b002ac8be5c65c2f485457b --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k18.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 18 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k18_pt.toml b/examples/grhd_tov/avmda_spherical1d_k18_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..8e7562b2804c0083143203bf7fef44722384cc70 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k18_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 18 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k19.toml b/examples/grhd_tov/avmda_spherical1d_k19.toml new file mode 100644 index 0000000000000000000000000000000000000000..42706682f1eeae1cca909c3ab3bfd5d39442b094 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k19.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 19 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k19_pt.toml b/examples/grhd_tov/avmda_spherical1d_k19_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..2699bca2b7b1be317857868e8fb1f444b688b55a --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k19_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 19 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k20.toml b/examples/grhd_tov/avmda_spherical1d_k20.toml new file mode 100644 index 0000000000000000000000000000000000000000..f5a5d217f90cf44ae544ac62a6ad37d189b4d295 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k20.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 20 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k20_pt.toml b/examples/grhd_tov/avmda_spherical1d_k20_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..4b9db6245a7b513d7a6631e4978b67a7f7e14134 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k20_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 20 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k21.toml b/examples/grhd_tov/avmda_spherical1d_k21.toml new file mode 100644 index 0000000000000000000000000000000000000000..5496d4905ca07384ede905c5839709215f03b51c --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k21.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 21 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k21_pt.toml b/examples/grhd_tov/avmda_spherical1d_k21_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..58ae54f0e0d0ade59a2e195087009abef5495434 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k21_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 21 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k22.toml b/examples/grhd_tov/avmda_spherical1d_k22.toml new file mode 100644 index 0000000000000000000000000000000000000000..4f68d6bab342f8ac4839025f5d394c4e699fbc8d --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k22.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 22 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k22_pt.toml b/examples/grhd_tov/avmda_spherical1d_k22_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..337493408a5e2bceeb611014878a0eb36f24a14a --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k22_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 22 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k23.toml b/examples/grhd_tov/avmda_spherical1d_k23.toml new file mode 100644 index 0000000000000000000000000000000000000000..919c14fd50ea95e00ea5c080fdef9dfaab2f9c44 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k23.toml @@ -0,0 +1,84 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 23 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k23_pt.toml b/examples/grhd_tov/avmda_spherical1d_k23_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..8826fa1d07c3fe8fe5dcbf32c17795e93c334684 --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k23_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 23 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/avmda_spherical1d_k24_pt.toml b/examples/grhd_tov/avmda_spherical1d_k24_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..8e1ea535fe41b666b2519d9bef81f0629861e01a --- /dev/null +++ b/examples/grhd_tov/avmda_spherical1d_k24_pt.toml @@ -0,0 +1,86 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +av_drag = 0.99 +perturbation_rho = true +# c2p_freeze_atm_reset = false +# training_wheels = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 24 +basis = "lgl" +periodic = false + +[Output] +# every_iteration = 1 +# # every_dt = 2.5 +# # every_dt = 5.0 +# # aligned_ts = "$(collect(range(0.1,1000.0,step=0.1)))" +# # aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# # aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +# variables = [ "D", "Sr", "Ï„", +# # "p", "Ï", "ϵ", "vr", +# "p", "Ï", "vr", "ÏhW2", +# # "p", "Ï", "vr", "v", +# # "src_D", "src_Sr", "src_Ï„", +# # "rhs_D", "rhs_Sr", "rhs_Ï„", +# # "flr_D", "flr_Sr", "flr_Ï„", +# # "smoothed_mu", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# # "E", #"flr_E" +# # "A", "∂Ar", "B", "∂Br" +# ] +# enable1d = true +# substep_every_dt = 5.0 + +# substep_every_iteration = 1 +# substep_variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "mu", "smoothed_mu", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] + +# every_iteration = 1 +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "mu", "smoothed_mu", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + ] + +[Evolution] +cfl = 0.8 +# tend = 500.0 +tend = 1000.0 +# tend = 12.0 + +[HRSC] +method = "av" +av_method = "mda" +mda_cmax = 100.0 +mda_n_low = 3.0 +mda_n_high = 4.0 +# mda_n_low = 2.0 +# mda_n_high = 3.0 +# av_method = "entropy" +# entropy_ce = 100.0 +# entropy_cmax = 100.0 +# av_smoother_order = 1 + +# [Log] +# progress_stdout = false diff --git a/examples/grhd_tov/dg_rescaled.toml b/examples/grhd_tov/dg_rescaled.toml new file mode 100644 index 0000000000000000000000000000000000000000..205e41b48405fb9fa3485892ed2f5389ca35fd07 --- /dev/null +++ b/examples/grhd_tov/dg_rescaled.toml @@ -0,0 +1,40 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "rescaled_spherical1d" +c2p_freeze_atm_reset = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 16 +basis = "lgl" +periodic = false + +[Output] +every_iteration = 1 +# every_dt = 2.5 +# every_dt = 5.0 +# aligned_ts = "$(collect(range(5.0,5000.0,step=5.0)))" +# aligned_ts = "$(collect(range(0.2,120.0,step=0.2)))" +variables = [ "D", "Sr", "Ï„", + # "p", "Ï", "ϵ", "vr", + "p", "Ï", "vr", + # "p", "Ï", "vr", "v", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "c2p_freeze_atm", + "max_v", + "E", + ] +enable1d = true + +[Evolution] +cfl = 0.5 +tend = 2000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k.toml b/examples/grhd_tov/dg_spherical1d_k.toml new file mode 100644 index 0000000000000000000000000000000000000000..36b6c358e8aa7fd6dd1422ed8db007e6abf38504 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 17 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k17.toml b/examples/grhd_tov/dg_spherical1d_k17.toml new file mode 100644 index 0000000000000000000000000000000000000000..2797e7be6748035c1d28c4f6d8e6b5242093abc5 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k17.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 17 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k17_pt.toml b/examples/grhd_tov/dg_spherical1d_k17_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..36b6c358e8aa7fd6dd1422ed8db007e6abf38504 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k17_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 17 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k18.toml b/examples/grhd_tov/dg_spherical1d_k18.toml new file mode 100644 index 0000000000000000000000000000000000000000..ed7b35b5c36f71ff5992a11bcd79ff3412dda498 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k18.toml @@ -0,0 +1,38 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 18 +basis = "lgl" +periodic = false + +[Output] + +# aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +# variables = [ "D", "Sr", "Ï„", +# "p", "Ï", "vr", "ÏhW2", +# "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", +# "max_v", +# ] +substep_every_iteration = 1 +substep_variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +# tend = 1000.0 +tend = 100.0 diff --git a/examples/grhd_tov/dg_spherical1d_k18_pt.toml b/examples/grhd_tov/dg_spherical1d_k18_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..cbe08744ea6e02f620b53548029cd082310fdde5 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k18_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 18 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k19.toml b/examples/grhd_tov/dg_spherical1d_k19.toml new file mode 100644 index 0000000000000000000000000000000000000000..ac9eaecfd2128d0152ddfd137fc173d437eed898 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k19.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 19 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k19_pt.toml b/examples/grhd_tov/dg_spherical1d_k19_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..9093fca4eba6936d92c11c45b89982e18eedaee3 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k19_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 19 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k20.toml b/examples/grhd_tov/dg_spherical1d_k20.toml new file mode 100644 index 0000000000000000000000000000000000000000..0d835a5832b092d156b5f1846202e65d61baab2d --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k20.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 20 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k20_pt.toml b/examples/grhd_tov/dg_spherical1d_k20_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..d0430cf86aa50a77d8fb5598b8e000cb8b3c1246 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k20_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 20 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k21.toml b/examples/grhd_tov/dg_spherical1d_k21.toml new file mode 100644 index 0000000000000000000000000000000000000000..9754dc7d46ca746af5f13101c050246a2e171f08 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k21.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 21 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k21_pt.toml b/examples/grhd_tov/dg_spherical1d_k21_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..f6090281aec1e98c77a41e50bc00817297bc9451 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k21_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 21 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k22.toml b/examples/grhd_tov/dg_spherical1d_k22.toml new file mode 100644 index 0000000000000000000000000000000000000000..f9b5de4ef30b1536f0b77fa3bc405b46e1723133 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k22.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 22 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k22_pt.toml b/examples/grhd_tov/dg_spherical1d_k22_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..0031befd03f5781b9e9e505a00806b3be28d3028 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k22_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 22 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k23.toml b/examples/grhd_tov/dg_spherical1d_k23.toml new file mode 100644 index 0000000000000000000000000000000000000000..8858323c768fb364957affb03ab1235f203f6e41 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k23.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 23 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k23_pt.toml b/examples/grhd_tov/dg_spherical1d_k23_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..dbf16ad1d1d8bb72506fb3c2765937df9c28c391 --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k23_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 23 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k24.toml b/examples/grhd_tov/dg_spherical1d_k24.toml new file mode 100644 index 0000000000000000000000000000000000000000..769eb5fa5e47159ad0985d7a3998d35bdd1d2efb --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k24.toml @@ -0,0 +1,31 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 24 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0 diff --git a/examples/grhd_tov/dg_spherical1d_k24_pt.toml b/examples/grhd_tov/dg_spherical1d_k24_pt.toml new file mode 100644 index 0000000000000000000000000000000000000000..d2268039841f3fbb54cd66a1d22931906431e7da --- /dev/null +++ b/examples/grhd_tov/dg_spherical1d_k24_pt.toml @@ -0,0 +1,32 @@ +[EquationOfState] +eos = "polytrope" +polytrope_k = 100.0 +polytrope_gamma = 2.0 + +[GRHD] +bc = "from_id" +id = "tov" +id_filename = "$(joinpath(ROOTDIR,\"initialdata\",\"TOV_stable.h5\"))" +atm_factor = 1e-8 +formulation = "spherical1d" +perturbation_rho = true + +[Mesh] +range = [ 0.0, 20.0 ] +n = 5 +k = 24 +basis = "lgl" +periodic = false + +[Output] + +aligned_ts = "$(collect(range(1.0,1000.0,step=1.0)))" +variables = [ "D", "Sr", "Ï„", + "p", "Ï", "vr", "ÏhW2", + "c2p_reset_ϵ", "c2p_reset_atm", "c2p_limit_vr", "c2p_init_admissible", + "max_v", + ] + +[Evolution] +cfl = 0.25 +tend = 1000.0