462 research outputs found
Stable finite element pair for Stokes problem and discrete Stokes complex on quadrilateral grids
In this paper, we first construct a nonconforming finite element pair for the
incompressible Stokes problem on quadrilateral grids, and then construct a
discrete Stokes complex associated with that finite element pair. The finite
element spaces involved consist of piecewise polynomials only, and the
divergence-free condition is imposed in a primal formulation. Combined with
some existing results, these constructions can be generated onto grids that
consist of both triangular and quadrilateral cells
Genetic Exponentially Fitted Method for Solving Multi-dimensional Drift-diffusion Equations
A general approach was proposed in this article to develop high-order
exponentially fitted basis functions for finite element approximations of
multi-dimensional drift-diffusion equations for modeling biomolecular
electrodiffusion processes. Such methods are highly desirable for achieving
numerical stability and efficiency. We found that by utilizing the one-one
correspondence between continuous piecewise polynomial space of degree
and the divergence-free vector space of degree , one can construct
high-order 2-D exponentially fitted basis functions that are strictly
interpolative at a selected node set but are discontinuous on edges in general,
spanning nonconforming finite element spaces. First order convergence was
proved for the methods constructed from divergence-free Raviart-Thomas space
at two different node set
Recovery-Based Error Estimators for Diffusion Problems: Explicit Formulas
We introduced and analyzed robust recovery-based a posteriori error
estimators for various lower order finite element approximations to interface
problems in [9, 10], where the recoveries of the flux and/or gradient are
implicit (i.e., requiring solutions of global problems with mass matrices). In
this paper, we develop fully explicit recovery-based error estimators for lower
order conforming, mixed, and non- conforming finite element approximations to
diffusion problems with full coefficient tensor. When the diffusion coefficient
is piecewise constant scalar and its distribution is local quasi-monotone, it
is shown theoretically that the estimators developed in this paper are robust
with respect to the size of jumps. Numerical experiments are also performed to
support the theoretical results
The TDNNS method for Reissner-Mindlin plates
A new family of locking-free finite elements for shear deformable
Reissner-Mindlin plates is presented. The elements are based on the
"tangential-displacement normal-normal-stress" formulation of elasticity. In
this formulation, the bending moments are treated as separate unknowns. The
degrees of freedom for the plate element are the nodal values of the
deflection, tangential components of the rotations and normal-normal components
of the bending strain. Contrary to other plate bending elements, no special
treatment for the shear term such as reduced integration is necessary. The
elements attain an optimal order of convergence
Convergence of an adaptive mixed finite element method for general second order linear elliptic problems
The convergence of an adaptive mixed finite element method for general second
order linear elliptic problems defined on simply connected bounded polygonal
domains is analyzed in this paper. The main difficulties in the analysis are
posed by the non-symmetric and indefinite form of the problem along with the
lack of the orthogonality property in mixed finite element methods. The
important tools in the analysis are a posteriori error estimators,
quasi-orthogonality property and quasi-discrete reliability established using
representation formula for the lowest-order Raviart-Thomas solution in terms of
the Crouzeix-Raviart solution of the problem. An adaptive marking in each step
for the local refinement is based on the edge residual and volume residual
terms of the a posteriori estimator. Numerical experiments confirm the
theoretical analysis.Comment: 24 pages, 8 figure
- …