An energy‐preserving Crank–Nicolson Galerkin method for Hamiltonian partial differential equations

A semidiscretization based method for solving Hamiltonian partial differential equations is proposed in this article. Our key idea consists of two approaches. First, the underlying equation is discretized in space via a selected finite element method and the Hamiltonian PDE can thus be casted to Hamiltonian ODEs based on the weak formulation of the system. Second, the resulting ordinary differential system is solved by an energy‐preserving integrator. The relay leads to a fully discretized and energy‐preserved scheme. This strategy is fully realized for solving a nonlinear Schrödinger equation through a combination of the Galerkin discretization in space and a Crank–Nicolson scheme in time. The order of convergence of our new method is if the discrete L²‐norm is employed. An error estimate is acquired and analyzed without grid ratio restrictions. Numerical examples are given to further illustrate the conservation and convergence of the energy‐preserving scheme constructed.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1485–1504, 2016     

A Crank–Nicolson collocation spectral method for the two‐dimensional viscoelastic wave equation

In this paper, we first establish the Crank–Nicolson collocation spectral (CNCS) method for two‐dimensional (2D) viscoelastic wave equation by means of the Chebyshev polynomials. And then, we analyze the existence, uniqueness, stability, and convergence of the CNCS solutions. Finally, we use some numerical experiments to verify the correctness of theoretical analysis. This implies that the CNCS model is very effective for solving the 2D viscoelastic wave equations.     

Numerical analysis of the Crank–Nicolson extrapolation time discrete scheme for magnetohydrodynamics flows

In this article, we consider the time‐discrete method for three‐dimensional incompressible magnetohydrodynamics (MHD) equations. The Crank–Nicolson extrapolation scheme is used for time discretization. From the previous articles, under the assumption that the solution has high regularity which cannot be realistically assumed, the convergence of this scheme is optimal two‐order. In this article, under modest assumptions of initial values and the body force, we prove some new regularity results of the MHD equations. In addition, we derive the unconditional convergence of our scheme, but the convergent order is not optimal. Furthermore, we provide another conditional convergence estimation to increase the order. It is shown that the convergent rate increase half order in ‐norm, and at least a quarter order increased in ‐norm than the uncondtional results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2169–2208, 2015     

Numerical study of Sivashinsky equation using a splitting scheme based on Crank‐Nicolson method

The mathematical modeling of a planar solid‐liquid interface in the solidification of a dilute binary alloy is formulating by one of nonintegrable, nonlinear evolution equation known as Sivashinsky equation. In the first part of this paper, the mathematical modeling of Sivashinsky equation is briefly discussed. Since, the exact solutions of this equation is yet unknown, obtaining its numerical solution plays an important role to simulate its behavior. Therefore, in the second part, a second‐order splitting finite difference scheme, based on Crank‐Nicolson method, is investigated to approximate the solution of the Sivashinsky equation with homogeneous boundary conditions. We prove the solvability of the present scheme and establish the error estimate of the numerical scheme.     

On the long‐time stability of the Crank–Nicolson scheme for the 2D Navier–Stokes equations

In this article we study the stability for all positive time of the Crank–Nicolson scheme for the two‐dimensional Navier–Stokes equations. More precisely, we consider the Crank–Nicolson time discretization together with a general spatial discretization, and with the aid of the discrete Gronwall lemma and of the discrete uniform Gronwall lemma we prove that the numerical scheme is stable, provided a CFL‐type condition is satisfied. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

An ETD Crank‐Nicolson method for reaction‐diffusion systems

A novel Exponential Time Differencing Crank‐Nicolson method is developed which is stable, second‐order convergent, and highly efficient. We prove stability and convergence for semilinear parabolic problems with smooth data. In the nonsmooth data case, we employ a positivity‐preserving initial damping scheme to recover the full rate of convergence. Numerical experiments are presented for a wide variety of examples, including chemotaxis and exotic options with transaction cost. © 2011Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

A modified Chebyshev ϑ‐weighted Crank–Nicolson method for analyzing fractional sub‐diffusion equations

This study presents a robust modification of Chebyshev ? ‐weighted Crank–Nicolson method for analyzing the sub‐diffusion equations in the Caputo fractional sense. In order to solve the problem, by discretization of the sub‐fractional diffusion equations using Taylor's expansion a linear system of algebraic equations that can be analyzed by numerical methods is presented. Furthermore, consistency, convergence, and stability analysis of the suggested method are discussed. In this framework, compact structures of sub‐diffusion equations are considered as prototype examples. The main advantage of the proposed method is that, it is more efficient in terms of CPU time, computational cost and accuracy in comparing with the existing ones in open literature.     

Partitioned penalty methods for the transport equation in the evolutionary Stokes–Darcy‐transport problem

There has been a surge of work on models for coupling surface‐water with groundwater flows which is at its core the Stokes–Darcy problem, as well as methods for uncoupling the problem into subdomain, subphysics solves. The resulting (Stokes–Darcy) fluid velocity is important because the flow transports contaminants. The numerical analysis and algorithm development for the evolutionary transport problem has, however, focused on a quasi‐static Stokes–Darcy model and a single domain (fully coupled) formulation of the transport equation. This report presents a numerical analysis of a partitioned method for contaminant transport for the fully evolutionary system. The algorithm studied is unconditionally stable with one subdomain solve per step. Numerical experiments are given using the proposed algorithm that investigates the effects of the penalty parameters on the convergence of the approximations.     

Convergence analysis of a linearized Crank–Nicolson scheme for the two‐dimensional complex Ginzburg–Landau equation

A linearized Crank–Nicolson‐type scheme is proposed for the two‐dimensional complex Ginzburg–Landau equation. The scheme is proved to be unconditionally convergent in the L² ‐norm by the discrete energy method. The convergence order is \begin{align*}\mathcal{O}(\tau^2+h_1^2+h^2_2)\end{align*}, where τ is the temporal grid size and h₁,h₂ are spatial grid sizes in the x ‐ and y ‐directions, respectively. A numerical example is presented to support the theoretical result. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

An adaptive algorithm for the Crank–Nicolson scheme applied to a time-dependent convection–diffusion problem

An a posteriori upper bound is derived for the nonstationary convection–diffusion problem using the Crank–Nicolson scheme and continuous, piecewise linear stabilized finite elements with large aspect ratio. Following Lozinski et al. (2009) [13], a quadratic time reconstruction is used.A space and time adaptive algorithm is developed to ensure the control of the relative error in the L²(H¹) norm. Numerical experiments illustrating the efficiency of this approach are reported; it is shown that the error indicator is of optimal order with respect to both the mesh size and the time step, even in the convection dominated regime and in the presence of boundary layers.     

Stability of the Crank–Nicolson–Adams–Bashforth scheme for the 2D Leray‐alpha model

We consider the stability of an efficient Crank–Nicolson–Adams–Bashforth method in time, finite element in space, discretization of the Leray‐α model. We prove finite‐time stability of the scheme in L², H¹, and H², as well as the long‐time L‐stability of the scheme under a Courant‐Freidrichs‐Lewy (CFL)‐type condition. Numerical experiments are given that are in agreement with the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1155–1183, 2016     

A decoupling method with different subdomain time steps for the nonstationary stokes–darcy model

This report analyzes a multirate, decoupling algorithm, which allows different time steps in the fluid region and the porous region for the nonstationary Stokes–Darcy problem. The method presented requires only one, uncoupled Stokes and Darcy subphysics and subdomain solve per time step. Under a time step restriction of the form △t ≤ C (physical parameters) we prove stability and convergence of the method. Numerical tests are given to show the convergence result and demonstrate the computational efficiency of the partitioned method. They also show that in the expected case of greater fluid velocities in the free‐flow region than in the porous media region, allowing smaller time steps in the subregion with the faster velocities increases both accuracy and efficiency. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Blow‐up of smooth solutions to the Navier–Stokes–Poisson equations

The main purpose of this paper is concerned with blow‐up smooth solutions to Navier–Stokes–Poisson (N‐S‐P) equations. First, we present a sufficient condition on the blow up of smooth solutions to the N‐S‐P system. Then we construct a family of analytical solutions that blow up in finite time. Copyright © 2010 John Wiley & Sons, Ltd.     

On the existence of solutions to the Navier–Stokes–Poisson equations of a two‐dimensional compressible flow

In this paper, we consider the Navier–Stokes–Poisson equations for compressible, barotropic flow in two space dimensions. We introduce useful tools from the theory of Orlicz spaces. Then we prove the existence of globally defined finite energy weak solutions for the pressure satisfying p(?)=a?log^d (?) for large ?. Here d>1 and a>0. Copyright © 2006 John Wiley & Sons, Ltd.     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

Zero‐viscosity‐capillarity limit to rarefaction waves for the 1D compressible Navier–Stokes–Korteweg equations

In this paper, we study the zero viscosity and capillarity limit problem for the one‐dimensional compressible isentropic Navier–Stokes–Korteweg equations when the corresponding Euler equations have rarefaction wave solutions. In the case that either the effects of initial layer are ignored or the rarefaction waves are smooth, we prove that the solutions of the Navier–Stokes–Korteweg equation with centered rarefaction wave data exist for all time and converge to the centered rarefaction waves as the viscosity and capillarity number vanish, and we also obtain a rate of convergence, which is valid uniformly for all time. These results are showed by a scaling argument and elementary energy analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

New exact solutions for seventh‐order Sawada–Kotera–Ito,Lax and Kaup–Kupershmidt equations using Exp‐function method

In this work, Exp‐function method is used to solve three different seventh‐order nonlinear partial differential KdV equations. Sawada–Kotera–Ito, Lax and Kaup–Kupershmidt equations are well known and considered for solve. Exp‐function method can be used as an alternative to obtain analytic and approximate solutions of different types of differential equations applied in engineering mathematics. Ultimately this method is implemented to solve these equations and convenient and effective solutions are obtained. Copyright © 2009 John Wiley & Sons, Ltd.     

On the pollution effect of quasi‐compressibility methods in magneto‐hydrodynamics and reactive flows

The study of the asymptotics of quasi‐compressibility methods is an essential tool to construct and improve numerical schemes in hydrodynamics, like e.g., stabilized finite element methods or time‐splitting projection methods. The goal of this work is to illustrate different asymptotical solution behaviour for equations in magneto‐hydrodynamics and those that describe reactive flows for standard quasi‐compressibility methods that influences the design of numerical algorithms. Copyright © 1999 John Wiley & Sons, Ltd.     

Decoupled energy‐law preserving numerical schemes for the Cahn–Hilliard–Darcy system

We study two novel decoupled energy‐law preserving and mass‐conservative numerical schemes for solving the Cahn‐Hilliard‐Darcy system which models two‐phase flow in porous medium or in a Hele–Shaw cell. In the first scheme, the velocity in the Cahn–Hilliard equation is treated explicitly so that the Darcy equation is completely decoupled from the Cahn–Hilliard equation. In the second scheme, an intermediate velocity is used in the Cahn–Hilliard equation which allows for the decoupling. We show that the first scheme preserves a discrete energy law with a time‐step constraint, while the second scheme satisfies an energy law without any constraint and is unconditionally stable. Ample numerical experiments are performed to gauge the efficiency and robustness of our scheme. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 936–954, 2016     

Application of He's variational iteration method and Adomian's decomposition method to Sawada–Kotera–Ito seventh‐order equation

In this article, the Sawada–Kotera–Ito seventh‐order equation is studied. He's variational iteration method and Adomian's decomposition method (ADM) are applied to obtain solution of this equation. We compare these methods together. The study highlights the significant features of the employed methods and its capability of handling completely integrable equations. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 887–897, 2011