摘要

The paper presents a numerical wave flume with porous media of variable porosity. A fully-Lagrangian computational method is developed for simulation of fluid flow interactions with porous media based on a projection based (or Incompressible) SPH (Smoothed Particle Hydrodynamics) method. The continuity and Navier-Stokes equations are reformulated with incorporation of volume fraction, based on two-phase mixture theory. The proposed method benefits from a set of enhanced schemes and is therefore referred to as Enhanced Incompressible SPH (ISPH). The effect of porous media is represented by considering linear and nonlinear resistance force terms as well as incorporation of a newly-derived source term of PPE (Poisson Pressure Equation). The results are characterized by smooth and continuous pressure fields, as well as regular particle distributions at the fluid porous media interface in absence of any artificial smoothing schemes. Since the local volume fraction of fluid phase particles is consistently calculated based on the volume fraction of surrounding porous phase particles, the method is capable of accurately reproducing fluid flow interactions with spatially varying porous media. Validations are made through a set of benchmark tests including fluid flow in a U-tube with uniform and multi-layered porous media, rapid seepage flows in rockfill dam, solitary wave attenuation over a porous bed, solitary wave interaction with a submerged porous structure and wave propagation over a submerged triangular permeable bar.