Shape-free polygonal hybrid displacement-function element method for analyses of Mindlin–Reissner plates

A high-performance shape-free polygonal hybrid displacement-function finite-element method is proposed for analyses of Mindlin–Reissner plates. The analytical solutions of displacement functions are employed to construct element resultant fields, and the three-node Timoshenko’s beam formulae are adopted to simulate the boundary displacements. Then, the element stiffness matrix is obtained by the modified principle of minimum complementary energy. With a simple division, the integration of all the necessary matrices can be performed within polygonal element region. Five new polygonal plate elements containing a mid-side node on each element edge are developed, in which element HDF-PE is for general case, while the other four, HDF-PE-SS1, HDF-PE-Free, IHDF-PE-SS1, and IHDF-PE-Free, are for the edge effects at different boundary types. Furthermore, the shapes of these new elements are quite free, i.e., there is almost no limitation on the element shape and the number of element sides. Numerical examples show that the new elements are insensitive to mesh distortions, possess excellent and much better performance and flexibility in dealing with challenging problems with edge effects, complicated loading, and material distributions.

     

A two-level finite element method for the Allen–Cahn equation

We consider the fully implicit treatment for the nonlinear term of the Allen–Cahn equation. To solve the nonlinear problem efficiently, the two-level scheme is employed. We obtain the discrete energy law of the fully implicit scheme and two-level scheme with finite element method. Also, the convergence of the two-level method is presented. Finally, some numerical experiments are provided to confirm the theoretical analysis.     

A mixed–primal finite element method for the Boussinesq problem with temperature-dependent viscosity

In this paper we focus on the analysis of a mixed finite element method for a class of natural convection problems in two dimensions. More precisely, we consider a system based on the coupling of the steady-state equations of momentum (Navier–Stokes) and thermal energy by means of the Boussinesq approximation (coined the Boussinesq problem), where we also take into account a temperature dependence of the viscosity of the fluid. The construction of this finite element method begins with the introduction of the pseudostress and vorticity tensors, and a mixed formulation for the momentum equations, which is augmented with Galerkin-type terms, in order to deal with the non-linearity of these equations and the convective term in the energy equation, where a primal formulation is considered. The prescribed temperature on the boundary becomes an essential condition, which is weakly imposed, leading us to the definition of the normal heat flux through the boundary as a Lagrange multiplier. We show that this highly coupled problem can be uncoupled and analysed as a fixed-point problem, where Banach and Brouwer theorems will help us to provide sufficient conditions to ensure well-posedness of the problems arising from the continuous and discrete formulations, along with several applications of continuous injections guaranteed by the Rellich–Kondrachov theorem. Finally, we show some numerical results to illustrate the performance of this finite element method, as well as to prove the associated rates of convergence.     

An edge-based smoothed finite element method (ES-FEM) with stabilized discrete shear gap technique for analysis of Reissner–Mindlin plates

An edge-based smoothed finite element method (ES-FEM) for static, free vibration and buckling analyses of Reissner–Mindlin plates using 3-node triangular elements is studied in this paper. The calculation of the system stiffness matrix is performed by using the strain smoothing technique over the smoothing domains associated with edges of elements. In order to avoid the transverse shear locking and to improve the accuracy of the present formulation, the ES-FEM is incorporated with the discrete shear gap (DSG) method together with a stabilization technique to give a so-called edge-based smoothed stabilized discrete shear gap method (ES-DSG). The numerical examples demonstrated that the present ES-DSG method is free of shear locking and achieves the high accuracy compared to the exact solutions and others existing elements in the literature.     

A combined scheme of edge-based and node-based smoothed finite element methods for Reissner–Mindlin flat shells

In this paper, a combined scheme of edge-based smoothed finite element method (ES-FEM) and node-based smoothed finite element method (NS-FEM) for triangular Reissner–Mindlin flat shells is developed to improve the accuracy of numerical results. The present method, named edge/node-based S-FEM (ENS-FEM), uses a gradient smoothing technique over smoothing domains based on a combination of ES-FEM and NS-FEM. A discrete shear gap technique is incorporated into ENS-FEM to avoid shear-locking phenomenon in Reissner–Mindlin flat shell elements. For all practical purpose, we propose an average combination (aENS-FEM) of ES-FEM and NS-FEM for shell structural problems. We compare numerical results obtained using aENS-FEM with other existing methods in the literature to show the effectiveness of the present method.     

An improved Reissner–Mindlin triangular element

A modified version of the low-order mixed finite element method proposed by Oñate et al. [Internat. J. Numer. Method. Engng. 37 (1994) 2569] for the Reissner–Mindlin plate model is analyzed. It is proved that the new method is optimal convergent, uniform in the thickness of the plate.     

Numerical results for mimetic discretization of Reissner–Mindlin plate problems

A low-order mimetic finite difference method for Reissner–Mindlin plate problems is considered. Together with the source problem, the free vibration and the buckling problems are investigated. Details about the scheme implementation are provided, and the numerical results on several different types of meshes are reported.     

A characteristic stabilized finite element method for the non-stationary Navier–Stokes equations

This paper is concerned with the analysis of a new stabilized method based on the local pressure projection. The proposed method has a number of attractive computational properties: parameter-free, avoiding higher-order derivatives or edge-based data structures. Error estimates of the velocity and the pressure are obtained for both the continuous and the fully discrete versions. Finally, some numerical experiments show that this method is highly efficient for the non-stationary Navier–Stokes equations.     

A stabilized MITC element for accurate wave response in Reissner–Mindlin plates

Residual based finite element methods are developed for accurate time-harmonic wave response of the Reissner–Mindlin plate model. The methods are obtained by appending a generalized least-squares term to the mixed variational form for the finite element approximation. Through judicious selection of the design parameters inherent in the least-squares modification, this formulation provides a consistent and general framework for enhancing the wave accuracy of mixed plate elements. In this paper, the mixed interpolation technique of the well-established MITC4 element is used to develop a new mixed least-squares (MLS4) four-node quadrilateral plate element with improved wave accuracy. Complex wave number dispersion analysis is used to design optimal mesh parameters, which for a given wave angle, match both propagating and evanescent analytical wave numbers for Reissner–Mindlin plates. Numerical results demonstrates the significantly improved accuracy of the new MLS4 plate element compared to the underlying MITC4 element.     

A priori and a posteriori analysis for a locking-free low order quadrilateral hybrid finite element for Reissner–Mindlin plates

This paper proposes a quadrilateral finite element method of the lowest order for Reissner–Mindlin (R–M) plates on the basis of Hellinger–Reissner variational principle, which includes variables of displacements, shear stresses and bending moments. This method uses continuous piecewise isoparametric bilinear interpolation for the approximation of transverse displacement and rotation. The piecewise-independent shear stress/bending moment approximation is constructed by following a self-equilibrium criterion and a shear-stress-enhanced condition. A priori and reliable a posteriori error estimates are derived and shown to be uniform with respect to the plate thickness t. Numerical experiments confirm the theoretical results.     

A stable space–time finite element method for parabolic evolution problems

This paper is concerned with the analysis of a new stable space–time finite element method (FEM) for the numerical solution of parabolic evolution problems in moving spatial computational domains. The discrete bilinear form is elliptic on the FEM space with respect to a discrete energy norm. This property together with a corresponding boundedness property, consistency and approximation results for the FEM spaces yield an a priori discretization error estimate with respect to the discrete norm. Finally, we confirm the theoretical results with numerical experiments in spatial moving domains.     

A stabilized finite element method for stochastic incompressible Navier–Stokes equations

The major emphasis of this work is the development of a stabilized finite element method for solving incompressible Navier–Stokes equations with stochastic input data. The polynomial chaos expansion is used to represent stochastic processes in the variational problem, resulting in a set of deterministic variational problems to be solved for each Wiener polynomial chaos. To obtain the chaos coefficients in the corresponding deterministic incompressible Navier–Stokes equations, we combine the modified method of characteristics with the finite element discretization. The obtained Stokes problem is solved using a robust conjugate-gradient algorithm. This algorithm avoids projection procedures and any special correction for the pressure. These numerical techniques associate the geometrical flexibility of the finite element method with the ability offered by the modified method of characteristics to solve convection-dominated problems using time steps larger than its Eulerian counterpart. Numerical results are shown for the benchmark problems of driven cavity flow and backward-facing step flow. We also present numerical results for a problem of stochastic natural convection. It is found that the proposed stabilized finite element method offers a robust and accurate approach for solving the stochastic incompressible Navier–Stokes equations, even when high Reynolds and Rayleigh numbers are used in the simulations.     

A finite element/Monte–Carlo method for polymer dilute solutions

A finite element/Monte–Carlo method is proposed for solving the flow of a polymer dilute solution. The mass and momentum equations are supplemented with stochastic differential equations which model the dynamics of the polymer chains. Finite elements and a Monte–Carlo method are used to solve the problem. A theoretical analysis is performed on a simplified, deterministic, Oldroyd-B problem. Numerical results are compared to experimental ones. Received: 22 December 2000 / Accepted: 30 May 2001     

Two-level stabilized finite element method for the transient Navier–Stokes equations

In this article, a two-level stabilized finite element method based on two local Gauss integrations for the two-dimensional transient Navier–Stokes equations is analysed. This new stabilized method presents attractive features such as being parameter-free, or being defined for nonedge-based data structures. Some new a priori bounds for the stabilized finite element solution are derived. The two-level stabilized method involves solving one small Navier–Stokes problem on a coarse mesh with mesh size 0<H<1, and a large linear Stokes problem on a fine mesh with mesh size 0<h?H. A H ¹-optimal velocity approximation and a L ²-optimal pressure approximation are obtained. If we choose h=O(H ²), the two-level method gives the same order of approximation as the standard stabilized finite element method.     

Characteristic stabilized finite element method for the transient Navier–Stokes equations

Based on the lowest equal-order conforming finite element subspace (X_h, M_h) (i.e. P₁–P₁ or Q₁–Q₁ elements), a characteristic stabilized finite element method for transient Navier–Stokes problem is proposed. The proposed method has a number of attractive computational properties: parameter-free, avoiding higher-order derivatives or edge-based data structures, and averting the difficulties caused by trilinear terms. Existence,uniqueness and error estimates of the approximate solution are proved by applying the technique of characteristic finite element method. Finally, a serious of numerical experiments are given to show that this method is highly efficient for transient Navier–Stokes problem.     

A stable nonconforming quadrilateral finite element method for the stationary Stokes and Navier–Stokes equations

Recently, Douglas et al. [4] introduced a new, low-order, nonconforming rectangular element for scalar elliptic equations. Here, we apply this element in the approximation of each component of the velocity in the stationary Stokes and Navier–Stokes equations, along with a piecewise-constant element for the pressure. We obtain a stable element in both cases for which optimal error estimates for the approximation of both the velocity and pressure in L ² can be established, as well as one in a broken H ¹-norm for the velocity. Received: January 1999 / Accepted: April 1999     

A Weak Galerkin Method for the Reissner–Mindlin Plate in Primary Form

A new finite element method is developed for the Reissner–Mindlin equations in its primary form by using the weak Galerkin approach. Like other weak Galerkin finite element methods, this one is highly flexible and robust by allowing the use of discontinuous approximating functions on arbitrary shape of polygons and, at the same time, is parameter independent on its stability and convergence. Error estimates of optimal order in mesh size h are established for the corresponding weak Galerkin approximations. Numerical experiments are conducted for verifying the convergence theory, as well as suggesting some superconvergence and a uniform convergence of the method with respect to the plate thickness.