Numerical Bifurcation and Spectral Stability of Wavetrains in Bidirectional Whitham Models

We consider several different bidirectional Whitham equations that have recently appeared in the literature. Each of these models combines the full two‐way dispersion relation from the incompressible Euler equations with a canonical shallow water nonlinearity, providing nonlocal model equations that may be expected to exhibit some of the interesting high‐frequency phenomena present in the Euler equations that standard “long‐wave” theories fail to capture. Of particular interest here is the existence and stability of periodic traveling wave solutions in such models. Using numerical bifurcation techniques, we construct global bifurcation diagrams for each system and compare the global structure of branches, together with the possibility of bifurcation branches terminating in a “highest” singular (peaked/cusped) wave. We also numerically approximate the stability spectrum along these bifurcation branches and compare the stability predictions of these models. Our results confirm a number of analytical results concerning the stability of asymptotically small waves in these models and provide new insights into the existence and stability of large amplitude waves.     

On the structure and growth rate of unstable modes to the Rossby‐Haurwitz wave

The normal mode instability study of a steady Rossby‐Haurwitz wave is considered both theoretically and numerically. This wave is exact solution of the nonlinear barotropic vorticity equation describing the dynamics of an ideal fluid on a rotating sphere, as well as the large‐scale barotropic dynamics of the atmosphere. In this connection, the stability of the Rossby‐Haurwitz wave is of considerable mathematical and meteorological interest. The structure of the spectrum of the linearized operator in case of an ideal fluid is studied. A conservation law for perturbations to the Rossby‐Haurwitz wave is obtained and used to get a necessary condition for its exponential instability. The maximum growth rate of unstable modes is estimated. The orthogonality of the amplitude of a non‐neutral or non‐stationary mode to the Rossby‐Haurwitz wave is shown in two different inner products. The analytical results obtained are used to test and discuss the accuracy of a numerical spectral method used for the normal mode stability study of arbitrary flow on a sphere. The comparison of the numerical and theoretical results shows that the numerical instability study method works well in case of such smooth solutions as the zonal flows and Rossby‐Haurwitz waves. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Linear stability of solitary waves of the Green‐Naghdi equations

We investigate the eigenvalue problem obtained from linearizing the Green‐Naghdi equations about solitary wave solutions. Unlike weakly nonlinear water wave models, the physical system considered here has nonlinearity in its highest derivative term. This results in more detailed asymptotic analysis of the eigenvalue problem in the presence of a large parameter. Combining the technique of singular perturbation with the Evans function, we show that for solitary waves of small amplitude, the problem has no eigenvalues of positive real part and the Evans function is nonvanishing everywhere except the origin. This fact then leads to the linear stability of these solitary waves. © 2001 John Wiley & Sons, Inc.     

Stability properties of steady water waves with vorticity

We present two stability analyses for exact periodic traveling water waves with vorticity. The first approach leads in particular to linear stability properties of water waves for which the vorticity decreases with depth. The second approach leads to a formal stability property for long water waves that have small vorticity and amplitude although we do not use a small‐amplitude or long‐wave approximation. © 2006 Wiley Periodicals, Inc.     

Time‐dependent scattering of generalized plane waves by a wedge

We obtain explicit formulas for the scattering of plane waves with arbitrary profile by a wedge under Dirichlet, Neumann and Dirichlet‐Neumann boundary conditions. The diffracted wave is given by a convolution of the profile function with a suitable kernel corresponding to the boundary conditions. We prove the existence and uniqueness of solutions in appropriate classes of distributions and establish the Sommerfeld type representation for the diffracted wave. As an application, we establish (i) stability of long‐time asymptotic local perturbations of the profile functions and (ii) the limiting amplitude principle in the case of a harmonic incident wave. Copyright © 2015 John Wiley & Sons, Ltd.     

Analyticity for Kuramoto–Sivashinsky‐type equations in two spatial dimensions

I. Stratis In this work, we investigate the analyticity properties of solutions of Kuramoto–Sivashinsky‐type equations in two spatial dimensions, with periodic initial data. In order to do this, we explore the applicability in three‐dimensional models of a spectral method, which was developed by the authors for the one‐dimensional Kuramoto–Sivashinsky equation. We introduce a criterion, which provides a sufficient condition for analyticity of a periodic function u∈C^∞, involving the rate of growth of ?ⁿu, in suitable norms, as n tends to infinity. This criterion allows us to establish spatial analyticity for the solutions of a variety of systems, including Topper–Kawahara, Frenkel–Indireshkumar, and Coward–Hall equations and their dispersively modified versions, once we assume that these systems possess global attractors. Copyright © 2015 John Wiley & Sons, Ltd.     

Complexiton solutions to the Hirota‐Satsuma‐Ito equation

The Hirota bilinear method is a powerful tool for solving nonlinear evolution equations. Together with the linear superposition principle, it can be used to find a special class of explicit solutions that correspond to complex eigenvalues of associated characteristic problems. These solutions are known as complexiton solutions or simply complexitons. In this article, we study complexiton solutions of the the Hirota‐Satsuma‐Ito equation which is a (2 + 1)‐dimensional extension of the Hirota‐Satsuma shallow water wave equation known to describe propagation of unidirectional shallow water waves. We first construct hyperbolic function solutions and consequently derive the so‐called complexitons via the Hirota bilinear method and the linear superposition principle. In particular, we find nonsingular complexiton solutions to the Hirota‐Satsuma‐Ito equation. Finally, we give some illustrative examples and a few concluding remarks.     

Instability of cnoidal‐peak for the NLS‐δ equation

We study the existence and stability of standing waves for the periodic cubic nonlinear Schrödinger equation with a point defect determined by the periodic Dirac distribution at the origin. We show that this model admits a smooth curve of periodic‐peak standing wave solutions with a profile determined by the Jacobi elliptic function of cnoidal type. Via a perturbation method and continuation argument, we obtain that in the repulsive defect, the cnoidal‐peak standing wave solutions are unstable in $H^1_{per}$ with respect to perturbations which have the same period as the wave itself. Global well‐posedness is verified for the Cauchy problem in $H^1_{per}$ .     

On resonant interactions of gravity‐capillary waves without energy exchange

We consider resonant triad interactions of gravity‐capillary waves and investigate in detail special resonant triads that exchange no energy during their interactions so that the wave amplitudes remain constant in time. After writing the resonance conditions in terms of two parameters (or two angles of wave propagation), we first identify a region in the two‐dimensional parameter space, where resonant triads can be always found, and then describe the variations of resonant wavenumbers and wave frequencies over the resonance region. Using the amplitude equations recovered from a Hamiltonian formulation for water waves, it is shown that any resonant triad inside the resonance region can interact without energy exchange if the initial wave amplitudes and relative phase satisfy the two conditions for fixed point solutions of the amplitude equations. Furthermore, it is shown that the symmetric resonant triad exchanging no energy forms a transversely modulated traveling wave field, which can be considered a two‐dimensional generalization of Wilton ripples.     

The gap lemma and geometric criteria for instability of viscous shock profiles

An obstacle in the use of Evans function theory for stability analysis of traveling waves occurs when the spectrum of the linearized operator about the wave accumulates at the imaginary axis, since the Evans function has in general been constructed only away from the essential spectrum. A notable case in which this difficulty occurs is in the stability analysis of viscous shock profiles. Here we prove a general theorem, the “gap lemma,” concerning the analytic continuation of the Evans function associated with the point spectrum of a traveling wave into the essential spectrum of the wave. This allows geometric stability theory to be applied in many cases where it could not be applied previously. We demonstrate the power of this method by analyzing the stability of certain undercompressive viscous shock waves. A necessary geometric condition for stability is determined in terms of the sign of a certain Melnikov integral of the associated viscous profile. This sign can easily be evaluated numerically. We also compute it analytically for solutions of several important classes of systems. In particular, we show for a wide class of systems that homoclinic (solitary) waves are linearly unstable, confirming these as the first known examples of unstable viscous shock waves. We also show that (strong) heteroclinic undercompressive waves are sometimes unstable. Similar stability conditions are also derived for Lax and overcompressive shocks and for n × n conservation laws, n ≥ 2. © 1998 John Wiley & Sons, Inc.     

Hamiltonian long‐wave expansions for free surfaces and interfaces

The theory of internal waves between two bodies of immiscible fluid is important both for its interest to ocean engineering and as a source of numerous interesting mathematical model equations that exhibit nonlinearity and dispersion. In this paper we derive a Hamiltonian formulation of the problem of a dynamic free interface (with rigid lid upper boundary conditions), and of a free surface and a free interface, this latter situation occurring more commonly in experiment and in nature. From the formulation, we develop a Hamiltonian perturbation theory for the long‐wave limits, and we carry out a systematic analysis of the principal long‐wave scaling regimes. This analysis provides a uniform treatment of the classical works of Peters and Stoker (28), Benjamin (3, 4), Ono (26), and many others. Our considerations include the Boussinesq and Korteweg–de Vries (KdV) regimes over finite‐depth fluids, the Benjamin‐Ono regimes in the situation in which one fluid layer is infinitely deep, and the intermediate long‐wave regimes. In addition, we describe a novel class of scaling regimes of the problem, in which the amplitude of the interface disturbance is of the same order as the mean fluid depth, and the characteristic small parameter corresponds to the slope of the interface. Our principal results are that we highlight the discrepancies between the case of rigid lid and of free surface upper boundary conditions, which in some circumstances can be significant. Motivated by the recent results of Choi and Camassa (6, 7), we also derive novel systems of nonlinear dispersive long‐wave equations in the large‐amplitude, small‐slope regime. Our formulation of the dynamical free‐surface, free‐interface problem is shown to be very effective for perturbation calculations; in addition, it holds promise as a basis for numerical simulations. © 2005 Wiley Periodicals, Inc.     

Rotational waves generated by current‐topography interaction

We study nonlinear free‐surface rotational waves generated through the interaction of a vertically sheared current with a topography. Equivalently, the waves may be generated by a pressure distribution along the free surface. A forced Korteweg–de Vries equation (fKdV) is deduced incorporating these features. The weakly nonlinear, weakly dispersive reduced model is valid for small amplitude topographies. To study the effect of gradually increasing the topography amplitude, the free surface Euler equations are formulated in the presence of a variable depth and a sheared current of constant vorticity. Under constant vorticity, the harmonic velocity component is formulated in a simplified canonical domain, through the use of a conformal mapping which flattens both the free surface as well as the bottom topography. Critical, supercritical, and subcritical Froude number regimes are considered, while the bottom amplitude is gradually increased in both the irrotational and rotational wave regimes. Solutions to the fKdV model are compared to those from the Euler equations. We show that for rotational waves the critical Froude number is shifted away from 1. New stationary solutions are found and their stability tested numerically.     

Existence and regularizing rate estimates of solutions to the 3‐D generalized micropolar system in Fourier‐Besov spaces

In this paper, we study global existence and asymptotic stability of solutions for the initial value problem of the three‐dimensional (3‐D) generalized incompressible micropolar system in Fourier‐Besov spaces. Besides, we also establish some regularizing rate estimates of the higher‐order spatial derivatives of solutions, which particularly imply the spatial analyticity and the temporal decay of global solutions.     

2‐D solitary waves in thermal media with nonsymmetric boundary conditions

Optical solitary waves and their stability in focusing thermal optical media, such as lead glasses, are studied numerically and theoretically in (2 + 1) dimensions. The optical medium is a square cell and mixed boundary conditions of Newton cooling and fixed temperature on different sides of the cell are used. Nonlinear thermal optical media have a refractive index which depends on temperature, so that heating from the optical beam and heat flow across the boundaries can change the refractive index of the medium. Solitary wave solutions are found numerically using the Newton conjugate‐gradient method, while their stability is studied using a linearized stability analysis and also via numerical simulations. It is found that the position of the solitary wave is dependent on the boundary conditions, with the center of the beam moving toward the warmer boundaries, as the parameters are varied. The stability of the solitary waves depends on the symmetry of the boundary conditions and the amplitude of the solitary waves.     

Stability of compact upwind‐biased finite‐difference schemes forfeedback‐coupled advection equations

A fifth‐order compact upwind‐biased finite‐difference scheme has been developed which is asymptotically stable when applied to linear 2 × 2 systems. This scheme was optimised with respect to wave propagation through the boundaries. In order to achieve asymptotic stability for such systems a sufficient stability condition, based on the Nyquist criterion of linear system theory, was included in the optimisation as additional constraint. The evaluation of the stability condition was done numerically. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

High‐order Padé and singly diagonally Runge‐Kutta schemes for linear ODEs,application to wave propagation problems

In this article, we address the problem of constructing high‐order implicit time schemes for wave equations. We consider two classes of one‐step A‐stable schemes adapted to linear Ordinary Differential Equation (ODE). The first class, which is not dissipative is based on the diagonal Padé approximant of exponential function. For this class, the obtained schemes have the same stability function as Gauss Runge‐Kutta (Gauss RK) schemes. They have the advantage to involve the solution of smaller linear systems at each time step compared to Gauss RK. The second class of schemes are constructed such that they require the inversion of a unique linear system several times at each time step like the Singly Diagonally Runge‐Kutta (SDIRK) schemes. While the first class of schemes is constructed for an arbitrary order of accuracy, the second‐class schemes is given up to order 12. The performance assessment we provide shows a very good level of accuracy for both classes of schemes, and the great interest of considering high‐order time schemes that are faster. The diagonal Padé schemes seem to be more accurate and more robust.     

Three‐dimensional internal gravity‐capillary waves in finite depth

We consider three‐dimensional inviscid‐irrotational flow in a two‐layer fluid under the effects of gravity and surface tension, where the upper fluid is bounded above by a rigid lid and the lower fluid is bounded below by a flat bottom. We use a spatial dynamics approach and formulate the steady Euler equations as an infinite‐dimensional Hamiltonian system, where an unbounded spatial direction x is considered as a time‐like coordinate. In addition, we consider wave motions that are periodic in another direction z. By analyzing the dispersion relation, we detect several bifurcation scenarios, two of which we study further: a type of 00(is)(iκ₀) resonance and a Hamiltonian Hopf bifurcation. The bifurcations are investigated by performing a center‐manifold reduction, which yields a finite‐dimensional Hamiltonian system. For this finite‐dimensional system, we establish the existence of periodic and homoclinic orbits, which correspond to, respectively, doubly periodic travelling waves and oblique travelling waves with a dark or bright solitary wave profile in the x direction. The former are obtained using a variational Lyapunov‐Schmidt reduction and the latter by first applying a normal form transformation and then studying the resulting canonical system of equations.     

THE LONG‐RUN INEFFICIENCY OF BLOCK‐RATE PRICING

ABSTRACT. In this paper we compare two regulation instruments, flat‐rate and increasing block‐rate pricing. The analysis applies to a competitive industry with free entry. Charge for irrigation water is a concrete example. It is shown that flat‐rate pricing leads to a first‐best social optimum, while with block‐rate pricing where the highest block set at the marginal cost, there is over production, firms are too small, and loss of economic surplus occurs. Moreover, first‐best is not implementable by increasing block‐rate pricing. This is in contrast to the commonly accepted view that block‐rate pricing is superior to flat‐rate pricing by allowing for income redistribution while preserving efficiency. Several second‐best situations are analyzed to show: 1) Block‐rate pricing with the highest block at the social marginal cost is optimal when the regulator must preserve the number of firms. 2) Water pricing alone cannot implement social optimum subject to a constant level of agricultural production. 3) Lobbying and political pressures, which force the regulator to sustain a constant average water price, result in optimal block‐rate pricing with the highest block below the social marginal cost.     

Regularity of Traveling Free Surface Water Waves with Vorticity

We prove real analyticity of all the streamlines, including the free surface, of a gravity- or capillary-gravity-driven steady flow of water over a flat bed, with a Hölder continuous vorticity function, provided that the propagating speed of the wave on the free surface exceeds the horizontal fluid velocity throughout the flow. Furthermore, if the vorticity possesses some Gevrey regularity of index s, then the stream function of class C ^2,μ admits the same Gevrey regularity throughout the fluid domain; in particular if the Gevrey index s equals 1, then we obtain analyticity of the stream function. The regularity results hold not only for periodic or solitary-water waves, but also for any solution to the hydrodynamic equations of class C ^2,μ .     

Neimark‐Sacker bifurcation of a fourth order difference equation

In this article, we study stability and bifurcation of a fourth order rational difference equation. We give condition for local stability, and we show that the equation undergoes a Neimark‐Sacker bifurcation. Moreover, we consider the direction of the Neimark‐Sacker bifurcation. Finally, we numerically validate our analytical results.