Fourth‐order compact difference schemes for 1D nonlinear Kuramoto–Tsuzuki equation

In this article, first, we establish some compact finite difference schemes of fourth‐order for 1D nonlinear Kuramoto–Tsuzuki equation with Neumann boundary conditions in two boundary points. Then, we provide numerical analysis for one nonlinear compact scheme by transforming the nonlinear compact scheme into matrix form. And using some novel techniques on the specific matrix emerged in this kind of boundary conditions, we obtain the priori estimates and prove the convergence in norm. Next, we analyze the convergence and stability for one of the linearized compact schemes. To obtain the maximum estimate of the numerical solutions of the linearized compact scheme, we use the mathematical induction method. The treatment is that the convergence in norm is obtained as well as the maximum estimate, further the convergence in norm. Finally, numerical experiments demonstrate the theoretical results and show that one of the linearized compact schemes is more accurate, efficient and robust than the others and the previous. It is worthwhile that the compact difference methods presented here can be extended to 2D case. As an example, we present one nonlinear compact scheme for 2D Ginzburg–Landau equation and numerical tests show that the method is accurate and effective. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2080–2109, 2015     

Finite element approximation of a phase field model arising in nanostructure patterning

We consider a fully practical finite element approximation of the nonlinear parabolic Cahn–Hilliard system subject to an initial condition on the conserved order parameter , and mixed boundary conditions. Here, is the interfacial parameter, is the field strength parameter, is the obstacle potential, is the diffusion coefficient, and denotes differentiation with respect to the second argument. Furthermore, w is the chemical potential and is the electrostatic potential. The system, in the context of nanostructure patterning, has been proposed to model the manipulation of morphologies in organic solar cells with the help of an applied electric field. In the limit , it reduces to a sharp interface problem that models the evolution of an unstable interface between two dielectric media in the presence of a quasistatic electric field. On introducing a finite element approximation for the above Cahn–Hilliard system, we prove existence and stability of a discrete solution. Moreover, in the case of two space dimensions, we are able to prove convergence and hence existence of a solution to the considered system of partial differential equations. We demonstrate the practicality of our finite element approximation with several numerical simulations in two and three space dimensions. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1890–1924, 2015     

A new cubic nonconforming finite element on rectangles

A new nonconforming rectangle element with cubic convergence for the energy norm is introduced. The degrees of freedom (DOFs) are defined by the 12 values at the three Gauss points on each of the four edges. Due to the existence of one linear relation among the above DOFs, it turns out the DOFs are 11. The nonconforming element consists of . We count the corresponding dimension for Dirichlet and Neumann boundary value problems of second‐order elliptic problems. We also present the optimal error estimates in both broken energy and norms. Finally, numerical examples match our theoretical results very well. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 691–705, 2015     

Numerical simulation and solutions of the two‐component second order KdV evolutionarysystem

In this study, with the aid of Wolfram Mathematica 11, the modified exp ‐expansion function method is used in constructing some new analytical solutions with novel structure such as the trigonometric and hyperbolic function solutions to the well‐known nonlinear evolutionary equation, namely; the two‐component second order KdV evolutionary system. Second, the finite forward difference method is used in analyzing the numerical behavior of this equation. We consider equation (6.5) and (6.6) for the numerical analysis. We examine the stability of the two‐component second order KdV evolutionary system with the finite forward difference method by using the Fourier‐Von Neumann analysis. We check the accuracy of the finite forward difference method with the help of and norm error. We present the comparison between the exact and numerical solutions of the two‐component second order KdV evolutionary system obtained in this article which and support with graphics plot. We observed that the modified exp ‐expansion function method is a powerful approach for finding abundant solutions to various nonlinear models and also finite forward difference method is efficient for examining numerical behavior of different nonlinear models.     

Decoupled schemes for unsteady MHD equations. I. time discretization

In this article, a time discretization decoupled scheme for two‐dimensional magnetohydrodynamics equations is proposed. The almost unconditional stability and convergence of this scheme are provided. The optimal error estimates for velocity and magnet are provided, and the optimal error estimate for pressure are deduced as well. Finite element spatial discretization and numerical implementation are considered in our article (Zhang and He, Comput Math Appl 69 (2015), 1390–1406). © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 956–973, 2017     

Analysis of expanded mixed finite element methods for the generalized forchheimer flows of slightly compressible fluids

The nonlinear Forchheimer equations are used to describe the dynamics of fluid flows in porous media when Darcy's law is not applicable. In this article, we consider the generalized Forchheimer flows for slightly compressible fluids, and then study the expanded mixed finite element method applied to the initial boundary value problem for the resulting degenerate parabolic equation for pressure. The bounds for the solutions, time derivative, and gradient of solutions are established. Utilizing the monotonicity properties of Forchheimer equation and boundedness of solutions, a priori error estimates for solution are obtained in ‐norm, ‐norm as well as for its gradient in ‐norm for all . Optimal ‐error estimates are shown for solutions under some additional regularity assumptions. Numerical results using the lowest order Raviart–Thomas mixed element confirm the theoretical analysis regarding convergence rates. Published 2015. Numer Methods Partial Differential Eq 32: 60–85, 2016     

A Fourier pseudospectral method for the “good” Boussinesq equation with second‐order temporal accuracy

In this article, we discuss the nonlinear stability and convergence of a fully discrete Fourier pseudospectral method coupled with a specially designed second‐order time‐stepping for the numerical solution of the “good” Boussinesq equation. Our analysis improves the existing results presented in earlier literature in two ways. First, a convergence for the solution and convergence for the time‐derivative of the solution are obtained in this article, instead of the convergence for the solution and the convergence for the time‐derivative, given in De Frutos, et al., Math Comput 57 (1991), 109–122. In addition, we prove that this method is unconditionally stable and convergent for the time step in terms of the spatial grid size, compared with a severe restriction time step restriction required by the proof in De Frutos, et al., Math Comput 57 (1991), 109–122.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 202–224, 2015     

Unconditional convergence and error estimates for bounded numerical solutions of the barotropic Navier–Stokes system

We consider a mixed finite‐volume finite‐element method applied to the Navier–Stokes system of equations describing the motion of a compressible, barotropic, viscous fluid. We show convergence as well as error estimates for the family of numerical solutions on condition that: (a) the underlying physical domain as well as the data are smooth; (b) the time step and the parameter of the spatial discretization are proportional, ; and (c) the family of numerical densities remains bounded for . No a priori smoothness is required for the limit (exact) solution. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1208–1223, 2017     

Two unconditionally stable and convergent difference schemes with the extrapolation method for the one‐dimensional distributed‐order differential equations

The Grünwald formula is used to solve the one‐dimensional distributed‐order differential equations. Two difference schemes are derived. It is proved that the schemes are unconditionally stable and convergent with the convergence orders and in maximum norm, respectively, where and are step sizes in time, space and distributed order. The extrapolation method is applied to improve the approximate accuracy to the orders and respectively. An illustrative numerical example is given to confirm the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 591–615, 2016     

Superconvergence of Ritz‐Galerkin finite element approximations for second order elliptic problems

In this paper, the author derives an ‐superconvergence for the piecewise linear Ritz‐Galerkin finite element approximations for the second‐order elliptic equation equipped with Dirichlet boundary conditions. This superconvergence error estimate is established between the finite element solution and the usual Lagrange nodal point interpolation of the exact solution, and thus the superconvergence at the nodal points of each element. The result is based on a condition for the finite element partition characterized by the coefficient tensor and the usual shape functions on each element, called ‐equilateral assumption in this paper. Several examples are presented for the coefficient tensor and finite element triangulations which satisfy the conditions necessary for superconvergence. Some numerical experiments are conducted to confirm this new theory of superconvergence.     

High‐order compact difference schemes for the modified anomalous subdiffusion equation

In this article, two kinds of high‐order compact finite difference schemes for second‐order derivative are developed. Then a second‐order numerical scheme for a Riemann–Liouvile derivative is established based on a fractional centered difference operator. We apply these methods to a fractional anomalous subdiffusion equation to construct two kinds of novel numerical schemes. The solvability, stability, and convergence analysis of these difference schemes are studied by using Fourier method. The convergence orders of these numerical schemes are and , respectively. Finally, numerical experiments are displayed which are in line with the theoretical analysis. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 213–242, 2016     

A new nonconforming finite element with a conforming finite element approximation for a coupled continuum pipe‐flow/Darcy model in Karst aquifers

Numerical method is considered for a coupled continuum pipe‐flow/Darcy model describing flow in porous media with an embedded conduit pipe. A new nonconforming element is constructed to solve the Darcy equation on porous matrix. The existence and uniqueness of the approximation solution are deduced. Optimal error estimates are obtained in and norms. Some numerical examples show the accuracy and efficiency of the presented method. With the same number of nodal‐points and the same amount of computation, the results using the new nonconforming element are much better than those by both conforming element and Wilson nonconforming element on the same mesh. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 778–798, 2016     

A fourth‐order Hermitian box‐scheme with fast solver for the Poisson problem in a cube

We present a fourth‐order Hermitian box‐scheme (HB‐scheme) for the Poisson problem in a cube. A single‐nonstaggered regular grid is used supporting the discrete unknowns u and . The scheme is fourth‐order accurate for u and in norm. The fast numerical resolution uses a matrix capacitance method, resulting in a computational complexity of . Numerical results are reported on several examples including nonseparable problems. The present scheme is the extension to the three‐dimensional case of the HB‐scheme presented in Abbas and Croisille [J Sci Comp 49 (2011), 239–267]. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 609–629, 2015     

Optimal error estimate of a conformal Fourier pseudo‐spectral method for the damped nonlinear Schrödinger equation

In this article, a Fourier pseudospectral method, which preserves the conforal conservation la, is proposed for solving the damped nonlinear Schrödinger equation. Based on the energy method and the semi‐norm equivalence between the Fourier pseudospectral method and the finite difference method, a priori estimate for the new method is established, which shows that the proposed method is unconditionally convergent with order of in the discrete ‐norm, where is the time step and is the number of collocation points used in the spectral method. Some numerical results are addressed to confirm our theoretical analysis.     

On the Euler implicit/explicit iterative scheme for the stationary Oldroyd fluid

In this article, we consider the stationary Oldroyd fluid equations from the large time behavior research of the nonstationary equations. Thus, to obtain its numerical solution, we first solve the nonstationary Oldroyd fluid equations via the Euler implicit/explicit finite element method with the integral term discretized by the right‐hand rectangle rule, then increase the total time (i.e., number of time steps) to approximate the solution of the original stationary equations. Under a new uniqueness condition (A2), we prove the exponential stability of the solution pair for the stationary equations and the almost unconditional stability of the numerical method. Furthermore, we also obtain the uniform optimal and error estimates in time integral . Finally, several numerical experiments are provided to verify our theoretical results.     

Convergence analysis of a new mixed finite element method for Biot's consolidation model

In this article, we propose a mixed finite element method for the two‐dimensional Biot's consolidation model of poroelasticity. The new mixed formulation presented herein uses the total stress tensor and fluid flux as primary unknown variables as well as the displacement and pore pressure. This method is based on coupling two mixed finite element methods for each subproblem: the standard mixed finite element method for the flow subproblem and the Hellinger–Reissner formulation for the mechanical subproblem. Optimal a‐priori error estimates are proved for both semidiscrete and fully discrete problems when the Raviart–Thomas space for the flow problem and the Arnold–Winther space for the elasticity problem are used. In particular, optimality in the stress, displacement, and pressure has been proved in when the constrained‐specific storage coefficient is strictly positive and in the weaker norm when is nonnegative. We also present some of our numerical results.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1189–1210, 2014     

Unconditional superconvergence analysis of an H1‐galerkin mixed finite element method for nonlinear Sobolev equations

An efficient H¹‐Galerkin mixed finite element method (MFEM) is presented with and zero order Raviart‐Thomas elements for the nonlinear Sobolev equations. On one hand, the existence and uniqueness of the solutions of the semidiscrete approximation scheme are proved and the super close results of order for the original variable u in a broken H¹ norm and the auxiliary variable in norm are deduced without the boundedness of the numerical solution in ‐norm. Conversely, a linearized Crank‐Nicolson fully discrete scheme with the unconditional super close property is also developed through a new approach, while previous literature always require certain time step conditions (see the references below). Finally, a numerical experiment is included to illustrate the feasibility of the proposed method. Here h is the subdivision parameter and τ is the time step.     

High‐order algorithms for Riesz derivative and their applications (V)

In this article, based on the idea of combing symmetrical fractional centred difference operator with compact technique, a series of even‐order numerical differential formulas (named the fractional‐compact formulas) are established for the Riesz derivatives with order . Properties of coefficients in the derived formulas are studied in details. Then applying the constructed fourth‐order formula, a difference scheme is proposed to solve the Riesz spatial telegraph equation. By the energy method, the constructed numerical algorithm is proved to be stable and convergent with order , where τ and h are the temporal and spatial stepsizes, respectively. Finally, several numerical examples are presented to verify the theoretical results.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1754–1794, 2017     

Least‐squares spectral method for velocity‐vorticity‐pressure form of the Stokes equations

The aim of this article is to present and analyze first‐order system least‐squares spectral method for the Stokes equations in two‐dimensional spaces. The Stokes equations are transformed into a first‐order system of equations by introducing vorticity as a new variable. The least‐squares functional is then defined by summing up the ‐ and ‐norms of the residual equations. The ‐norm in the least‐squares functional is replaced by suitable operator. Continuous and discrete homogeneous least‐squares functionals are shown to be equivalent to ‐norm of velocity and ‐norm of vorticity and pressure for spectral Galerkin and pseudospectral method. The spectral convergence of the proposed methods are given and the theory is validated by numerical experiment. Mass conservation is also briefly investigated. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 661–680, 2016     

A finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors

An implicit Euler finite‐volume scheme for a spinorial matrix drift‐diffusion model for semiconductors is analyzed. The model consists of strongly coupled parabolic equations for the electron density matrix or, alternatively, of weakly coupled equations for the charge and spin‐vector densities, coupled to the Poisson equation for the electric potential. The equations are solved in a bounded domain with mixed Dirichlet–Neumann boundary conditions. The charge and spin‐vector fluxes are approximated by a Scharfetter–Gummel discretization. The main features of the numerical scheme are the preservation of nonnegativity and bounds of the densities and the dissipation of the discrete free energy. The existence of a bounded discrete solution and the monotonicity of the discrete free energy are proved. For undoped semiconductor materials, the numerical scheme is unconditionally stable. The fundamental ideas are reformulations using spin‐up and spin‐down densities and certain projections of the spin‐vector density, free energy estimates, and a discrete Moser iteration. Furthermore, numerical simulations of a simple ferromagnetic‐layer field‐effect transistor in two space dimensions are presented. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 819–846, 2016