Explore the vast range of titles available.

Rogue waves (RWs) can form on the ocean surface due to the well-known quasi-four-wave resonant interaction or superposition principle. The first is known as the nonlinear focusing mechanism and leads to an increased probability of RWs when unidirectionality and narrowband energy of the wave field are satisfied. This work delves into the dynamics of extreme wave focusing in crossing seas, revealing a distinct type of nonlinear RWs, characterised by a decisive longevity compared with those generated by the dispersive focusing (superposition) mechanism. In fact, through fully nonlinear hydrodynamic numerical simulations, we show that the interactions between two crossing unidirectional wave beams can trigger fully localised and robust development of RWs. These coherent structures, characterised by a typical spectral broadening then spreading in the form of dual bimodality and recurrent wave group focusing, not only defy the weakening expectation of quasi-four-wave resonant interaction in directionally spreading wave fields, but also differ from classical focusing mechanisms already mentioned. This has been determined following a rigorous lifespan-based statistical analysis of extreme wave events in our fully nonlinear simulations. Utilising the coupled nonlinear Schrödinger framework, we also show that such intrinsic focusing dynamics can be captured by weakly nonlinear wave evolution equations. This opens new research avenues for further explorations of these complex and intriguing wave phenomena in hydrodynamics as well as other nonlinear and dispersive multi-wave systems.

Interaction with an opposing current amplifies wave modulation and accelerates nonlinear wave focusing in regular wavepackets. This results in large-amplitude waves, usually known as rogue waves, even if the wave conditions are less prone to extremes. Laboratory experiments in three independent facilities are presented here to assess the role of opposing currents in changing the statistical properties of unidirectional and directional mechanically generated random wavefields. The results demonstrate in a consistent and robust manner that opposing currents induce a sharp and rapid transition from weakly to strongly non-Gaussian properties. This is associated with a substantial increase in the probability of occurrence of rogue waves for unidirectional and directional sea states, for which the occurrence of extreme and rogue waves is normally the least expected.

We simulate focusing surface gravity wave groups with directional spreading using the modified nonlinear Schrödinger (MNLS) equation and compare the results with a fully-nonlinear potential flow code, OceanWave3D. We alter the direction and characteristic wavenumber of the MNLS carrier wave, to assess the impact on the simulation results. Both a truncated (fifth-order) and exact version of the linear dispersion operator are used for the MNLS equation. The wave groups are based on the theory of quasi-determinism and a narrow-banded Gaussian spectrum. We find that the truncated and exact dispersion operators both perform well if: (1) the direction of the carrier wave aligns with the direction of wave group propagation; (2) the characteristic wavenumber of the carrier wave coincides with the initial spectral peak. However, the MNLS simulations based on the exact linear dispersion operator perform significantly better if the direction of the carrier wave does not align with the wave group direction or if the characteristic wavenumber does not coincide with the initial spectral peak. We also perform finite-depth simulations with the MNLS equation for dimensionless depths ( k p d ) between 1.36 and 5.59, incorporating depth into the boundary conditions as well as the dispersion operator, and compare the results with those of fully-nonlinear potential flow code to assess the finite-depth limitations of the MNLS.

Theoretical researches have established that the energy dynamics of a mature tropical cyclone may be idealized to be very similar to a theoretical Carnot heat engine. Assuming the dissipative heating of the atmospheric boundary layer and the net production of mechanical energy in the cyclone dominate the energy budget of the storm, the potential maximum wind speed of the cyclone can be approximated as a function of the air–sea temperature difference (Ts − T0) and specific enthalpy (k0*−k) difference: |Vmax|2≈CkCDTs−T0T0(k0*−k). Although this theory gives a straighforward estimate of maximum tropical cyclone intensity, studies found that few real storms achieve this theoretical maximum estimated using climatological atmospheric conditions and sea surface temperatures. The discrepancies were attributed to a lack of knowledge of the values of the drag coefficient (CD) and surface exchange coefficient for enthalpy (Ck), and on insufficient upper ocean thermal measurements under hurricanes. Recent observational and numerical studies have unearthed another possible factor for these discrepancies by showing that the energy flux into surface gravity waves under extreme wind conditions can be an order of magnitude greater than formerly believed, and thus may play an important role in the energy budget of tropical cyclones. In this study, numerical experiments are performed to investigate the effect of surface gravity waves under a range of idealized tropical cyclone winds. The wave fields are simulated using the WAVEWATCH III model. Our results demonstrate that by considering the energy flux to surface gravity waves, the potential maximum wind speed can be reduced by up to 12% and this ratio varies with the storm size, intensity, and translation speed.

We report laboratory experiments of long-crested irregular water surface waves propagating over a shoal. For a sufficiently shallow shoal we find that the surface elevation can have a local maximum of skewness and kurtosis above the shallower part of the shoal close to the edge on the incoming side, and a local minimum of skewness over the downward slope on the lee side of the shoal. We find that the horizontal fluid velocity can have a local maximum and minimum of skewness at the same locations as those for the surface elevation. However, the kurtosis of the horizontal fluid velocity can have a local maximum over the downward slope on the lee side of the shoal, different from the location of the maximum of kurtosis of the surface elevation.

We revisit the classical but as yet unresolved problem of predicting the breaking onset of 2D and 3D irrotational gravity water waves. Based on a fully nonlinear 3D boundary element model, our numerical simulations investigate geometric, kinematic and energetic differences between maximally tall non-breaking waves and marginally breaking waves in focusing wave groups. Our study focuses initially on unidirectional domains with flat bottom topography and conditions ranging from deep to intermediate depth (depth to wavelength ratio from 1 to 0.2). Maximally tall non-breaking (maximally recurrent) waves are clearly separated from marginally breaking waves by their normalised energy fluxes localised near the crest tip region. The initial breaking instability occurs within a very compact region centred on the wave crest. On the surface, this reduces to the local ratio of the energy flux velocity (here the fluid velocity) to the crest point velocity for the tallest wave in the evolving group. This provides a robust threshold parameter for breaking onset for 2D wave packets propagating in uniform water depths from deep to intermediate. Further targeted study of representative cases of the most severe laterally focused 3D wave packets in deep and intermediate depth water shows that the threshold remains robust. These numerical findings for 2D and 3D cases are closely supported by our companion observational results. Warning of imminent breaking onset is detectable up to a fifth of a carrier wave period prior to a breaking event.

In previous computational fluid dynamics studies of breaking waves, there has been a marked tendency to severely over-estimate turbulence levels, both pre- and post-breaking. This problem is most likely related to the previously described (though not sufficiently well recognized) conditional instability of widely used turbulence models when used to close Reynolds-averaged Navier–Stokes (RANS) equations in regions of nearly potential flow with finite strain, resulting in exponential growth of the turbulent kinetic energy and eddy viscosity. While this problem has been known for nearly 20 years, a suitable and fundamentally sound solution has yet to be developed. In this work it is demonstrated that virtually all commonly used two-equation turbulence closure models are unconditionally, rather than conditionally, unstable in such regions. A new formulation of the $k$ – $\\unicode[STIX]{x1D714}$ closure is developed which elegantly stabilizes the model in nearly potential flow regions, with modifications remaining passive in sheared flow regions, thus solving this long-standing problem. Computed results involving non-breaking waves demonstrate that the new stabilized closure enables nearly constant form wave propagation over long durations, avoiding the exponential growth of the eddy viscosity and inevitable wave decay exhibited by standard closures. Additional applications on breaking waves demonstrate that the new stabilized model avoids the unphysical generation of pre-breaking turbulence which widely plagues existing closures. The new model is demonstrated to be capable of predicting accurate pre- and post-breaking surface elevations, as well as turbulence and undertow velocity profiles, especially during transition from pre-breaking to the outer surf zone. Results in the inner surf zone are similar to standard closures. Similar methods for formally stabilizing other widely used closure models ( $k$ – $\\unicode[STIX]{x1D714}$ and $k$ – $\\unicode[STIX]{x1D700}$ variants) are likewise developed, and it is recommended that these be utilized in future RANS simulations of surface waves. (In the above $k$ is the turbulent kinetic energy density, $\\unicode[STIX]{x1D714}$ is the specific dissipation rate, and $\\unicode[STIX]{x1D700}$ is the dissipation.)

Abrupt depth transitions (ADTs) have recently been identified as potential causes of ‘rogue’ ocean waves. When stationary and (close-to-) normally distributed waves travel into shallower water over an ADT, distinct spatially localized peaks in the probability of extreme waves occur. These peaks have been predicted numerically, observed experimentally, but not explained theoretically. Providing this theoretical explanation using a leading-order-physics-based statistical model, we show, by comparing to new experiments and numerical simulations, that the peaks arise from the interaction between linear free and second-order bound waves, also present in the absence of the ADT, and new second-order free waves generated due to the ADT.

Freak or rogue waves are so called because of their unexpectedly large size relative to the population of smaller waves in which they occur. The 25.6 m high Draupner wave, observed in a sea state with a significant wave height of 12 m, was one of the first confirmed field measurements of a freak wave. The physical mechanisms that give rise to freak waves such as the Draupner wave are still contentious. Through physical experiments carried out in a circular wave tank, we attempt to recreate the freak wave measured at the Draupner platform and gain an understanding of the directional conditions capable of supporting such a large and steep wave. Herein, we recreate the full scaled crest amplitude and profile of the Draupner wave, including bound set-up. We find that the onset and type of wave breaking play a significant role and differ significantly for crossing and non-crossing waves. Crucially, breaking becomes less crest-amplitude limiting for sufficiently large crossing angles and involves the formation of near-vertical jets. In our experiments, we were only able to reproduce the scaled crest and total wave height of the wave measured at the Draupner platform for conditions where two wave systems cross at a large angle.

Bragg scattering of nonlinear surface waves over a wavy bottom is studied using two-dimensional fully nonlinear numerical wave tanks (NWTs). In particular, we consider cases of high nonlinearity which lead to complex wave generation and transformations, hence possible multiple Bragg resonances. The performance of the NWTs is well verified by benchmarking experiments. Classic Bragg resonances associated with second-order triad interactions among two surface (linear incident and reflected waves) and one bottom wave components (class I), and third-order quartet interactions among three surface (linear incident and reflected waves, and second-order reflected/transmitted waves) and one bottom wave components (class III) are observed. In addition, class I Bragg resonance occurring for the second-order (rather than linear) transmitted waves, and Bragg resonance arising from quintet interactions among three surface and two bottom wave components, are newly captured. The latter is denoted class IV Bragg resonance which magnifies bottom nonlinearity. It is also found that wave reflection and transmission at class III Bragg resonance have a quadratic rather than a linear relation with the bottom slope if the bottom size increases to a certain level. The surface wave and bottom nonlinearities are found to play opposite roles in shifting the Bragg resonance conditions. Finally, the results indicate that Bragg resonances are responsible for the phenomena of beating and parasitic beating, leading to a significantly large local free surface motion in front of the depth transition.