Oscillatory integrals generated by the Ostrovsky equation

The Ostrovsky equation governs the propagation of long nonlinear surface waves in the presence of rotation. It is related to the Korteweg-de Vries (KdV) and the Kadomtsev-Petviashvili models. KdV can be obtained from the equation in question when the rotation parameter γ equals zero. A fundamental solution of the Cauchy problem for the linear Ostrovsky equation is presented in the form of an oscillatory Fourier integral. Another integral representation involving Airy and Bessel functions is derived for it. It is shown that its asymptotic expansion as γ → 0 contains the KdV fundamental solution as the zero term. The Airy transform is used to establish some of its properties. Higher-order asymptotics for γ → 0 on a bounded time interval are obtained for both the fundamental solution and the solution of the linear Cauchy problem for the Ostrovsky equation.     

Asymptotic properties of solutions to two discrete airy equations

We construct two finite difference models for the Airy differential equation. In one model, the form of the complete asymptotic representation of the solution can be found. However, this is not the case for the second model which is based on the use of a nonstandard difference scheme. This scheme leads to a second-order, linear difference equation that is not of a form for which the theorems of Poincaré and Perron can be directly applied to obtain the asymptotic behavior of the solutions.     

Asymptotic solutions to a discrete airy equation

A nonstandard finite difference scheme for the Airy equation leads to a linear, second—order difference equation for which the theorems of Poincaré and Perron do not apply if asymptotic representations of the solutions are desired. Using the method of dominant balance, a suggested form for the asymptotic solutions is obtained. This relation is then used to construct the required asymptotic representations for the two linearly independent solutions     

Oscillatory integrals generated by the Ostrovsky equation

The Ostrovsky equation governs the propagation of long nonlinear surface waves in the presence of rotation. It is related to the Korteweg-de Vries (KdV) and the Kadomtsev-Petviashvili models. KdV can be obtained from the equation in question when the rotation parameter γ equals zero. A fundamental solution of the Cauchy problem for the linear Ostrovsky equation is presented in the form of an oscillatory Fourier integral. Another integral representation involving Airy and Bessel functions is derived for it. It is shown that its asymptotic expansion as γ → 0 contains the KdV fundamental solution as the zero term. The Airy transform is used to establish some of its properties. Higher-order asymptotics for γ → 0 on a bounded time interval are obtained for both the fundamental solution and the solution of the linear Cauchy problem for the Ostrovsky equation. Received: November 23, 2004; revised: March 13, 2005 Research is supported by US Department of Defense, under grant No. DAAD19-03-1-0204     

Maslov’s canonical operator,Hörmander’s formula,and localization of the Berry-Balazs solution in the theory of wave beams

For three-dimensional Schrödinger equations, we study how to localize exact solutions represented as the product of an Airy function (Berry-Balazs solutions) and a Bessel function and known as Airy-Bessel beams in the paraxial approximation in optics. For this, we represent such solutions in the form of Maslov’s canonical operator acting on compactly supported functions on special Lagrangian manifolds. We then use a result due to Hörmander, which permits using the formula for the commutation of a pseudodifferential operator with Maslov’s canonical operator to “move” the compactly supported amplitudes outside the canonical operator and thus obtain effective formulas preserving the structure based on the Airy and Bessel functions. We discuss the influence of dispersion effects on the obtained solutions.     

Eigenfunctions of the Airyʼs integral transform: Properties,numerical computations and asymptotic behaviors

This paper is devoted to the study of the spectrum of the integral operator with Airy?s kernel. We provide the reader with some efficient methods of computing the eigenfunctions and the eigenvalues of this operator. The first method is based on the use of a differential operator which commutes with this integral operator. The second method is based on a Gaussian quadrature method applied to the finite Airy?s transform. The asymptotic behavior of the eigenfunctions is studied by the use of a WKB method. Some numerical examples are given to illustrate results of this work.     

The Generating Function for the Airy Point Process and a System of Coupled Painlevé II Equations

For a wide class of Hermitian random matrices, the limit distribution of the eigenvalues close to the largest one is governed by the Airy point process. In such ensembles, the limit distribution of the k th largest eigenvalue is given in terms of the Airy kernel Fredholm determinant or in terms of Tracy–Widom formulas involving solutions of the Painlevé II equation. Limit distributions for quantities involving two or more near‐extreme eigenvalues, such as the gap between the k th and the ℓth largest eigenvalue or the sum of the k largest eigenvalues, can be expressed in terms of Fredholm determinants of an Airy kernel with several discontinuities. We establish simple Tracy–Widom type expressions for these Fredholm determinants, which involve solutions to systems of coupled Painlevé II equations, and we investigate the asymptotic behavior of these solutions.     

Lie group symmetry analysis for granular media stress equations

The Airy stress function, although frequently employed in classical linear elasticity, does not receive similar usage for granular media problems. For plane strain quasi-static deformations of a cohesionless Coulomb-Mohr granular solid, a single nonlinear partial differential equation is formulated for the Airy stress function by combining the equilibrium equations with the yield condition. This has certain advantages from the usual approach, in which two stress invariants and a stress angle are introduced, and a system of two partial differential equations is needed to describe the flow. In the present study, the symmetry analysis of differential equations is utilised for our single partial differential equation, and by computing an optimal system of one-dimensional Lie algebras, a complete set of group-invariant solutions is derived. By this it is meant that any group-invariant solution of the governing partial differential equation (provided it can be derived via the classical symmetries method) may be obtained as a member of this set by a suitable group transformation. For general values of the parameters (angle of internal friction ? and gravity g) it is found there are three distinct classes of solutions which correspond to granular flows considered previously in the literature. For the two limiting cases of high angle of internal friction and zero gravity, the governing partial differential equation admit larger families of Lie point symmetries, and from these symmetries, further solutions are derived, many of which are new. Furthermore, the majority of these solutions are exact, which is rare for granular flow, especially in the case of gravity driven flows.     

Asymptotic Solutions of Two-Dimensional Hartree-Type Equations Localized in the Neighborhood of Line Segments

We consider the eigenvalue problem for the two-dimensional Schrödinger equation containing an integral Hartree-type nonlinearity with an interaction potential having a logarithmic singularity. Global asymptotic solutions localized in the neighborhood of a line segment in the plane are constructed using the matching method for asymptotic expansions. The Bogoliubov and Airy polarons are used as model functions in these solutions. An analogue of the Bohr–Sommerfeld quantization rule is established to find the related series of eigenvalues.     

A New and Elegant Approach for Solving $$\varvec{n \times n}$$-Order Linear Fractional Differential Equations

The application of fractional differential equations (FDEs) in the fields of science and engineering are gradually increasing day by day during the last two decades. The solutions of linear systems of FDEs are of great importance. Several investigations are carried out on such systems using eigenvalue analysis or Laplace transform method. But both the methods have limitations, and as of now there are no methods for solving \(n \times n\)-order linear FDEs. In the present investigation, the issues of such difficulties are addressed, and the exact solutions of linear \(2 \times 2\)-order linear FDEs are presented by Laplace transform. We are unable to provide the exact solutions of such system of order \(n \times n\) by Laplace transform. To overcome this, we provide a new and elegant approach to find the approximate solutions of \(n \times n\)-order linear FDEs with the help of residual power series (RPS) method. The results thus obtained are verified by providing numerous examples.     

A commuting-vector-field approach to some dispersive estimates

We prove the pointwise decay of solutions to three linear equations: (1) the transport equation in phase space generalizing the classical Vlasov equation, (2) the linear Schrödinger equation, (3) the Airy (linear KdV) equation. The usual proofs use explicit representation formulae, and either obtain \(L^1\)—\(L^\infty \) decay through directly estimating the fundamental solution in physical space or by studying oscillatory integrals coming from the representation in Fourier space. Our proof instead combines “vector field” commutators that capture the inherent symmetries of the relevant equations with conservation laws for mass and energy to get space–time weighted energy estimates. Combined with a simple version of Sobolev’s inequality this gives pointwise decay as desired. In the case of the Vlasov and Schrödinger equations, we can recover sharp pointwise decay; in the Schrödinger case we also show how to obtain local energy decay as well as Strichartz-type estimates. For the Airy equation we obtain a local energy decay that is almost sharp from the scaling point of view, but nonetheless misses the classical estimates by a gap. This work is inspired by the work of Klainerman on \(L^2\)—\(L^\infty \) decay of wave equations, as well as the recent work of Fajman, Joudioux, and Smulevici on decay of mass distributions for the relativistic Vlasov equation.     

Nonlinear Transition Layers—The Second Painleve Transcendent

We investigate a large class of weakly nonlinear second-order ordinary differential equations with slowly varying coefficients. We show that the standard two-timing perturbation solution is not valid during the transition from oscillatory to exponentially decaying behavior. In all cases this difficulty is remedied by a nonlinear transition layer, whose leading-order character is described by one special nonlinear differential equation known as the second Painlevé transcendent (in essence a nonlinear Airy equation). The method of matched asymptotic expansions yields the desired connection formula. The second Painlevé transcendent also provides two other types of transitions: (1) between weakly nonlinear solutions (either oscillatory or exponentially decaying) and special fully nonlinear solutions, and (2) between two of these special nonlinear solutions. These special solutions are of three: different kinds: (a) slowly varying stable equilibrium solutions, (b) “exploding” solutions, and (c) solutions depending on both the fast and slow scales (which emerge from the unstable zero equilibrium solution).     

Global Asymptotics of the Second Painlevé Equation in Okamoto’s Space

We study the solutions of the second Painlevé equation (P II) in the space of initial conditions first constructed by Okamoto, in the limit as the independent variable, x, goes to infinity. Simultaneously, we study solutions of the related equation known as the thirty-fourth Painlevé equation (P 34). By considering degenerate cases of the autonomous flow, we recover the known special solutions, which are either rational functions or expressible in terms of Airy functions. We show that the solutions that do not vanish at infinity possess an infinite number of poles. An essential element of our construction is the proof that the union of exceptional lines is a repeller for the dynamics in Okamoto’s space. Moreover, we show that the limit set of the solutions exists and is compact and connected.     

Integrals containing products of Airy functions and planar photon propagation in a magnetic field

We derive a number of new results on integrals of products of Airy functions, using various integral transform techniques. As an application, we deduce a compact integral representation, suitable for numerical integration, of the one‐loop photon propagator in a magnetic field in 2 + 1 dimensional scalar quantum electrodynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Homogenization and dispersion effects in the problem of propagation of waves generated by a localized source

We construct asymptotic solutions to the wave equation with velocity rapidly oscillating against a smoothly varying background and with localized initial perturbations. First, using adiabatic approximation in the operator form, we perform homogenization that leads to a linearized Boussinesq-type equation with smooth coefficients and weak “anomalous” dispersion. Then, asymptotic solutions to this and, as a consequence, to the original equations are constructed by means of a modified Maslov canonical operator; for initial perturbations of special form, these solutions are expressed in terms of combinations of products of the Airy functions of a complex argument. On the basis of explicit formulas obtained, we analyze the effect of fast oscillations of the velocity on the solution fronts and solution profiles near the front.     

Computing Complex Airy Functions by Numerical Quadrature

Integral representations are considered of solutions of the Airy differential equation w –zw=0 for computing Airy functions for complex values of z. In a first method contour integral representations of the Airy functions are written as non-oscillating integrals for obtaining stable representations, which are evaluated by the trapezoidal rule. In a second method an integral representation is evaluated by using generalized Gauss–Laguerre quadrature; this approach provides a fast method for computing Airy functions to a predetermined accuracy. Comparisons are made with well-known algorithms of Amos, designed for computing Bessel functions of complex argument. Several discrepancies with Amos' code are detected, and it is pointed out for which regions of the complex plane Amos' code is less accurate than the quadrature algorithms. Hints are given in order to build reliable software for complex Airy functions.     

Exact and analytical solutions for a nonlinear sigma model

We consider wave solutions to nonlinear sigma models in n dimensions. First, we reduce the system of governing PDEs into a system of ODEs through a traveling wave assumption. Under a new transform, we then reduce this system into a single nonlinear ODE. Making use of the method of homotopy analysis, we are able to construct approximate analytical solutions to this nonlinear ODE. We apply two distinct auxiliary linear operators and show that one of these permits solutions with lower residual error than the other. This demonstrates the effectiveness of properly selecting the auxiliary linear operator when performing homotopy analysis of a nonlinear problem. From here, we then obtain residual error‐minimizing values of the convergence control parameter. We find that properly selecting the convergence control parameter makes a drastic difference in the magnitude of the residual error. Together, appropriate selection of the auxiliary linear operator and of the convergence control parameter is shown to allow approximate solutions that quickly converge to the true solution, which means that few terms are needed in the construction of such solution. This, in turn, greatly improves computational efficiency. Copyright © 2013 John Wiley & Sons, Ltd.     

Uniformization of WKB functions by Wigner transform

We develop a technique for uniformizing WKB functions which fail to correctly represent wave fields on caustics due to geometric singularities of ray fields. The uniformization technique is based on appropriate asymptotic surgery of the Wigner transform of the WKB functions, in different regions of the phase space. We present the details of the computation for the model example of the semiclassical Airy equation and we explain how the method can be applied to higher dimensional WKB functions for fold caustics     

A study on the convergence conditions of generalized differential transform method

This paper deals with constructing generalized ‘fractional’ power series representation for solutions of fractional order differential equations. We present a brief review of generalized Taylor's series and generalized differential transform methods. Then, we study the convergence of fractional power series. Our emphasis is to address the sufficient condition for convergence and to estimate the truncated error. Numerical simulations are performed to estimate maximum absolute truncated error when the generalized differential transform method is used to solve non‐linear differential equations of fractional order. The study highlights the power of the generalized differential transform method as a tool in obtaining fractional power series solutions for differential equations of fractional order. Copyright © 2016 John Wiley & Sons, Ltd.     

The numerical solutions of differential transform method and the Laplace transform method for a system of differential equations

The differential transform method is one of the approximate methods which can be easily applied to many linear and nonlinear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. In this paper, we present the definition and operation of the one-dimensional differential transform and investigate the particular exact solutions of system of ordinary differential equations that usually arise in mathematical biology by a one-dimensional differential transform method. The numerical results of the present method are presented and compared with the exact solutions that are calculated by the Laplace transform method.