Numerical analysis for a system coupling curve evolution attached orthogonally to a fixed boundary,to a reaction–diffusion equation on the curve

We consider a semidiscrete finite element approximation for a system consisting of the evolution of a planar curve evolving by forced curve shortening flow inside a given bounded domain , such that the curve meets the boundary orthogonally, and the forcing is a function of the solution of a reaction–diffusion equation that holds on the evolving curve. We prove optimal order error bounds for the resulting approximation and present numerical experiments.     

Superconvergence analysis of an energy stable scheme with three step backward differential formula-finite element method for nonlinear reaction–diffusion equation

A three step backward differential formula scheme is proposed for nonlinear reaction–diffusion equation and superconvergence results are studied with Galerkin finite element method unconditionally. Energy stability is testified for the constructed scheme with an artificial term. Splitting technique is utilized to get rid of the ratio between the time step size and the subdivision parameter . Temporal error estimate in H²-norm is derived, which leads to the boundedness of the solutions of the time-discrete equations. Unconditional spatial error estimate in L²-norm is deduced which help bound the numerical solutions in L^∞-norm. Superconvergent property of in H¹-norm with order is obtained by taking difference between two time levels of the error equations unconditionally. The global superconvergent property is deduced through the above results. Two numerical examples show the validity of the theoretical analysis.     

A numerical algorithm for the nonlinear Timoshenko beam system

An initial boundary value problem is considered for the dynamic beam system Its solution is found by means of an algorithm, the constituent parts of which are the finite element method, the implicit symmetric difference scheme used to approximate the solution with respect to the spatial and time variables, and also a Picard type iteration process for solving the system of nonlinear equations obtained by discretization. Errors of three parts of the algorithm are estimated and, as a result, its total error estimate is obtained. A numerical example is solved.     

On the numerical approximation of p‐biharmonic and ∞‐biharmonic functions

The ∞ ‐Bilaplacian is a third‐order fully nonlinear PDE given by (1) In this work, we build a numerical method aimed at quantifying the nature of solutions to this problem, which we call ∞ ‐biharmonic functions. For fixed p we design a mixed finite element scheme for the prelimiting equation, the p‐Bilaplacian (2) We prove convergence of the numerical solution to the weak solution of and show that we are able to pass to the limit p → ∞ . We perform various tests aimed at understanding the nature of solutions of and we prove convergence of our discretization to an appropriate weak solution concept of this problem that of ‐solutions.     

A posteriori error estimator based on derivative recovery for the discontinuous Galerkin method for nonlinear hyperbolic conservation laws on Cartesian grids

In this article, we develop and analyze a new recovery‐based a posteriori error estimator for the discontinuous Galerkin (DG) method for nonlinear hyperbolic conservation laws on Cartesian grids, when the upwind flux is used. We prove, under some suitable initial and boundary discretizations, that the ‐norm of the solution is of order , when tensor product polynomials of degree at most are used. We further propose a very simple derivative recovery formula which gives a superconvergent approximation to the directional derivative. The order of convergence is showed to be . We use our derivative recovery result to develop a robust recovery‐type a posteriori error estimator for the directional derivative approximation which is based on an enhanced recovery technique. The proposed error estimators of the recovery‐type are easy to implement, computationally simple, asymptotically exact, and are useful in adaptive computations. Finally, we show that the proposed recovery‐type a posteriori error estimates, at a fixed time, converge to the true errors in the ‐norm under mesh refinement. The order of convergence is proved to be . Our theoretical results are valid for piecewise polynomials of degree and under the condition that each component, , of the flux function possesses a uniform positive lower bound. Several numerical examples are provided to support our theoretical results and to show the effectiveness of our recovery‐based a posteriori error estimator. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1224–1265, 2017     

Asymptotically exact a posteriori local discontinuous Galerkin error estimates for the one‐dimensional second‐order wave equation

In this article, we analyze a residual‐based a posteriori error estimates of the spatial errors for the semidiscrete local discontinuous Galerkin (LDG) method applied to the one‐dimensional second‐order wave equation. These error estimates are computationally simple and are obtained by solving a local steady problem with no boundary condition on each element. We apply the optimal L² error estimates and the superconvergence results of Part I of this work [Baccouch, Numer Methods Partial Differential Equations 30 (2014), 862–901] to prove that, for smooth solutions, these a posteriori LDG error estimates for the solution and its spatial derivative, at a fixed time, converge to the true spatial errors in the L²‐norm under mesh refinement. The order of convergence is proved to be , when p‐degree piecewise polynomials with are used. As a consequence, we prove that the LDG method combined with the a posteriori error estimation procedure yields both accurate error estimates and superconvergent solutions. Our computational results show higher convergence rate. We further prove that the global effectivity indices, for both the solution and its derivative, in the L²‐norm converge to unity at rate while numerically they exhibit and rates, respectively. Numerical experiments are shown to validate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1461–1491, 2015     

Ultraconvergence of finite element method by Richardson extrapolation for elliptic problems with inhomogeneous boundary conditions

In this article, Richardson extrapolation technique is employed to investigate the local ultraconvergence properties of Lagrange finite element method using piecewise polynomials of degrees () for the second order elliptic problem with inhomogeneous boundary. A sequence of special graded partition are proposed and a new interpolation operator is introduced to achieve order local ultraconvergence for the displacement and derivative.     

Fast high order difference schemes for the time fractional telegraph equation

In this paper, a fast high order difference scheme is first proposed to solve the time fractional telegraph equation based on the ℱℒ 2-1_σ formula for the Caputo fractional derivative, which reduces the storage and computational cost for calculation. A compact scheme is then presented to improve the convergence order in space. The unconditional stability and convergence in maximum norm are proved for both schemes, with the accuracy order and , respectively. Difficulty arising from the two Caputo fractional derivatives is overcome by some detailed analysis. Finally, we carry out numerical experiments to show the efficiency and accuracy, by comparing with the ℒ 2-1_σ method.     

Error estimates of fully discrete finite element solutions for the 2D Cahn–Hilliard equation with infinite time horizon

In this article, we deal with a rigorous error analysis for the finite element solutions of the two‐dimensional Cahn–Hilliard equation with infinite time. The error estimates with respect to are proven for the fully discrete conforming piecewise linear element solution under Assumption (A1) on the initial value and Assumption (A2) on the discrete spectrum estimate in the finite element space. The analysis is based on sharp a‐priori estimates for the solutions, particularly reflecting their behavior as . Numerical experiments are carried out to support the theoretical analysis and demonstrate the efficiency of the fully discrete mixed finite element methods. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 742–762, 2017     

An adaptive least‐squares mixed finite element method for the Signorini problem

We present and analyze a least squares formulation for contact problems in linear elasticity which employs both, displacements and stresses, as independent variables. As a consequence, we obtain stability and high accuracy of our discretization also in the incompressible limit. Moreover, our formulation gives rise to a reliable and efficient a posteriori error estimator. To incorporate the contact constraints, the first‐order system least squares functional is augmented by a contact boundary functional which implements the associated complementarity condition. The bilinear form related to the augmented functional is shown to be coercive and therefore constitutes an upper bound, up to a constant, for the error in displacements and stresses in . This implies the reliability of the functional to be used as an a posteriori error estimator in an adaptive framework. The efficiency of the use of the functional as an a posteriori error estimator is monitored by the local proportion of the boundary functional term with respect to the overall functional. Computational results using standard conforming linear finite elements for the displacement approximation combined with lowest‐order Raviart‐Thomas elements for the stress tensor show the effectiveness of our approach in an adaptive framework for two‐dimensional and three‐dimensional Hertzian contact problems. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 276–289, 2017     

Maximum norm stability and error estimates for the evolving surface finite element method

We show convergence in the natural and norm for a semidiscretization with linear finite elements of a linear parabolic partial differential equations on evolving surfaces. To prove this, we show error estimates for a Ritz map, error estimates for the material derivative of a Ritz map and a weak discrete maximum principle.     

H1-superconvergence of finite difference method based on Q1-element on quasi-uniform mesh for the 3D Poisson equation

In this paper, the finite difference (FD) method is considered for the 3D Poisson equation by using the Q₁-element on a quasi-uniform mesh. First, under the regularity assumption of , the H¹-superconvergence of the FD solution u_h based on the Q₁-element to the first-order interpolation function is obtained. Next, the H¹-superconvergence of the second-order interpolation postprocessing function based on the FD solution u_h to u is provided. Finally, numerical tests are presented to show the H¹-superconvergence result of the FD postprocessing function to u if .     

A robust scheme based on novel-operational matrices for some classes of time-fractional nonlinear problems arising in mechanics and mathematical physics

In this paper, we present a novel approach based on shifted Gegenbauer wavelets to attain approximate solutions of some classed of time-fractional nonlinear problems. First, we present the approximation of a function of two variables u(x,t) with help of shifted Gegenbauer wavelets and then some novel operational matrices are proposed with the help of piecewise functions to investigate the positive integer derivative ( D _x and D _t), fractional-order derivative ( and ), fractional-order integration ( and ) and delay terms ( and ) of approximated function u(x,t). In order to transform the discussed nonlinear problem into linear problem Picard iterative scheme has been adopt. The current scheme converts the discussed highly nonlinear time-fractional problem into system of linear algebraic equation the help of developed operational matrices and Picard idea. Analysis on the error bound and convergence to authenticate the mathematical formulation of the computational algorithm. We solve various test problems, such as the van der Pol oscillator model, generalized Burger–Huxley, neutral delay parabolic differential equations, sine-Gordon, parabolic integro-differential equation and nonlinear Schrödinger equations to show the efficiency and accuracy of linearized shifted Gegenbauer wavelets method. A comprehensive comparative examination shows the credibility, accuracy, and reliability of the presently proposed computational approach. Also, this scheme can be extended conveniently to other multi-dimensional physical problems of highly nonlinear fractional or variable order of complex nature.     

Local error estimates of the finite element method for an elliptic problem with a Dirac source term

The solutions of elliptic problems with a Dirac measure right‐hand side are not in dimension and therefore the convergence of the finite element solutions is suboptimal in the ‐norm. In this article, we address the numerical analysis of the finite element method for the Laplace equation with Dirac source term: we consider, in dimension 3, the Dirac measure along a curve and, in dimension 2, the punctual Dirac measure. The study of this problem is motivated by the use of the Dirac measure as a reduced model in physical problems, for which high accuracy of the finite element method at the singularity is not required. We show a quasioptimal convergence in the ‐norm, for on subdomains which exclude the singularity; in the particular case of Lagrange finite elements, an optimal convergence in ‐norm is shown on a family of quasiuniform meshes. Our results are obtained using local Nitsche and Schatz‐type error estimates, a weak version of Aubin‐Nitsche duality lemma and a discrete inf‐sup condition. These theoretical results are confirmed by numerical illustrations.     

Error estimates to smooth solutions of semi‐discrete discontinuous Galerkin methods with quadrature rules for scalar conservation laws

In this article, we focus on error estimates to smooth solutions of semi‐discrete discontinuous Galerkin (DG) methods with quadrature rules for scalar conservation laws. The main techniques we use are energy estimate and Taylor expansion first introduced by Zhang and Shu in (Zhang and Shu, SIAM J Num Anal 42 (2004), 641–666). We show that, with (piecewise polynomials of degree k) finite elements in 1D problems, if the quadrature over elements is exact for polynomials of degree , error estimates of are obtained for general monotone fluxes, and optimal estimates of are obtained for upwind fluxes. For multidimensional problems, if in addition quadrature over edges is exact for polynomials of degree , error estimates of are obtained for general monotone fluxes, and are obtained for monotone and sufficiently smooth numerical fluxes. Numerical results validate our analysis. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 467–488, 2017     

On concentration of solution to a Schrödinger logarithmic equation with deepening potential well

In this work, we prove the existence of positive solution for the following class of problems where λ>0 and is a potential satisfying some conditions. Using the variational method developed by Szulkin for functionals, which are the sum of a C¹ functional with a convex lower semicontinuous functional, we prove that for each large enough λ>0, there exists a positive solution for the problem, and that, as λ→+∞, such solutions converge to a positive solution of the limit problem defined on the domain Ω=int(V^?1({0})).     

Half‐inverse problems for the quadratic pencil of the Sturm–Liouville equations with impulse

In this paper, we consider the inverse spectral problem for the impulsive Sturm–Liouville differential pencils on [0, π] with the Robin boundary conditions and the jump conditions at the point . We prove that two potentials functions on the whole interval and the parameters in the boundary and jump conditions can be determined from a set of eigenvalues for two cases: (i) the potentials given on and (ii) the potentials given on , where 0 < α < 1 , respectively. Inverse spectral problems, Sturm–Liouville operator, spectrum, uniqueness.     

Weak Galerkin finite element methods for linear parabolic integro‐differential equations

The semidiscrete and fully discrete weak Galerkin finite element schemes for the linear parabolic integro‐differential equations are proposed. Optimal order error estimates are established for the corresponding numerical approximations in both and norms. Numerical experiments illustrating the error behaviors are provided.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1357–1377, 2016     

Numerical study of a conservative weighted compact difference scheme for the symmetric regularized long wave equations

In this article, a new weighted and compact conservative difference scheme for the symmetric regularized long wave (SRLW) equations is considered. The new scheme is decoupled and linearized in practical computation, that is, at each time step only two tridiagonal systems of linear algebraic equations need to be solved. It is proved by the discrete energy method that the compact scheme is uniquely solvable, the convergence and stability of the difference scheme are obtained, and its numerical convergence order is in the ‐norm. Numerical experiment results show that the scheme is efficient and reliable.     

A steady barycentric Lagrange interpolation method for the 2D higher‐order time‐fractional telegraph equation with nonlocal boundary condition with error analysis

In this paper, we consider the numerical solution of the time‐fractional telegraph equation with a nonlocal boundary condition. A novel barycentric Lagrange interpolation collocation method is developed to solve this equation. Two difficulties have been sorted: the singularity of the integration and the higher accuracy. At the same, we put forward a steady barycentric Lagrange interpolation technique to overcome the new “Runge” phenomenon in computation. Error estimates of the barycentric Lagrange interpolation and the time‐fractional telegraph system for the present method are presented in Sobolev spaces. High convergence rates of the proposed method are obtained and are consisted with the numerical values. Especially in the time dimension, we get the error bound, for h‐refinement and for n_t‐density in the L₂ norms. The numerical results obtained show that the proposed numerical algorithm is accurate and computationally efficient for solving time‐fractional telegraph equation. Experiments demonstrate the high convergence rates of the proposed method are consisted with the theoretical values.