A non-conformal eXtended Finite Element approach: Integral matching Xfem

This work is dedicated to the mathematical and numerical analysis of a new Xfem approach: the integral maching Xfem. It is known that the quality of the approximation and the convergence rate of Xfem type methods is broadly influenced by the transition layer between the singular enrichment area and the rest of the domain. In the presented method, this transition layer is replaced by an interface associated with an integral matching condition of mortar type. We prove an optimal convergence result for such a non-conformal approximation method and we perform some numerical experiments showing the advantages of the integral matching Xfem with respect to former Xfem approaches.     

Numerical Model for Electro-Mechanical Analysis Based on Finite Element Method

In this paper, we present a numerical procedure that can be used to model the electro-mechanical coupled behavior of the dielectric actuator domain. The equation describing the electrostatical part is given by the reduced form of the Maxwell equation and the electrostatic potential [1]. The mechanical problem is described by the constitutive equations and equilibrium equations. Using the finite element method, this technique is to divide a whole problem into sub-problems. The complexity of the original problem is therefore reduced by focusing only on the most relevant areas. A finite element analysis is then performed by applying the electrostatic Maxwell pressure as Neumann boundary conditions to compute the displacements. Once the displacement is computed, the electrostatic domain or the conductor is updated. Electrostatic analysis is performed on the updated geometry and the finite element method is then used to determine the change in potential due to geometric perturbations. Once the surface charge densities are known, the new electrostatic Maxwell pressure is computed. The mechanical and electrostatic analysis is repeated until an equilibrium state is computed. The procedure is demonstrated in the paper by the solution of some two-dimensional and three-dimensional problems. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Cascadic Multigrid Method for Isoparametric Finite Element with Numerical Integration

The purpose of this paper is to study the cascadic multigrid method for the secondorder elliptic problems with curved boundary in two-dimension which are discretized by the isoparametric finite element method with numerical integration. We show that the CCG method is accurate with optimal complexity and traditional multigrid smoother (likesymmetric Gauss-Seidel, SSOR or damped Jacobi iteration) is accurate with suboptimal complexity.     

角质层渗透性质的有限元数值模拟计算

角质层是皮肤屏障作用的最主要部分，它决定了外界物质对皮肤的渗透情况．在假设角质层细胞为一种三维的十四面体（物理学经典的tetrakaidecahedron体）的情况下，利用有限元法对角质层渗透性质进行了数值模拟研究．为此，首先完成了对角质层空间结构的网格拆分，拆分过程分两步进行：1．对角蛋白细胞的网格拆分；2．对角蛋白细胞周围的网状脂质体的网格拆分．在数值模拟过程中，则用有限元法得到方程离散的格式，用多重网格算法降低高频误差，提高计算精度．最后，给出了数值模拟结果的可视化效果图．     

水污染问题特征有限元方法的数值计算及理论分析

本文研究了水污染二维对流占优数学模型特征有限元方法的计算问题，导出的计算格式对时间变量用特征线方法离散，对空间变量用Galerkin有限元方法离散，得到的H^1－模和L^2－模误差估计是最优阶的。     

Generalizations of the Finite Element Method

This paper is concerned with the generalization of the finite element method via the use of non-polynomial enrichment functions. Several methods employ this general approach, e.g. the extended finite element method and the generalized finite element method. We review these approaches and interpret them in the more general framework of the partition of unity method. Here we focus on fundamental construction principles, approximation properties and stability of the respective numerical method. To this end, we consider meshbased and meshfree generalizations of the finite element method and the use of smooth, discontinuous, singular and numerical enrichment functions.     

Numerical Analysis of the Behaviour of Tooth Intrarradial Posts by the Finite Element Method

The aim of this work is to study the behaviour of a carbon/epoxi post housed in a canine tooth after endodontic treatment in order to support the typical loads present during mastication. The three-D basic design of the dental piece consisting of tooth + post was carried out with a three-dimensional parametric design program. We study the stresses and displacements of the different elements of the dental piece under normal load conditions, and present the results and conclusions.     

精确有限元法

本文提出构造有限单元的新方法——精确有限元法.它可以求解在任意边界条件下任意变系数正定或非正定偏微分方程。文中给出它的收敛性证明和计算偏微分方程的一般格式。用精确元法所得到的单元是一个非协调元,单元之间的相容条件容易处理.与相同自由度普通有限元相比,由精确元法所得到的解的高阶导数具有较高的收敛精度.文末给出数值算例,所得到的结果均收敛于精确解,并有较好的数值精度.     

Semilinear Finite Element Method

In the Ritz-Galerkin method the linear subspace of the trial solution is extended to a closed subset. Some results, such as orthogonalization and minimum property of the error function, are obtained. A second order scheme is developed for solving a linear singular perturbation elliptic problem and error estimates are given for a uniform mesh size. Numerical results for linear and semilinear singular perturbation problems are included.     

基于浸入边界-有限元法的流固耦合碰撞数值模拟方法

针对流固耦合碰撞问题,建立了流体中固体与固体碰撞界面解析直接模拟方法,采用清晰界面浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时网格出现负体积的问题,采用基于罚函数的有限元方法对固体的运动和碰撞进行求解,以分域耦合方式实现流体域和固体域的耦合求解.通过与静止流体中球形颗粒与壁面正碰撞和斜碰撞的实验数据对比,验证了建立的数值模拟方法对流体中固体与固体碰撞数值模拟的正确性,获得了流体域流场在碰撞前后随时间的变化,同时通过该文建立的数值模拟方法也获得了固体域中固体的碰撞力和应力.未来,将把该数值模拟方法应用到流体流动环境中,如固体颗粒对管道的冲蚀、流体诱导海洋立管之间的碰撞、坠物对海底管道的撞击等.     

弹性波方程有限元模拟的曲边地表处理

本文对比研究了关于弹性波模拟中的曲边地表形状处理的两种方法,一种是用给定的实际介质数值划定的地表形状,另一种是用样条插值逼近地表形状.本文采用有限元方法进行弹性波数值模拟,给出了基于这两种方法计算的数值例子,并对结果进行了分析比较.结果表明使用后一种方法对地表进行处理时,地表人工离散产生的干扰明显减少,优于前一种方法.     

A New Finite Element Space for Expanded Mixed Finite Element Method

In this article, we propose a new finite element space Λ$_h$ for the expanded mixed finite element method (EMFEM) for second-order elliptic problems to guarantee its computing capability and reduce the computation cost. The new finite element space Λ$_h$ is designed in such a way that the strong requirement V$_h\subset$Λ$_h$ in [9] is weakened to {v$_h\in$V$_h$; divv$_h$=0}$\subset$Λ$_h$ so that it needs fewer degrees of freedom than its classical counterpart. Furthermore, the new Λ$_h$ coupled with the Raviart-Thomas space satisfies the inf-sup condition, which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix, and thus the existence, uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in $\mathbb{R}^d, d=2,3$ and for triangular partitions in $\mathbb{R}^2$. Also, the solvability of the EMFEM for triangular partition in $\mathbb{R}^3$ can be directly proved without the inf-sup condition. Numerical experiments are conducted to confirm these theoretical findings.     

Forbidden Shapes in the Finite Element Method

In two dimensional finite element analysis, increased accuracyis obtained by the addition of a number of nodal points alongthe sides of the elements. Curvilinear coordinates (p, q) wereintroduced to enable the interpolating functions to be obtainedin these complicated situations. Unfortunately as was pointedout by Jordan (1970) and is further outlined here, many usefulelement shapes with side nodes are "forbidden" by this proceduredue to the vanishing of the Jacobian of the transformation fromthe (p, q) system to the fixed (x, y) system. A new technique is introduced in the present paper, based ongeometrical considerations, whereby interpolating functionsare obtained directly in terms of x and y for the triangle andquadrilateral with arbitrarily placed side points. These localfunctions can be used to construct piecewise smooth global interpolatingfunctions with C^° continuity over regions possessing curvedboundaries and composed of elements which are triangles andparallelograms with arbitrarily positioned side points.     

An Approximation of Three-Dimensional Semiconductor Devices by Mixed Finite Element Method and Characteristics-Mixed Finite Element Method

The mathematical model for semiconductor devices in three space dimensions are numerically discretized. The system consists of three quasi-linear partial differential equations about three physical variables: the electrostatic potential, the electron concentration and the hole concentration. We use standard mixed finite element method to approximate the elliptic electrostatic potential equation. For the two convection-dominated concentration equations, a characteristics-mixed finite element method is presented. The scheme is locally conservative. The optimal $L^2$-norm error estimates are derived by the aid of a post-processing step. Finally, numerical experiments are presented to validate the theoretical analysis.     

Approximation properties of the Generalized Finite Element Method

In this paper, we have obtained an approximation result in the Generalized Finite Element Method (GFEM) that reflects the global approximation property of the Partition of Unity (PU) as well as the approximability of the local approximation spaces. We have considered a GFEM, where the underlying PU functions reproduce polynomials of degree l. With the space of polynomials of degree k serving as the local approximation spaces of the GFEM, we have shown, in particular, that the energy norm of the GFEM approximation error of a smooth function is O(h ^l + k ). This result cannot be obtained from the classical approximation result of GFEM, which does not reflect the global approximation property of the PU.     

The P-version Penalty Finite Element Method

A p-version penalty finite element method is used to solve themodel problem –u=f in , u=g on . Error estimates are derivedin H¹-norm. The p-version penalty method with extrapolationyields an approximate solution which converges at the optimalrate. Numerical results show the effectiveness of the p-versionpenalty method with extrapolation.     

Nonconforming Elements in the Mixed Finite Element Method

In this paper, an abstract error estimate of mixed finite element methods using nonconforming elements is presented. In addition, a class of nonconforming rectangular elements is proposed, and applied to Stokes equations. The optimal error estimate is given.     

On Exact Results in the Finite Element Method

We prove that the finite element method for one-dimensional problems yields no discretization error at nodal points provided the shape functions are appropriately chosen. Then we consider a biharmonic problem with mixed boundary conditions and the weak solution u. We show that the Galerkin approximation of u based on the so-called biharmonic finite elements is independent of the values of u in the interior of any subelement.