supersede trixi-framework#852

Author: LasNikas

src/schemes/boundary/open_boundary/boundary_zones.jl

@@ -196,6 +196,19 @@ function BoundaryZone(; plane, plane_normal, density, particle_spacing,
     reference_values = (reference_velocity=velocity_ref, reference_pressure=pressure_ref,
                         reference_density=density_ref)
 
+    coordinates_svector = reinterpret(reshape, SVector{NDIMS, ELTYPE}, ic.coordinates)
+
+    if prescribed_pressure
+        ic.pressure .= pressure_ref.(coordinates_svector, 0)
+    end
+    if prescribed_density
+        ic.density .= density_ref.(coordinates_svector, 0)
+        ic.mass .= ic.density * ic.particle_spacing^NDIMS
+    end
+    if prescribed_velocity
+        ic.velocity .= stack(velocity_ref.(coordinates_svector, 0))
+    end
+
     return BoundaryZone(ic, spanning_set_, zone_origin, zone_width,
                         flow_direction, plane_normal_, reference_values,
                         average_inflow_velocity, prescribed_density, prescribed_pressure,

test/schemes/boundary/open_boundary/characteristic_variables.jl

@@ -15,7 +15,8 @@
     # Prescribed quantities
     reference_velocity = (pos, t) -> SVector(t, 0.0)
     reference_pressure = (pos, t) -> 50_000.0 * t
-    reference_density = (pos, t) -> 1000.0 * t
+    # Add small offset to avoid "ArgumentError: density must be positive and larger than `eps()`"
+    reference_density = (pos, t) -> 1000.0 * (t + sqrt(eps()))
 
     # Plane points of open boundary
     plane_points_1 = [[0.0, 0.0], [0.5, -0.5], [1.0, 0.5]]