--- title: 'NLD-HONLS: Nonlinear Mean-Flow Damping Model' url: https://www.emergentmind.com/topics/nonlinear-mean-flow-damping-model-nld-honls type: topic --- # NLD-HONLS: Nonlinear Mean-Flow Damping Model The nonlinear mean-flow damping model, usually denoted NLD-HONLS, is a dissipative higher-order nonlinear Schrödinger formulation in which the Dysthe-type nonlocal mean-flow feedback is assigned a complex coefficient, so that dissipation enters through the same carrier–sideband interaction channels as the conservative mean-flow term rather than through uniform modewise decay. In the periodic deep-water setting studied in recent work, NLD-HONLS is the \(\Gamma=0,\beta>0\) specialization of a unified damped HONLS family, and its distinctive signature is the appearance of interaction-dependent dissipative phase corrections that preserve organized Floquet-band dynamics during recurrent focusing while also producing permanent downshifting and strongly localized damping near steep crests [2605.26811] [2507.20375]. ## 1. Governing equation and model class In the formulation used for spatially periodic complex envelopes \(u(x,t)\) of period \(L\), with \(\mu=2\pi/L\), the damped higher-order NLS family is $$ i u_t + u_{xx} + 2 |u|^2 u + \epsilon \left[ 2 u (1 + i \beta) \mathscr{H}\!\left((|u|^2)_x\right) - 8 |u|^2 u_x + \frac{1}{2} u_{xxx} \right] + i \Gamma u + 2 i \epsilon \Gamma u_x = 0, $$ where \(\epsilon>0\) is a small higher-order parameter, \(\Gamma\ge 0\) is the viscous-damping coefficient, \(\beta\ge 0\) is the nonlinear mean-flow damping coefficient, and \(\mathscr H\) denotes the Hilbert transform [2605.26811]. The leading \(u_{xx}\) and \(2|u|^2u\) terms are the focusing NLS dispersion and cubic nonlinearity; \(\frac12 u_{xxx}\) is the third-order dispersive correction; \(-8|u|^2u_x\) is the self-steepening correction; and the nonlocal mean-flow feedback is carried by the Hilbert-transform term. Three closely related regimes are distinguished in the literature. | Regime | Parameters | Dominant damping mechanism | |---|---|---| | HONLS | \(\beta=\Gamma=0\) | None | | NLD-HONLS | \(\Gamma=0,\ \beta>0\) | Nonlocal interaction-dependent damping | | V-HONLS | \(\beta=0,\ \Gamma>0\) | Linear modewise viscous decay | The defining feature of NLD-HONLS is the factor \((1+i\beta)\) multiplying the conservative mean-flow interaction. Its real part yields the usual HONLS mean-flow feedback, while the imaginary part introduces nonlinear dissipation through the same nonlocal structure. In the physical interpretation adopted in later Floquet studies, this damping is steepness-selective and localized near strongly modulated crests, and it is associated with dissipative mean-flow response induced by viscosity, turbulence, or micro-breaking [2507.20375]. By contrast, the viscous terms \(i\Gamma u+2i\epsilon \Gamma u_x\) are diagonal in Fourier space and primarily produce amplitude decay without directly entering the leading interaction-phase equations [2605.26811]. The energy–momentum budgets make this distinction explicit. For the unified model, $$ \frac{dE}{dt} = -\big(2\Gamma E + 4\beta B + 2\Gamma P\big),\qquad \frac{dP}{dt} = -\big(2\Gamma P - 4\beta S + 4\Gamma Q\big), $$ with \(E,P,B,Q,S\) defined by spatial averages of \(|u|^2\), \(u_x\), and mean-flow couplings. In the NLD-HONLS regime \((\Gamma=0)\), one has \(dE/dt=-4\beta B<0\) and \(dP/dt=+4\beta S\), showing that energy decay and momentum evolution are governed by the nonlocal damping rather than by uniform linear attenuation [2507.20375]. ## 2. Five-mode reduction and carrier–sideband structure To isolate the dominant coherent exchange between the carrier and its first few sidebands, the 2026 interaction-phase study adopts a closed five-mode truncation, $$ u(x,t)=\sum_{n=-2}^{2} A_n(t)e^{in\mu x},\qquad A_n(t)\in\mathbb C, $$ retaining the carrier \(n=0\), first sidebands \(n=\pm1\), and second sidebands \(n=\pm2\) [2605.26811]. Projecting the PDE onto these Fourier modes yields a finite-dimensional nonlinear system in which linear dispersion, cubic interactions, mean-flow couplings, self-steepening, and dissipation remain explicitly separated. The linear dispersive coefficient is $$ \Omega_n=\mu^2 n^2-\frac{\epsilon}{2}\mu^3 n^3, $$ while the projected quadratic spectral field \(|u|^2=\sum_{p=-4}^{4}Q_p e^{ip\mu x}\) determines the cubic and mean-flow interaction coefficients \(C_n\), \(M_n\), and \(D_n\). The resulting modal equations are $$ i \dot A_n = \Omega_n A_n - 2 C_n + 2 \epsilon \mu (1+i\beta) M_n - 8\epsilon\mu D_n - i\Gamma A_n - 2 i \epsilon \Gamma \mu n A_n,\qquad n=-2,-1,0,1,2. $$ This decomposition makes the structural difference between the two dissipative mechanisms transparent. In V-HONLS, dissipation stays diagonal through \(-i\Gamma A_n-2i\epsilon\Gamma\mu nA_n\), so the nonlinear interaction tensors \(C_n\), \(M_n\), and \(D_n\) remain unchanged. In NLD-HONLS, the dissipative contribution is \(2i\epsilon\beta\mu M_n\), so the damping is injected directly into the same projected mean-flow interaction operator that governs conservative carrier–sideband exchange [2605.26811]. Writing \(A_n=r_n e^{i\phi_n}\) yields amplitude and phase equations. The amplitude dynamics contain both conservative and dissipative projections of \(C_n\), \(M_n\), and \(D_n\), together with the viscous decay terms. The phase dynamics contain the modal detuning \(-\Omega_n\) and the real or imaginary projections of the same interaction coefficients divided by \(r_n\). This representation is central because the carrier–sideband dynamics are governed less by the individual modal phases \(\phi_n\) than by specific four-wave phase combinations. ## 3. Interaction-phase dynamics and the \(-\kappa_j\sin(\psi_j)\) mechanism The dominant carrier–sideband interaction phases are $$ \psi_1=\phi_1+\phi_{-1}-2\phi_0,\qquad \psi_2=\phi_2+\phi_{-2}-2\phi_0, $$ the arguments of the principal four-wave products \(A_1A_{-1}\overline{A_0}^2\) and \(A_2A_{-2}\overline{A_0}^2\) [2605.26811]. Differentiation yields exact phase-balance equations of the form $$ \dot\psi_j=\Delta_j+\mathcal C_j+\mathcal M_j+\mathcal D_j,\qquad j\in\{1,2\}, $$ where \(\Delta_j\) is the linear phase mismatch and \(\mathcal C_j,\mathcal M_j,\mathcal D_j\) are the projected cubic, mean-flow, and self-steepening contributions. The decisive simplification occurs in the carrier–sideband regime $$ r_0 \gg r_{\pm1} \gg r_{\pm2}. $$ In that regime, the dissipative part of the mean-flow projection contributes leading phase-dependent terms proportional to \(\sin(\psi_j)\), and the interaction-phase equations acquire the explicit corrections $$ \dot\psi_j \supset -\kappa_j\sin(\psi_j),\qquad j=1,2, $$ with $$ \kappa_j=2\beta\epsilon\mu\left[jr_0^2\left(\frac{r_{-j}}{r_j}+\frac{r_j}{r_{-j}}\right)+4jr_jr_{-j}\right]. $$ Equivalently, $$ \kappa_1=2\beta\epsilon\mu\left[r_0^2\left(\frac{r_{-1}}{r_1}+\frac{r_1}{r_{-1}}\right)+4r_1r_{-1}\right], $$ $$ \kappa_2=2\beta\epsilon\mu\left[2r_0^2\left(\frac{r_{-2}}{r_2}+\frac{r_2}{r_{-2}}\right)+8r_2r_{-2}\right]. $$ These terms are interaction-dependent dissipative corrections produced by the imaginary part of the mean-flow projection, not by any diagonal damping channel [2605.26811]. The interpretation given in the 2026 analysis is precise. The \(-\kappa_j\sin(\psi_j)\) terms provide restoring-type feedback within the dominant carrier–sideband interaction dynamics. They do not imply rigid phase locking, but they oppose sustained monotone drift away from phase configurations favorable to recurrent focusing. This explains why NLD-HONLS can exhibit substantial interaction-phase restructuring without losing spectral organization. A common misconception is that strong phase evolution necessarily signals diffuse or disordered modulation; the finite-gap benchmarks used in the same study show that localized restructuring and cumulative drift in one interaction phase can coexist with organized quasiperiodic Floquet structure [2605.26811]. By contrast, viscous damping does not produce any direct \(-\kappa_j\sin(\psi_j)\) contribution at leading order. It alters \(\psi_j\) only indirectly through its effect on the amplitudes \(r_n\). This distinction is central to the observed contrast between recurrent carrier–sideband focusing under NLD-HONLS and progressively diffuse multimode evolution under V-HONLS [2605.26811]. ## 4. Floquet spectrum, spectral organization, and finite-gap interpretation The Floquet spectral framework used throughout this literature is inherited from the focusing NLS Zakharov–Shabat problem. For an \(L\)-periodic potential \(u(x,t)\), the spatial operator is $$ \mathcal L^{(x)}v= \begin{pmatrix} \partial_x+i\lambda & -u\\ \overline u & \partial_x-i\lambda \end{pmatrix}v=0, $$ with spectral parameter \(\lambda\in\mathbb C\). The monodromy over one spatial period defines the Floquet discriminant $$ \Delta(u,\lambda)=\mathrm{Trace}\big(\Psi(x+L;\lambda)\Psi^{-1}(x;\lambda)\big), $$ and the Floquet spectrum is $$ \sigma(u)=\{\lambda\in\mathbb C:\Delta(u,\lambda)\in\mathbb R,\ |\Delta(u,\lambda)|\le 2\}. $$ Simple periodic or antiperiodic points are zeros of \(\Delta\pm2\) with nonzero derivative; double points satisfy \(\partial_\lambda\Delta=0\) as well; and critical points satisfy \(\partial_\lambda\Delta=0\) without necessarily obeying \(\Delta=\pm2\) [2605.26811]. The 2026 comparison between nonlinear mean-flow damping and viscous damping distinguishes two types of nonintegrable spectral evolution. In an “organized” evolution, dominant bands remain separated and localized in the upper half-plane, band lengths vary smoothly, and no critical-point crossings occur. In a “reconnection” regime, bands repeatedly cross through critical points, connectivity changes, band identities swap, and the spectrum becomes diffuse and multimode [2605.26811]. The localization criterion used both in the 2026 interaction-phase paper and in the earlier rogue-wave studies is $$ |\gamma(t;\lambda_m,\lambda_n)|<0.025, $$ so that a band of length below \(0.025\) is treated as localized or “soliton-like” [2605.26811] [2203.13488]. Finite-gap NLS benchmarks play an interpretive role. A one-mode finite-gap benchmark exhibits bounded oscillatory interaction phases and negligible spectral drift, while a symmetric two-mode quasiperiodic benchmark shows recurrent focusing–defocusing cycles in which one interaction phase remains relatively bounded and the other undergoes cumulative drift with localized restructuring near focusing events. In both cases the Floquet spectrum stays organized in the integrable limit, demonstrating that pronounced interaction-phase evolution does not by itself imply spectral disorder [2605.26811]. Spatially periodic breather initial data are then used to launch damped HONLS simulations. These SPBs are relatives of the Akhmediev breather and Kuznetsov–Ma soliton and provide controlled access to modulational-instability-driven focusing cycles [2203.13488]. This spectral viewpoint also underlies the “soliton-like rogue wave” classification. Under NLD-HONLS, tiny upper-half-plane bands can pinch off, generating one- or two-soliton-like spectral states from which strongly localized rogue events emerge. In early-to-middle modulational instability, all rogue waves in the 2022 SPB study occur while the spectrum is one- or two-soliton-like, whereas near instability saturation rogue waves can also arise after the spectrum has left the soliton-like state [2203.13488]. ## 5. Numerical framework and observed dynamics Across the recent NLD-HONLS studies, the computational setting is a periodic domain, typically \(L=4\sqrt2\,\pi\), with Fourier pseudo-spectral discretization in space, \(N=256\) modes, and ETDRK4 time stepping with \(\Delta t=10^{-3}\). Representative higher-order parameter values are \(\epsilon=0.05\); NLD-HONLS simulations use \(\beta\in[0.1,0.8]\), often with \(\beta=0.1\), while V-HONLS comparisons use values such as \(\Gamma=0.002\) [2605.26811] [2507.20375]. Initial data include coalesced two-mode SPBs \(U(x,0;0,0)\), other SPB snapshots indexed by \(T_0\in[-5,0]\), and perturbed Stokes waves of the form \(u(x,0)=a(1+\alpha\cos(\mu x))\) with \(a=0.45\) and \(\alpha=10^{-2}\) [2203.13488] [2507.20375]. The diagnostics are correspondingly multi-layered. Floquet spectra are computed by solving the Zakharov–Shabat system numerically, constructing \(\Delta(u,\lambda)\), locating zeros of \(\Delta\pm2\) with Müller’s method, and tracking critical points via \(\partial_\lambda\Delta=0\) [2605.26811]. Interaction phases are extracted from the four-wave products \(I_j=A_jA_{-j}\overline{A_0}^2\), spectral drift is monitored through the dominant Fourier mode \(k_{\mathrm{peak}}(t)\), and focusing is quantified either by the maximum envelope amplitude \(S(t)=\|u\|_\infty\) or by the strength function \(S(t)=U_{\max}/H_s\), with rogue-wave threshold \(S\ge 2.2\) in the SPB and SRW studies [2605.26811] [2507.20375]. The numerical findings are unusually consistent across studies. In the interaction-phase analysis, NLD-HONLS exhibits persistent recurrent carrier–sideband focusing together with organized Floquet evolution, even when one interaction phase undergoes substantial restructuring. The spectrum avoids critical-point crossings, and the recurrent carrier–sideband exchange remains coherent. V-HONLS instead shows progressively diffuse modulation dynamics, repeated band reconnection through both real and complex critical points, weakening persistence of the recurrent focusing structure, and loss of clear carrier–sideband organization [2605.26811]. The 2025 Floquet comparison sharpens this picture. For steep SPB data at \(\beta=0.1\), NLD-HONLS supports a two-mode soliton-like regime up to approximately \(t\approx 40\), then a one-mode soliton-like regime up to approximately \(t\approx 90\), with no real or complex critical-point crossings observed on \(0