In this paper we use the Computational Fluid Dynamics (CFD) toolbox OpenFOAM to perform numerical simulations of multiple floating point absorber wave energy converters (WECs) arranged in a ...geometrical array configuration inside a numerical wave tank (NWT). The two-phase Navier-Stokes fluid solver is coupled with a motion solver to simulate the hydrodynamic flow field around the WECs and the wave-induced rigid body heave motion of each WEC within the array. In this study, the numerical simulations of a single WEC unit are extended to multiple WECs and the complexity of modelling individual floating objects close to each other in an array layout is tackled. The NWT is validated for fluid-structure interaction (FSI) simulations by using experimental measurements for an array of two, five and up to nine heaving WECs subjected to regular waves. The validation is achieved by using mathematical models to include frictional forces observed during the experimental tests. For all the simulations presented, a good agreement is found between the numerical and the experimental results for the WECs’ heave motions, the surge forces on the WECs and the perturbed wave field around the WECs. As a result, our coupled CFD–motion solver proves to be a suitable and accurate toolbox for the study of fluid-structure interaction problems of WEC arrays.
Long-term changes in dune volume at the Belgian coast are analyzed based on measured data by airborne surveys available from 1979. For most of the 65 km long coastal stretch, dune volume increases ...linearly in time at a constant rate. Dune growth varies between 0–12.3 m3/m/year with an average dune growth of 6.2 m3/m/year, featuring large variations in longshore directions. Based on a wind data set from 2000–2017, it is found that potential aeolian sediment transport has its main drift from the west to southwest direction (onshore to oblique onshore). Based on a modified Bagnold model, onshore potential aeolian sediment transport ranges to a maximum of 9 m3/m/year, while longshore potential aeolian sediment transport could reach up to 20 m3/m/year. We found an important correlation between observed and predicted dune development at decadal timescales when zones with dune management activities are excluded. Most of the predicted data are within a factor of two of the measured values. The variability in potential transport is well related to the variability in dune volume changes at the considered spatial–temporal scale, suggesting that natural dune growth is primarily caused by aeolian sediment transport from the beach. It also suggests that annual differences in forcing and transport limiting conditions (wind and moisture) only have a modest effect on the overall variability of dune volume trends.
The objective of the present work is to investigate wave run-up around a monopile subjected to regular waves inside a numerical wave flume using the Computational Fluid Dynamics (CFD) toolbox ...OpenFOAM®. Reynolds-Averaged Navier-Stokes (RANS) turbulence modelling is performed by applying the k-ω SST model. Boundary conditions for wave generation and absorption are adopted from the IHFOAM toolbox. Simulations of propagating water waves show sometimes excessive wave damping (i.e. a significant decrease in wave height over the length of the numerical wave flume) based on RANS turbulence modelling. This anomaly is prevented by implementing a buoyancy term in the turbulent kinetic energy equation. The additional term suppresses the turbulence level at the interface between water and air. The proposed buoyancy-modified k-ω SST turbulence model results in an overall stable wave propagation model without significant wave damping over the length of the flume. Firstly, the necessity of a buoyancy-modified k-ω SST turbulence model is demonstrated for the case of propagating water waves in an empty wave flume. Secondly, numerical results of wave run-up around a monopile under regular waves using the buoyancy-modified k-ω SST turbulence model are validated by using experimental data measured in a wave flume by De Vos et al. (2007). Furthermore, time-dependent high spatial resolutions of the numerically obtained wave run-up around the monopile are presented. These results are in line with the experimental data and available analytical formulations.
•Simulations of wave run-up around a monopile in a wave flume using OpenFOAM®.•RANS turbulence modelling causes wave damping over the length of the flume.•A buoyancy-modified k-ω SST turbulence model is proposed to avoid wave damping.•The buoyancy-modified turbulence model is validated with experimental results.
In this work, the performance of a buoyancy-modified turbulence model is shown for simulating wave breaking in a numerical wave flume. Reynolds-Averaged Navier-Stokes (RANS) modelling is performed by ...applying both a k-ω and a k-ω SST turbulence model using the Computational Fluid Dynamics (CFD) toolbox OpenFOAM. In previous work of the authors (Devolder et al., 2017), the observed significant decrease in wave height over the length of the numerical wave flume based on RANS turbulence modelling for the case of propagating waves has been avoided by developing a buoyancy-modified k-ω SST model in which (i) the density is explicitly included in the turbulence transport equations and (ii) a buoyancy term is added to the turbulent kinetic energy (TKE) equation. In this paper, two buoyancy-modified turbulence models are applied for the case of wave breaking simulations: k-ω and k-ω SST. Numerical results of wave breaking under regular waves are validated with experimental data measured in a wave flume by Ting and Kirby (1994). The numerical results show a good agreement with the experimental measurements for the surface elevations, undertow profiles of the horizontal velocity and turbulent kinetic energy profiles. Moreover, the underlying motivations for the concept of a buoyancy-modified turbulence model are demonstrated by the numerical results and confirmed by the experimental observations. Firstly, the buoyancy term forces the solution of the flow field near the free water surface to a laminar solution in case of wave propagation. Secondly in the surf zone where waves break, the buoyancy term goes to zero and a fully turbulent solution of the flow field is calculated. Finally and most importantly, the buoyancy-modified turbulence models significantly reduce the common overestimation of TKE in the flow field.
•Simulations of wave breaking in a numerical wave flume using OpenFOAM®.•Performance of a buoyancy-modified k-ω and k-ω SST turbulence model are evaluated.•The buoyancy-modified turbulence models are validated with experimental results.•Accurate predictions of the turbulent flow field at the wave breaking zone.
Nature-based solutions to mitigate the impact of future climate change depend on restoring biological diversity and natural processes. Coastal foredunes represent the most important natural flood ...barriers along coastlines worldwide, but their area has been squeezed dramatically because of a continuing urbanization of coastlines, especially in Europe. Dune development is steered by the development of vegetation in interaction with sand fluxes from the beach. Marram grass (
Calamagrostis arenaria
, formerly
Ammophila arenaria
) is the main dune building species along most European coasts, but also in other continents where the species was introduced. Engineering of coastal dunes, for instance by building dunes in front of dikes, needs to be based on a solid understanding of the species’ interactions with the environment. Only quantitative approaches enable the further development of mechanistic models and coastal management strategies that encapsulate these biomorphogenic interactions. We here provide a quantitative review of the main biotic and physical interactions that affect marram grass performance, their interactions with sand fluxes and how they eventually shape dune development. Our review highlights that the species’ spatial organization is central to dune development. We further demonstrate this importance by means of remote sensing and a mechanistic model and provide an outlook for further research on the use of coastal dunes as a nature-based solution for coastal protection.
Abstract Coastal regions globally face escalating challenges from climate change, including rising sea levels and intensified storm events. To address these threats, coastal resilience emerges as a ...critical paradigm advocating the integration of nature‐based solutions with traditional engineering approaches. Coastal dunes, acting as protective barriers, offer a promising avenue. This three‐year study assesses the efficacy of an artificial dune system to address local sand‐related nuisances on the adjacent seawall featuring planted marram grass in Oosteroever, Belgium. The focus is on understanding sediment accumulation, dune morphology and vegetation development. The results demonstrate a significant increase in dune height, reaching up to 2 m in the zones planted with marram grass, surpassing the height of the adjacent seawall. Comprehensive profiles and drone surveys revealed consistent dune growth of 27 m 3 /m, which contrasted with the substantial erosion in the adjacent unvegetated beach areas of up to 30 m 3 /m. One storm event caused dune toe erosion of 1.5 m 3 /m, but the dune demonstrated rapid recovery through natural aeolian processes. Marram grass development was not impacted by the initial planting configuration and density and was more pronounced at the perimeter edges of the dune. This study highlights the success of the ‘dune‐in‐front‐of‐a‐dike’ approach, offering insights for sustainable coastal resilience strategies.
This paper presents a study on the coupling between a fluid solver and a motion solver to perform fluid–structure interaction (FSI) simulations of floating bodies such as point absorber wave energy ...converters heaving under wave loading. The two-phase fluid solver with dynamic mesh handling, interDyMFoam, is a part of the Computational Fluid Dynamics (CFD) toolbox OpenFOAM. The incompressible Navier–Stokes (NS) equations are solved together with a conservation equation for the Volume of Fluid (VoF). The motion solver is computing the kinematic body motion induced by the fluid flow. A coupling algorithm is needed between the fluid solver and the motion solver to obtain a converged solution between the hydrodynamic flow field around and the kinematic motion of the body during each time step in the transient simulation. For body geometries with a significant added mass effect, simple coupling algorithms show slow convergence or even instabilities. In this paper, we identify the mechanism for the numerical instability and we derive an accelerated coupling algorithm (based on a Jacobian) to enhance the convergence speed between the fluid and motion solver. Secondly, we illustrate the coupling algorithm by presenting a free decay test of a heaving wave energy converter. Thirdly and most challenging, a water impact test of a free falling wedge with a significant added mass effect is successfully simulated. For both test cases, the numerical results obtained by using the accelerated coupling algorithm are in a very good agreement with the experimental measurements.
In preparation for the use of computational fluid dynamics (CFD) simulation results as ‘numerical experiments’ in fire research, the agreement with experimental data for two different small-scale ...set-ups is discussed. The first configuration concerns the position of smoke-free height in case of fire with vertical ventilation in an atrium. The second set-up deals with the critical velocity for smoke backlayering in case of fire in a horizontally ventilated tunnel. An
N-percent rule is introduced for the determination of the presence of smoke in the simulation results, based on the local temperature rise. The CFD package FDS is used for the numerical simulations. The paper does not scrutinize the detailed accuracy of the results, as this is hardly possible with any state-of-the-art experimental data at hand. Rather, the global accuracy is discussed with current numerical implementation and models in FDS, considering continuous evolution over different version releases with time. The agreement between the experiments and numerical simulations is very promising. Even when quantitative agreement with experimental data is not perfect, the trends are very well reproduced in the simulations. While much additional work is required, both in CFD as in ‘real’ experiments, the results are encouraging for the potential of state-of-the-art CFD to be used as numerical experiments.
This study examines the effect of fetch distance on aeolian sediment transport and the corresponding growth rate of aeolian transport rate on a sandy beach under varying wind speeds. Eight ...experiments were conducted on both dry and intertidal beach areas, with wind speeds ranging from 9 to 17 m/s. Vertically stacked sand traps were used to measure transport rates over distances up to 100 m, with experiments lasting 20–40 minutes. Results indicate an initial increase in transport rate until the critical fetch distance, followed by a subsequent decrease in sediment transport with further fetch distance. Critical fetch distances, corresponding to maximum transport rates, ranged from 23.8 to 103.5 m. A strong linear correlation was observed between wind speed and critical fetch distance, suggesting that higher wind speeds necessitate longer fetch distances for maximum transport. Furthermore, established formulations by van Rijn and Strypsteen are tested to estimate maximum sediment transport and critical fetch distance based on grain‐related shear velocity. Combining all experiments showed that the growth rate in transport up to the point of critical fetch can be best described by a trigonometric function. These findings contribute to a better understanding of the fetch effect on aeolian sediment dynamics in coastal sandy environments.
This study investigates the fetch effect on aeolian sediment transport in coastal sandy environments. Results reveal an initial increase in transport rate up to a critical fetch distance, followed by a decline. A new linear correlation with wind speed suggests higher speeds necessitate longer fetch for maximum transport. Grain‐related shear velocity can also be used to predict the critical fetch distance. The findings enhance our understanding of aeolian dynamics on sandy beaches.
Abstract
The dynamics of the dune toe along the urbanised 65 km macro‐tidal coast of Belgium has been examined based on (bi‐)annual cross‐shore profiles derived from airborne LiDAR surveys conducted ...between 2000 and 2019. Results indicate that the average dune toe level is located at +5.9 m Tweede Algemene Waterpassing (TAW, Belgian Ordnance Datum), which is 1 m lower than the conventional dune toe level. However, this level is not static but rather increasing over time with an average rate of 2.3 cm/year, making it comparable with coastal areas along western Europe. Both landward retreat and seaward progradation of the dune toes are observed, with rates up to 2 m/year. The analysis revealed that dune growth is primarily facilitated by the development of incipient or embryonic dunes, new foredunes and modifications to the stoss slope of the original dune profile. Half of the sea‐fronting dunes were observed to be developed by brushwood fences, which had a positive effect on dune growth. At some locations, incipient dune development and shoreline progradation were observed seaward of old foredunes. Other dune regions were characterised by natural dune blowouts and management activities for recreational purposes. Dune toe reinforcements were implemented in the west to stabilise the dune toe and prevent erosion by waves. However, it was noted that if the adjacent beach accreted because of natural growth or periodic nourishment, incipient dunes could form in front of the reinforcement.