Jo/support 1d

Project.toml

@@ -12,6 +12,8 @@ Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Profile = "9abbd945-dff8-562f-b5e8-e1ebf5ef1b79"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
+Reduce = "93e0c654-6965-5f22-aba9-9c1ae6b3c259"
+Revise = "295af30f-e4ad-537b-8983-00126c2a3abe"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
 Tullio = "bc48ee85-29a4-5162-ae0b-a64e1601d4bc"

examples/1d/parameters.toml

@@ -0,0 +1,47 @@
+####################################################################################################
+# Simulation
+ndims = 1
+# equations = "CompressibleEulerEquations"
+equations = "LinearScalarAdvectionEquation"
+advectionvelocity = [1.0] # Need only for linarscalaradvection
+
+initial_conditions = "initial_conditions_convergence_test"
+#initial_conditions = "initial_conditions_constant"
+#initial_conditions = "initial_conditions_linear_x_y"
+surface_flux = "flux_lax_friedrichs"
+# source_terms =
+t_start = 0.0
+t_end   = 1.0
+
+# restart = true
+# restart_filename = "out/restart_000100.h5"
+
+
+####################################################################################################
+# Solver
+solver = "dg"
+polydeg = 4
+cfl = 0.2
+n_steps_max = 10000
+analysis_interval = 100
+extra_analysis_quantities = ["entropy", "energy_total"]
+
+
+####################################################################################################
+# Mesh
+n_cells_max = 30000
+coordinates_min = [-1]
+coordinates_max =  [1]
+initial_refinement_level = 4
+# refinement_patches = [
+#   {type = "box", coordinates_min = [-0.5, -0.5], coordinates_max = [0.5, 0.5]},
+#   {type = "box", coordinates_min = [-0.1, -0.1], coordinates_max = [0.1, 0.1]},
+# ]
+
+
+####################################################################################################
+# I/O
+# save_initial_solution = false
+solution_interval = 100
+solution_variables = "primitive"
+restart_interval = 10

examples/1d/parameters_density_wave.toml

@@ -0,0 +1,41 @@
+####################################################################################################
+# Simulation
+ndims = 1
+equations = "CompressibleEulerEquations"
+initial_conditions = "initial_conditions_density_wave"
+surface_flux = "flux_central"
+volume_flux = "flux_central"
+t_start = 0.0
+t_end   = 2.0
+
+# restart = true
+# restart_filename = "out/restart_000100.h5"
+
+
+####################################################################################################
+# Solver
+solver = "dg"
+polydeg = 5
+cfl = 0.8
+n_steps_max = 10000
+analysis_interval = 100
+
+
+####################################################################################################
+# Mesh
+n_cells_max = 30000
+coordinates_min = [-1]
+coordinates_max =  [1]
+initial_refinement_level = 2
+# refinement_patches = [
+#   {type = "box", coordinates_min = [-0.5, -0.5], coordinates_max = [0.5, 0.5]},
+#   {type = "box", coordinates_min = [-0.1, -0.1], coordinates_max = [0.1, 0.1]},
+# ]
+
+
+####################################################################################################
+# I/O
+# save_initial_solution = false
+solution_interval = 100
+solution_variables = "primitive"
+restart_interval = 10

examples/1d/parameters_source_terms.toml

@@ -0,0 +1,42 @@
+####################################################################################################
+# Simulation
+ndims = 1
+equations = "CompressibleEulerEquations"
+
+initial_conditions = "initial_conditions_convergence_test"
+# initial_conditions = "initial_conditions_constant"
+surface_flux = "flux_lax_friedrichs"
+source_terms = "source_terms_convergence_test"
+t_start = 0.0
+t_end   = 2.0
+
+# restart = true
+# restart_filename = "out/restart_000100.h5"
+
+
+####################################################################################################
+# Solver
+solver = "dg"
+polydeg = 3
+cfl = 0.6
+n_steps_max = 10000
+analysis_interval = 100
+
+
+####################################################################################################
+# Mesh
+n_cells_max = 10000
+coordinates_min = [0]
+coordinates_max =  [2]
+initial_refinement_level = 4
+# refinement_patches = [
+#   {type = "box", coordinates_min = 0, coordinates_max = 16},
+# ]
+
+
+####################################################################################################
+# I/O
+# save_initial_solution = false
+solution_interval = 100
+solution_variables = "primitive"
+restart_interval = 100

examples/2d/parameters.toml

@@ -6,8 +6,8 @@ equations = "LinearScalarAdvectionEquation"
 advectionvelocity = [1.0, 1.0] # Need only for linarscalaradvection
 
 initial_conditions = "initial_conditions_convergence_test"
-# initial_conditions = "initial_conditions_constant"
-# initial_conditions = "initial_conditions_linear_x_y"
+#initial_conditions = "initial_conditions_constant"
+#initial_conditions = "initial_conditions_linear_x_y"
 surface_flux = "flux_lax_friedrichs"
 # source_terms =
 t_start = 0.0
@@ -21,7 +21,7 @@ t_end   = 1.0
 # Solver
 solver = "dg"
 polydeg = 3
-cfl = 0.8
+cfl = 0.2
 n_steps_max = 10000
 analysis_interval = 100
 extra_analysis_quantities = ["entropy", "energy_total"]