Comments

trixi-framework · Mar 28, 2024 · ec0f3b0 · ec0f3b0
1 parent aa2d8e3
commit ec0f3b0
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl b/examples/tree_1d_dgsem/elixir_linearizedeuler_gauss_wall.jl
@@ -18,6 +18,9 @@ mesh = TreeMesh(coordinates_min, coordinates_max,
                 n_cells_max = 100_000,
                 periodicity = false)
 
+# Initialize density and pressure perturbation with a Gaussian bump
+# that is advected to left with v - c and to the right with v + c.
+# Correspondigly, the bump splits in half.
 function initial_condition_gauss_wall(x, t, equations::LinearizedEulerEquations1D)
     v1_prime = 0.0
     rho_prime = p_prime = 2 * exp(-(x[1] - 45)^2 / 25)