From 24977275e691c0793484d70948e7f44dbacfa403 Mon Sep 17 00:00:00 2001 From: Daniel_Doehring Date: Sun, 7 Jan 2024 11:17:57 +0100 Subject: [PATCH] two boundaries --- .../elixir_euler_free_stream_boundaries.jl | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/examples/p4est_3d_dgsem/elixir_euler_free_stream_boundaries.jl b/examples/p4est_3d_dgsem/elixir_euler_free_stream_boundaries.jl index 20d374da0f3..bdc4da26c1f 100644 --- a/examples/p4est_3d_dgsem/elixir_euler_free_stream_boundaries.jl +++ b/examples/p4est_3d_dgsem/elixir_euler_free_stream_boundaries.jl @@ -16,17 +16,18 @@ solver = DGSEM(polydeg = polydeg, surface_flux = flux_lax_friedrichs) ############################################################################### # Get the uncurved mesh from a file (downloads the file if not available locally) -default_mesh_file = joinpath(@__DIR__, "cube_boundaries.inp") +default_mesh_file = joinpath(@__DIR__, "mesh_cube_with_boundaries.inp") isfile(default_mesh_file) || - download("https://gist.github.com/DanielDoehring/eefe73ae5d113bc3944a518b6e88e663/raw/359a58a808790f3c3efc050273270eb1cc8ee353/cube_boundaries.inp", + download("https://gist.githubusercontent.com/DanielDoehring/710eab379fe3042dc08af6f2d1076e49/raw/38e9803bc0dab9b32a61d9542feac5343c3e6f4b/mesh_cube_with_boundaries.inp", default_mesh_file) mesh_file = default_mesh_file -boundary_symbols = [:PhysicalSurface1] +boundary_symbols = [:PhysicalSurface1, :PhysicalSurface2] mesh = P4estMesh{3}(mesh_file, polydeg = polydeg, boundary_symbols = boundary_symbols) -boundary_conditions = Dict(:PhysicalSurface1 => BoundaryConditionDirichlet(initial_condition)) +boundary_conditions = Dict(:PhysicalSurface1 => BoundaryConditionDirichlet(initial_condition), + :PhysicalSurface2 => BoundaryConditionDirichlet(initial_condition)) semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver, boundary_conditions = boundary_conditions)