From e20025c3ec519ac486cc55f19b326232fcc356a9 Mon Sep 17 00:00:00 2001 From: Florian Atteneder <florian.atteneder@uni-jena.de> Date: Fri, 16 Feb 2024 12:07:02 +0100 Subject: [PATCH] wip: stash some tov examples --- examples/grhd_tov/avmda_spherical1d.toml | 4 +- examples/grhd_tov/avmda_spherical1d_k17.toml | 85 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k17_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k18.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k18_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k19.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k19_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k20.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k20_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k21.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k21_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k22.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k22_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/avmda_spherical1d_k23.toml | 84 ++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k23_pt.toml | 86 +++++++++++++++++++ .../grhd_tov/avmda_spherical1d_k24_pt.toml | 86 +++++++++++++++++++ examples/grhd_tov/dg_rescaled.toml | 40 +++++++++ examples/grhd_tov/dg_spherical1d_k.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k17.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k17_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k18.toml | 38 ++++++++ examples/grhd_tov/dg_spherical1d_k18_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k19.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k19_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k20.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k20_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k21.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k21_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k22.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k22_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k23.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k23_pt.toml | 32 +++++++ examples/grhd_tov/dg_spherical1d_k24.toml | 31 +++++++ examples/grhd_tov/dg_spherical1d_k24_pt.toml | 32 +++++++ 34 files changed, 1862 insertions(+), 2 deletions(-) create mode 100644 examples/grhd_tov/avmda_spherical1d_k17.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k17_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k18.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k18_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k19.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k19_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k20.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k20_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k21.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k21_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k22.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k22_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k23.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k23_pt.toml create mode 100644 examples/grhd_tov/avmda_spherical1d_k24_pt.toml create mode 100644 examples/grhd_tov/dg_rescaled.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k17.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k17_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k18.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k18_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k19.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k19_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k20.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k20_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k21.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k21_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k22.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k22_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k23.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k23_pt.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k24.toml create mode 100644 examples/grhd_tov/dg_spherical1d_k24_pt.toml diff --git a/examples/grhd_tov/avmda_spherical1d.toml b/examples/grhd_tov/avmda_spherical1d.toml index d9bfc250..805f9fbe 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 00000000..3ab403f8 --- /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 00000000..8bcc3747 --- /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 00000000..5aa66c8a --- /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 00000000..8e7562b2 --- /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 00000000..42706682 --- /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 00000000..2699bca2 --- /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 00000000..f5a5d217 --- /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 00000000..4b9db624 --- /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 00000000..5496d490 --- /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 00000000..58ae54f0 --- /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 00000000..4f68d6ba --- /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 00000000..33749340 --- /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 00000000..919c14fd --- /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 00000000..8826fa1d --- /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 00000000..8e1ea535 --- /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 00000000..205e41b4 --- /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 00000000..36b6c358 --- /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 00000000..2797e7be --- /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 00000000..36b6c358 --- /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 00000000..ed7b35b5 --- /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 00000000..cbe08744 --- /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 00000000..ac9eaecf --- /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 00000000..9093fca4 --- /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 00000000..0d835a58 --- /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 00000000..d0430cf8 --- /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 00000000..9754dc7d --- /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 00000000..f6090281 --- /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 00000000..f9b5de4e --- /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 00000000..0031befd --- /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 00000000..8858323c --- /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 00000000..dbf16ad1 --- /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 00000000..769eb5fa --- /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 00000000..d2268039 --- /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 -- GitLab