Least‐squares spectral collocation method for the Stokes equations

First‐order system least‐squares spectral collocation methods are presented for the Stokes equations by adopting the first‐order system and modifying the least‐squares functionals in 2 . Then homogeneous Legendre and Chebyshev (continuous and discrete) functionals are shown to be elliptic and continuous with respect to appropriate product weighted norms. The spectral convergence is analyzed for the proposed methods with some numerical experiments. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 128–139, 2004     

Negative norm least‐squares methods for the velocity‐vorticity‐pressure Navier–Stokes equations

We develop and analyze a least‐squares finite element method for the steady state, incompressible Navier–Stokes equations, written as a first‐order system involving vorticity as new dependent variable. In contrast to standard L² least‐squares methods for this system, our approach utilizes discrete negative norms in the least‐squares functional. This allows us to devise efficient preconditioners for the discrete equations, and to establish optimal error estimates under relaxed regularity assumptions. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 237–256, 1999     

Least‐squares mixed finite element methods for the incompressible Navier‐Stokes equations

Least‐squares mixed finite element schemes are formulated to solve the evolutionary Navier‐Stokes equations and the convergence is analyzed. We recast the Navier‐Stokes equations as a first‐order system by introducing a vorticity flux variable, and show that a least‐squares principle based on L² norms applied to this system yields optimal discretization error estimates. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 441–453, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10015     

Analysis and computation of least‐squares methods for a compressible Stokes problem

We develop and analyze a negative norm least‐squares method for the compressible Stokes equations with an inflow boundary condition. Least‐squares principles are derived for a first‐order form of the equations obtained by using ω = ?×u and φ = ? · u as new dependent variables. The resulting problem is incompletely elliptic, i.e., it combines features of elliptic and hyperbolic equations. As a result, well‐posedness of least‐squares functionals cannot be established using the ADN elliptic theory and so we use direct approaches to prove their norm‐equivalence. The article concludes with numerical examples that illustrate the theoretical convergence rates. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

An analysis for the compressible Stokes equations by first‐order system of least‐squares finite element method

This article applies the first‐order system least‐squares (fosls) finite element method developed by Cai, Manteuffel and McCormick to the compressible Stokes equations. By introducing a new dependent velocity flux variable, we recast the compressible Stokes equations as a first‐order system. Then it is shown that the ellipticity and continuity hold for the least‐squares functionals employing the mixture of H^?1 and L², so that the fosls finite element methods yield best approximations for the velocity flux and velocity. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:689–699, 2001     

Stability and error analysis for a spectral Galerkin method for the Navier‐Stokes equations with H2 or H1 initial data

In this article we consider a spectral Galerkin method with a semi‐implicit Euler scheme for the two‐dimensional Navier‐Stokes equations with H² or H¹ initial data. The H²‐stability analysis of this spectral Galerkin method shows that for the smooth initial data the semi‐implicit Euler scheme admits a large time step. The L²‐error analysis of the spectral Galerkin method shows that for the smoother initial data the numerical solution u exhibits faster convergence on the time interval [0, 1] and retains the same convergence rate on the time interval [1, ∞). © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Least‐squares method for the Oseen equation

This article studies the least‐squares finite element method for the linearized, stationary Navier–Stokes equation based on the stress‐velocity‐pressure formulation in d dimensions (d = 2 or 3). The least‐squares functional is simply defined as the sum of the squares of the L² norm of the residuals. It is shown that the homogeneous least‐squares functional is elliptic and continuous in the norm. This immediately implies that the a priori error estimate of the conforming least‐squares finite element approximation is optimal in the energy norm. The L² norm error estimate for the velocity is also established through a refined duality argument. Moreover, when the right‐hand side f belongs only to , we derive an a priori error bound in a weaker norm, that is, the norm. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1289–1303, 2016     

Least‐squares mixed finite element methods for the RLW equations

A least‐squares mixed finite element (LSMFE) schemes are formulated to solve the 1D regularized long wave (RLW) equations and the convergence is discussed. The L² error estimates of LSMFE methods for RLW equations under the standard regularity assumption on the finite element partition are given.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Spectral method for vorticity‐stream function form of Navier–Stokes equations with slip boundary conditions

In this paper, we propose a spectral method for the vorticity‐stream function form of the Navier–Stokes equations with slip boundary conditions. The numerical solutions fulfill the incompressibility and the physical boundary conditions automatically. The stability and convergence of the proposed methods are proven. Numeric results demonstrate the efficiency of suggested algorithm. Copyright © 2011 John Wiley & Sons, Ltd.     

Enhanced algorithms for implementing the spectral discretization of the vorticity‐velocity‐pressure formulation of Stokes problem

In this paper, we deal with the numerical resolution of spectral discretization of the vorticity‐velocity‐pressure formulation of Stokes problem in a square or a cube provided with nonstandard boundary conditions, which involve the normal component of the velocity and the tangential components of the vorticity. Therefore, we propose two algorithms: the Uzawa algorithm and the global resolution. We implemented the two algorithms and compared their results. With global resolution, we obtained a very good accuracy with a small number of iteration.     

A two‐level method in space and time for the Navier‐Stokes equations

A two‐level method in space and time for the time‐dependent Navier‐Stokes equations is considered in this article. The approximate solution u_M∈H_M is decomposed into the large eddy component v∈H_m(m < M) and the small eddy component w∈H. We obtain the large eddy component v by solving a standard Galerkin equation in a coarse‐level subspace H_m with a time step length k, whereas the small eddy component w is derived by solving a linear equation in an orthogonal complement subspace H with a time step length pk, where p is a positive integer. The analysis shows that our two‐level scheme has long‐time stability and can reach the same accuracy as the standard Galerkin method in fine‐level subspace H_M for an appropriate configuration of p and m. Moreover, some numerical examples are provided to complement our theoretical analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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     

Spectral discretization of Darcy equations coupled with Stokes equations by vorticity–velocity–pressure formulations

We propose to make the numerical analysis of a model coupling the Darcy equations in a porous medium with the Stokes equations in the cracks. The coupling is provided by a pressure continuity on the interface. We describe a discretization by spectral element methods. We derive a priori optimal error estimates and we present some numerical experiments which confirm the results of the analysis.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1628–1651, 2017     

An optimal nonlinear Galerkin method with mixed finite elements for the steady Navier‐Stokes equations

An optimal nonlinear Galerkin method with mixed finite elements is developed for solving the two‐dimensional steady incompressible Navier‐Stokes equations. This method is based on two finite element spaces X_H and X_h for the approximation of velocity, defined on a coarse grid with grid size H and a fine grid with grid size h ? H, respectively, and a finite element space M_h for the approximation of pressure. We prove that the difference in appropriate norms between the solutions of the nonlinear Galerkin method and a classical Galerkin method is of the order of H⁵. If we choose H = O(h^2/5), these two methods have a convergence rate of the same order. We numerically demonstrate that the optimal nonlinear Galerkin method is efficient and can save a large amount of computational time. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 762–775, 2003.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

H least-squares method for the velocity–pressure–stress formulation of Stokes equations

This paper studies the discrete H⁻¹-norm least-squares method for the incompressible Stokes equations based on the velocity–pressure–stress formulation by the least-squares functional defined as the sum of L²-norms and H⁻¹-norm of the residual equations. Some computational experiments by multigrid method and preconditioning conjugate gradient method (PCGM) on this method are shown by taking efficient and β in the discrete solution operator T_h=h²I+βB_h corresponding to the minus one norm. We also propose a new method and compare it with PCGM and multigrid method through the analysis of numerical experiments depending on the choice of β.     

A new defect‐correction method for the stationary Navier‐Stokes equations based on pressure projection

A new defect‐correction method based on the pressure projection for the stationary Navier‐Stokes equations is proposed in this paper. A local stabilized technique based on the pressure projection is used in both defect step and correction step. The stability and convergence of this new method is analyzed detailedly. Finally, numerical examples confirm our theory analysis and validate high efficiency and good stability of this new method.     

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     

A two‐grid stabilization method for solving the steady‐state Navier‐Stokes equations

We formulate a subgrid eddy viscosity method for solving the steady‐state incompressible flow problem. The eddy viscosity does not act on the large flow structures. Optimal error estimates are obtained for velocity and pressure. The numerical illustrations agree completely with the theoretical results. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005