Solution of the scalar wave equation over very long distances using nonlinear solitary waves: Relation to finite difference methods

The linear wave equation represents the basis of many linear electromagnetic and acoustic propagation problems. Features that a computational model must have, to capture large scale realistic effects (for over the horizon or “OTH” radar communication, for example), include propagation of short waves with scattering and partial absorption by complex topography. For these reasons, it is not feasible to use Green’s Function or any simple integral method, which neglects these intermediate effects and requires a known propagation function between source and observer. In this paper, we describe a new method for propagating such short waves over long distances, including intersecting scattered waves. The new method appears to be much simpler than conventional high frequency schemes: Lagrangian “particle” based approaches, such as “ray tracing” become very complex in 3-D, especially for waves that may be expanding, or even intersecting. The other high frequency scheme in common use, the Eikonal, also has difficulty with intersecting waves.Our approach, based on nonlinear solitary waves concentrated about centroid surfaces of physical wave features, is related to that of Whitham [1], which involves solving wave fronts propagating on characteristics. Then, the evolving electromagnetic (or acoustic) field can be approximated as a collection of propagating co-dimension one surfaces (for example, 2-D surfaces in three dimensions). This approach involves solving propagation equations discretely on an Eulerian grid to approximate the linear wave equation. However, to propagate short waves over long distances, conventional Eulerian numerical methods, which attempt to resolve the structure of each wave, require far too many grid cells and are not feasible on current or foreseeable computers. Instead, we employ an “extended” wave equation that captures the important features of the propagating waves. This method is first formulated at the partial differential equation (PDE) level, as a wave equation with an added “confining” term that involves both a positive and a negative dissipation. Once we have the stable PDE, the discrete formulation is simply a multidimensional PDE with (stable) perturbations caused by the discretization. The resulting discrete solution can then be low order and very simple and yet remain stable over arbitrarily long times. When discretized and solved on an Eulerian grid, this new method allows far coarser grids than required by conventional resolution considerations, while still accounting for the effects of varying atmospheric and topographic features. An important point is that the new method is in the same form as conventional discrete wave equation methods. However, the conventional solution eventually decays, and only the “intermediate asymptotic” solution can be used. Simply by adding an extra term, we show that a nontrivial true asymptotic solution can be obtained. A similar solitary wave based approach has been used successfully in a different problem (involving “Vorticity Confinement”), for a number of years.     

Acoustic-Gravity Waves in the Atmosphere with a Piecewise-Linear Temperature Profile

We consider the problem of propagation of acoustic-gravity waves in the atmosphere with a constant temperature gradient in the near-surface layer. The assumption of linear temperature dependence on height allowed us to reduce the wave equation to the hypergeometric form, regardless of the compressibility of the medium. The solution of this equation is represented in terms of degenerate hypergeometric functions. To analyze the obtained solution, we consider a two-layer model of a half-bounded atmosphere with a height-independent background temperature in the upper layer. The results are studied in detail under the approximation of an incompressible medium. For the model specified above, we find analytical expressions for the perturbation fields and obtain a characteristic equation whose solution allows us to calculate wave dispersion characteristics at frequencies close to the Brunt-Väisälä frequency for large horizontal scales as compared to the layer thickness.     

A Method for Solving Initial-value Problem and Finding Exact Solutions of Nonlinear Partial Differential Equations

Homogeneous balance method for solving nonlinear partial differential equation(s) is extended to solving initial-value problem and getting new solution(s) from a known solution of the equation(s) under consideration. The approximate equations for long water waves are chosen to illustrate the method, infinitely many simple-solitary-wave solutions and infinitely many rational function solutions, especially the closed form of the solution for initial-value problem, are obtained by using the extended homogeneous balance method given here.     

Dynamical systems in the theory of solitons in the presence of nonlocal interactions

The behavior of solitons in models which take into account complex dispersion or nonlocal interaction of nonlinear waves is examined. A method is proposed to reduce this problem to one involving special trajectories (homoclinic and heteroclinic) of the dynamic system. This method involves replacing the nonlinear integrodifferential equation with the differential equations which link the original nonlinear field with the auxiliary linear fields. The interaction of fields in such a model is a local interaction. The number of introduced linear fields is determined by the Laplace transform of the integral operator kernel of the basic integrodifferential equation. The problem involving topological solitons for the nonlocal generalization of the Klein-Gordon equation is considered. Nonlocal interactions are found to lead to a number of singularities (unrestricted increase in the slope of the topological soliton front, break in the solutions, and other singularities).     

Solutions of scalar and electromagnetic wave equations in the metric of gravitational and electromagnetic waves

The wave equation for a scalar field ? and vector potentialA* are solved in the background metric of a gravitational wave. The corresponding solutions when the metric is generated by a plane electromagnetic wave, is obtained from these solutions. The solution for the scalar wave is discussed in detail. It is found that because of the interaction, two new waves are generated in the lower order approximations. One of them has the same phase dependence as the original wave while the other shows a transient character. There is no interaction when the waves are along the same direction.     

Propagation of thickness-twist waves in elastic plates with periodically varying thickness and phononic crystals

We study the propagation of thickness-twist (TT) waves in a crystal plate of AT-cut quartz with periodically varying, piecewise constant thickness. The scalar differential equation by Tiersten and Smythe is employed. The problem is found to be mathematically equivalent to the motion of an electron in a periodic potential field governed by Schrodinger’s equation. An analytical solution is obtained. Numerical results show that the eigenvalue (frequency) spectrum of the waves has a band structure with allowed and forbidden bands. Therefore, for TT waves, plates with periodically varying thickness can be considered as phononic crystals. The effects of various parameters on the frequency spectrum are examined.     

Fractional diffusion-wave problem in cylindrical coordinates

In this Letter, we present analytical and numerical solutions for an axis-symmetric diffusion-wave equation. For problem formulation, the fractional time derivative is described in the sense of Riemann-Liouville. The analytical solution of the problem is determined by using the method of separation of variables. Eigenfunctions whose linear combination constitute the closed form of the solution are obtained. For numerical computation, the fractional derivative is approximated using the Grünwald-Letnikov scheme. Simulation results are given for different values of order of fractional derivative. We indicate the effectiveness of numerical scheme by comparing the numerical and the analytical results for α=1 which represents the order of derivative.     

Extreme waves statistics for the Ablowitz-Ladik system

We examine statistics of waves for the problem of modulation instability development in the framework of discrete integrable Ablowitz-Ladik (AL) system. Modulation instability depends on one free parameter h that has the meaning of the coupling between the nodes on the lattice. For strong coupling h ? 1, the probability density functions (PDFs) for waves amplitudes coincide with that for the continuous classical nonlinear Schrödinger equation; the PDFs for both systems are very close to Rayleigh ones. When the coupling is weak h ～ 1, there appear highly localized waves with very large amplitudes, that drastically change the PDFs to significantly non-Rayleigh ones, with so-called “fat tails” when the probability of a large wave occurrence is by several orders of magnitude higher than that predicted by the linear theory. Evolution of amplitudes for such rogue waves with time is similar to that of the Peregrine solution for the classical nonlinear Schrödinger equation.     

Efficient simulation of non-linear effects in 2D optical nanostructures to TM waves

We develop a theory to study stationary TM-type waves propagating in a nanostructured layer of 2D non-linear optical metamaterial or plasmonic device. It is assumed that the layer is inhomogeneous and contains non-linear isotropic elemental materials with non-linearity and loss mechanisms, including both linear and non-linear losses. While modeling of the non-linear propagation of the TE-type scalar waves is straightforward, the TM-type waves within the standard E-field formulations of non-linear optics cannot be treated in a purely scalar H-field context since an implicit equation for the non-linear dielectric functions should be resolved otherwise. A new formulation, which is built on the solution of the implicit equation for the non-linear dielectric function, is proposed. We use a general cubic non-linearity to illustrate all of the important features of the proposed approach. The general solution for scalar H-field waves is validated versus our previously tested particular cases, and important differences are shown between those cases and the general solution. These details, for example, include the link between linear and non-linear loss mechanisms, and connection between the linear and non-linear dielectric functions. The proposed approach is used for modeling a non-linear focusing device with optically controlled isotropic Kerr-type non-linearity; the simulation results prove the predicted functioning of the device.     

Analytic solutions of the optical bistability equations for a standing wave cavity

The problem of optical bistability in a standing wave cavity in the steady state leads to a pair of coupled, nonlinear, ordinary differential equations for the forward and backward waves. Only numerical solutions have so far been presented for these equations. We give their exact analytic solutions and find good agreement with the numerical results. The exact solutions are shown to reduce to the mean field equation for the input and output fields in the double limits T → 0 and αL → 0 for the mirror transmission and the linear absorption absorption, respectively.     

Electromagnetic waves propagation in real multimode waveguide structures

A power flow equation describing the power flow of electromagnetic waves in a real multimode waveguide represented by a linear multichannel trasmission system exhibiting an attenuation, dispersion, and an interchannel interaction is solved in three iterative ways. In the case of constant attenuation and constant propagation velocities in all channels and for the convolution interchannel interaction, a closed analytical solution of the problem is presented. Interesting forms of the solution, transfer function, and impulse response of the system which enable us to separate and compare the attenuation, dispersion, and coupling effects independently of each other are derived. The transfer function, the impulse response, and excitation conditions for their measurements are further discussed. Relations which make it possible to compare different experimental data obtained under specific excitation conditions in various laboratories are determined. Finally, a linear system with memory in time is studied and it is shown that the real multimode waveguide as a linear multichannel transmission system can be considered as a general linear system with memory in space.     

Folded novel accurate analytical and semi-analytical solutions of a generalized Calogero–Bogoyavlenskii–Schiff equation

This paper studies the analytical and semi-analytic solutions of the generalized Calogero–Bogoyavlenskii–Schiff(CBS) equation. This model describes the(2 + 1)–dimensional interaction between Riemann-wave propagation along the y-axis and the x-axis wave. The extended simplest equation(ESE) method is applied to the model, and a variety of novel solitarywave solutions is given. These solitary-wave solutions prove the dynamic behavior of soliton waves in plasma. The accuracy of the obtained solution is verified using a variational iteration(VI) semi-analytical scheme. The analysis and the match between the constructed analytical solution and the semi-analytical solution are sketched using various diagrams to show the accuracy of the solution we obtained. The adopted scheme's performance shows the effectiveness of the method and its ability to be applied to various nonlinear evolution equations.     

场方法的改进及其在积分Riemann-Cartan空间运动方程中的应用

和Hamilton-Jacobi方法类似,Vujanovi?场方法把求解常微分方程组特解的问题转化为寻找一个一阶拟线性偏微分方程(基本偏微分方程)完全解的问题,但Vujanovi?场方法依赖于求出基本偏微分方程的完全解,而这通常是困难的,这就极大地限制了场方法的应用.本文将求解常微分方程组特解的Vujanovi?场方法改进为寻找动力学系统运动方程第一积分的场方法,并将这种方法应用于一阶线性非完整约束系统Riemann-Cartan位形空间运动方程的积分问题中.改进后的场方法指出,只要找到基本偏微分方程的包含m(m≤ n,n为基本偏微分方程中自变量的数目)个任意常数的解,就可以由此找到系统m个第一积分.特殊情况下,如果能够求出基本偏微分方程的完全解(完全解是m=n时的特例),那么就可以由此找到≤系统全部第一积分,从而完全确定系统的运动.Vujanovi?场方法等价于这种特殊情况.     

On MHD transient flow of a Maxwell fluid in a porous medium and rotating frame

Analytic solution for unsteady magnetohydrodynamic (MHD) flow is constructed in a rotating non-Newtonian fluid through a porous medium. Constitutive equations for a Maxwell fluid have been taken into consideration. The hydromagnetic flow in the uniformly rotating fluid is generated by a suddenly moved infinite plate in its own plane. Analytic solution of the governing flow problem is obtained by means of the Fourier sine transform. It is shown that the obtained solution satisfies both the associate partial differential equation and the initial and boundary conditions. The solution for a Navier-Stokes fluid is recovered if λ→0. The steady state solution is also obtained for t→∞.     

Numerical solution of fractional differential equations via a Volterra integral equation approach

The main focus of this paper is to present a numerical method for the solution of fractional differential equations. In this method, the properties of the Caputo derivative are used to reduce the given fractional differential equation into a Volterra integral equation. The entire domain is divided into several small domains, and by collocating the integral equation at two adjacent points a system of two algebraic equations in two unknowns is obtained. The method is applied to solve linear and nonlinear fractional differential equations. Also the error analysis is presented. Some examples are given and the numerical simulations are also provided to illustrate the effectiveness of the new method.     

Functional quantum statistics of light propagation in a two-level system

The functional Fokker-Planck formalism developed in a preceding paper is applied to the problem of a radiation field propagating in a medium, which contains resonant two-level atoms. Besides the electromagnetic field also the medium is described by continuous space dependent fields. We give the masterequation and transform it into ac-number functional differential equation for a characteristic functional. This equation is reduced considerably by the projection onto one dimension and the introduction of the diffusion approximation. It forms a solid basis for the study of all types of light propagation in resonant media including classical and quantum noise. We give an approximate solution of this equation by considering the problem of an externally pumped optical transmission line, in the case that saturation effects are absent. The spectral function of the electric field strength is obtained which describes a statistical mixture of photons with the quasiparticles of the polarization field. It shows the onset of a condensation of the quasiparticles into a single state. Self excitation of the transmission line is obtained at a certain threshold of the atomic inversion. This threshold is characterized by a finite occupation number of one single quasiparticle state. The influence of a finite length of the transmission line is briefly considered.     

Eigenstate expansions in a soluble model of identical particle scattering

Solutions to the Schrödinger equation and the inhomogeneous equation for the case of two identical particles interacting with a center of force are studied. Eigenstate expansions for solving each equation are explicitly introduced and their properties discussed. The case when the interparticle interaction v₁₂ is zero is then examined; this is a completely soluble problem. The eigenstate expansion solutions for the Schrödinger and inhomogeneous equations are used to explore the means by which the correct solution is obtained. Finally, approximate solutions, obtained by truncating the eigenfunction expansions, are introduced. It is seen that both methods lead to the correct amplitude when τ₁₂ = 0, even though the approximate solution to the inhomogeneous equation does not lead, in the end, to an antisymmetric solution.     

Solutions to Class of Linear and Nonlinear Fractional Differential Equations

In this paper, the fractional auxiliary sub-equation expansion method is proposed to solve nonlinear fractional differential equations. To illustrate the effectiveness of the method, we discuss the space-time fractional Kd V equation, the space-time fractional RLW equation, the space-time fractional Boussinesq equation, and the(3+1)-spacetime fractional ZK equation. The solutions are expressed in terms of fractional hyperbolic and fractional trigonometric functions. These solutions are useful to understand the mechanisms of the complicated nonlinear physical phenomena and fractional differential equations. Among these solutions, some are found for the first time. The analytical solution of homogenous linear FDEs with constant coefficients are obtained by using the series and the Mittag–Leffler function methods. The obtained results recover the well-know solutions when α = 1.     

充水粘弹性管道的频散曲线计算分析*

针对谱方法分析计算充水粘弹性管道的广义特征值问题,根据Chebyshev多项式及微分矩阵、位移和应力连续条件,将波动方程离散为相应的线性方程。利用MATLAB数值编程计算充水弹性和粘弹性管道对应频率下的轴对称纵向导波频散曲线和衰减曲线。分析表明,波传播在粘弹性管道中不仅具有衰减特性,而且由于水和粘弹性壳体交叉耦合作用,在一定频率范围内产生两种截断模态。     

Solitary wave dynamics in shallow water over periodic topography

The problem of long-wave scattering by piecewise-constant periodic topography is studied both for a linear solitary-like wave pulse, and for a weakly nonlinear solitary wave [Korteweg-de Vries (KdV) soliton]. If the characteristic length of the topographic irregularities is larger than the pulse length, the solution of the scattering problem is obtained analytically for a leading wave in the framework of linear shallow-water theory. The wave decrement in the case of the small height of the topographic irregularities is proportional to delta2, where delta is the relative height of the topographic obstacles. An analytical approximate solution is also obtained for the weakly nonlinear problem when the length of the irregularities is larger than the characteristic nonlinear length scale. In this case, the Korteweg-de Vries equation is solved for each piece of constant depth by using the inverse scattering technique; the solutions are matched at each step by using linear shallow-water theory. The weakly nonlinear solitary wave decays more significantly than the linear solitary pulse. Solitary wave dynamics above a random seabed is also discussed, and the results obtained for random topography (including experimental data) are in reasonable agreement with the calculations for piecewise topography.