GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

1．引言本文的工作主要是讨论非定常的热传导一对流问题的向后一步的Euler全离散化的非线性Galerkin混合元解的存在性及其误差估计．该工作是对山中的同一问题研究的第二部分．在第一部分［1］，我们已经讨论了此问题的半离散化的情形．由于所研究的目标都是非定常的热传导一对流问题，其背景是相同的，在此将不重复了，请参考［1］．本文的安排如下，52先回顾非定常的热传导一对流问题的混合元解的经典性质．53回顾半离散化的非线性Galerkin混合元解的性质，并导出后续讨论需要的一些关于时间导数的估计．54讨论向后一步的Euler全离散化…     

非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式

In this paper, a fully discrete format of nonlinear Galerkin mixed element method with backward one-step Euler discretization of time for the non stationary conduction-convection problems is presented. The scheme is based on two finite element spaces XH and Xh for the approximation of the velocity, defined respectively on a coarse grid with grids size H and another fine grid with grid size h<< H, a finite element space Mh for the approximation of the pressure and two finite element spaces AH and Wh, for the approximation of the temperature,also defined respectivply on the coarse grid with grid size H and another fine grid with grid size h. The existence and the convergence of the fully discrete mixed element solution are shown. The scheme consists in using standard backward one step Euler-Galerkin fully discrete format at first L0 steps (L0 2) on fine grid with grid size h, but using nonlinear Galerkin mixed element method of backward one step Euler-Galerkin fully discrete format through L0 + 1 step to end step. We have proved that the fully discrete nonlinear Galerkin mixed element procedure with respect to the coarse grid spaces with grid size H holds superconvergence.     

Navier-Stokes方程的最佳有限元非线性Galerkin算法

1．引言对于Navier－Stokes方程有限元数值求解方面的研究已有很多的文章和专著,多数是采用有限元Galerkin算法,例见文献[1-4]．然而,由于Navier-Stokes方程在大雷诺数时有其强的非线性性和对时间土的长期依赖性,用计算机求解Navier－Stokes方程在速度和容量方面是难以承受的．为了克服这些困难,最近人们提出了有限元非线性Galerkin算法,见文献卜8］,然而这种算法只是在某一有限时刻之后具有好的收敛速度,在初始时刻的某一区间不能达到好的收敛速度．本文应用Taylor展开技术导出了数值求解二维非定常Navier－Stokes方程的最佳…     

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.     

Nonlinear Galerkin methods and mixed finite elements: two-grid algorithms for the Navier-Stokes equations

Summary. A nonlinear Galerkin method using mixed finite elements is presented for the two-dimensional incompressible Navier-Stokes equations. The scheme is based on two finite element spaces and for the approximation of the velocity, defined respectively on one coarse grid with grid size and one fine grid with grid size and one finite element space for the approximation of the pressure. Nonlinearity and time dependence are both treated on the coarse space. We prove that the difference between the new nonlinear Galerkin method and the standard Galerkin solution is of the order of $H^2$, both in velocity ( and pressure norm). We also discuss a penalized version of our algorithm which enjoys similar properties. Received October 5, 1993 / Revised version received November 29, 1993     

Numerical analysis of a modified finite element nonlinear Galerkin method

Summary. A fully discrete modified finite element nonlinear Galerkin method is presented for the two-dimensional equation of Navier-Stokes type. The spatial discretization is based on two finite element spaces X_H and X_h defined on a coarse grid with grid size H and a fine grid with grid size h << H, respectively; the time discretization is based on the Euler explicit scheme with respect to the nonlinear term. We analyze the stability and convergence rate of the method. Comparing with the standard finite element Galerkin method and the nonlinear Galerkin method, this method can admit a larger time step under the same convergence rate of same order. Hence this method can save a large amount of computational time. Finally, we provide some numerical tests on this method, the standard finite element Galerkin method, and the nonlinear Galerkin method, which are in a good agreement with the theoretical analysis.Mathematics Subject Classification (2000): 35Q30, 65M60, 65N30, 76D05     

Convergence and stability of finite element nonlinear Galerkin method for the Navier-Stokes equations

In this article, we present a new fully discrete finite element nonlinear Galerkin method, which are well suited to the long time integration of the Navier-Stokes equations. Spatial discretization is based on two-grid finite element technique; time discretization is based on Euler explicit scheme with variable time step size. Moreover, we analyse the boundedness, convergence and stability condition of the finite element nonlinear Galerkin method. Our discussion shows that the time step constraints of the method depend only on the coarse grid parameter and the time step constraints of the finite element Galerkin method depend on the fine grid parameter under the same convergence accuracy. Received February 2, 1994 / Revised version received December 6, 1996     

Two grid discretizations with backtracking of the stream function form of the Navier-Stokes equations

We analyze a two grid finite element method with backtracking for the stream function formulation of the stationary Navier—Stokes equations. This two grid method involves solving one small, nonlinear coarse mesh system, one linearized system on the fine mesh and one linear correction problem on the coarse mesh. The algorithm and error analysis are presented.     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

Comparison of generalized differential quadrature and Galerkin methods for the analysis of micro-electro-mechanical coupled systems

This paper presents a comprehensive comparison study between the generalized differential quadrature (GDQ) and the well-known global Galerkin method for analysis of pull-in behavior of nonlinear micro-electro-mechanical coupled systems. The nonlinear governing integro-differential equation for double clamped MEMS devices which was derived using variational principle by the authors [Sadeghian H, Rezazadeh G, Osterberg PM. Application of the generalized differential quadrature method to the study of pull-in phenomena of MEMS switches. J Microelectromech Syst 2007;16(6):1334–40] is discretized by applying Galerkin and GDQ methods. The divergence instability or pull-in phenomenon is analyzed. Obtained results are compared with the results of the pervious works. The Galerkin method is implemented with effect of number of used shape functions. Different types of trail functions on calculated pull-in voltage are examined.Furthermore, compare to one term and two terms truncation Galerkin method, it is observed that the GDQ with small number of grid points (non-uniform) performs accurate results for nonlinear micro-electro-mechanical coupled behavior which requires a large number of grid points at high-order approximation.     

Comparison of the finite difference and Galerkin methods as applied to the solution of the hydrodynamic equations

The three-dimensional nonlinear hydrodynamic equations which describe wind induced flow in a homogeneous sea are transformed from Cartesian coordinates into sigma coordinates. The solution of these equations in the horizontal is accomplished using a standard finite difference grid and established finite difference methods.The accuracy and computational efficiency, in terms of both computer time and main memory requirements, of using either the Galerkin method or a finite difference grid through the vertical is considered. Calculations, using the same number of functions in the Galerkin method as grid bases through the vertical shows that the Galerkin method has superior accuracy over the grid box method. Hence, for a given accuracy a smaller number of functions than grid boxes may be used, with associated saving in computational resources.For the case in which the vertical variation of eddy viscosity is fixed, an eigenvalue problem can be solved to yield a set of eigenfunctions. Using these eigenfunctions as a basis set with the Galerkin approach, a Galerkin-eigenfunction method is developed. Calculations show that the Galerkin-eigenfunction technique is accurate and in a linear model is clearly computationally more economic than the use of grid boxes through the vertical.     

外部非定常Navier-Stokes方程的非线性Galerkin算法

提出了求解外部非定常Navier-stokes方程的有限元边界元耦合的非线性Galerkin算法,证明了相应变分问题的正则性和数值解的收敛速度。收敛性分析表明如果选取粗网格尺度H是细网格尺度h的开平方数量级,则该算法提供了与古典Galerkin算法同阶的收敛速度。然而非线性Galerkin算法仅仅需要在粗网格解非线性问题,在细网格上解线性问题。因此,该算法可以节省计算工作量。     

加罚Navier—Stokes方程的最佳非线性Galerkin算法

该文提出了求解二维加罚Navier-Stokes方程的最佳非线性Galerkin算法．这个算法在于在粗网格有限元空间上求解一非线性子问题，在细网格增量有限元空间Wh上求解一线性子问题．如果线性有限元被使用及，则该算法具有和有限元Galerkin算法同阶的收敛速度．然而该文提出的算法可以节省可观的计算时间．     

Model reduction on inertial manifolds for N–S equations approached by multilevel finite element method

Approximate Inertial Manifolds (AIMs) is approached by multilevel finite element method, which can be referred to as a Post-processed nonlinear Galerkin finite element method, and is applied to the model reduction for fluid dynamics, a typical kind of nonlinear continuous dynamic system from viewpoint of nonlinear dynamics. By this method, each unknown variable, namely, velocity and pressure, is divided into two components, that is the large eddy and small eddy components. The interaction between large eddy and small eddy components, which is negligible if standard Galerkin algorithm is used to approach the original governing equations, is considered essentially by AIMs, and consequently a coarse grid finite element space and a fine grid incremental finite element space are introduced to approach the two components. As an example, the flow field of incompressible flows around airfoil is simulated numerically and discussed, and velocity and pressure distributions of the flow field are obtained accurately. The results show that there exists less essential degrees-of-freedom which can dominate the dynamic behaviors of the discretized system in comparison with the traditional methods, and large computing time can be saved by this efficient method. In a sense, the small eddy component can be captured by AIMs with fewer grids, and an accurate result can also be obtained.     

A two-grid method for finite volume element approximations of second-order nonlinear hyperbolic equations

The two-grid method is studied for solving a two-dimensional second-order nonlinear hyperbolic equation using finite volume element method. The method is based on two different finite element spaces defined on one coarse grid with grid size H and one fine grid with grid size h, respectively. The nonsymmetric and nonlinear iterations are only executed on the coarse grid and the fine grid solution can be obtained in a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. A prior error estimate in the H¹-norm is proved to be O(h+H³|lnH|) for the two-grid semidiscrete finite volume element method. With these proposed techniques, solving such a large class of second-order nonlinear hyperbolic equations will not be much more difficult than solving one single linearized equation. Finally, a numerical example is presented to validate the usefulness and efficiency of the method.     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS I： SPATIAL DISCRETIZATION

In this article we consider a two-level finite element Galerkin method using mixed finite elements for the two-dimensional nonstationary incompressible Navier-Stokes equations. The method yields a $H^1$-optimal velocity approximation and a $L_2$-optimal pressure approximation. The two-level finite element Galerkin method involves solving one small, nonlinear Navier-Stokes problem on the coarse mesh with mesh size $H$, one linear Stokes problem on the fine mesh with mesh size $h << H$. The algorithm we study produces an approximate solution with the optimal, asymptotic in $h$, accuracy.     

Decoupled two‐grid finite element method for the time‐dependent natural convection problem I: Spatial discretization

In this article, a decoupled two grid finite element method (FEM) is proposed and analyzed for the nonsteady natural convection problem using the coarse grid numerical solutions to decouple the nonlinear coupled terms, and the corresponding optimal error estimates are derived. Compared with the standard Galerkin FEM and the usual two‐grid FEM, our algorithm not only keeps good accuracy but also saves a lot of computational cost. Some numerical examples are provided to verify the performances of the decoupled two‐grid FEM. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the decoupled two‐grid FEM for the nonsteady natural convection problem. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2135–2168, 2015     

TWO-LEVEL METHOD FOR UNSTEADY NAVIER-STOKES EQUATIONS IN STREAM FUNCTION FORM

Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied. This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid. Moreover,the scaling between these two grid sizes is super-linear. Approximation,stability and convergence aspects of a fully discrete scheme are analyzed. At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.     

A Two-Level Method for Nonsymmetric Eigenvalue Problems

A two-level discretization method for eigenvalue problems is studied. Compared to the standard Galerkin finite element discretization technique performed on a fine grid this method discretizes the eigenvalue problem on a coarse grid and obtains an improved eigenvector (eigenvalue) approximation by solving only a linear problem on the fine grid (or two linear problems for the case of eigenvalue approximation of nonsymmetric problems). The improved solution has the asymptotic accuracy of the Galerkin discretization solution. The link between the method and the iterated Galerkin method is established. Error estimates for the general nonsymmetric case are derived.