Performance of interFoam on the simulation of progressive waves

The performance of interFoam (a widely-used solver within the popular open source CFD package OpenFOAM) in simulating the propagation of a nonlinear (stream function solution) regular wave is investigated in this work, with the aim of systematically documenting its accuracy. It is demonstrated that over time there is a tendency for surface elevations to increase, wiggles to appear in the free surface, and crest velocities to become (severely) overestimated. It is shown that increasing the temporal and spatial resolution can mitigate these undesirable effects, but that a relatively small Courant number is required. It is further demonstrated that the choice of discretization schemes and solver settings (often treated as a "black box" by users) can have a major impact on the results. This impact is documented, and it is shown that obtaining a "diffusive balance" is crucial to accurately propagate a surface wave over long distances without requiring exceedingly high temporal and spatial resolutions. Finally, the new code isoAdvector is compared to interFoam, which is demonstrated to produce comparably accurate results, while maintaining a sharper surface. It is hoped that the systematic documentation of the performance of the interFoam solver will enable its more accurate and optimal use, as well as increase awareness of potential shortcomings, by CFD researchers interested in the general CFD simulation of free surface waves.Comment: 18 pages and 23 figure

Tsunami Induced Scour Around Monopile Foundations

While the run-up, inundation, and destructive potential of tsunami events has received considerable attention in the literature, the associated interaction with the sea bed i.e. boundary layer dynamics, induced sediment transport, and resultant sea bed morphology, has received relatively little specific attention. The present paper aims to further the understanding of tsunami-induced scour, by numerically investigating tsunami-induced flow and scour processes around a monopile structure, representative of those commonly utilized as offshore wind turbine foundations. The simulations are based on a model solving Reynolds averaged Navier-Stokes (RANS) equations, fully coupled with turbulence closure, bed and suspended load sediment transport descriptions, and a seabed morphological model. The model was developed and utilized in simulating breaker bar morphology by Jacobsen et al. (2014); It has been additionally been used in simulating wave induced scour beneath pipelines by Fuhrman et al. (2014) and Larsen et al. (2016) as well as scour around a monopile by Baykal et al. (2015).Due to the large computational expenses it is presently only feasible to simulate such scour processes around a monopile at model (laboratory) spatial and temporal scales. Therefore, prior to conducting such numerical simulations involving tsunami-induced scour, it is necessary to first establish a methodology for maintaining similarity of model and full field scales. To achieve hydrodynamic similarity we will select the flow parameters such that we maintain similarity in terms of the diameter-based Froude number, as well as the boundary layer thickness-to-diameter ratio δ/D. Equating the Froude number ensures that the adverse pressure gradient induced by the presence of the structure itself will be similar at both model and field scales, i.e. that the ratio of the excess stagnation pressure head in front of the monopile will be maintained. Similarly, by maintaining similarity in δ/D, we ensure that the relative size of the horseshoe vortex, which is expected to largely govern the scouring process, will be similar at both model and full scales. This strategy also yields reasonable similarity in the expected tsunami period-to-scour time scale ratio.As an example, three full tsunami periods have been simulated in succession, taking a full scale period of 13 min. Snapshots of the computed scour hole at selected times when the flow is rightward (left sub-plots) as well as leftward (right sub-plots) are depicted Figure 1. These snapshots illustrate the generally stepwise build-up of scour on the two opposing sides of the monopile during each successive half-cycle of the simulated tsunami. A complementary, simple and practical engineering method for predicting tsunami-induced scour is likewise developed, founded upon existing experimentally-based expressions for use in steady current scour, but invoking the boundary layer thickness and Shields parameter expected from tsunami wave events i.e. it effectively combines both current-like and wave-like properties of tsunamis (see e.g. Williams and Fuhrman, 2016).<br/