Vegetation meadows in coastal waters are a key constituent of a future green defense package due to the ecosystem services they provide and the potential to attenuate wave energy. To numerically describe the vegetation dynamics under wave action, this paper presents a novel application of a numerical coupling for solving fluid–elastic structure interactions (FSI) problems involving ultra-thin elements in a 3-D environment. The extended two-way coupling employed in this work combines the mesh-free Smoothed Particle Hydrodynamics (SPH) method in the DualSPHysics code to solve the fluid flow, and the Finite Element Analysis (FEA) structural solver in Project Chrono to solve the structural dynamics. To represent the vegetation, a flexible structure based on the Euler–Bernoulli beam model is used. The beam element is embedded into the SPH domain using an envelope subdomain that is discretized using dummy boundary particles. As such, this dummy envelope serves as a decoupling interface for the geometrical properties of the structure, allowing for ultra-thin structures smaller than the initial inter-particle distance (dp). The numerical approach is validated against an experimental setup including a flexible blade swaying under the action of an oscillatory flow. The results demonstrate that the numerical model is able to resolve the wave–vegetation interaction problem. Furthermore, additional insights into the blade dynamics reveal that the swaying velocity increases linearly along the length, with the upper part swaying at a speed comparable to the fluid velocity while the stem remains relatively stationary. Additionally, the findings indicate that rigid vegetation experiences higher forces per unit length, and in systems with substantial swaying motion, energy dissipation predominantly occurs around the lower base of the vegetation.

As waves interact with the slopes of coral reefs and other steep bathymetry profiles, plunging breaking usually occurs where the free surface overturns and violent water motion is triggered. Resolving these surf zone processes pose significant challenges for conventional mesh-based hydrodynamic models, due to the rapidly-deforming nature of the free surface and associated flows. Yet the accurate prediction of these surf zone hydrodynamics is critical for predicting a wide range of nearshore processes driven by wave breaking (e.g., wave dissipation and energy transfers; mean water levels and currents; and wave runup). In this study we assess the ability of the mesh-free, Lagrangian particle-based numerical modelling approach Smoothed Particle Hydrodynamics (SPH) based on DualSPHysics, to simulate the fine-scale hydrodynamic processes driven by irregular wave transformation over a fringing reef profile, by comparing results against detailed experimental observations from a physical modelling study. To greatly improve the computational efficiency, the SPH model was coupled to the mesh-based multi-layer nonhydrostatic wave-flow model SWASH. With this coupled approach, SWASH was used to efficiently simulate the evolution of non-breaking waves from the wavemaker up to the fore reef slope, with the SPH model then used to simulate the detailed hydrodynamic processes over the reef from just offshore of the breakpoint to the shoreline. The SPH model was able to accurately reproduce the complex free surface deformations during plunging breaking, the spectral evolution of waves across the reef flat (including nonlinear wave shape), the mean water levels and currents, and wave runup at the shoreline. Using the long duration simulations (>400 wave periods), the model was able to reproduce the full range of wave motions over the reef (from sea-swell to infragravity frequencies), including the increasing dominance of low frequency waves towards the shoreline and the large cross-reef standing wave motions excited by the reef geometry.