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     

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     

The Crank‐Nicolson/Adams‐Bashforth scheme for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

In this article, we study the stability and convergence of the Crank‐Nicolson/Adams‐Bashforth scheme for the two‐dimensional nonstationary Navier‐Stokes equations with a nonsmooth initial data. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the implicit Crank‐Nicolson scheme for the linear terms and the explicit Adams‐Bashforth scheme for the nonlinear term. Moreover, we prove that the scheme is almost unconditionally stable for a nonsmooth initial data u₀ with div u₀ = 0, i.e., the time step τ satisfies: τ ≤ C₀ if u₀ ∈ H¹ ∩ L^∞; τ |log h| ≤ C₀ if u₀ ∈ H¹ for the mesh size h and some positive constant C₀. Finally, we obtain some error estimates for the discrete velocity and pressure under the above stability condition. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 155‐187, 2012     

A two-level finite element method for the Navier–Stokes equations based on a new projection

We consider a fully discrete two-level approximation for the time-dependent Navier–Stokes equations in two dimension based on a time-dependent projection. By defining this new projection, the iteration between large and small eddy components can be reflected by its associated space splitting. Hence, we can get a weakly coupled system of large and small eddy components. This two-level method applies the finite element method in space and Crank–Nicolson scheme in time. Moreover,the analysis and some numerical examples are shown that the proposed two-level scheme can reach the same accuracy as the classical one-level Crank–Nicolson method with a very fine mesh size h by choosing a proper coarse mesh size H. However, the two-level method will involve much less work.     

Stokes and Navier–Stokes equations with nonhomogeneous boundary conditions

In this paper, we study the existence and regularity of solutions to the Stokes and Oseen equations with nonhomogeneous Dirichlet boundary conditions with low regularity. We consider boundary conditions for which the normal component is not equal to zero. We rewrite the Stokes and the Oseen equations in the form of a system of two equations. The first one is an evolution equation satisfied by Pu, the projection of the solution on the Stokes space – the space of divergence free vector fields with a normal trace equal to zero – and the second one is a quasi-stationary elliptic equation satisfied by (I−P)u, the projection of the solution on the orthogonal complement of the Stokes space. We establish optimal regularity results for Pu and (I−P)u. We also study the existence of weak solutions to the three-dimensional instationary Navier–Stokes equations for more regular data, but without any smallness assumption on the initial and boundary conditions.     

The large time behavior of solutions to 3D Navier–Stokes equations with nonlinear damping

This paper studies the Cauchy problem of the 3D Navier–Stokes equations with nonlinear damping term | u | ^β?1u (β ≥ 1). For β ≥ 3, we derive a decay rate of the L²‐norm of the solutions. Then, the large time behavior is given by comparing the equation with the classic 3D Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.     

A remark on the decay of solutions to the 3‐D Navier–Stokes equations

In this paper we derive a decay rate of the L²‐norm of the solution to the 3‐D Navier–Stokes equations. Although the result which is proved by Fourier splitting method is well known, our method is new, concise and direct. Moreover, it turns out that the new method established here has a wide application on other equations. Copyright © 2007 John Wiley & Sons, Ltd.     

Moment estimates for weak solutions to the Navier–Stokes equations in half‐space

We establish the moment estimates for a class of global weak solutions to the Navier–Stokes equations in the half‐space. Copyright © 2009 John Wiley & Sons, Ltd.     

Stabilized spectral element approximation for the Navier–Stokes equations

The conforming spectral element methods are applied to solve the linearized Navier–Stokes equations by the help of stabilization techniques like those applied for finite elements. The stability and convergence analysis is carried out and essential numerical results are presented demonstrating the high accuracy of the method as well as its robustness. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 115–141, 1998     

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.     

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.     

On the asymptotic stability of steady solutions of the Navier–Stokes equations in unbounded domains

We consider the problem of the asymptotic behaviour in the L²‐norm of solutions of the Navier–Stokes equations. We consider perturbations to the rest state and to stationary motions. In both cases we study the initial‐boundary value problem in unbounded domains with non‐compact boundary. In particular, we deal with domains with varying and possibly divergent exits to infinity and aperture domains. Copyright © 2007 John Wiley & Sons, Ltd.     

Sufficient conditions for the regularity of the solutions of the Navier–Stokes equations

In this paper we find sufficient conditions, involving only the pressure, that ensure the regularity of weak solutions to the Navier–Stokes equations. Conditions involving only the pressure were previously obtained in [7,4]. Following a remark in this last reference we improve, in particular, Kaniel's result [7]. Our condition can be seen at the light of the heuristic idea that the pressure behaves similarly to the modulus squared of the velocity. Copyright © 1999 John Wiley & Sons, Ltd.     

On the regularity criterion for the solutions of 3D Navier–Stokes equations in weak multiplier spaces

In this work, we improved the regularity criterion on the Cauchy problem for the Navier–Stokes equations in multiplier space in terms of the two partial derivatives of velocity fields, ?₁u₁ and ?₂u₂. Copyright © 2011 John Wiley & Sons, Ltd.     

Well‐posedness for fractional Navier–Stokes equations in the largest critical spaces

This note studies the well‐posedness of the fractional Navier–Stokes equations in some supercritical Besov spaces as well as in the largest critical spaces for β ∈ (1/2,1). Meanwhile, the well‐posedness for fractional magnetohydrodynamics equations in these Besov spaces is also studied. Copyright © 2012 John Wiley & Sons, Ltd.     

Two-level Galerkin–Lagrange multipliers method for the stationary Navier–Stokes equations

In this paper, we combine the Galerkin–Lagrange multiplier (GLM) method with the two-level method to solve the stationary Navier–Stokes equations in order to avoid the time-consuming process and the construction of zero-divergence elements. Different quadrilateral partitions are used for approximating the velocity and the pressure. Then some error estimates are obtained and some numerical results of the GLM method and the two-level GLM method are given. The results show that the two-level method based on the GLM method is more efficient than the GLM method under the convergence rate of same order.     

Regularity criteria for the Navier–Stokes equations involving the ratio pressure‐gradient of velocity

A number of bounds upon the pressure are known to guarantee regularity of the solutions of the Navier–Stokes equations. Since the pressure is the potential whose source is the product of the velocity and its gradient, it is worth to consider bounds depending on the quotient of the pressure and some quantity measuring the size of this source. Estimates involving the ratio pressure–velocity are known. Our result includes the velocity gradient: if the ratio remains bounded for some r<1, so does the velocity and therefore it retains its regularity. Copyright © 2009 John Wiley & Sons, Ltd.     

Strong solutions for the nonhomogeneous Navier–Stokes equations in unbounded domains

We show the existence of strong solutions for the nonhomogeneous Navier–Stokes equations in three‐dimensional domains with boundary uniformly of class C³. Under suitable assumptions, uniqueness is also proved. Copyright © 2009 John Wiley & Sons, Ltd.     

A two-grid method based on Newton iteration for the Navier–Stokes equations

In this paper, we consider a two-grid method for resolving the nonlinearity in finite element approximations of the equilibrium Navier–Stokes equations. We prove the convergence rate of the approximation obtained by this method. The two-grid method involves solving one small, nonlinear coarse mesh system and two linear problems on the fine mesh which have the same stiffness matrix with only different right-hand side. The algorithm we study produces an approximate solution with the optimal asymptotic in h and accuracy for any Reynolds number. Numerical example is given to show the convergence of the method.     

A fully discrete nonlinear Galerkin method for the 3D Navier–Stokes equations

The purpose of this paper is twofold: (i) We show that the Fourier‐based Nonlinear Galerkin Method (NLGM) constructs suitable weak solutions to the periodic Navier–Stokes equations in three space dimensions provided the large scale/small scale cutoff is appropriately chosen. (ii) If smoothness is assumed, NLGM always outperforms the Galerkin method by a factor equal to 1 in the convergence order of the H ¹‐norm for the velocity and the L²‐norm for the pressure. This is a purely linear superconvergence effect resulting from standard elliptic regularity and holds independently of the nature of the boundary conditions (whether periodicity or no‐slip BC is enforced). © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008