A mixed finite element method for the unilateral contact problem in elasticity

In this paper, we provide a new mixed finite element approximation of the varia-tional inequality resulting from the unilateral contact problem in elasticity. We use the continuous piecewise P2-P1 finite element to approximate the displacement field and the normal stress component on the contact region. Optimal convergence rates are obtained under the reasonable regularity hypotheses. Numerical example verifies our results.     

A coupling of nonconforming and mixed finite element methods for Biot's consolidation model

In this article, we develop a nonconforming mixed finite element method to solve Biot's consolidation model. In particular, this work has been motivated to overcome nonphysical oscillations in the pressure variable, which is known as locking in poroelasticity. The method is based on a coupling of a nonconforming finite element method for the displacement of the solid phase with a standard mixed finite element method for the pressure and velocity of the fluid phase. The discrete Korn's inequality has been achieved by adding a jump term to the discrete variational formulation. We prove a rigorous proof of a‐priori error estimates for both semidiscrete and fully‐discrete schemes. Optimal error estimates have been derived. In particular, optimality in the pressure, measured in different norms, has been proved for both cases when the constrained specific storage coefficient c₀ is strictly positive and when c₀ is nonnegative. Numerical results illustrate the accuracy of the method and also show the effectiveness of the method to overcome the nonphysical pressure oscillations. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A stabilized mixed formulation for finite strain deformation for low-order tetrahedral solid elements

A stabilized mixed finite element formulation for four-noded tetrahedral elements is introduced for robustly solving small and large deformation problems. The uniqueness of the formulation lies within the fact that it is general in that it can be applied to any type of material model without requiring specialized geometric or material parameters. To overcome the problem of volumetric locking, a mixed element formulation that utilizes linear displacement and pressure fields was implemented. The stabilization is provided by enhancing the rate of deformation tensor with a term derived using a bubble function approach. The element was implemented through a user-programmable element of the commercial finite element software ANSYS. Using the ANSYS platform, the performance of the element was evaluated by comparing the predicted results with those obtained using mixed quadratic tetrahedral elements and hexahedral elements with a B-bar formulation. Based on the quality of the results, the new element formulation shows significant potential for use in simulating complex engineering processes.     

A mixed nonconforming finite element for linear elasticity

This article considers a mixed finite element method for linear elasticity. It is based on a modified mixed formulation that enforces the continuity of the stress weakly by adding a jump term of the approximated stress on interior edges. The symmetric stress are approximated by nonconforming linear elements and the displacement by piecewise constants. We establish ??(h) error bound in the (broken) L² norm for the divergence of the stress and ??(h) error bound in the L² norm for both the displacement and the stress tensor. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

A finite element method for contact/impact

Ideas from the analysis of differential-algebraic equations are applied to the numerical solution of frictionless contact/impact problems in solid mechanics. Index-one and two formulations for dynamic contact–impact within the context of the finite element method are derived. The resulting equations are shown to stabilize the kinematic fields at the contact interface, at the expense of a small energy loss, which is shown to decrease consistently with mesh refinement. This energy dissipation is shown to be necessary for the establishment of persistent contact. A Newmark-type time integration scheme is derived from the proposed formulation, and shown to yield excellent results in modeling the transition to contact/impact.     

A posteriori error estimate for the mixed finite element method

A computable error bound for mixed finite element methods is established in the model case of the Poisson-problem to control the error in the H(div,) -norm. The reliable and efficient a posteriori error estimate applies, e.g., to Raviart-Thomas, Brezzi-Douglas-Marini, and Brezzi-Douglas-Fortin-Marini elements.

     

A two‐grid method for characteristic expanded mixed finite element solution of miscible displacement problem

A combined method consisting of mixed finite element method (MFEM) for the pressure equation and expanded mixed finite element method with characteristics(CEMFEM) for the concentration equation is presented to solve the coupled system of incompressible miscible displacement problem. To solve the resulting nonlinear system of equations efficiently, the two‐grid algorithm relegates all of the Newton‐like iterations to grids much coarser than the original one, with no loss in order of accuracy. It is shown that coarse space can be extremely coarse and our algorithm achieve asymptotically optimal approximation when the mesh sizes satisfy $H=O(h^{\frac{1}{4}})$. Numerical experiment is provided to confirm our theoretical results.     

A new family of stable mixed finite elements for the 3D Stokes equations

A natural mixed-element approach for the Stokes equations in the velocity-pressure formulation would approximate the velocity by continuous piecewise-polynomials and would approximate the pressure by discontinuous piecewise-polynomials of one degree lower. However, many such elements are unstable in 2D and 3D. This paper is devoted to proving that the mixed finite elements of this - type when satisfy the stability condition--the Babuska-Brezzi inequality on macro-tetrahedra meshes where each big tetrahedron is subdivided into four subtetrahedra. This type of mesh simplifies the implementation since it has no restrictions on the initial mesh. The new element also suits the multigrid method.

     

A new stabilized finite element method for optimal control for a Ladyzhenskaya model for unsteady flows

This article considers the time‐dependent optimal control problem of tracking the velocity for the viscous incompressible flows which is governed by a Ladyzhenskaya equations with distributed control. The existence of the optimal solution is shown and the first‐order optimality condition is established. The semidiscrete‐in‐time approximation of the optimal control problem is also given. The spatial discretization of the optimal control problem is accomplished by using a new stabilized finite element method which does not need a stabilization parameter or calculation of high order derivatives. Finally a gradient algorithm for the fully discrete optimal control problem is effectively proposed and implemented with some numerical examples. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 263–287, 2012     

A mixed finite element method for fourth order eigenvalue problems

A mixed finite element method for approximating eigenpairs of IV order elliptic eigenvalue problems with Dirichlet boundary conditions has been given. The method can be applied to the vibration analysis of anisotropic/orthotropic/isotropic/biharmonic plates. Computer implementation procedures for this mixed method are given along with the results of numerical experiments.     

A finite element model for nonlinear behaviour of piezoceramics under weak electric fields

It has been experimentally observed that piezoceramic materials exhibit different types of nonlinearities under different combinations of electric and mechanical fields. When excited near resonance in the presence of weak e to a Duffinor such as jump phenomena and presence of superharmonics in the response spectra. There has not been much work in the litrature to model these types of nonlinearities. Some authors have developed one-dimensional models for the above phenomenon and derived closed-form solutions for the displacement response of piezo-actuators. However, the generalized three-dimensional (3-D) formulation of electric enthalpy, the variational formulation and the FEM implementation have not yet been addressed, which are the focus of this paper. In this work, these nonlinearities have been modelled in a 3-D piezoelectric continuum using higher order quadratic and cubic terms in the generalized electric enthalpy density function. The coupled nonlinear finite element equations have been derived using variational formulation. A special linearization technique for assembling the nonlinear matrices and solution of the resulting nonlinear equations has been developed. The method has been used for simulating the nonlinear frequency response of a lead zirconate titanate plate excited near its first in-plane vibration resonance frequency with sinusoidal excitations of different electric field strengths. The results have been compared with those of the experiment.     

A finite element method for interface problems in domains with smooth boundaries and interfaces

This paper is concerned with the analysis of a finite element method for nonhomogeneous second order elliptic interface problems on smooth domains. The method consists in approximating the domains by polygonal domains, transferring the boundary data in a natural way, and then applying a finite element method to the perturbed problem on the approximate polygonal domains. It is shown that the error in the finite element approximation is of optimal order for linear elements on a quasiuniform triangulation. As such the method is robust in the regularity of the data in the original problem.     

A two‐grid stabilized mixed finite element method for Darcy‐Forchheimer model

A two‐grid stabilized mixed finite element method based on pressure projection stabilization is proposed for the two‐dimensional Darcy‐Forchheimer model. We use the derivative of a smooth function, , to approximate the derivative of in constructing the two‐grid algorithm. The two‐grid method consists of solving a small nonlinear system on the coarse mesh and then solving a linear system on the fine mesh. There are a substantial reduction in computational cost. We prove the existence and uniqueness of solution of the discrete schemes on the coarse grid and the fine grid and obtain error estimates for the two‐grid algorithm. Finally, some numerical experiments are carried out to verify the accuracy and efficiency of the method.     

A combined mixed finite element method and local discontinuous Galerkin method for miscible displacement problem in porous media

A combined method consisting of the mixed finite element method for flow and the local discontinuous Galerkin method for transport is introduced for the one-dimensional coupled system of incompressible miscible displacement problem. Optimal error estimates in L∞(0,T;L2) for concentration c,in L2(0,T;L2)for cxand L∞(0,T;L2) for velocity u are derived. The main technical difficulties in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method,the nonlinearity,and the coupling of the models. Numerical experiments are performed to verify the theoretical results. Finally,we apply this method to the one-dimensional compressible miscible displacement problem and give the numerical experiments to confirm the efficiency of the scheme.     

A robust optimal preconditioner for the mixed finite element discretization of elliptic optimal control problems

In this paper, we consider the efficient solving of the resulting algebraic system for elliptic optimal control problems with mixed finite element discretization. We propose a block‐diagonal preconditioner for the symmetric and indefinite algebraic system solved with minimum residual method, which is proved to be robust and optimal with respect to both the mesh size and the regularization parameter. The block‐diagonal preconditioner is constructed based on an isomorphism between appropriately chosen solution space and its dual for a general control problem with both state and gradient state observations in the objective functional. Numerical experiments confirm the efficiency of our proposed preconditioner.     

A cut finite element method for a Stokes interface problem

We present a finite element method for the Stokes equations involving two immiscible incompressible fluids with different viscosities and with surface tension. The interface separating the two fluids does not need to align with the mesh. We propose a Nitsche formulation which allows for discontinuities along the interface with optimal a priori error estimates. A stabilization procedure is included which ensures that the method produces a well conditioned stiffness matrix independent of the location of the interface.     

Sobolev方程一个新的H~1-Galerkin混合有限元分析

研究了Sobolev方程的H~1-Galerkin混合有限元方法.利用不完全双二次元Q_2~-和一阶BDFM元,建立了一个新的混合元模式,通过Bramble-Hilbert引理,证明了单元对应的插值算子具有的高精度结果.进一步,对于半离散和向后欧拉全离散格式,分别导出了原始变量u在H~1-模和中间变量p在H(div)-模意义下的超逼近性质.     

A numerical method based on extended Raviart–Thomas (ER‐T) mixed finite element method for solving damped Boussinesq equation

In this paper, we present the approximate solution of damped Boussinesq equation using extended Raviart–Thomas mixed finite element method. In this method, the numerical solution of this equation is obtained using triangular meshes. Also, for discretization in time direction, we use an implicit finite difference scheme. In addition, error estimation and stability analysis of both methods are shown. Finally, some numerical examples are considered to confirm the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

Hybridized weak Galerkin finite element method for linear elasticity problem in mixed form

A hybridization technique is applied to the weak Galerkin finite element method (WGFEM) for solving the linear elasticity problem in mixed form. An auxiliary function, the Lagrange multiplier defined on the boundary of elements, is introduced in this method. Consequently, the computational costs are much lower than the standard WGFEM. Optimal order error estimates are presented for the approximation scheme. Numerical results are provided to verify the theoretical results.