Search

Nearshore rocky reefs with scales of order 10–100 m are common along the world's coastline and often shape wave-driven hydrodynamics and shoreline morphology in their lee. The interaction of waves with these reefs generally results in either two or four-cell mean circulation systems (2CC and 4CC, respectively), with diverging flows behind the reefs and at the shoreline in the 2CC case and flows that diverge in the lee and converge at the shoreline in the 4CC case. By applying a phase-resolving wave-flow model to conduct a detailed analysis of mean momentum balances for waves interacting with nearshore reefs, we develop an understanding of the drivers of 2CC and 4CC flow dynamics and how they vary for different reef geometries and wave and water level conditions. The 2CC or 4CC patterns were primarily driven by alongshore pressure gradients toward the exposed (nonreef fronted) or reef-fronted beach. These alongshore pressure gradients were dependent on the cross-shore setup dynamics governed by the balance between pressure (i.e., related to the setup) and radiation stress gradients, and mean bottom stresses exerted on the water column. If shoreline wave setup in the lee of the reef was less than the exposed beach, a 4CC pattern developed with convergent flow at the shoreline in the lee of the reef; otherwise, a 2CC emerged with divergent flow at the shoreline. Across the parameter space investigated, reef roughness, distance to the shoreline, and beach slope were the three parameters most likely to change the flow patterns between 2CC and 4CC., Environmental Fluid Mechanics

Wave breaking on the steep fore-reef slopes of shallow fringing reefs can be effective at dissipating incident sea-swell waves prior to reaching reef shorelines. However, wave setup and free infragravity waves generated during the sea-swell breaking process are often the largest contributors to wave-driven water levels (wave runup) at the shoreline. Laboratory flume experiments and a two-dimensional vertical phase-resolving nonhydrostatic wave-flow model, which includes a canopy model to predict drag forces generated by roughness elements, were used to investigate wave-driven water levels for along-shore uniform fringing reefs. In contrast to many previous studies, both the laboratory experiment and the numerical model account for the effects of large bottom roughness. The numerical model reproduced the observations of the wave transformation and runup over both smooth and rough reef profiles. The numerical model was then extended to quantify the influence of reef geometry and compared to simulations of plane beaches lacking a reef. For a fixed offshore forcing condition, the fore-reef slope controlled wave runup on reef-fronted beaches, whereas the beach slope controlled wave runup on plane beaches. As a result, the coastal protection utility of reefs is dependent on these slopes. For our examples, with a fore-reef slope of 1/5 and a 500 m prototype reef flat length, a beach slope of ∼1/30 marked the transition between the reef providing runup reduction for steeper beach slopes and enhancing wave runup for milder slopes. Roughness coverage, spacing, dimensions, and drag coefficient were investigated, with results indicating the greatest runup reductions were due to tall roughness elements on the reef flat., Green Open Access added to TU Delft Institutional Repository 'You share, we take care!' - Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

Engineered and natural submerged coastal structures (e.g., submerged breakwaters and reefs) modify incident wave fields and thus can alter hydrodynamic processes adjacent to coastlines. Although submerged structures are generally assumed to promote beach protection by dissipating waves offshore and creating sheltered conditions in their lee, their interaction with waves can result in mean wave-driven circulation patterns that may either promote shoreline accretion or erosion. Here, we analyse the mean flow patterns and shoreline water levels (wave runup) in the lee of idealised impermeable submerged structures with a phase-resolved nonhydrostatic numerical model. Waves propagating over submerged structures can drive either a 2-cell mean (wave-averaged) circulation, which is characterised by diverging flows behind the structure and at the shoreline, or 4-cell circulation, with converging flows at the shoreline and diverging flows in the immediate lee of the structure. The numerical results show that the mode of circulation can be predicted with a set of relationships depending on the incoming wave heights, the structure crest level, and distance to the shoreline (or structure depth). Qualitative agreement between the mean flow and proxies for the sediment transport using an energetics approach suggest that the mean flow can be a robust proxy for inferring sediment transport patterns. For the cases considered, the submerged structures had a minimal influence on shoreline wave setup and wave runup despite the wave energy dissipation by the structures due to alongshore wave energy fluxes in the lee. Consequently, these results suggest that the coastal protection provided by the range of impermeable submerged structures we modelled is primarily due to their capacity to promote beach accretion., Green Open Access added to TU Delft Institutional Repository ‘You share, we take care!’ – Taverne project https://www.openaccess.nl/en/you-share-we-take-care Otherwise as indicated in the copyright section: the publisher is the copyright holder of this work and the author uses the Dutch legislation to make this work public., Environmental Fluid Mechanics

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-ree, Peer Reviewed, Postprint (author's final draft)

Two-dimensional mean wave-driven flow and setup dynamics were investigated at a reef-lagoon system at Ningaloo Reef, Western Australia, using the numerical wave-flow model, SWASH. Phase-resolved numerical simulations of the wave and flow fields, validated with highly detailed field observations (including >10 sensors through the energetic surf zone), were used to quantify the main mechanisms that govern the mean momentum balances and resulting mean current and setup patterns, with particular attention to the role of nonlinear wave shapes. Momentum balances from the phase-resolved model indicated that onshore flows near the reef crest were primarily driven by the wave force (dominated by radiation stress gradients) due to intense breaking, whereas the flow over the reef flat and inside the lagoon and channels was primarily driven by a pressure gradient. Wave setup inside the lagoon was primarily controlled by the wave force and bottom stress. The bottom stress reduced the setup on the reef flat and inside the lagoon. Excluding the bottom stress contribution in the setup balance resulted in an over prediction of the wave-setup inside the lagoon by up to 200–370%. The bottom stress was found to be caused by the combined presence of onshore directed wave-driven currents and (nonlinear) waves. Exclusion of the bottom stress contribution from nonlinear wave shapes led to an over prediction of the setup inside the lagoon by approximately 20–40%. The inclusion of the nonlinear wave shape contribution to the bottom stress term was found to be particularly relevant in reef regions that experience a net onshore mass flux over the reef crest., Environmental Fluid Mechanics

In this study we assess the capabilities of the mesh-free, Lagrangian particle method (Smooth Particle Hydrodynamics, SPH) method to simulate the detailed hydrodynamic processes generated by both spilling and plunging breaking waves within the surf zone, and present how the approach can be used to predict a broad spectrum of hydrodynamic processes relevant to coastal applications where wave breaking is important. The weakly-compressible SPH code DualSPHysics was applied to simulate wave breaking over two bathymetric profiles (a plane beach and fringing reef) and compared to experimental flume measurements of waves, currents, and mean water levels. We demonstrate how the model can accurately reproduce a broad range of relevant hydrodynamic processes, ranging from the nonlinear evolution of wave shapes across the surfzone, wave setup distributions, mean current profiles and wave runup. We compare the surfzone predictions with results from other classes of wave models, and illustrate some of the advantages of the SPH approach (particularly in resolving the hydrodynamics above the wave trough). Overall, the results reveal how the mesh-free SPH approach can accurately reproduce the detailed wave breaking processes with comparable skill to state-of-the-art mesh-based Computational Fluid Dynamic models, and how it can be used as a valuable tool to develop new physical insight into surf zone processes., Accepted Author Manuscript, Hydraulic Structures and Flood Risk

In this study we investigated the capabilities of the mesh-free, Lagrangian particle method (Smoothed Particle Hydrodynamics, SPH) to simulate the detailed hydrodynamic processes generated by both spilling and plunging breaking waves within the surf zone. The weakly-compressible SPH code DualSPHysics was applied to simulate wave breaking over two distinct bathymetric profiles (a plane beach and fringing reef) and compared to experimental flume measurements of waves, flows, and mean water levels. Despite the simulations spanning very different wave breaking conditions (including an extreme case with violently plunging waves on an effectively dry reef slope), the model was able to reproduce a wide range of relevant surf zone hydrodynamic processes using a fixed set of numerical parameters. This included accurate predictions of the nonlinear evolution of wave shapes (e.g., asymmetry and skewness properties), rates of wave dissipation within the surf zone, and wave setup distributions. By using this mesh-free approach, the model was able to resolve the critical crest region within the breaking waves, which provided robust predictions of the wave-induced mass fluxes within the surf zone responsible for the undertow. Within this breaking crest region, the model results capture how the potential energy of the organized wave motion is initially converted to kinetic energy and then dissipated, which reproduces the distribution of wave forces responsible for wave setup generation across the surf zone. Overall, the results reveal how the mesh-free SPH approach can accurately reproduce the detailed wave breaking processes with comparable skill to state-of-the-art mesh-based Computational Fluid Dynamics (CFD) models, and thus can be applied to provide valuable new physical insight into surf zone dynamics., Peer Reviewed, Postprint (author's final draft)

A combination of numerical hydrodynamic models, a large-clast inverse sediment-transport model, and extensive field measurements were used to discriminate between a tsunami and a storm striking Anegada, BVI a few centuries ago. In total 161 cobbles and boulders were measured, ranging from 1.5 to 830 kg at distances of up to 1 km from the shoreline and 2 km from the crest of a fringing coral reef. Transported limestone clasts were derived from outcrops in the low lying interior of Anegada. Estimates of the near-bed flow velocities required to transport the observed boulders were calculated using a simple sediment-transport model, which accounts for fluid drag, inertia, buoyancy and lift forces on boulders and includes both sliding and overturning transport mechanisms. Estimated near-bed flow velocities are converted to depth-averaged velocities using a linear eddy viscosity model and compared with flow depth, acceleration, and depth-averaged velocity time series from high resolution coastal inundation models. Such models simulate overwash by the storm surge and waves of a category 5 hurricane and tsunami inundation from a Lisbon earthquake of M 9.0 and two hypothetical earthquakes along the North America Caribbean Plate boundary. A modeled category 5 hurricane and three simulated tsunamis were all capable of inundating the boulder fields and transporting a portion of the observed clasts, but only an earthquake of M 8.0 on a normal fault of the outer-rise along the Puerto Rico Trench was found to be capable of transporting the largest clasts at their current locations. Model results show that, while both storm waves and tsunamis are capable of generating velocities and accelerations necessary to transport large boulders near the reef crest, attenuation of wave energy due to wave breaking and bottom friction limits the capacity of storm waves to transport large clasts at great inland distances. Through sensitivity analysis, it is shown that even when using coefficients