Implement Tensile Instability Control (TIC)

src/schemes/fluid/pressure_acceleration.jl

@@ -58,6 +58,14 @@ end
     return -volume_term * pressure_tilde * W_a
 end
 
+# Same as `pressure_acceleration_continuity_density`, but using the minus formulation
+# when pressures are negative to avoid tensile instability.
+# See Lyu et al. (2021). https://doi.org/10.1016/j.apor.2021.102938
+@inline function tensile_instability_correction(m_a, m_b, rho_a, rho_b, p_a, p_b, W_a)
+    p = p_b + abs(p_a)
+    return -m_b * p / (rho_a * rho_b) * W_a
+end
+
 function choose_pressure_acceleration_formulation(pressure_acceleration,
                                                   density_calculator, NDIMS, ELTYPE,
                                                   correction)