cheaper structured mesh RTI with wall BCs

trixi-framework · Sep 6, 2021 · 6921de2 · 6921de2
1 parent ea44d7b
commit 6921de2
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 14 deletions.
diff --git a/examples/structured_2d_dgsem/elixir_euler_rayleigh_taylor_instability.jl b/examples/structured_2d_dgsem/elixir_euler_rayleigh_taylor_instability.jl
@@ -70,23 +70,23 @@ num_elements_per_dimension = 32
 cells_per_dimension = (num_elements_per_dimension, num_elements_per_dimension * 4)
 mesh = StructuredMesh(cells_per_dimension, mapping)
 
-# boundary_conditions = (
-#                        x_neg=boundary_condition_slip_wall,
-#                        x_pos=boundary_condition_slip_wall,
-#                        y_neg=boundary_condition_slip_wall,
-#                        y_pos=boundary_condition_slip_wall,
-#                       )
-
 initial_condition = initial_condition_rayleigh_taylor_instability
 
-# Alternative setup: left/right periodic BCs and Dirichlet BCs on the top/bottom.
 boundary_conditions = (
-                       x_neg=boundary_condition_periodic,
-                       x_pos=boundary_condition_periodic,
-                       y_neg=BoundaryConditionDirichlet(initial_condition),
-                       y_pos=BoundaryConditionDirichlet(initial_condition),
+                       x_neg=boundary_condition_slip_wall,
+                       x_pos=boundary_condition_slip_wall,
+                       y_neg=boundary_condition_slip_wall,
+                       y_pos=boundary_condition_slip_wall,
                       )
 
+# # Alternative setup: left/right periodic BCs and Dirichlet BCs on the top/bottom.
+# boundary_conditions = (
+#                        x_neg=boundary_condition_periodic,
+#                        x_pos=boundary_condition_periodic,
+#                        y_neg=BoundaryConditionDirichlet(initial_condition),
+#                        y_pos=BoundaryConditionDirichlet(initial_condition),
+#                       )
+
 semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver;
                                     source_terms=source_terms_rayleigh_taylor_instability,
                                     boundary_conditions=boundary_conditions)

diff --git a/test/test_structured_2d.jl b/test/test_structured_2d.jl
@@ -182,8 +182,9 @@ isdir(outdir) && rm(outdir, recursive=true)
 
   @trixi_testset "elixir_euler_rayleigh_taylor_instability.jl" begin
     @test_trixi_include(joinpath(EXAMPLES_DIR, "elixir_euler_rayleigh_taylor_instability.jl"),
-      l2   = [0.07623624040606676, 0.007463787908070064, 0.0035904828287556386, 0.010477848481477177],
-      linf = [1.1657733821734124, 0.04325677476747364, 0.04922542769288251, 0.0381568152843168],
+      l2   = [0.06365630381017849, 0.007166887387738937, 0.002878708825497772, 0.010247678114070121],
+      linf = [0.4799214336153155, 0.024595483032220266, 0.02059808120543466, 0.03190756362943725],
+      cells_per_dimension = (8,8),
       tspan = (0.0, 0.3))
   end