Stability and convergence of staggered Runge-Kutta schemes for semilinear wave equations

A staggered Runge-Kutta (staggered RK) scheme is a Runge-Kutta type scheme using a time staggered grid, as proposed by Ghrist et al. in 2000 [6]. Afterwards, Verwer in two papers investigated the efficiency of a scheme proposed by Ghrist et al. [6] for linear wave equations. We study stability and convergence properties of this scheme for semilinear wave equations. In particular, we prove convergence of a fully discrete scheme obtained by applying the staggered RK scheme to the MOL approximation of the equation.     

Mathematical analysis of a wave energy converter model

In a context where sustainable development should be a priority, Orazov et al. have proposed in 2010, an excitation scheme for buoy-type ocean wave energy converters. The simplest model for this scheme is a non autonomous piecewise linear second order differential equation. The goal of the present paper is to give a mathematical framework for this model and to highlight some properties of its solutions. In particular, we will look at bounded and periodic solutions, and compare the energy-harvesting capabilities of this novel WEC with respect to that of a wave energy converter without mass modulation.     

p-Adic Analogue of the Wave Equation

In this paper, a p-adic analogue of the wave equation with Lipschitz source is considered. Since it is hard to arrive the solution of the problem, we propose a regularized method to solve the problem from a modified p-adic integral equation. Moreover, we give an iterative scheme for numerical computation of the regularlized solution.

     

A global extrapolated procedure for the boussinesq equation

In this paper we use the global extrapolation procedure to study the Boussinesq equation in one dimension. The application of a parametric finite-difference method leads to a three-time level nonlinear scheme, which is analyzed for local truncation error, stability and convergence. Then, the nonlinear term of the equation is properly linearized so, that the scheme becomes linear. The results of a number of numerical experiments for the single-soliton wave are given.     

Integrability,existence of global solutions,and wave breaking criteria for a generalization of the Camassa–Holm equation

Recent generalizations of the Camassa–Holm equation are studied from the point of view of existence of global solutions, criteria for wave breaking phenomena and integrability. We provide conditions, based on lower bounds for the first spatial derivative of local solutions, for global well-posedness in Sobolev spaces for the family under consideration. Moreover, we prove that wave breaking phenomena occurs under certain mild hypothesis. Based on the machinery developed by Dubrovin [Commun. Math. Phys. 267, 117–139 (2006)] regarding bi-Hamiltonian deformations, we introduce the notion of quasi-integrability and prove that there exists a unique bi-Hamiltonian structure for the equation only when it is reduced to the Dullin–Gotwald–Holm equation. Our results suggest that a recent shallow water model incorporating Coriollis effects is integrable only in specific situations. Finally, to finish the scheme of geometric integrability of the family of equations initiated in a previous work, we prove that the Dullin–Gotwald–Holm equation describes pseudo-spherical surfaces.     

Multisymplecticity of the centred box discretization for hamiltonian PDEs with m ≥ 2 space dimensions

In this letter, it is shown that the centred box discretization for Hamiltonian PDEs with m ≥ 2 space dimensions is multisymplectic in the sense of Bridges and Reich in [1–6]. Multisymplectic discretizations for the generalized KP equation and the wave equation with 2 space dimensions, respectively, are given. A multisymplectically numerical scheme of the wave equation is derived.     

Convergent numerical schemes for the compressible hyperelastic rod wave equation

We propose a fully discretised numerical scheme for the hyperelastic rod wave equation on the line. The convergence of the method is established. Moreover, the scheme can handle the blow-up of the derivative which naturally occurs for this equation. By using a time splitting integrator which preserves the invariants of the problem, we can also show that the scheme preserves the positivity of the energy density.     

A Conservative Space-time Mesh Refinement Method for the 1-D Wave Equation. Part II: Analysis

Summary. We introduced in [2] a new method for space-time refinement for the 1-D wave equation. This method is based on the conservation of a discrete energy through two different discretization grids which guarantees the stability of the scheme. In this second part, we analyse the accuracy of this scheme in a detailed way by means of a plane wave analysis and numerical experiments that permit us to point out spurious numerical phenomena and explain how to control them. Mathematics Subject Classification (2000):65N12     

Marching schemes for inverse acoustic scattering problems

Summary For the numerical solution of inverse Helmholtz problems the boundary value problem for a Helmholtz equation with spatially variable wave number has to be solved repeatedly. For large wave numbers this is a challenge. In the paper we reformulate the inverse problem as an initial value problem, and describe a marching scheme for the numerical computation that needs only n² log n operations on an n × n grid. We derive stability and error estimates for the marching scheme. We show that the marching solution is close to the low-pass filtered true solution. We present numerical examples that demonstrate the efficacy of the marching scheme.     

Solitary wave solutions of the modified equal width wave equation

In this paper we use a linearized numerical scheme based on finite difference method to obtain solitary wave solutions of the one-dimensional modified equal width (MEW) equation. Two test problems including the motion of a single solitary wave and the interaction of two solitary waves are solved to demonstrate the efficiency of the proposed numerical scheme. The obtained results show that the proposed scheme is an accurate and efficient numerical technique in the case of small space and time steps. A stability analysis of the scheme is also investigated.     

Unconditionally optimal convergence of an energy-conserving and linearly implicit scheme for nonlinear wave equations

In this paper, we present and analyze an energy-conserving and linearly implicit scheme for solving the nonlinear wave equations. Optimal error estimates in time and superconvergent error estimates in space are established without certain time-step restrictions. The key is to estimate directly the solution bounds in the H²-norm for both the nonlinear wave equation and the corresponding fully discrete scheme, while the previous investigations rely on the temporal-spatial error splitting approach. Numerical examples are presented to confirm energy-conserving properties, unconditional convergence and optimal error estimates, respectively, of the proposed fully discrete schemes.

     

Remarks on the asymptotic behavior of scalar auxiliary variable (SAV) schemes for gradient-like flows

We introduce a time semi-discretization of a damped wave equation by a SAV scheme with second order accuracy. The energy dissipation law is shown to hold without any restriction on the time step. We prove that any sequence generated by the scheme converges to a steady state (up to a subsequence). We notice that the steady state equation associated to the SAV scheme is a modified version of the steady state equation associated to the damped wave equation. We show that a similar result holds for a SAV fully discrete version of the Cahn-Hilliard equation and we compare numerically the two steady state equations.     

无界域上半线性强阻尼波动方程的全离散有理谱逼近

本文运用Chebyshev有理谱方法来讨论半线性强阻尼波动方程.通过建立时间、空间方向全离散的Chebyshev有理谱格式,证明了由此格式所确定的离散算子半群存在整体吸引子,并从理论上建立了在有限时间上近似解的误差估计.     

A finite difference scheme for traveling wave solutions to Burgers equation

We construct a finite difference scheme for the ordinary differential equation describing the traveling wave solutions to the Burgers equation. This difference equation has the property that its solution can be calculated. Our procedure for determining this solution follows closely the analysis used to obtain the traveling wave solutions to the original ordinary differential equation. The finite difference scheme follows directly from application of the nonstandard rules proposed by Mickens. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 815–820, 1998     

On the Numerical Solution to a Nonlinear Wave Equation Associated with the First Painlev′e Equation: an Operator-Splitting Approach

The main goal of this article is to discuss the numerical solution to a nonlinear wave equation associated with the first of the celebrated Painlevé transcendent ordinary differential equations. In order to solve numerically the above equation, whose solutions blow up in finite time, the authors advocate a numerical methodology based on the Strang’s symmetrized operator-splitting scheme. With this approach, one can decouple nonlinearity and differential operators, leading to the alternate solution at every time step of the equation as follows: (i) The first Painlevé ordinary differential equation, (ii) a linear wave equation with a constant coefficient. Assuming that the space dimension is two, the authors consider a fully discrete variant of the above scheme, where the space-time discretization of the linear wave equation sub-steps is achieved via a Galerkin/finite element space approximation combined with a second order accurate centered time discretization scheme. To handle the nonlinear sub-steps, a second order accurate centered explicit time discretization scheme with adaptively variable time step is used, in order to follow accurately the fast dynamic of the solution before it blows up. The results of numerical experiments are presented for different coefficients and boundary conditions. They show that the above methodology is robust and describes fairly accurately the evolution of a rather “violent” phenomenon.     

Instability of solitary waves for a generalized Benney–Luke equation

We study linear instability of solitary wave solutions of a one-dimensional generalized Benney–Luke equation, which is a formally valid approximation for describing two-way water wave propagation in the presence of surface tension. Further, we implement a finite difference numerical scheme which combines an explicit predictor and an implicit corrector step to compute solutions of the model equation which is used to validate the theory presented.     

Development of a numerical phase optimized upwinding combined compact difference scheme for solving the Camassa–Holm equation with different initial solitary waves

In this article, the solution of Camassa–Holm (CH) equation is solved by the proposed two‐step method. In the first step, the sixth‐order spatially accurate upwinding combined compact difference scheme with minimized phase error is developed in a stencil of four points to approximate the first‐order derivative term. For the purpose of retaining both of the long‐term accurate Hamiltonian property and the geometric structure inherited in the CH equation, the time integrator used in this study should be able to conserve symplecticity. In the second step, the Helmholtz equation governing the pressure‐like variable is approximated by the sixth‐order accurate three‐point centered compact difference scheme. Through the fundamental and numerical verification studies, the integrity of the proposed high‐order scheme is demonstrated. Another aim of this study is to reveal the wave propagation nature for the investigated shallow water equation subject to different initial wave profiles, whose peaks take the smooth, peakon, and cuspon forms. The transport phenomena for the cases with/without inclusion of the linear first‐order advection term κu_x in the CH equation will be addressed. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1645–1664, 2015     

Numerical solution of the two-dimensional Poincaré equation

This paper deals with numerical approximation of the two-dimensional Poincaré equation that arises as a model for internal wave motion in enclosed containers. Inspired by the hyperbolicity of the equation we propose a discretisation particularly suited for this problem, which results in matrices whose size varies linearly with the number of grid points along the coordinate axes. Exact solutions are obtained, defined on a perturbed boundary. Furthermore, the problem is seen to be ill-posed and there is need for a regularisation scheme, which we base on a minimal-energy approach.     

l 1-error estimates on the immersed interface upwind scheme for linear convection equations with piecewise constant coefficients: A simple proof

Abstract A linear convection equation with discontinuous coefcients arises in wave propagation through interfaces.An interface condition is needed at the interface to select a unique solution.An upwind scheme that builds this interface condition into its numerical flux is called the immersed interface upwind scheme.An l1-error estimate of such a scheme was frst established by Wen et al.(2008).In this paper,we provide a simple analysis on the l1-error estimate.The main idea is to formulate the solution to the underline initial-value problem into the sum of solutions to two convection equations with constant coefcients,which can then be estimated using classical methods for the initial or boundary value problems.     

Multisymplectic Structure and Multisymplectic Scheme for the Nonlinear Ware Equation

Abstract The multisvmplectic structure of the nonlinear wave equation is derived directly from the variationalprinciple. In the numerical aspect,we present a multisymplectic nine points scheme which is equivalent to themultisymplectic Preissman scheme.A series of numerical results are reported to illustrate the effectiveness ofthe scheme.