An alternative smooth particle hydrodynamics formulation to simulate chemotaxis in porous media

Towards a High Order Convergent ALE-SPH Scheme with Efficient WENO Spatial Reconstruction

This paper studies the convergence properties of an arbitrary Lagrangian–Eulerian (ALE) Riemann-based SPH algorithm in conjunction with a Weighted Essentially Non-Oscillatory (WENO) high-order spatial reconstruction, in the framework of the DualSPHysics open-source code. A convergence analysis is carried out for Lagrangian and Eulerian simulations and the numerical results demonstrate that, in absence of particle disorder, the overall convergence of the scheme is close to the one guaranteed by the WENO spatial reconstruction. Moreover, an alternative method for the WENO spatial reconstruction is introduced which guarantees a speed-up of 3.5, in comparison with the classical Moving Least-Squares (MLS) approach.

SPH Simulation of Interior and Exterior Flow Field Characteristics of Porous Media

At the present time, one of the most relevant challenges in marine and ocean engineering and practice is the development of a mathematical modeling that can accurately replicate the interaction of water waves with porous coastal structures. Over the last 60 years, multiple techniques and solutions have been identified, from linearized solutions based on wave theories and constant friction coefficients to very sophisticated Eulerian or Lagrangian solvers of the Navier-Stokes (NS) equations. In order to explore the flow field interior and exterior of the porous media under different working conditions, the Smooth Particle Hydrodynamics (SPH) numerical simulation method was used to simulate the flow distribution inside and outside a porous media applied to interact with the wave propagation. The flow behavior is described avoiding Euler’s description of the interface problem between the Euler mesh and the material selected. Considering the velocity boundary conditions and the cyclical circulation boundary conditions at the junction of the porous media and the water flow, the SPH numerical simulation is used to analyze the flow field characteristics, as well as the longitudinal and vertical velocity distribution of the back vortex flow field and the law of eddy current motion. This study provides innovative insights on the mathematical modelling of the interaction between porous structures and flow propagation. Furthermore, there is a good agreement (within 10%) between the numerical results and the experimental ones collected for scenarios with porosity of 0.349 and 0.475, demonstrating that SPH can simulate the flow patterns of the porous media, the flow through the inner and outer areas of the porous media, and the flow field of the back vortex region. Results obtained and the new mathematical approach used can help to effectively simulate with high-precision the changes along the water depth, for a better design of marine and ocean engineering solutions adopted to protect coastal areas.