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.     

Existence of global attractors for the three‐dimensional Brinkman–Forchheimer equation

In this paper, we investigate the asymptotic behavior of solutions of the three‐dimensional Brinkman–Forchheimer equation. We first prove the existence and uniqueness of solutions of the equation in L², and then show that the equation has a global attractor in H² when the external forcing term belongs to L². Copyright © 2008 John Wiley & Sons, Ltd.     

Stability and error analysis for spectral Galerkin method for the Navier–Stokes equations with L2 initial data

In this article we consider the spectral Galerkin method with the implicit/explicit Euler scheme for the two‐dimensional Navier–Stokes equations with the L² initial data. Due to the poor smoothness of the solution on [0,1), we use the the spectral Galerkin method based on high‐dimensional spectral space H_M and small time step Δt² on this interval. While on [1,∞), we use the spectral Galerkin method based on low‐dimensional spectral space H_m(m = O(M^1/2)) and large time step Δt. For the spectral Galerkin method, we provide the standard H²‐stability and the L²‐error analysis. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A fully discrete spectral method for fractional Cattaneo equation based on Caputo–Fabrizo derivative

Recently Caputo and Fabrizio introduced a new derivative with fractional order without singular kernel. The derivative can be used to describe the material heterogeneities and the fluctuations of different scales. In this article, we derived a new discretization of Caputo–Fabrizio derivative of order α (1 < α < 2) and applied it into the Cattaneo equation. A fully discrete scheme based on finite difference method in time and Legendre spectral approximation in space is proposed. The stability and convergence of the fully discrete scheme are rigorously established. The convergence rate of the fully discrete scheme in H¹ norm is O(τ² + N^1?m), where τ, N and m are the time‐step size, polynomial degree and regularity in the space variable of the exact solution, respectively. Furthermore, the accuracy and applicability of the scheme are confirmed by numerical examples to support the theoretical results.     

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.     

A linearized Crank–Nicolson–Galerkin FEM for the time‐dependent Ginzburg–Landau equations under the temporal gauge

We propose a decoupled and linearized fully discrete finite element method (FEM) for the time‐dependent Ginzburg–Landau equations under the temporal gauge, where a Crank–Nicolson scheme is used for the time discretization. By carefully designing the time‐discretization scheme, we manage to prove the convergence rate , where τ is the time‐step size and r is the degree of the finite element space. Due to the degeneracy of the problem, the convergence rate in the spatial direction is one order lower than the optimal convergence rate of FEMs for parabolic equations. Numerical tests are provided to support our error analysis. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1279–1290, 2014     

Analysis and computational method based on quadratic B‐spline FEM for the Rosenau–Burgers equation

L_∞‐error estimates for B‐spline Galerkin finite element solution of the Rosenau–Burgers equation are considered. The semidiscrete B‐spline Galerkin scheme is studied using appropriate projections. For fully discrete B‐spline Galerkin scheme, we consider the Crank–Nicolson method and analyze the corresponding error estimates in time. Numerical experiments are given to demonstrate validity and order of accuracy of the proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 877–895, 2016     

A Petrov‐Galerkin method with quadrature for semi‐linear second‐order hyperbolic problems

We propose and analyze a Crank–Nicolson quadrature Petrov–Galerkin (CNQPG) ‐spline method for solving semi‐linear second‐order hyperbolic initial‐boundary value problems. We prove second‐order convergence in time and optimal order H² norm convergence in space for the CNQPG scheme that requires only linear algebraic solvers. We demonstrate numerically optimal order H^k, k = 0,1,2, norm convergence of the scheme for some test problems with smooth and nonsmooth nonlinearities. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Numerical approximation of the generalized regularized long wave equation using Petrov–Galerkin finite element method

The generalized regularized long wave (GRLW) equation has been developed to model a variety of physical phenomena such as ion‐acoustic and magnetohydrodynamic waves in plasma, nonlinear transverse waves in shallow water and phonon packets in nonlinear crystals. This paper aims to develop and analyze a powerful numerical scheme for the nonlinear GRLW equation by Petrov–Galerkin method in which the element shape functions are cubic and weight functions are quadratic B‐splines. The proposed method is implemented to three reference problems involving propagation of the single solitary wave, interaction of two solitary waves and evolution of solitons with the Maxwellian initial condition. The variational formulation and semi‐discrete Galerkin scheme of the equation are firstly constituted. We estimate rate of convergence of such an approximation. Using Fourier stability analysis of the linearized scheme we show that the scheme is unconditionally stable. To verify practicality and robustness of the new scheme error norms L₂, L_∞ and three invariants I₁, I₂, and I₃ are calculated. The computed numerical results are compared with other published results and confirmed to be precise and effective.     

A three‐level linearized finite difference scheme for the camassa–holm equation

The Camassa–Holm (CH) system is a strong nonlinear third‐order evolution equation. So far, the numerical methods for solving this problem are only a few. This article deals with the finite difference solution to the CH equation. A three‐level linearized finite difference scheme is derived. The scheme is proved to be conservative, uniquely solvable, and conditionally second‐order convergent in both time and space in the discrete L_∞ norm. Several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 451–471, 2014     

A superconvergent nonconforming mixed finite element method for the Navier–Stokes equations

The superconvergence for a nonconforming mixed finite element approximation of the Navier–Stokes equations is analyzed in this article. The velocity field is approximated by the constrained nonconforming rotated Q₁ (CNRQ₁) element, and the pressure is approximated by the piecewise constant functions. Under some regularity assumptions, the superconvergence estimates for both the velocity in broken H¹‐norm and the pressure in L²‐norm are obtained. Some numerical examples are presented to demonstrate our theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 646–660, 2016     

Ergodicity of 2D stochastic Ginzburg–Landau–Newell equations driven by degenerate noise

The current paper is devoted to stochastic Ginzburg–Landau–Newell equation with degenerate random forcing. The existence and pathwise uniqueness of strong solutions with H¹‐initial data is established, and then the existence of an invariant measure for the Feller semigroup is shown by Krylov–Bogoliubov theorem. Because of the coupled items in the stochastic Ginzburg–Landau–Newell equations, the higher order momentum estimates can be only obtained in the L²‐norm. We show the ergodicity of invariant measure for the transition semigroup by asymptotically strong Feller property and the support property. Copyright © 2017 John Wiley & Sons, Ltd.     

A Priori Error Estimates of Crank–Nicolson Finite Volume Element Method for a Hyperbolic Optimal Control Problem

In this article, a Crank–Nicolson linear finite volume element scheme is developed to solve a hyperbolic optimal control problem. We use the variational discretization technique for the approximation of the control variable. The optimal convergent order O(h² + k²) is proved for the numerical solution of the control, state and adjoint‐state in a discrete L²‐norm. To derive this result, we also get the error estimate (convergent order O(h² + k²)) of Crank–Nicolson finite volume element approximation for the second‐order hyperbolic initial boundary value problem. Numerical experiments are presented to verify the theoretical results.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1331–1356, 2016     

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     

Unconditional superconvergent analysis of a new mixed finite element method for Ginzburg–Landau equation

In this article, unconditional superconvergent analysis of a linearized fully discrete mixed finite element method is presented for a class of Ginzburg–Landau equation based on the bilinear element and zero‐order Nédélec's element pair (Q₁₁/Q₀₁ × Q₁₀). First, a time‐discrete system is introduced to split the error into temporal error and spatial error, and the corresponding error estimates are deduced rigorously. Second, the unconditional superclose and optimal estimate of order O(h² + τ) for u in H¹‐norm and p = ?u in L²‐norm are derived respectively without the restrictions on the ratio between h and τ, where h is the subdivision parameter and τ, the time step. Third, the global superconvergent results are obtained by interpolated postprocessing technique. Finally, some numerical results are carried out to confirm the theoretical analysis.     

A Petrov–Galerkin spectral method for fourth‐order problems

In this paper, we consider the Petrov–Galerkin spectral method for fourth‐order elliptic problems on rectangular domains subject to non‐homogeneous Dirichlet boundary conditions. We derive some sharp results on the orthogonal approximations in one and two dimensions, which play important roles in numerical solutions of higher‐order problems. By applying these results to a fourth‐order problem, we establish the H²‐error and L²‐error bounds of the Petrov–Galerkin spectral method. Numerical experiments are provided to illustrate the high accuracy of the proposed method and coincide well with the theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Existence and exponential stability in Lr‐spaces of stationary Navier–Stokes flows with prescribed flux in infinite cylindrical domains

We prove existence, uniqueness and exponential stability of stationary Navier–Stokes flows with prescribed flux in an unbounded cylinder of ?ⁿ,n?3, with several exits to infinity provided the total flux and external force are sufficiently small. The proofs are based on analytic semigroup theory, perturbation theory and L^r ? L^q‐estimates of a perturbation of the Stokes operator in L^q‐spaces. Copyright © 2006 John Wiley & Sons, Ltd.     

A high‐order compact scheme for the nonlinear fractional Klein–Gordon equation

In this article, a high‐order finite difference scheme for a kind of nonlinear fractional Klein–Gordon equation is derived. The time fractional derivative is described in the Caputo sense. The solvability of the difference system is discussed by the Leray–Schauder fixed point theorem, while the stability and L_∞ convergence of the finite difference scheme are proved by the energy method. Numerical examples are provided to demonstrate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 706–722, 2015     

Approximation of the long‐term dynamics of the dynamical system generated by the multilayer quasigeostrophic equations of the ocean

In this article, we study the multilayer quasigeostrophic equations of the ocean. More precisely, we discretize these equations in time using the implicit Euler scheme and using the classical and uniform discrete Gronwall lemmas we prove that the approximate solution is uniformly bounded in H^?1, L² and H¹. Using the uniform stability of the scheme and the theory of the multivalued attractors, we then prove that the discrete attractors generated by the numerical scheme converge to the global attractor of the continuous system as the time‐step approaches zero. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1041–1065, 2016     

Numerical solution of the Rosenau–KdV–RLW equation by operator splitting techniques based on B‐spline collocation method

In the present study, the operator splitting techniques based on the quintic B‐spline collocation finite element method are presented for calculating the numerical solutions of the Rosenau–KdV–RLW equation. Two test problems having exact solutions have been considered. To demonstrate the efficiency and accuracy of the present methods, the error norms L₂ and L_∞ with the discrete mass Q and energy E conservative properties have been calculated. The results obtained by the method have been compared with the exact solution of each problem and other numerical results in the literature, and also found to be in good agreement with each other. A Fourier stability analysis of each presented method is also investigated.