We propose a new two-parameter family of hybrid traveling-standing (TS) water waves in infinite depth that evolve to a spatial translation of their initial condition at a later time. We use the square root of the energy as an amplitude parameter and introduce a traveling parameter that naturally interpolates between pure traveling waves moving in either direction and pure standing waves in one of four natural phase configurations. The problem is formulated as a two-point boundary value problem and a quasi-periodic torus representation is presented that exhibits TS-waves as nonlinear superpositions of counter-propagating traveling waves. We use an overdetermined shooting method to compute nearly 50,000 TS-wave solutions and explore their properties. Examples of waves that periodically form sharp crests with high curvature or dimpled crests with negative curvature are presented. We find that pure traveling waves maximize the magnitude of the horizontal momentum among TS-waves of a given energy. Numerical evidence suggests that the two-parameter family of TS-waves contains many gaps and disconnections where solutions with the given parameters do not exist. Some of these gaps are shown to persist to zero-amplitude in a fourth-order perturbation expansion of the solutions in powers of the amplitude parameter. Analytic formulas for the coefficients of this perturbation expansion are identified using Chebyshev interpolation of solutions computed in quadruple-precision.
- Award ID(s):
- 1800241
- NSF-PAR ID:
- 10358787
- Date Published:
- Journal Name:
- Nonlinearity
- Volume:
- 35
- Issue:
- 5
- ISSN:
- 0951-7715
- Page Range / eLocation ID:
- 2424 to 2473
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Analytic perturbation theory for matrices and operators is an immensely useful mathematical technique. Most elementary introductions to this method have their background in the physics literature, and quantum mechanics in particular. In this note, we give an introduction to this method that is independent of any physics notions, and relies purely on concepts from linear algebra. An additional feature of this presentation is that matrix notation and methods are used throughout. In particular, we formulate the equations for each term of the analytic expansions of eigenvalues and eigenvectors as {\em matrix equations}, namely Sylvester equations in particular. Solvability conditions and explicit expressions for solutions of such matrix equations are given, and expressions for each term in the analytic expansions are given in terms of those solutions. This unified treatment simplifies somewhat the complex notation that is commonly seen in the literature, and in particular, provides relatively compact expressions for the non-Hermitian and degenerate cases, as well as for higher order terms.more » « less
-
Tobias Ekholm (Ed.)We prove nonlinear asymptotic stability of a large class of monotonic shear flows among solutions of the 2D Euler equations in the channel $\mathbb{T}\times[0,1]$. More precisely, we consider shear flows $(b(y),0)$ given by a function $b$ which is Gevrey smooth, strictly increasing, and linear outside a compact subset of the interval $(0,1)$ (to avoid boundary contributions which are incompatible with inviscid damping). We also assume that the associated linearized operator satisfies a suitable spectral condition, which is needed to prove linear inviscid damping. Under these assumptions, we show that if $u$ is a solution which is a small and Gevrey smooth perturbation of such a shear flow $(b(y),0)$ at time $t=0$, then the velocity field $u$ converges strongly to a nearby shear flow as the time goes to infinity. This is the first nonlinear asymptotic stability result for Euler equations around general steady solutions for which the linearized flow cannot be explicitly solved.more » « less
-
more » « less
Abstract We consider a Celestial Mechanics model: the spin–orbit problem with a dissipative tidal torque, which is a singular perturbation of a conservative system. The goal of this paper is to show that it is possible to maintain the accuracy and reliability of the computation of quasi-periodic attractors for parameter values extremely close to the breakdown and, therefore, it is possible to obtain information on the breakdown mechanism of these quasi-periodic attractors. The method uses at the same time numerical and rigorous improvements to provide (i) a very accurate computation of the time-1 map of the spin–orbit problem (which reduces the dimensionality of the problem); (ii) a very efficient KAM method for maps which computes the attractor and its tangent spaces (by quadratically convergent, low storage requirements, and low operation count); (iii) explicit algorithms backed by a rigorous a posteriori KAM theorem, which establishes that if the algorithm is successful and produces a small residual, then there is a true solution nearby; and (iv) guaranteed algorithms to reach arbitrarily close to the border of existence as long as there are enough computer resources. As a by-product of the accuracy that we maintain till breakdown, we study several scale-invariant observables of the tori used in the renormalization group of infinite-dimensional spaces. In contrast with previously studied simple models, the behavior at breakdown of the spin–orbit problem does not satisfy standard scaling relations which implies that the spin–orbit problem is not described by a hyperbolic fixed point of the renormalization operator.
-
more » « less
SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1361-6544/ac5cd5
