Symplectic wavelet collocation method for Hamiltonian wave equations

This paper introduces a novel symplectic wavelet collocation method for solving nonlinear Hamiltonian wave equations. Based on the autocorrelation functions of Daubechies compactly supported scaling functions, collocation method is conducted for the spatial discretization, which leads to a finite-dimensional Hamiltonian system. Then, appropriate symplectic scheme is employed for the integration of the Hamiltonian system. Under the hypothesis of periodicity, the properties of the resulted space differentiation matrix are analyzed in detail. Conservation of energy and momentum is also investigated. Various numerical experiments show the effectiveness of the proposed method.     

A spectral collocation algorithm for two-point boundary value problem in fiber Raman amplifier equations

A novel algorithm implementing Chebyshev spectral collocation (pseudospectral) method in combination with Newton’s method is proposed for the nonlinear two-point boundary value problem (BVP) arising in solving propagation equations in fiber Raman amplifier. Moreover, an algorithm to train the known linear solution for use as a starting solution for the Newton iteration is proposed and successfully implemented. The exponential accuracy obtained by the proposed Chebyshev pseudospectral method is demonstrated on a case of the Raman propagation equations with strong nonlinearities. This is in contrast to algebraic accuracy obtained by typical solvers used in the literature. The resolving power and the efficiency of the underlying Chebyshev grid are demonstrated in comparison to a known BVP solver.     

非线性波动方程的Jacobi椭圆函数包络周期解

应用Jacobi椭圆函数展开法求得了一类非线性波方程的包络周期解,而且用这种方法得到的周期解在一定条件下可以退化为包络冲击波解或包络孤立波解 关键词： Jacobi椭圆函数 非线性方程 包络周期解 包络孤立波解     

A surface wave dispersion relation for non-local media

Conditions are obtained for the existence of surface waves at the interface of vacuum and a semi-infinite non-local medium whose inverse dielectric function is assumed to be symmetric in spatial coordinate in the direction perpendicular to the interface. It is shown that these conditions reduce to those obtained by Maradudin in the appropriate limits: for the isotropic case; and for no retardation.     

长短波相互作用方程的Jacobi椭圆函数求解

推广了Jacobi椭圆函数展开方法,研究了复非线性演化方程组的求解问题,得到了长短波相互作用方程的准确包络周期解.该结果在一定条件下包含了相应的孤波解. 关键词： Jacobi椭圆函数方法 长短波相互作用方程 孤波解     

Quasiclassical trajectory-coherent states for equations of wave type

Asymptotic solutions are constructed for equations of second order with respect to the time which contain a small parameter. The solutions are constructed on the basis of the Maslov complex-growth method.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 5, pp. 43–48, May, 1989.     

An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations

In recent years, there has been a growing interest in analyzing and quantifying the effects of random inputs in the solution of ordinary/partial differential equations. To this end, the spectral stochastic finite element method (SSFEM) is the most popular method due to its fast convergence rate. Recently, the stochastic sparse grid collocation method has emerged as an attractive alternative to SSFEM. It approximates the solution in the stochastic space using Lagrange polynomial interpolation. The collocation method requires only repetitive calls to an existing deterministic solver, similar to the Monte Carlo method. However, both the SSFEM and current sparse grid collocation methods utilize global polynomials in the stochastic space. Thus when there are steep gradients or finite discontinuities in the stochastic space, these methods converge very slowly or even fail to converge. In this paper, we develop an adaptive sparse grid collocation strategy using piecewise multi-linear hierarchical basis functions. Hierarchical surplus is used as an error indicator to automatically detect the discontinuity region in the stochastic space and adaptively refine the collocation points in this region. Numerical examples, especially for problems related to long-term integration and stochastic discontinuity, are presented. Comparisons with Monte Carlo and multi-element based random domain decomposition methods are also given to show the efficiency and accuracy of the proposed method.     

A wave automaton for Maxwell's equations

This paper presents an extension to electromagnetic fields of the wave automaton, which was introduced in recent years for describing wave propagation in inhomogeneous media. Using elementary processes obeying a discrete Huygens' principle and satisfying fundamental symmetries such as time reversal and reciprocity, this new wave automaton is capable of modeling Maxwell's equations in 3+1 dimensions. It supplements the methods that were developed early for scalar and spinor fields. Received 19 July 2001     

A new construction for spinor wave equations

The construction of discrete scalar wave propagation equations in arbitrary inhomogeneous media was recently achieved by using elementary dynamical processes realizing a discrete counterpart of the Huygens principle. In this paper, we generalize this approach to spinor wave propagation. Although the construction can be formulated on a discrete lattice of any dimension, for simplicity we focus on spinors living in 1+1 space-time dimensions. The Dirac equation in the Majorana-Weyl representation is directly recovered by incorporating appropriate symmetries of the elementary processes. The Dirac equation in the standard representation is also obtained by using its relationship with the Majorana-Weyl representation. Received: 3 November 1997 / Received in final form: 9 February 1998 / Accepted: 16 February 1998     

Driven nonequilibrium lattice systems with shifted periodic boundary conditions

We present the first study of a driven nonequilibrium lattice system in the two-phase region, withshifted periodic boundary conditions, forcing steps into the interface. When the shift corresponds to small angles with respect to the driving field, we find nonanalytic behavior in the (internal) energy of the system, supporting numerical evidence that interface roughness is suppressed by the field. For larger shifts, the competition between the driving field and the boundary induces the breakup of a single strip with tilted interfaces into many narrower strips with aligned interfaces. The size and temperature dependences of the critical angles of such breakup transitions are studied.     

Well-posedness and the energy and charge conservation for nonlinear wave equations in discrete space-time

We consider the discretization problem for U(1)-invariant nonlinear wave equations in any dimension. We show that the classical finite-difference scheme used by Strauss and Vazquez (in J. Comput. Phys. 28, 271–278 (1978)) conserves the positive-definite discrete analog of the energy if the grid ratio satisfies \(dt/dx \leqslant 1/\sqrt n \), where dt and dx are the mesh sizes of the time and space variables and n is the spatial dimension. We also show that, if the grid ratio is \(dt/dx \leqslant 1/\sqrt n \), then there is a discrete analog of charge, and this discrete analog is conserved.We prove the existence and uniqueness of solutions to the discrete Cauchy problem. We use energy conservation to obtain a priori bounds for finite energy solutions, thus showing that the Strauss-Vazquez finite-difference scheme for the nonlinear Klein-Gordon equation with positive nonlinear term in the Hamiltonian is conditionally stable.     

Thermal momentum distribution from path integrals with shifted boundary conditions

For a thermal field theory formulated in the grand canonical ensemble, the distribution of the total momentum is an observable characterizing the thermal state. We show that its cumulants are related to thermodynamic potentials. In a relativistic system, for instance, the thermal variance of the total momentum is a direct measure of the enthalpy. We relate the generating function of the cumulants to the ratio of (a) a partition function expressed as a Matsubara path integral with shifted boundary conditions in the compact direction and (b) the ordinary partition function. In this form the generating function is well suited for Monte Carlo evaluation, and the cumulants can be extracted straightforwardly. We test the method in the SU(3) Yang-Mills theory and obtain the entropy density at three different temperatures.     

求sine-Gordon 型方程孤波解的一种统一方法

借助于一个辅助常微分方程的解，提出了一种求sine-Gordon 型方程孤波解的统一方法，并用该法简洁地求得了三个著名的sine-Gordon 型方程，即单sine-Gordon 方程、双sine-Gordon 方程和mKdV-sine-Gordon方程的精确孤波解。     

A staggered-grid finite-difference method with perfectly matched layers for poroelastic wave equations.

A particle velocity-strain, finite-difference (FD) method with a perfectly matched layer (PML) absorbing boundary condition is developed for the simulation of elastic wave propagation in multidimensional heterogeneous poroelastic media. Instead of the widely used second-order differential equations, a first-order hyperbolic leap-frog system is obtained from Biot's equations. To achieve a high accuracy, the first-order hyperbolic system is discretized on a staggered grid both in time and space. The perfectly matched layer is used at the computational edge to absorb the outgoing waves. The performance of the PML is investigated by calculating the reflection from the boundary. The numerical method is validated by analytical solutions. This FD algorithm is used to study the interaction of elastic waves with a buried land mine. Three cases are simulated for a mine-like object buried in "sand," in purely dry "sand" and in "mud." The results show that the wave responses are significantly different in these cases. The target can be detected by using acoustic measurements after processing.     

Bäcklund transformation,local and non-local conservation laws for super-chiral fields

New non-local conservation laws, parametric Bäcklund transformation and local conservation laws are constructed for super-chiral fields in general, using similar methods for ordinary chiral fields. We thus have a unified view of these field theories.     

A dynamical monte carlo algorithm for master equations with time-dependent transition rates

A Monte Carlo algorithm for simulating master equations with time-dependent transition rates is described. It is based on a waiting time image, and takes into account that the system can become frozen when the transition rates tend to zero fast enough in time. An analytical justification is provided. The algorithm reduces to the Bortz-Kalos-Lebowitz one when the transition rates are constant. Since the exact evaluation of waiting times is rather involved in general, a simple and efficient iterative method for accurately calculating them is introduced. As an example, the algorithm is applied to a one-dimensional Ising system with Glauber dynamics. It is shown that it reproduces the exact analytical results, being more efficient than the direct implementation of the Metropolis algorithm     

计算相位响应曲线的方波扰动直接算法

通过相位响应曲线可对具有极限环周期运动的动力系统的性质有更为深入的理解.神经元是一个典型的动力系统,因此相位响应曲线提供了一种研究神经元重复周期放电行为的新思路.本文提出一种求解相位响应曲线的方法,即方波扰动的直接算法,通过Hodgkin-Huxley,Fitz Hugh-Nagumo,Morris-Lecar和Hindmarsh-Rose神经元模型验证该算法可计算周期峰放电、周期簇放电的相位响应曲线.该算法克服了其他算法在运用过程中的局限性.利用该算法计算结果表明:周期峰放电的相位响应曲线类型是由其分岔类型所决定;在Morris-Lecar模型中发现一种开始于Hopf分岔终止于鞍点同宿轨道分岔的阈上周期振荡,其相位响应曲线属于第二类型.通过大量的相位响应曲线的计算发现相位响应的相对大小及正负性仅取决于扰动所施加的时间,而且周期簇放电的相位响应曲线比周期峰放电的相位响应曲线更为复杂.     

A class of nondispersive evolution equations with solitary wave solutions

We present a class of nonlinear evolution equations possessing stable solitary wave solutions with a sech² profile. These equations are related to the Korteweg-de Vries (KdV) and regularised longwave (RLW) equations, but, unlike the latter, are dispersion free in the linear limit.