Grain-scale investigation of swash zone sediment transport on a gravel beach using DEM-MPS coupled scheme

Sediment transport in the swash zone directly affects beach changes such as shoreline recession; thus, detailed understandings of sediment transport mechanisms are necessary to accurately estimate the short-time scales sediment transport rate. However, these detailed mechanisms under runup waves have not been elucidated because of the complex solid-gas-liquid multiphase turbulence flow. In this study, we attempt to numerically investigate the sediment grain-scale mechanism to overcome the shortcomings of experimental measurements and the free surface treatment in many numerical simulations. The gravel transport process on a sloped beach under regular waves was simulated using a 2D coupled model of the discrete element method (DEM) and a modified moving particle semi-implicit (MPS) method; a sub-model was built into the DEM-MPS model to improve fluid volume conservation. After validating the simulated performance by comparing it to a previous experiment, the gravel motions were investigated for turbulence and inner beach structure. The Shields number, estimated using the drag force distribution, revealed that significant turbulence contributed to onshore gravel transport near the rundown limit. The inter-gravel contact structure inside the beach explained the decrease in offshore sediment transport during backwash as increased resistance to gravel motions resulting from beach compaction.