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.     

Numerical method for generalized time fractional KdV-type equation

In this article, an efficient numerical method for linearized and nonlinear generalized time-fractional KdV-type equations is proposed by combining the finite difference scheme and Petrov–Galerkin spectral method. The scale and weight functions involved in generalized fractional derivative cause too much difficulty in discretization and numerical analysis. Fortunately, motivated by finite difference method for fractional differential equation on graded mesh, the stability and convergence of the constructed method are established rigorously. It is proved that the full discretization schemes of generalized time-fractional KdV-type equation is unconditionally stable in linear case. While for nonlinear case, it is stable under a CFL condition and for not small ϵ, coefficient of the high-order spatial differential term. In addition, the full discretization schemes with respect to linear and nonlinear cases respectively converge to the associated exact solutions with orders and , where τ, α, N and m accordingly indicate the time step size, the order of the fractional derivative, polynomial degree, and regularity of the exact solution. Numerical experiments are carried out to support the theoretical results.     

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.     

Analysis and efficient implementation of alternating direction implicit finite volume method for Riesz space‐fractional diffusion equations in two space dimensions

In this article, we develop a Crank–Nicolson alternating direction implicit finite volume method for time‐dependent Riesz space‐fractional diffusion equation in two space dimensions. Norm‐based stability and convergence analysis are given to show that the developed method is unconditionally stable and of second‐order accuracy both in space and time. Furthermore, we develop a lossless matrix‐free fast conjugate gradient method for the implementation of the numerical scheme, which only has memory requirement and computational complexity per iteration with N being the total number of spatial unknowns. Several numerical experiments are presented to demonstrate the effectiveness and efficiency of the proposed scheme for large‐scale modeling and simulations.     

Regularization of a final value problem for a nonlinear biharmonic equation

In this paper, we consider the nonlinear biharmonic equation. The problem is ill‐posed in the sense of Hadamard. To obtain a stable numerical solution, we consider a regularization method. We show rigourously, with error estimates provided, that the corresponding regularized solutions converge to the true solution strongly in uniformly with respect to the space coordinate under some a priori assumptions on the solution.     

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.     

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.     

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.     

Optimal error estimate of two linear and momentum-preserving Fourier pseudo-spectral schemes for the RLW equation

In this paper, two novel linear-implicit and momentum-preserving Fourier pseudo-spectral schemes are proposed and analyzed for the regularized long-wave equation. The numerical methods are based on the blend of the Fourier pseudo-spectral method in space and the linear-implicit Crank–Nicolson method or the leap-frog scheme in time. The two fully discrete linear schemes are shown to possess the discrete momentum conservation law, and the linear systems resulting from the schemes are proved uniquely solvable. Due to the momentum conservative property of the proposed schemes, the Fourier pseudo-spectral solution is proved to be bounded in the discrete L^∞ norm. Then by using the standard energy method, both the linear-implicit Crank–Nicolson momentum-preserving scheme and the linear-implicit leap-frog momentum-preserving scheme are shown to have the accuracy of in the discrete L^∞ norm without any restrictions on the grid ratio, where N is the number of nodes and τ is the time step size. Numerical examples are carried out to verify the correction of the theory analysis and the efficiency of the proposed schemes.     

Two‐grid method for two‐dimensional nonlinear Schrödinger equation by finite element method

A conservative two‐grid finite element scheme is presented for the two‐dimensional nonlinear Schrödinger equation. One Newton iteration is applied on the fine grid to linearize the fully discrete problem using the coarse‐grid solution as the initial guess. Moreover, error estimates are conducted for the two‐grid method. It is shown that the coarse space can be extremely coarse, with no loss in the order of accuracy, and still achieve the asymptotically optimal approximation as long as the mesh sizes satisfy in the two‐grid method. The numerical results show that this method is very effective.     

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 .     

An efficient numerical approach to solve Schrödinger equations with space fractional derivative

We design and analyze an efficient numerical approach to solve the coupled Schrödinger equations with space‐fractional derivative. The numerical scheme is based on leap‐frog in time direction and Fourier method in spatial direction. The advantage of the numerical scheme is that only a linear equation needs to be solved for each time step size, and we proved that the energy and mass of space‐fractional coupled Schrödinger equations (SFCSEs) are conserved in the case of full‐discrete scheme. Moreover, we also analyze the error estimate of the numerical scheme, and numerical solutions converge with the order in L² norm. Numerical examples are illustrated to verify the theoretical results.     

Orthogonal spline collocation scheme for multiterm fractional convection‐diffusion equation with variable coefficients

The orthogonal spline collocation (OSC) technique is an efficient way to solve a wide variety of problems that are modeled by ordinary and partial differential equations. In this article, by using OSC method in spatial direction and classical L1 approximation in temporal direction, a fully discrete scheme is established for a class of two‐dimensional multiterm fractional convection‐diffusion reaction equation with variable coefficients. The optimal estimates in H^j (j = 0, 1, 2) norms at each time step are derived. Also, estimate in space is provided. At last, we provide some numerical results to verify the accuracy and efficiency of the proposed algorithm.     

Finite analytic numerical method for three‐dimensional quasi‐laplace equation with conductivity in tensor form

The finite analytic numerical method for 3D quasi‐Laplace equation with conductivity in full tensor form is constructed in this article. For cubic grid system, the gradient of the potential variable will diverge when tending to the common edge joining the four grids with different conductivities. However, the potential gradient along the tangential direction is of limited value. As a consequence, the 3D quasi‐Laplace equations will behave as a quasi‐2D one. An approximate analytical solution of the 3D quasi‐Laplace equation can be found around the common edge, which is expressed as a combination of a power‐law function and a linear function. With the help of this approximate analytical solution, a 3D finite analytical numerical scheme is then constructed. Numerical examples show that the proposed numerical scheme can provide rather accurate solutions only with or subdivisions. More important, the convergent speed of the numerical scheme is independent of the conductivity heterogeneity. In contrast, when using the traditional numerical schemes, typically such as the MPFA method, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result for the strong heterogeneous case.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1475–1492, 2017     

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.     

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     

Heinz‐type inequality and bi‐Lipschitz continuity for quasiconformal mappings satisfying inhomogeneous biharmonic equations

Let be the class of all sense‐preserving homeomorphic self‐mappings of . The aim of this paper is twofold. First, we obtain Heinz‐type inequality for (K,K′)‐quasiconformal mappings satisfying inhomogeneous biharmonic equation Δ(Δω) = g in unit disk with associated boundary value conditions and . Second, we establish biLipschitz continuity for (K,K′)‐quasiconformal mappings satisfying aforementioned inhomogeneous biharmonic equation when and are small enough.     

Convergence analysis for second‐order accurate schemes for the periodic nonlocal Allen‐Cahn and Cahn‐Hilliard equations

In this paper, we provide a detailed convergence analysis for fully discrete second‐order (in both time and space) numerical schemes for nonlocal Allen‐Cahn and nonlocal Cahn‐Hilliard equations. The unconditional unique solvability and energy stability ensures ? ⁴ stability. The convergence analysis for the nonlocal Allen‐Cahn equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nonlocal Cahn‐Hilliard equation, because of the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken, and an H ^?1 inner‐product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a priori bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (eg, 1964), 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 O (s ³+h ⁴) convergence in norm, in which s and h denote the time step and spatial mesh sizes, respectively. This in turn leads to the necessary bound under a standard constraint s ≤C h . Here, we also prove convergence of the scheme in the maximum norm under the same constraint.     

Numerical solutions to a BBM‐Burgers model with a nonlocal viscous term

In this paper, we numerically investigate the BBM‐Burgers equation with a nonlocal viscous term (1) where is the Riemann‐Liouville half derivative. In particular, we implement different numerical schemes to approximate the solution and its asymptotical behavior. Also, we compare our numerical results with those given in 2013, 2014 for similar models.     

An efficient difference scheme for the coupled nonlinear fractional Ginzburg–Landau equations with the fractional Laplacian

In this article, an efficient difference scheme for the coupled fractional Ginzburg–Landau equations with the fractional Laplacian is studied. We construct the discrete scheme based on the implicit midpoint method in time and a weighted and shifted Grünwald difference method in space. Then, we prove that the scheme is uniquely solvable, and the numerical solutions are bounded and unconditionally convergent in the norm. Finally, numerical tests are given to confirm the theoretical results and show the effectiveness of the scheme.