add TLSPH support (#954)

Co-authored-by: LasNikas <[email protected]>
Co-authored-by: Erik Faulhaber <[email protected]>

Author: 3 people

src/general/semidiscretization.jl

@@ -825,6 +825,26 @@ function update_nhs!(neighborhood_search,
             semi, points_moving=(true, true), eachindex_y=each_active_particle(neighbor))
 end
 
+function update_nhs!(neighborhood_search,
+                     system::OpenBoundarySystem{<:BoundaryModelDynamicalPressureZhang},
+                     neighbor::TotalLagrangianSPHSystem,
+                     u_system, u_neighbor, semi)
+    update!(neighborhood_search,
+            current_coordinates(u_system, system),
+            current_coordinates(u_neighbor, neighbor),
+            semi, points_moving=(true, true), eachindex_y=each_active_particle(neighbor))
+end
+
+function update_nhs!(neighborhood_search,
+                     system::TotalLagrangianSPHSystem,
+                     neighbor::OpenBoundarySystem{<:BoundaryModelDynamicalPressureZhang},
+                     u_system, u_neighbor, semi)
+    update!(neighborhood_search,
+            current_coordinates(u_system, system),
+            current_coordinates(u_neighbor, neighbor),
+            semi, points_moving=(true, true), eachindex_y=each_active_particle(neighbor))
+end
+
 function update_nhs!(neighborhood_search,
                      system::OpenBoundarySystem, neighbor::TotalLagrangianSPHSystem,
                      u_system, u_neighbor, semi)

src/schemes/structure/total_lagrangian_sph/system.jl

@@ -496,11 +496,12 @@ end
 
 # To account for boundary effects in the viscosity term of the RHS, use the viscosity model
 # of the neighboring particle systems.
-@inline function viscosity_model(system::TotalLagrangianSPHSystem, neighbor_system)
+@inline function viscosity_model(system::TotalLagrangianSPHSystem,
+                                 neighbor_system::AbstractFluidSystem)
     return neighbor_system.viscosity
 end
 
-@inline function viscosity_model(system::AbstractFluidSystem,
+@inline function viscosity_model(system::Union{AbstractFluidSystem, OpenBoundarySystem},
                                  neighbor_system::TotalLagrangianSPHSystem)
     return neighbor_system.boundary_model.viscosity
 end