Source author record

Francesco Fedele

Francesco Fedele appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2026arXiv

Fold catastrophe in breaking waves

We present a dynamical-systems perspective on wave breaking for ideal incompressible free-surface flows. By tracking the most energetic hotspot on the wave surface, we find that near breaking the surface slope m evolves on a fast timescale governed by the small parameter epsilon = (partial_z u)^(-1), the inverse vertical velocity gradient at the hotspot, while the focusing parameter A = (U - Ce)/(U - Creq) varies slowly and adiabatically. Here U is the horizontal fluid velocity at the energetic point, Ce its propagation speed, and Creq the equivalent crest speed. This slow-fast structure reveals a fold catastrophe in the (m, A) space whose boundary forms the geometric skeleton organizing the dynamics near breaking. Finite-time blowup occurs when the trajectory crosses this boundary, marking the onset of breaking. The inception of breaking is further characterized by crossing the slope threshold theta* = arctan(sqrt(2) - 1) = 22.5 degrees. This critical angle marks the maximum anisotropy that can be sustained between the Hessians of the velocity and pressure fields, reflecting an imbalance between kinetic and potential energy fluxes. The anisotropy of the velocity Hessian also gives rise to the classical 30-degree slope observed at the inflection point of steep waves near breaking inception. The crest height is limited by the maximum excess of kinetic over potential energy that the flow can sustain, beyond which breaking becomes inevitable. Wave breaking can also be interpreted as a gravity analogue of a collapsing black hole, with apparent and event horizons representing the onset and inception of breaking.

preprint2022arXiv

On the stability of the laminar boundary layer beneath a Stokes wave

The linear stability of the laminar boundary layer flow of a Stokes wave in deep waters is investigated by means of a 'momentary' criterion of instability for unsteady flows (Blondeaux and Seminara, 1979). In the parameter range investigated, it is found that the flow is stable to 2-D perturbations. The least stable eigenmode of the resulting Orr-Sommerfield spectrum attains its maximum beneath the boundary layer of the Stokes wave. Moreover, an analysis of the associated pseudospectrum indicates that the laminar flow when modified by imperfections is unstable due to the non-normality of the Orr-Sommerfeld operator, and the unstable pseudo-eigenmodes tend to peak within the boundary layer. The laminar flow of the Stokes wave is also stable to 3-D streamwise-independent perturbations. Instability is observed for the laminar flow with imperfections. The associated unstable pseudo-eigenmodes are streamwise vortical rolls similar to Langmuir cells. The laminar boundary layer flow of a Stokes wave appears to be stable to infinitesimal perturbations, but it may likely be unstable to finite perturbations, as in Poiseuille pipe flows. The present results are supportive of the recent experimental evidence of spontaneous occurrence of turbulence beneath unforced non-breaking surface waves.

preprint2020arXiv

Crest speeds of unsteady surface water waves

Intuitively, crest speeds of water waves are assumed to match their phase velocities. However, this is generally not the case for natural waves within unsteady wave groups. This motivates our study, which presents new insights into the generic behavior of crest speeds of linear to highly nonlinear unsteady waves. While our major focus is on gravity waves where a generic crest slowdown occurs cyclically, results for capillary-dominated waves are also discussed, for which crests cyclically speed up. This curious phenomenon arises when the theoretical constraint of steadiness is relaxed, allowing waves to change their form, or shape. In particular, a kinematic analysis of both simulated and observed open ocean gravity waves reveals a forward-to-backward leaning cycle for each individual crest within a wave group. This is clearly manifest during the focusing of dominant wave groups essentially due to the dispersive nature of waves. It occurs routinely for focusing linear (vanishingly small steepness) wave groups, and it is enhanced as the wave spectrum broadens. It is found to be relatively insensitive to the degree of phase coherence and focusing of wave groups. The nonlinear nature of waves limits the crest slowdown. This reduces when gravity waves become less dispersive, either as they steepen or as they propagate over finite water depths. This is demonstrated by numerical simulations of the unsteady evolution of 2D and 3D dispersive gravity wave packets in both deep and intermediate water depths, and by open ocean space-time measurements.

preprint2017arXiv

Some special solutions to the Hyperbolic NLS equation

The Hyperbolic Nonlinear Schrodinger equation (HypNLS) arises as a model for the dynamics of three-dimensional narrowband deep water gravity waves. In this study, the Petviashvili method is exploited to numerically compute bi-periodic time-harmonic solutions of the HypNLS equation. In physical space they represent non-localized standing waves. Non-trivial spatial patterns are revealed and an attempt is made to describe them using symbolic dynamics and the language of substitutions. Finally, the dynamics of a slightly perturbed standing wave is numerically investigated by means a highly acccurate Fourier solver.

preprint2016arXiv

Are rogue waves really unexpected?

An unexpected wave is defined by Gemmrich & Garrett (2008) as a wave that is much taller than a set of neighboring waves. Their definition of "unexpected" refers to a wave that is not anticipated by a casual observer. Clearly, unexpected waves defined in this way are predictable in a statistical sense. They can occur relatively often with a small or moderate crest height, but large unexpected waves that are rogue are rare. Here, this concept is elaborated and statistically described based on a third-order nonlinear model. In particular, the conditional return period of an unexpected wave whose crest exceeds a given threshold is developed. This definition leads to greater return periods or on average less frequent occurrences of unexpected waves than those implied by the conventional return periods not conditioned on a reference threshold. Ultimately, it appears that a rogue wave that is also unexpected would have a lower occurrence frequency than that of a usual rogue wave. As specific applications, the Andrea and WACSIS rogue wave events are examined in detail. Both waves appeared without warning and their crests were nearly $2$-times larger than the surrounding $O(10)$ wave crests, and thus unexpected. The two crest heights are nearly the same as the threshold~$h_{0.3\cdot10^{6}}\sim1.6H_{s}$ exceeded on average once every~$0.3\cdot 10^{6}$ waves, where $H_s$ is the significant wave height. In contrast, the Andrea and WACSIS events, as both rogue and unexpected, would occur slightly less often and on average once every~$3\cdot10^{6}$ and~$0.6\cdot10^6$ waves respectively.

preprint2016arXiv

Kinematics of fluid particles on the sea surface. Hamiltonian theory

We derive the John-Sclavounos equations describing the motion of a fluid particle on the sea surface from first principles using Lagrangian and Hamiltonian formalisms applied to the motion of a frictionless particle constrained on an unsteady surface. The main result is that vorticity generated on a stress-free surface vanishes at a wave crest when the horizontal particle velocity equals the crest propagation speed, which is the kinematic criterion for wave breaking. If this holds for the largest crest, then the symplectic two-form associated with the Hamiltonian dynamics reduces instantaneously to that associated with the motion of a particle in free flight, as if the surface did not exist. Further, exploiting the conservation of the Hamiltonian function for steady surfaces and traveling waves we show that particle velocities remain bounded at all times, ruling out the possibility of the finite-time blowup of solutions.

preprint2015arXiv

On Oceanic Rogue Waves

We propose a new conceptual framework for the prediction of rogue waves and third-order space-time extremes of wind seas that relies on the Tayfun (1980) and Janssen (2009) models coupled with Adler-Taylor (2009) theory on the Euler characteristics of random fields. Extreme statistics of the Andrea rogue wave event are examined capitalizing on European Reanalysis (ERA)-interim data. A refinement of Janssen's (2003) theory suggests that in realistic oceanic seas characterized by short-crested multidirectional waves, homogeneous and Gaussian initial conditions become irrelevant as the wave field adjusts to a non-Gaussian state dominated by bound nonlinearities over time scales $t\gg t_{c}\approx0.13T_{0}/νσ_θ$, where $T_{0}$, $ν$ and $σ_θ$ denote mean wave period, spectral bandwidth and angular spreading of dominant waves. For the Andrea storm, ERA-interim predictions yield $t_{c}/T_{0}\sim O(1)$ indicating that quasi-resonant interactions are negligible. Further, the mean maximum sea surface height expected over the Ekofisk platform's area is higher than that expected at a fixed point. However, both of these statistics underestimate the actual crest height $h_{obs}\sim1.63H_s$ observed at a point near the Ekofisk site, where $H_s$ is the significant wave height. To explain the nature of such extreme, we account for both skewness and kurtosis effects and consider the threshold $h_{q}$ exceeded with probability $q$ by the maximum surface height of a sea state over an area in time. We find that $h_{obs}$ nearly coincides with the threshold $h_{1/1000}\sim1.62H_s$ estimated at a point for a typical $3$-hour sea state, suggesting that the Andrea rogue wave is likely to be a rare occurrence in quasi-Gaussian seas.

preprint2015arXiv

On the kurtosis of deep-water gravity waves

In this paper, we revisit Janssen's (2003) formulation for the dynamic excess kurtosis of weakly nonlinear gravity waves at deep water. For narrowband directional spectra, the formulation is given by a sixfold integral that depends upon the Benjamin-Feir index and the parameter $R=σ_θ^{2}/2ν^{2}$, a measure of short-crestedness for the dominant waves with $ν$ and $σ_θ$} denoting spectral bandwidth and angular spreading. Our refinement leads to a new analytical solution for the dynamic kurtosis of narrowband directional waves described with a Gaussian type spectrum. For multidirectional or short-crested seas initially homogenous and Gaussian, in a focusing (defocusing) regime dynamic kurtosis grows initially, attaining a positive maximum (negative minimum) at the intrinsic time scale \[ τ_{c}=ν^{2}ω_{0}t_{c}=1/\sqrt{3R},\qquad\mathrm{or}\qquad t_{c}/T_{0}\approx0.13/νσ_θ, \] where $ω_{0}=2π/T_{0}$ denotes the dominant angular frequency. Eventually the dynamic excess kurtosis tends monotonically to zero as the wave field reaches a quasi-equilibrium state characterized with nonlinearities mainly due to bound harmonics. Quasi-resonant interactions are dominant only in unidirectional or long-crested seas where the longer-time dynamic kurtosis can be larger than that induced by bound harmonics, especially as the Benjamin-Feir index increases. Finally, we discuss the implication of these results on the prediction of rogue waves.

preprint2014arXiv

Geometric phases of water waves

Recently, Banner et al. (2014) highlighted a new fundamental property of open ocean wave groups, the so-called crest slowdown. For linear narrowband waves, this is related to the geometric and dynamical phase velocities $U_d$ and $U_g$ associated with the parallel transport through the principal fiber bundle of the wave motion with $\mathit{U}(1)$ symmetry. The theoretical predictions are shown to be in fair agreement with ocean field observations, from which the average crest speed $c=U_d+U_g$ with $c/U_d\approx0.8$ and $U_{g}/U_d\approx-0.2$.

preprint2014arXiv

Interval-based parameter identification for structural static problems

We present an interval-based approach for parameter identification in structural static inverse problems. The proposed inverse formulation exploits the Interval Finite Element Method (IFEM) combined with adjoint-based optimization. The inversion consists of a two-step algorithm: first, an estimate of the parameters is obtained by means of a deterministic iterative solver. Then, the algorithm switches to the interval extension of the previous solver, using the deterministic estimate of the parameters as an initial guess. The iterations are terminated based on a new containment-stopping criterion, which is intrinsic to intervals. Various numerical examples show that the proposed method provides guaranteed interval enclosures of the parameters.

preprint2013arXiv

Linking reduced breaking crest speeds to unsteady nonlinear water wave group behavior

Observations show that maximally-steep breaking water wave crest speeds are much slower than expected. We report a wave-crest slowdown mechanism generic to unsteady propagating deep water wave groups. Our fully nonlinear computations show that just prior to reaching its maximum height, each wave crest slows down significantly and either breaks at this reduced speed, or accelerates forward unbroken. This finding is validated in our extensive laboratory and field observations. This behavior appears to be generic to unsteady dispersive wave groups in other natural systems.

preprint2013arXiv

On certain properties of the compact Zakharov equation

Long-time evolution of a weakly perturbed wavetrain near the modulational instability threshold is investigated within the framework of the compact Zakharov equation for unidirectional deep-water waves, recently derived by Zakharov & Dyachenko (2011). Multiple-scale solutions reveal that a perturbation to a slightly unstable uniform wavetrain of steepness μslowly evolves according to a Nonlinear Schrodinger equation. In particular, for small carrier wave steepness μ<μ_1~0.27 the perturbation dynamics is of focusing type and the long-time behavior is characterized by Fermi-Pasta-Ulam recurrence, the signature of breather interactions. However, the amplitude of breathers and their likelihood of occurrence tend to diminish as μincreases while the Benjamin-Feir index decreases and becomes nil at \mu1. Thus, homoclinic orbits persist only for small values of wave steepness μ<<μ_1, in agreement with recent experimental and numerical observations of breathers. When the Zakharov equation is already beyond its range of nominal validity for μ>μ_1, predictions seem to foreshadow a dynamical trend to wave breaking. In particular, the perturbation dynamics becomes of defocussing type, and nonlinearities tend to stabilize a linearly unstable wavetrain as Fermi-Pasta-Ulam recurrence is suppressed. At μ=μ_c~0.577, subharmonic perturbations restabilize and superharmonic instability appears, possibly indicating that wave dynamic behavior changes at large steepness, in qualitative agreement with the numerical simulations of Longuet-Higgins and Cokelet (1978) for steep waves. Indeed, for μ>μ_c a multiple-scale perturbation analysis reveals that a weak narrowband perturbation to a uniform wavetrain evolves in accord with a modified Korteweg-de Vries/Camassa-Holm type equation, again implying a possible mechanism conducive to initiating wave breaking.

preprint2012arXiv

Geometric numerical schemes for the KdV equation

Geometric discretizations that preserve certain Hamiltonian structures at the discrete level has been proven to enhance the accuracy of numerical schemes. In particular, numerous symplectic and multi-symplectic schemes have been proposed to solve numerically the celebrated Korteweg-de Vries (KdV) equation. In this work, we show that geometrical schemes are as much robust and accurate as Fourier-type pseudo-spectral methods for computing the long-time KdV dynamics, and thus more suitable to model complex nonlinear wave phenomena.

preprint2012arXiv

Hamiltonian description and traveling waves of the spatial Dysthe equations

The spatial version of the fourth-order Dysthe equations describe the evolution of weakly nonlinear narrowband wave trains in deep waters. For unidirectional waves, the hidden Hamiltonian structure and new invariants are unveiled by means of a gauge transformation to a new canonical form of the evolution equations. A highly accurate Fourier-type spectral scheme is developed to solve for the equations and validate the new conservation laws, which are satisfied up to machine precision. Further, traveling waves are numerically investigated using the Petviashvili method. It is found that their collision appears inelastic, suggesting the non-integrability of the Dysthe equations.

preprint2012arXiv

Solitary wave interaction in a compact equation for deep-water gravity waves

In this study we compute numerical traveling wave solutions to a compact version of the Zakharov equation for unidirectional deep-water waves recently derived by Dyachenko & Zakharov (2011) Furthermore, by means of an accurate Fourier-type spectral scheme we find that solitary waves appear to collide elastically, suggesting the integrability of the Zakharov equation.

preprint2012arXiv

Special solutions to a compact equation for deep-water gravity waves

Recently, Dyachenko & Zakharov (2011) have derived a compact form of the well known Zakharov integro-differential equation for the third order Hamiltonian dynamics of a potential flow of an incompressible, infinitely deep fluid with a free surface. Special traveling wave solutions of this compact equation are numerically constructed using the Petviashvili method. Their stability properties are also investigated. Further, unstable traveling waves with wedge-type singularities, viz. peakons, are numerically discovered. To gain insights into the properties of singular traveling waves, we consider the academic case of a perturbed version of the compact equation, for which analytical peakons with exponential shape are derived. Finally, by means of an accurate Fourier-type spectral scheme it is found that smooth solitary waves appear to collide elastically, suggesting the integrability of the Zakharov equation.

preprint2012arXiv

Vortexons in axisymmetric Poiseuille pipe flows

We present a study on the nonlinear dynamics of small long-wave disturbances to the laminar state in non-rotating axisymmetric Poiseuille pipe flows. At high Reynolds numbers, the associated Navier-Stokes equations can be reduced to a set of coupled Korteweg-de Vries-type (KdV) equations that support inviscid and smooth travelling waves numerically computed using the Petviashvili method. In physical space they correspond to localized toroidal vortices concentrated near the pipe boundaries (wall vortexons) or that wrap around the pipe axis (centre vortexons), in agreement with the analytical soliton solutions derived by Fedele (2012). The KdV dynamics of a perturbation is also investigated by means of an high accurate Fourier-based numerical scheme. We observe that an initial vortical patch splits into a centre vortexon radiating patches of vorticity near the wall. These can undergo further splitting leading to a proliferation of centre vortexons that eventually decay due to viscous effects. The splitting process originates from a radial flux of azimuthal vorticity from the wall to the pipe axis in agreement with the inverse cascade of cross-stream vorticity identified in channel flows by Eyink (2008). The inviscid vortexon most likely is unstable to non-axisymmetric disturbances and may be a precursor to puffs and slug flow formation.

preprint2011arXiv

Hamiltonian form and solitary waves of the spatial Dysthe equations

The spatial Dysthe equations describe the envelope evolution of the free-surface and potential of gravity waves in deep waters. Their Hamiltonian structure and new invariants are unveiled by means of a gauge transformation to a new canonical form of the evolution equations. An accurate Fourier-type spectral scheme is used to solve for the wave dynamics and validate the new conservation laws, which are satisfied up to machine precision. Traveling waves are numerically constructed using the Petviashvili method. It is shown that their collision appears inelastic, suggesting the non-integrability of the Dysthe equations.