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.     

On general form of Navier-Stokes equations and implicit factored scheme

A general weak conservative form of Navier-Stokes equations expressed with respect to non-orthogonal Curvilinear coordinates and with primitive variables was obtained by using tensor analysis technique, where the contravariant and covariant velocity components were employed. Compared with the current coordinate transformation method, the established equations are concise and forthright, and they are more convenient to be used for solving problems in body-fitted curvilinear coordinate system. An implicit factored scheme for solving the equations is presented with detailed discussions in this paper. For n-dimensional flow the algorithm requires n-steps and for each step only a block tridiagonal matrix equation needs to be solved. It avoids inverting the matrix for large systems of equations and enhances the speed of arithmetic. In this study, the Beam- Warming’s implicit factored schceme is extended and developed in non-orthogonal curvilinear coordinate system.     

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.     

A nonlinear Galerkin/Petrov-least squares mixed element method for the stationary Navier-Stokes equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｏｎｌｉｎｅａｒＧａｌｅｒｋｉｎｍｅｔｈｏｄｉｓａｍｕｌｔｉ＿ｌｅｖｅｌｓｃｈｅｍｅｔｏｆｉｎｄｔｈｅａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｆｏｒｔｈｅｄｉｓｓｉｐａｔｉｖｅＰＤＥ .ＴｈｉｓｍｅｔｈｏｄｈａｓｆｉｒｓｔｍａｉｎｌｙｂｅｅｎａｄｄｒｅｓｓｅｄｂｙＦｏｉａｓ＿Ｍａｎｌｅｙ＿Ｔｅｍａｍ[1],Ｍａｒｉｏｎ＿Ｔｅｍａｍ[2 ],Ｆｏｉａｓ＿Ｊｏｌｌｙ＿Ｋｅｖｒｅｋｉｄｉｓ＿Ｔｉｔｉ[3]ａｎｄＤｅｖｕｌｄｅｒ＿Ｍａｒｉｏｎ＿Ｔｉｔｉ[4 ]ｉｎｔｈｅｃａｓｅｏｆｓｐｅｃｔ…     

Streamline upwind/full Galerkin method for solution of convection dominated solute transport problems

A new streamline method for the solution of convection-dominated transport problems is introduced. First, an analysis is made of the nature of classic streamline upwinding from the perspective of space and spacetime solution domains, as they pertain to the nature of the problems. From this analysis emerges a rigorous logic for upwinding, which can easily be implemented in the full (spacetime) Galerkin formulation of the transport equation. Comparative performance testing of this technique, in solving a number of examples, proves it to be robust and versatile. The advantage of this method resides in its applicability to a wider range of Courant numbers.     

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

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

Difference scheme and numerical simulation based on mixed finite element method for natural convection problem

ntroductionLetΩ R2 beaboundeddomain .Weconsiderthefollowingnon_stationarynaturalconvectionproblem :Problem (Ⅰ )　Findu =(u1,u2 ) ,p ,andTsuchthat,foranyt1>0 ,ut- μΔu +(u· )u + p=λjT　　 ((x ,y ,t) ∈Ω× (0 ,t1) ) ,divu =0　　　　　　　　　 ((x ,y,t) ∈Ω× (0 ,t1) ) ,Tt-ΔT +λu· T =0　　 ((x,y,t) ∈Ω× (0 ,t1) ) ,u =0 ,T =0　　　　　　 ((x,y,t)∈ Ω× (0 ,t1) ) ,u(x ,y ,0 ) =0 ,　T(x,y,0 ) =f(x,y)　　 ((x,y) ∈Ω) ,whereuisthefluidvelocityvectorfield ,pthepressurefield ,Tthet…     

A mixed finite element method for the generalized Stokes problem

We present and analyse a new mixed finite element method for the generalized Stokes problem. The approach, which is a natural extension of a previous procedure applied to quasi‐Newtonian Stokes flows, is based on the introduction of the flux and the tensor gradient of the velocity as further unknowns. This yields a two‐fold saddle point operator equation as the resulting variational formulation. Then, applying a slight generalization of the well known Babu?ka–Brezzi theory, we prove that the continuous and discrete formulations are well posed, and derive the associated a priori error analysis. In particular, the finite element subspaces providing stability coincide with those employed for the usual Stokes flows except for one of them that needs to be suitably enriched. We also develop an a posteriori error estimate (based on local problems) and propose the associated adaptive algorithm to compute the finite element solutions. Several numerical results illustrate the performance of the method and its capability to localize boundary layers, inner layers, and singularities. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.     

不同材料契合弹性力学问题的虚边界元最小二乘法

本文针对不同材料契合弹性力学问题,由虚边界元方法出发,建立了引入拉氏乘子的最小二科解法,对该类问题避免了采用Ｈｅｔｅｎｙｉ’ｓ基本解的麻烦和限制。数值结果表明本文算法的有效性和计算精度高。     

Least-squares mixed finite element method for a class of stokes equation

ＩｎｔｒｏｄｕｃｔｉｏｎＡｇｅｎｅｒａｌｔｈｅｏｒｙｏｆｔｈｅｌｅａｓｔ＿ｓｑｕａｒｅｓｍｅｔｈｏｄｈａｓｂｅｅｎｄｅｖｅｌｏｐｅｄｂｙＡＫＡｚｉｚ,ＲＢＫｅｌｌｏｇｇａｎｄＡＢＳｔｅｐｈｅｎｓｉｎ[1].Ｔｈｅｍｏｓｔｉｍｐｏｒｔａｎｔａｄｖａｎｔａｇｅｌｅａｄｓｔｏａｓｙｍｍｅｔｒｉｃｐｏｓｉｔｉｖｅｄｅｆｉｎｉｔｅｐｒｏｂｌｅｍ.ＪＨＢｒａｍｂｌｅａｎｄＪＡＮｉｔｓｈｅｐｒｅｓｅｎｔｅｄａｌｅａｓｔ＿ｓｑｕａｒｅｓｍｅｔｈｏｄｆｏｒＤｉｒｉｃｈｌｅｔｐｒｏｂｌｅｍｓｉｎ[2].Ｔｈｅｍｅｔｈｏｄｇｅ…     

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.     

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.     

A least squares finite element method for viscoelastic fluid flow problems

Viscoelastic flows remain a demanding class of problems for approximate analysis, particularly at increasing Weissenberg numbers. Part of the difficulty stems from the convective behavior and in the treatment of the stress field as a primary unknown. This latter aspect has led to the use of higher-order piecewise approximations for the stress approximation spaces in recent finite element research. The computational complexity of the discretized problem is increased significantly by this approach but at present it appears the most viable technique for solving these problems. Motivated by recent success in treating mixed systems and convective problems, we formulate here a least squares finite element method for the viscoelastic flow problem. Numerical experiments are conducted to test the method and examine its strengths and limitations. Some difficulties and open issues are identified through the numerical experiments. We consider the use of high degree elements (p refinement) to improve performance and accuracy.     

NEW ALGORITHM OF COUPLING ELEMENT-FREE GALERKIN WITH FINITE ELEMENT METHOD

IntroductionMeshless methods, as a special numerical method, originated from1970s. Since thediffuse element method was proposed by Nayroleset al.[1]in1992, the meshless methodshave received wide attentions in the mechanics area, and have shown obvious adv…     

流体饱和两相多孔介质拟静态问题的混合有限元方法

针对基于混合物理论的两相多孔介质模型,采用Galerkin加权残值有限元法,导出求解所静态问题的基于us-uF-P变量的混合有限元方程,由于系统方程的系数矩阵非定,进而针对该方程组提出了一种失代求解方法,并由分片试验得出节点压力插值函数的阶须低于固体相节点的位移插值函数的阶的结论,算例结果表明,采用基于u2-uF-p变量的混合法计算所得的固体相和流体相速度以及固体相的有效应力与罚方法一致,而压力值的粗度高于罚方法。     

偶应力问题的杂交／混合元分析

将弹性力学中Hellinger—Reissner交分原理推广到偶应力理论中，并以罚函数的形式引入其约束条件，提出了一种有效的杂交／混合单元。文中分别分析了带中心小孔平板在轴向均匀加载时的应力集中情况，以及含中问裂纹的无限平板单轴拉伸时的位移场和应力场。算例表明，该单元计算效率高，精度好，即使在材料本征长度很小时，仍然能够得到相当理想的结果。