Shear dilatation/contraction of granular materials has long been recognized as an important process in granular flows but a comprehensive theoretical description of this process for a wide range of shear rates is not yet available. In this paper, a theoretical formulation of dilatation/contraction is proposed for continuum modelling of granular flows, in which the dilatation/contraction effects consist of a frictional component, which results from the rearrangement of enduring-contact force chains among particles, and a collisional component, which arises from inter-grain collisions. In this formulation, a frictional solid pressure, which considers the rearrangement of contact force chains under shear deformation, is proposed for the frictional dilatation/contraction, while well-established rheological laws are adopted for the collisional inter-grain pressure to account for the collisional dilatancy effect. The proposed formulation is first verified analytically by describing the shear-weakening behaviour of granular samples in a torsional shear rheometer and by capturing the incipient failure of both dry and immersed granular slopes. The proposed dilatation/contraction formulation is then further validated numerically by integrating it into a two-fluid continuum model and applying the model to study the collapse of submerged granular columns, in which the dilatation/contraction plays a critical role.

When dealing with the complex deformation of free surface such as wave breaking, traditional mesh-based Computational Fluid Dynamics (CFD) methods often face problems arising alongside grid distortion and re-meshing. Therefore, the meshless method became robust for treating large displaced free surface and other boundaries caused by moving structures. The particle method, which is an important branch of meshless method, is mainly divided into the Smoothed Particle Hydrodynamics (SPH) and Moving Particle Semi-implicit (MPS) methods. Different from the SPH method, which involves continuity and treat density as a variable when building kernel functions, the kernel function in the MPS method is a weight function which treats density as a constant, and the spatial derivatives are discretized by establishing the gradient operator and Laplace operator separately. In other words, the first- or second-order continuity of the kernel functions in the MPS method is not a necessity as in SPH, though it might be desirable. At present, the MPS method has been successfully applied to various violent-free surface flow problems in ocean engineering and diverse applications have been comprehensively demonstrated in a number of review papers. This work will focus on algorithm developments of the MPS method and to provide all perspectives in terms of numerical algorithms along with their pros and cons.

Meshless methods have become an essential numerical tool for simulating a wide range of flow-structure interaction problems. However, the way by which the impermeable solid boundary condition is implemented can significantly affect the accuracy of the results and computational cost. This paper develops an improved boundary scheme through a weak formulation for the boundary particles based on Pressure Poisson’s Equation (PPE). In this scheme, the wall boundary particles simultaneously satisfy the PPE in the local integration domain by adopting the Meshless Local Petrov–Galerkin method with the Rankine source solution (MLPG_R) integration scheme (Ma, 2005b) and the Neumann boundary condition, i.e., normal pressure gradient condition, on the wall boundary which truncates the local integration domain. The new weak formulation vanishes the derivatives of the unknown pressure at wall particles and is discretized in the truncated support domain without extra artificial treatment. This improved boundary scheme is validated by analytical solutions, numerical benchmarks, and experimental data in the cases of patch tests, lid-driven cavity, flow over a cylinder and monochromic wave generation. Second-order convergent rate is achieved even for disordered particle distributions. The results show higher accuracy in pressure and velocity, especially near the boundary, compared to the existing boundary treatment methods that directly discretize the pressure Neumann boundary condition.

Piles are common support elements for marine and coastal structures. The scour around pile foundations caused by currents is a major threat to the stability and safety of these structures. The empirical equations commonly used for estimating the equilibrium scour depth around pile groups are limited in their predicative capability, especially when the current approaches the pile group at an angle. This study applies a Multi-Layer Perceptron Backpropagation (MLP/BP) neural network to develop a general model for predicting the local maximum equilibrium scour depth around pile groups in steady currents. The input parameters for the model include all relevant non-dimensional hydrodynamic and structural variables taking full account of the effects of the pile group arrangement and its orientation relative to the approaching current. The model’s performance was evaluated by comparing its predictions against those generated by multiple other machine learning methods, as well as against results from widely used empirical formulas. A comprehensive sensitivity analysis is carried out to determine the importance ranking of the input parameters on model accuracy.