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.     

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.     

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

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

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.     

Three dimensional finite element analysis for Rayleight-Benard convection in rectangular enclosures

The pattern of Rayleigh-Benard convection of air in a rectangular box heated-from-below is studied by numerically solving the three-dimensional time-dependent Navier-Stokes equations under the Boussinesq approximation. Slightly supercritical Rayleigh number was adopted to track the evolutions of flow structure as a function of enclosure's aspect ratio (A=L/H). The flow will asymptotically evolve to different patterns, among which, two possible types of flow pattern are found. One consists of the pair of straight vortex rolls and the other appears as closed vortex rings. The transition between the flow patterns indicates that there exists a flow bifurcation with the variation of container's aspect ratio. In addition, both steady and oscillatory flows have been observed, corresponding to the pair of straight vortex rolls and the vortex ring, respectively. The complexity of flow structure tends to increase with the increasing aspect ratio of the rectangular enclosure. The project supported by the National Natural Science Foundation of China (19889210), the National Distinguished Young Fund (10125210), the Hundred Talents Program of CAS, and the Training Program for the Trans-Century Outstanding Young of MOE     

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).     

Introduction of turbulent model in a mixed finite volume/finite element method

The aim of this work is to present a new numerical method to compute turbulent flows in complex configurations. With this in view, a k-? model with wall functions has been introduced in a mixed finite volume/finite element method. The numerical method has been developed to deal with compressible flows but is also able to compute nearly incompressible flows. The physical model and the numerical method are first described, then validation results for an incompressible flow over a backward-facing step and for a supersonic flow over a compression ramp are presented. Comparisons are performed with experimental data and with other numerical results. These simulations show the ability of the present method to predict turbulent flows, and this method will be applied to simulate complex industrial flows (flow inside the combustion chamber of gas turbine engines). The main goal of this paper is not to test turbulence models, but to show that this numerical method is a solid base to introduce more sophisticated turbulence model.     

Finite cell method compared to h-version finite element method for elasto-plastic problems

The finite cell method （FCM） combines the high-order finite element method （FEM） with the fictitious domain approach for the purpose of simple meshing. In the present study, the FCM is used to the Prandtl-Reuss flow theory of plasticity, and the results are compared with the h-version finite element method （h-FEM）. The numerical results show that the FCM is more efficient compared to the h-FEM for elasto-plastic problems, although the mesh does not conform to the boundary. It is also demonstrated that the FCM performs well for elasto-plastic loading and unloading.     

p-version least squares finite element formulation for two-dimensional,incompressible fluid flow

A p-version least squares finite element formulation for non-linear problems is applied to the problem of steady, two-dimensional, incompressible fluid flow. The Navier-Stokes equations are cast as a set of first-order equations involving viscous stresses as auxiliary variables. Both the primary and auxiliary variables are interpolated using equal-order C⁰ continuity, p-version hierarchical approximation functions. The least squares functional (or error functional) is constructed using the system of coupled first-order non-linear partial differential equations without linearization, approximations or assumptions. The minimization of this least squares error functional results in finding a solution vector {δ} for which the partial derivative of the error functional (integrated sum of squares of the errors resulting from individual equations for the entire discretization) with respect to the nodal degrees of freedom {δ} becomes zero. This is accomplished by using Newton's method with a line search. Numerical examples are presented to demonstrate the convergence characteristics and accuracy of the method.     

A finite element method for the one‐dimensional extended Boussinesq equations

A new finite element method for Nwogu's (O. Nwogu, ASCE J. Waterw., Port, Coast., Ocean Eng., 119 , 618–638 (1993)) one‐dimensional extended Boussinesq equations is presented using a linear element spatial discretisation method coupled with a sophisticated adaptive time integration package. The accuracy of the scheme is compared to that of an existing finite difference method (G. Wei and J.T. Kirby, ASCE J. Waterw., Port, Coast., Ocean Eng., 121 , 251–261 (1995)) by considering the truncation error at a node. Numerical tests with solitary and regular waves propagating in variable depth environments are compared with theoretical and experimental data. The accuracy of the results confirms the analytical prediction and shows that the new approach competes well with existing finite difference methods. The finite element formulation is shown to enable the method to be extended to irregular meshes in one dimension and has the potential to allow for extension to the important practical case of unstructured triangular meshes in two dimensions. This latter case is discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

The three-dimensional prolonged adaptive unstructured finite element multigrid method for the Navier–Stokes equations

This paper presents the development of the three- dimensional prolonged adaptive finite element equation solver for the Navier–Stokes equations. The finite element used is the tetrahedron with quadratic approximation of the velocities and linear approximation of the pressure. The equation system is formulated in the basic variables. The grid is adapted to the solution by the element Reynolds number. An element in the grid is refined when the Reynolds number of the element exceeds a preset limit. The global Reynolds number in the investigation is increased by scaling the solution for a lower Reynolds number. The grid is refined according to the scaled solution and the prolonged solution for the lower Reynolds number constitutes the start vector for the higher Reynolds number. Since the Reynolds number is the ratio of convection to diffusion, the grid refinements act as linearization and symmetrization of the equation system. The linear equation system of the Newton formulation is solved by CGSTAB with coupled node fill-in preconditioner. The test problem considered is the three-dimensional driven cavity flow. © 1997 John Wiley & Sons, Ltd.     

Multiphase flow in heterogeneous porous media: A classical finite element method versus an implicit pressure–explicit saturation‐based mixed finite element–finite volume approach

Various discretization methods exist for the numerical simulation of multiphase flow in porous media. In this paper, two methods are introduced and analyzed—a full‐upwind Galerkin method which belongs to the classical finite element methods, and a mixed‐hybrid finite element method based on an implicit pressure–explicit saturation (IMPES) approach. Both methods are derived from the governing equations of two‐phase flow. Their discretization concepts are compared in detail. Their efficiency is discussed using several examples. Copyright © 1999 John Wiley & Sons, Ltd.     

A finite element model for numerical simulation of thermo-mechanical frictional contact problems

Two kinds of variational principles for numerical simulation of heat transfer and contact analysis are respectively presented. A finite element model for numerical simulation of the thermal contact problems is developed with a pressure dependent heat transfer constitutive model across the contact surface. The numerical algorithm for the finite element analysis of the thermomechanical contact problems is thus developed. Numerical examples are computed and the results demonstrate the validity of the model and algorithm developed. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (50178016, 10225212) and the Foundation for University Key Teacher by the Ministry of Education of China     

P-version hybrid analytical finite element method for plate bending problems

Two sets of trial functions with different variables are constructed for the admissible space of the finite element analysis. The trial functions satisfy the equilibrium differential equation inside elements, while the deflections and rotations on the edges of the elements are approximated by the Peano hierarchical interpolation functions. Then, a generalized variational principle is applied to set up the p-version hybrid analytical finite element method for plate bending problems. The accuracy of finite element computation can be improved by increasing the order of the interpolation polynomials with fixed mesh. In the finite element formulation, to obtain the stiffness matrices and the load vectors, it is only necessary to perform quadrature over the edges of the elements. These matrices and vectors possess an embedding structure. The conformability between the elements can be controlled automatically.This work is supported by the Natural Science Foundation of China and the Aeronautical Science Foundation of China.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

自然单元法的自适应研究

自然单元法(NEM)是一种新兴的无网格数值计算方法,具有前处理简单和易于准确施加本质边界条件等优点.本文基于Z-Z后验误差估计方法,给出了一种自然单元法的误差估计因子和自适应分析细化方案.采用结点处的光滑应变计算相应的恢复应力,并用于构造全域上的恢复应力场.通过结点Voronoi单胞内的能量范数相对误差指示需要进行结点...     

Discrete formulation of mixed finite element methods for vapor deposition chemical reaction equations

The vapor deposition chemical reaction processes, which are of extremely extensive applications, can be classified as a mathematical model by the following governing nonlinear partial differential equations containing velocity vector, temperature field, pressure field, and gas mass field. The mixed finite element (MFE) method is employed to study the system of equations for the vapor deposition chemical reaction processes. The semidiscrete and fully discrete MFE formulations are derived. And the existence and convergence (error estimate) of the semidiscrete and fully discrete MFE solutions are demonstrated. By employing MFE method to treat the system of equations for the vapor deposition chemical reaction processes, the numerical solutions of the velocity vector, the temperature field, the pressure field, and the gas mass field can be found out simultaneously. Thus, these researches are not only of important theoretical means, but also of extremely extensive applied vistas.     

p-version least squares finite element formulation for two-dimensional,incompressible, non-Newtonian isothermal and non-isothermal fluid flow

This paper presents a p- version least squares finite element formulation (LSFEF) for two-dimensional, incompressible, non-Newtonian fluid flow under isothermal and non-isothermal conditions. The dimensionless forms of the diffential equations describing the fluid motion and heat transfer are cast into a set of first-order differential equations using non-Newtonian stresses and heat fluxes as auxiliary variables. The velocities, pressure and temperature as well as the stresses and heat fluxes are interpolated using equal-order, C⁰-continuous, p-version hierarchical approximation functions. The application of least squares minimization to the set of coupled first-order non-linear partial differential equations results in finding a solution vector {δ} which makes the partial derivatives of the error functional with respect to {δ} a null vector. This is accomplished by using Newton's method with a line search. The paper presents the implementation of a power-law model for the non-Newtonian Viscosity. For the non-isothermal case the fluid properties are considered to be a function of temperature. Three numerical examples (fully developed flow between parallel plates, symmetric sudden expansion and lid-driven cavity) are presented for isothermal power-law fluid flow. The Couette shear flow problem and the 4:1 symmetric sudden expansion are used to present numerical results for non-isothermal power-law fluid flow. The numerical examples demonstrate the convergence characteristics and accuracy of the formulation.