An energy-conserving second order numerical scheme for nonlinear hyperbolic equation with an exponential nonlinear term

摘要

We present a second order accurate numerical scheme for a nonlinear hyperbolic equation with an exponential nonlinear term. The solution to such an equation is proven to have a conservative nonlinear energy. Due to the special nature of the nonlinear term, the positivity is proven to be preserved under a periodic boundary condition for the solution. For the numerical scheme, a highly nonlinear fractional term is involved, for the theoretical justification of the energy stability. We propose a linear iteration algorithm to solve this nonlinear numerical scheme. A theoretical analysis shows a contraction mapping property of such a linear iteration under a trivial constraint for the time step. We also provide a detailed convergence analysis for the second order scheme, in the l(infinity) (0, T; l(infinity)) norm. Such an error estimate in the maximum norm can be obtained by performing a higher order consistency analysis using asymptotic expansions for the numerical solution. As a result, instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an 0(Delta t(3) + h(4)) convergence in l(infinity) (0, T; l(2)) norm, which leads to the necessary l(infinity) error estimate using the inverse inequality, under a standard constraint Delta t <= Ch. A numerical accuracy check is given and some numerical simulation results are also presented. (C) 2014 Elsevier B.V. All rights reserved.