add for fluid and boundary

Author: LasNikas

examples/fluid/hydrostatic_water_column_2d.jl

@@ -73,7 +73,7 @@ semi = Semidiscretization(fluid_system, boundary_system,
 ode = semidiscretize(semi, tspan)
 
 info_callback = InfoCallback(interval=50)
-saving_callback = SolutionSavingCallback(dt=0.02, prefix="", dv=TrixiParticles.acceleration)
+saving_callback = SolutionSavingCallback(dt=0.02, prefix="")
 
 # This is to easily add a new callback with `trixi_include`
 extra_callback = nothing

src/general/custom_quantities.jl

@@ -123,7 +123,3 @@ end
 function avg_density(system, v_ode, u_ode, semi, t)
     return NaN
 end
-
-function acceleration(system, dv_ode, du_ode, v_ode, u_ode, semi, t)
-    return wrap_v(dv_ode, system, semi)
-end

src/schemes/boundary/system.jl

@@ -456,19 +456,20 @@ function system_correction(system::BoundarySPHSystem{<:BoundaryModelDummyParticl
 end
 
 function system_data(system::BoundarySPHSystem, dv_ode, du_ode, v_ode, u_ode, semi)
+    dv = [current_acceleration(system, particle) for particle in eachparticle(system)]
     v = wrap_v(v_ode, system, semi)
     u = wrap_u(u_ode, system, semi)
 
     coordinates = current_coordinates(u, system)
-    velocity = 0 # TODO this is broken
+    velocity = [current_velocity(v, system, particle) for particle in eachparticle(system)]
     density = current_density(v, system)
     pressure = current_pressure(v, system)
 
-    return (; coordinates, velocity, density, pressure)
+    return (; coordinates, velocity, density, pressure, acceleration=dv)
 end
 
 function available_data(::BoundarySPHSystem)
-    return (:coordinates, :velocity, :density, :pressure)
+    return (:coordinates, :velocity, :density, :pressure, :acceleration)
 end
 
 function system_data(system::BoundaryDEMSystem, dv_ode, du_ode, v_ode, u_ode, semi)

src/schemes/fluid/fluid.jl

@@ -177,6 +177,7 @@ end
 function system_data(system::FluidSystem, dv_ode, du_ode, v_ode, u_ode, semi)
     (; mass) = system
 
+    dv = wrap_v(dv_ode, system, semi)
     v = wrap_v(v_ode, system, semi)
     u = wrap_u(u_ode, system, semi)
 
@@ -185,11 +186,11 @@ function system_data(system::FluidSystem, dv_ode, du_ode, v_ode, u_ode, semi)
     density = current_density(v, system)
     pressure = current_pressure(v, system)
 
-    return (; coordinates, velocity, mass, density, pressure)
+    return (; coordinates, velocity, mass, density, pressure, acceleration=dv)
 end
 
 function available_data(::FluidSystem)
-    return (:coordinates, :velocity, :mass, :density, :pressure)
+    return (:coordinates, :velocity, :mass, :density, :pressure, :acceleration)
 end
 
 include("pressure_acceleration.jl")