A Nonlinear Galerkin/Petrov-Least Squares Mixed Element Method for the Stationary Navier-Stokes Equation

A nonlinear Galerkin/ Petrov- least squares mixed element (NGPLSME) method for the stationary Navier-Stokes equations is presented and analyzed. The scheme is that Petrov-least squares forms of residuals are added to the nonlinear Galerkin mixed element method so that it is stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. The existence, uniqueness and convergence ( at optimal rate ) of the NGPLSME solution is proved in the case of sufficient viscosity ( or small data).     

Adaptive mixed least squares Galerkin/Petrov finite element method for stationary conduction convection problems

An adaptive mixed least squares Galerkin/Petrov finite element method (FEM) is developed for stationary conduction convection problems. The mixed least squares Galerkin/Petrov FEM is consistent and stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. Using the general theory of Verfürth, the posteriori error estimates of the residual type are derived. Finally, numerical tests are presented to illustrate the effectiveness of the method.     

NONLINEAR GALERKIN MIXED ELEMENT METHODS FOR STATIONARY INCOMPRESSIBLE MAGNETOHYDRODYNAMICS

A nonlinear Galerkin mixed element (NGME) method for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of the NGME solution are derived.     

A NONLINEAR GALERKIN MIXED ELEMENT METHOD AND A POSTERIORI ERROR ESTIMATOR FOR THE STATIONARY NAVIER-STOKES EQUATIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｏｎｌｉｎｅａｒＧａｌｅｒｋｉｎｍｅｔｈｏｄｉｓａｍｕｌｔｉ＿ｌｅｖｅｌｓｃｈｅｍｅｔｏｆｉｎｄｔｈｅａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｆｏｒｔｈｅｄｉｓｓｉｐａｔｉｖｅＰＤＥ (ｐａｒｔｉａｌｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎ) .Ｔｈｉｓｍｅｔｈｏｄｃｏｎｓｉｓｔｓｉｎｓｐｌｉｔｔｉｎｇｔｈｅｕｎｋｎｏｗｎｉｎｔｏｔｗｏ (ｏｒｍｏｒｅ)ｔｅｒｍｓ ,ｗｈｉｃｈｂｅｌｏｎｇｔｏｔｈｅｄｉｓｃｒｅｔｅｓｐａｃｅｓｗｉｔｈｄｉｆｆｅｒｅｎｔｍｅｓｈｓｉｚｅ .Ｔｈｅ…     

线性定常对流占优对流扩散问题的无网格解法

应用无网格Galerkin方法求解对流占优对流扩散问题时会出现非物理现象的数值伪振荡,本文将SUPG方法、GLS方法、SGS方法与无网格Galerkin方法相耦合,成功解决了对流扩散方程中对流项占优时的数值伪振荡问题。运用本文构造的方法,采用线性基和具有C2连续的权函数,应用移动最小二乘法可容易地构造高阶导数连续的形函数,从而避免了有限元方法中当采用线性元插值时,因忽略稳定项中二阶导数项而降低计算精度和稳定性的问题。数值实验表明:本文构造的方法具有计算精度高、稳定性好、计算算法实施简单、前后处理方便的优点,这些方法不仅能适用于对流项占优问题,而且也能很好地消除反应项占优时的数值伪振荡问题。     

动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径.自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法.自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势.在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程.然后,采用预校正形式的Newmark法在时间域上进行求解.最后给出了数值算例,并验证了所提方法的有效性和正确性.     

Numerical simulation of natural and mixed convection flows by Galerkin‐characteristic method

A numerical investigation is performed to study the solution of natural and mixed convection flows by Galerkin‐characteristic method. The method is based on combining the modified method of characteristics with a Galerkin finite element discretization in primitive variables. It can be interpreted as a fractional step technique where convective part and Stokes/Boussinesq part are treated separately. The main feature of the proposed method is that, due to the Lagrangian treatment of convection, the Courant–Friedrichs–Lewy (CFL) restriction is relaxed and the time truncation errors are reduced in the Stokes/Boussinesq part. Numerical simulations are carried out for a natural convection in squared cavity and for a mixed convection flow past a circular cylinder. The computed results are compared with those obtained using other Eulerian‐based Galerkin finite element solvers, which are used for solving many convective flow models. The Galerkin‐characteristic method has been found to be feasible and satisfactory. Copyright © 2006 John Wiley & Sons, Ltd.     

反平面断裂问题的无单元伽辽金比例边界法

将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。     

Galerkin-Petrov least squares mixed element method for stationary incompressible magnetohydrodynamics

The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.     

An optimizing reduced PLSMFE formulation for non‐stationary conduction–convection problems

In this paper, proper orthogonal decomposition (POD) is combined with the Petrov–Galerkin least squares mixed finite element (PLSMFE) method to derive an optimizing reduced PLSMFE formulation for the non‐stationary conduction–convection problems. Error estimates between the optimizing reduced PLSMFE solutions based on POD and classical PLSMFE solutions are presented. The optimizing reduced PLSMFE formulation can circumvent the constraint of Babu?ka–Brezzi condition so that the combination of finite element subspaces can be chosen freely and allow optimal‐order error estimates to be obtained. Numerical simulation examples have shown that the errors between the optimizing reduced PLSMFE solutions and the classical PLSMFE solutions are consistent with theoretical results. Moreover, they have also shown the feasibility and efficiency of the POD method. Copyright © 2008 John Wiley & Sons, Ltd.     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

Nonlinear Galerkin method for the exterior nonstationary Navier-Stokes equations

ＩｎｔｒｏｄｕｃｔｉｏｎＬｅｔΩｃｏｎｔａｉｎｉｎｇｚｅｒｏｐｏｉｎｔｂｅａｓｉｍｐｌｙ＿ｃｏｎｎｅｃｔｅｄｂｏｕｎｄｅｄｏｐｅｎｓｅｔｏｆＲ2 ｗｉｔｈｓｍｏｏｔｈｂｏｕｎｄａｒｙΓａｎｄｌｅｔΩ′ｄｅｎｏｔｅｔｈｅｃｏｍｐｌｅｍｅｎｔｏｆΩ ∪Γ .ＴｈｅｅｘｔｅｒｉｏｒｎｏｎｓｔａｔｉｏｎａｒｙＮａｖｉｅｒ＿ＳｔｏｋｅｓｐｒｏｂｌｅｍｆｏｒａｆｌｕｉｄｏｃｃｕｐｙｉｎｇΩ′ｃｏｎｓｉｓｔｓｉｎｆｉｎｄｉｎｇｔｈｅｖｅｌｏｃｉｔｙ ｕ(ｘ,ｔ)ｏｆｔｈｅｆｌｕｉｄａｎｄｉｔｓｐｒｅｓｓｕｒｅ ｐ(ｘ ,…     

A comparative study of GLS finite elements with velocity and pressure equally interpolated for solving incompressible viscous flows

A comparative study of the bi‐linear and bi‐quadratic quadrilateral elements and the quadratic triangular element for solving incompressible viscous flows is presented. These elements make use of the stabilized finite element formulation of the Galerkin/least‐squares method to simulate the flows, with the pressure and velocity fields interpolated with equal orders. The tangent matrices are explicitly derived and the Newton–Raphson algorithm is employed to solve the resulting nonlinear equations. The numerical solutions of the classical lid‐driven cavity flow problem are obtained for Reynolds numbers between 1000 and 20 000 and the accuracy and converging rate of the different elements are compared. The influence on the numerical solution of the least square of incompressible condition is also studied. The numerical example shows that the quadratic triangular element exhibits a better compromise between accuracy and converging rate than the other two elements. Copyright © 2008 John Wiley & Sons, Ltd.     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

Augmented mixed finite element methods for a vorticity‐based velocity–pressure–stress formulation of the Stokes problem in 2D

In this paper, we consider an augmented velocity–pressure–stress formulation of the 2D Stokes problem, in which the stress is defined in terms of the vorticity and the pressure, and then we introduce and analyze stable mixed finite element methods to solve the associated Galerkin scheme. In this way, we further extend similar procedures applied recently to linear elasticity and to other mixed formulations for incompressible fluid flows. Indeed, our approach is based on the introduction of the Galerkin least‐squares‐type terms arising from the corresponding constitutive and equilibrium equations, and from the Dirichlet boundary condition for the velocity, all of them multiplied by stabilization parameters. Then, we show that these parameters can be suitably chosen so that the resulting operator equation induces a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well posed for any choice of finite element subspaces. In particular, we can use continuous piecewise linear velocities, piecewise constant pressures, and rotated Raviart–Thomas elements for the stresses. Next, we derive reliable and efficient residual‐based a posteriori error estimators for the augmented mixed finite element schemes. In addition, several numerical experiments illustrating the performance of the augmented mixed finite element methods, confirming the properties of the a posteriori estimators, and showing the behavior of the associated adaptive algorithms are reported. The present work should be considered as a first step aiming finally to derive augmented mixed finite element methods for vorticity‐based formulations of the 3D Stokes problem. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of piezoelectric ceramic multilayer actuators based on an electro-mechanical coupled meshless method

This paper presents an efficient meshless method for analyzing cracked piezoelectric structures subjected to mechanical and electrical loading. In this method, an element free Galerkin (EFG) formulation, an enriched basic function and some special shape functions that contain discontinuous derivatives are employed. Based on the moving least squares (MLS) interpolation approach, the EFG method is one of the promising methods for dealing with problems involving progressive crack growth. Since the method is meshless and no element connectivity data are needed, the burdensome remeshing procedure required in the conventional finite element method (FEM) is avoided. The numerical results show that the proposed method can yield an accurate near-tip stress field in an infinite piezoelectric plate containing an interior hole. In another example studying a ceramic multilayer actuator, the proposed model was found to be accurate in the simulation of stress and electric field concentrations arround the abrupt end of an internal electrode. The project supported by the National Natural Science Foundation of China (10025209, 10132010, and 90208002), and the Research Grants Council of the Hong Kong Special Administrative Region, China (HKU 7203/03E). The English text was polished by Yunming Chen.     

极限下限分析的正交基无单元Galerkin法

基于极限分析的下限定理,建立了用正交基无单元Galerkin法进行理想弹塑性结构极 限分析的整套求解算法．下限分析所需的虚拟弹性应力场可由正交基无单元Galerkin法直接 得到,所需的自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模 拟．这些自平衡应力场基矢量可由弹塑性增量分析中的平衡迭代得到．通过对自平衡应力场 子空间的不断修正,整个问题的求解将化为一系列非线性数学规划子问题,并通过复合形法 进行求解．算例表明该方法有效地克服了维数障碍问题,使计算效率得到了充分的提高,是 切实可行的．     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

Parallel domain decomposition method for finite element approximation of 3D steady state non‐Newtonian fluids

We introduce a stabilized finite element method for the 3D non‐Newtonian Navier–Stokes equations and a parallel domain decomposition method for solving the sparse system of nonlinear equations arising from the discretization. Non‐Newtonian flow problems are, generally speaking, more challenging than Newtonian flows because the nonlinearities are not only in the convection term but also in the viscosity term, which depends on the shear rate. Many good iterative methods and preconditioning techniques that work well for the Newtonian flows do not work well for the non‐Newtonian flows. We employ a Galerkin/least squares finite element method, with stabilization parameters adjusted to count the non‐Newtonian effect, to discretize the equations, and the resulting highly nonlinear system of equations is solved by a Newton–Krylov–Schwarz algorithm. In this study, we apply the proposed method to some inelastic power‐law fluid flows through the eccentric annuli with inner cylinder rotation and investigate the robustness of the method with respect to some physical parameters, including the power‐law index and the Reynolds number ratios. We then report the superlinear speedup achieved by the domain decomposition algorithm on a computer with up to 512 processors. Copyright © 2015 John Wiley & Sons, Ltd.