The effects of smoothing length on the onset of wave breaking in smoothed particle hydrodynamics (SPH) simulations of highly directionally spread waves

• Published in: Computational particle mechanics 2022-09, Vol.9 (5), p.1031-1047
• Main Authors: Kanehira, Taiga, McAllister, Mark L., Draycott, Samuel, Nakashima, Takuji, Ingram, David M., van den Bremer, Ton S., Mutsuda, Hidemi
• Format: Article
• Language: English
• Subjects: Breaking waves; Classical and Continuum Physics; Computational Science and Engineering; Engineering; Experiments; Extreme values; Floating bodies; Fluid mechanics; Free surfaces; Mathematical models; Numerical methods; Offshore structures; Simulation; Smooth particle hydrodynamics; Smoothing; Temperature; Theoretical and Applied Mechanics; Wave breaking; Wave crest; Wave power
• ISSN: 2196-4378; 2196-4386
• DOI: 10.1007/s40571-022-00463-z

Ocean wave breaking is a difficult-to-model oceanographic process, which has implications for extreme wave statistics, the dissipation of wave energy, and air–sea interaction. Numerical methods capable of reliably simulating real-world directionally spread breaking waves are useful for investigating the physics of wave breaking and for the design of offshore structures and floating bodies. Smoothed particle hydrodynamics is capable of modelling highly steep and overturning free surfaces, which makes it a promising method for simulating breaking waves. This paper investigates the effect of smoothing length on simulated wave breaking in both following and crossing seas. To do so, we reproduce numerically the experiments of highly directionally spread breaking waves in McAllister et al. (J Fluid Mech 860:767–786, 2019. https://doi.org/10.1017/jfm.2018.886 ) using a range of normalised smoothing lengths: h / d p = 1.4 , 1.7, 2.0, 2.3, with h smoothing length and d p particle spacing. The smallest smoothing length we use appears to adversely affect the fidelity of the simulated surface elevation, so that the tallest wave crest observed in experiments is not fully reproduced (coefficient of determination r 2 ≈ 0.7 ). For smoothing lengths h / d p = 1.7 , 2.0, and 2.3, the experiments are well reproduced ( r 2 ≥ 0.88 ); in these simulations smoothing length predominantly affects the spatial extent and duration of breaking. Qualitative and quantitative comparison of our simulations shows that values of h / d p in the range 1.7 - 2 best reproduce the wave breaking phenomena observed in experiments.