面向有限元分析的三角网格布尔运算方法

为加速有限元分析流程,提出了一种面向有限元分析的自适应三角网格模型布尔运算方法。将ADT（alternating digital tree）数据结构应用于三角面片的判交计算,提升了布尔运算中三角面片的判交效率;借鉴结合球填充算法和插点/去点算法的网格重生成方法来重剖分局部区域,保证了算法的效率和生成网格的质量;通过改进的八叉树背景网格来记录和光滑尺寸场,生成了自适应的网格;通过尺寸场合理地控制重剖分区域,使算法效率和生成网格的质量得以很好地平衡。实验结果表明,所提方法可以高效地生成满足有限元分析需求的网格。     

ACCURACY ANALYSIS FOR QUASI-CAREY ELEMENT*

In this paper, a new triangular element （Quasi-Carey element） is constructed by the idea of Specht element. It is shown that this Quasi-Carey element possesses a very special property, i.e., the consistency error is of order O（h^2）, one order higher than its interpolation error when the exact solution belongs to H^3（Ω）. However, the interpolation error and consistency error of Carey element are of order O（h）. It seems that the above special property has never been seen for other triangular elements for the second order problems.     

Nonconforming finite element method for nonlinear parabolic equations

<正> A nonconforming finite element method for the nonlinear parabolic equations is studied inthis paper.The convergence analysis is presented and the optimal error estimate in L~2(‖·‖_h)norm isobtained through Ritz projection technique,where ‖·‖_h is a norm over the finite element space.     

Two-scale finite element Green’s function approximations with applications to electrostatic potential computation

<正> In this paper,a two-scale finite element approach is proposed and analyzed for approximationsof Green's function in three-dimensions.This approach is based on a two-scale finite elementspace defined,respectively,on the whole domain with size H and on some subdomain containing singularpoints with size h (h H).It is shown that this two-scale discretization approach is very efficient.In particular,the two-scale discretization approach is applied to solve Poisson-Boltzmann equationssuccessfully.     

L ∞-estimates of mixed finite element methods for general nonlinear optimal control problems

This paper investigates L~∞-estimates for the general optimal control problems governed by two-dimensional nonlinear elliptic equations with pointwise control constraints using mixed finite element methods.The state and the co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control is approximated by piecewise constant functions.The authors derive L~∞-estimates for the mixed finite element approximation of nonlinear optimal control problems.Finally,the numerical examples are given.     

High accuracy analysis of tensor-product linear pentahedral finite elements for variable coefficient elliptic equations

For a general second-order variable coefficient elliptic boundary value problem in three dimensions,the authors derive the weak estimate of the first type for tensor-product linear pentahedral finite elements.In addition,the estimate for the W^1,1 -seminorm of the discrete derivative Green’s function is given.Finally,the authors show that the derivatives of the finite element solution u_h and the corresponding interpolantΠ_u are superclose in the pointwise sense of the L^∞-norm.     

基于B样条小波有限元IMU结构的模态分析

给出了一种惯性测量单元(IMU)系统结构,建立了整个模型B样条小波有限元模型;研究了一类新的有限元空间,它以B样条小波函数作为有限等参元的形状函数;建立了IMU模型的B样条小波有限元序列,利用B样条小波函数的变尺度特性在不改变网格的剖分下提高分辨率,因此在处理局部应力集中,曲率分布突变边界等应用中具有一定的优势,其结果更接近于实际;最后对系统模型进行了仿真,给出了各个阶次模态值,通过与传统有限元比较发现,小波有限元具有收敛性好、求解迅速和网格划分灵活等优点。     

Upwind finite volume element methods for one-dimensional semiconductor device

A fully discrete implicit Euler upwind finite volume element method is derived and studied for one-dimensional semiconductor device. Upwind scheme is introduced to deal with the convection-dominated diffusion equations in the semiconductor model. With different time steps for the electrostatic potential and the other unknown quantities, the computational procedure of the method is obtained. The local mass conservation laws are preserved under the framework of the upwind finite volume element schemes. A first-order accuracy in the L ²-norm is proved. Numerical experiments are given to validate the usefulness and efficiency of the method.     

Nonhomogeneous boundary value problems for the Korteweg-de Vries equation on a bounded domain

<正> This paper studies an initial-boundary-value problem (IBVP) of the Korteweg-de Vriesequation posed on a finite interval with general nonhomogeneous boundary conditions.Using thestrong Kato smoothing property of the associated linear problem,the IBVP is shown to be locallywell-posed in the space H~s(0,1) for any s≥0 via the contraction mapping principle.     

A discontinuous Galerkin method combined with mixed finite element for seawater intrusion problem

<正> Seawater intrusion problem is considered in this paper.Its mathematical model is anonlinear coupled system of partial differential equations with initial boundary problem.It consistsof the water head equation and the salt concentration equation.A combined method is developedto approximate the water head equation by mixed finite element method and concentration equationby discontinuous Galerkin method.The scheme is continuous in time and optimal order estimates inH~1-norm and L~2-norm are derived for the errors.     

Least-squares Galerkin procedure for second-order hyperbolic equations

This paper proposes the least-squares Galerkin finite element scheme to solve second-order hyperbolic equations. The convergence analysis shows that the method yields the approximate solutions with optimal accuracy in (L ²(Ω))² × L ²(Ω) norms. Moreover, the method gets the approximate solutions with second-order accuracy in time increment. A numerical example testifies the efficiency of the novel scheme.     

Minimizing a linear fractional function subject to a system of sup-T equations with a continuous Archimedean triangular norm

This paper shows that the problem of minimizing a linear fractional function subject to a system of sup-T equations with a continuous Archimedean triangular norm T can be reduced to a 0-1 linear fractional optimization problem in polynomial time. Consequently, parametrization techniques, e.g., Dinkelbach’s algorithm, can be applied by solving a classical set covering problem in each iteration. Similar reduction can also be performed on the sup-T equation constrained optimization problems with an objective function being monotone in each variable separately. This method could be extended as well to the case in which the triangular norm is non-Archimedean. The work is supported by the National Science Foundation of the United States under Grant No. #DMI-0553310.     

THE SUPERCONVERGENCE ANALYSIS OF AN ANISOTROPIC FINITE ELEMENT

This paper deals with the high accuracy analysis of bilinear finite element on the class of anisotropic rectangular meshes. The inverse inequalities on anisotropic meshes are established. The superclose and the superconvergence are obtained for the second order elliptic problem. A numerical test is given, which coincides with our theoretical analysis.     

High accuracy analysis of the finite element method for nonlinear viscoelastic wave equations with nonlinear boundary conditions

The standard finite elements of degree p over the rectangular meshes are applied to solve a kind of nonlinear viscoelastic wave equations with nonlinear boundary conditions, and the superclose property of the continuous Galerkin approximation is derived without using the nonclassical elliptic projection of the exact solution of the model problem. The global superconvergence of one order higher than the traditional error estimate is also obtained through the postprocessing technique.     

An anisotropic nonconforming finite element scheme with moving grids for parabolic integro-differential equations

A Crank-Nicolson scheme based on nonconforming finite element with moving grids is discussed for a class of parabolic integro-differential equations under anisotropic meshes. The corresponding convergence analysis is presented and the error estimates are obtained by using the interpolation operator instead of the conventional elliptic projection which is an indispensable tool in the convergence analysis of traditional finite element methods in previous literature.     

多分辨率网格的数据压缩

针对任一个三角形网格模型，按其简化顺序，对其三角形表数据进行重新排列，将其直接构造成一个多分辨率表示，构造的多分辨率表示与单分辨率表示完全一样，其存储空间是原模型的存储空间100％，达到了多分辨率表示存储空间的下限，实现了多分辨率表示与单分辨率表示的统一。进一步地，对多分辨率模型的三角形表数据的存储空间再进行压缩。实验表明：压缩后的三角形表数据的存储空间是压缩前对应的存储空间的70％，比单分辨率表示方式还简单、省空间。本文的多分辨率表示的构造方法及数据压缩方法非常简单，适合一般的三角形网格模型。     

RESIDUAL A POSTERIORI ERROR ESTIMATE OF A NEW TWO-LEVEL METHOD FOR STEADY NAVIER-STOKES EQUATIONS

Residual-based a posteriori error estimate for conforming finite element solutions of incom-pressible Navier-Stokes equations,which is computed with a new two-level method that is differentfrom Volker John,is derived.A posteriori error estimate contains additional terms in comparison tothe estimate for the solution obtained by the standard finite element method.The importance of theadditional terms in the error estimates is investigated by studying their asymptotic behavior.For opti-mal scaled meshes,these bounds are not of higher order than the convergence of discrete solution.Thetwo-level method aims to solve the nonlinear problem on a coarse grid with less computational work,then to solve the linear problem on a fine grid,which is superior to the usual finite element methodsolving a similar nonlinear problem on the fine grid.     

An FE-Inexact Heterogeneous ADMM for Elliptic Optimal Control Problems with <Emphasis Type="Italic">L</Emphasis><Superscript>1</Superscript>-Control Cost

Elliptic PDE-constrained optimal control problems with L¹-control cost (L¹-EOCP) are considered. To solve L¹-EOCP, the primal-dual active set (PDAS) method, which is a special semismooth Newton (SSN) method, used to be a priority. However, in general solving Newton equations is expensive. Motivated by the success of alternating direction method of multipliers (ADMM), we consider extending the ADMM to L¹-EOCP. To discretize L¹-EOCP, the piecewise linear finite element (FE) is considered. However, different from the finite dimensional l1-norm, the discretized L¹-norm does not have a decoupled form. To overcome this difficulty, an effective approach is utilizing nodal quadrature formulas to approximately discretize the L¹-norm and L²-norm. It is proved that these approximation steps will not change the order of error estimates. To solve the discretized problem, an inexact heterogeneous ADMM (ihADMM) is proposed. Different from the classical ADMM, the ihADMM adopts two different weighted inner products to define the augmented Lagrangian function in two subproblems, respectively. Benefiting from such different weighted techniques, two subproblems of ihADMM can be efficiently implemented. Furthermore, theoretical results on the global convergence as well as the iteration complexity results o(1/k) for ihADMM are given. In order to obtain more accurate solution, a two-phase strategy is also presented, in which the primal-dual active set (PDAS) method is used as a postprocessor of the ihADMM. Numerical results not only confirm error estimates, but also show that the ihADMM and the two-phase strategy are highly efficient.     

Superconvergence of mixed covolume method on quadrilateral grids for elliptic problems

This paper considers the mixed covolume method for the second-order elliptic equations over quadrilaterals.Superconvergence results are established in this paper on quadrilateral grids satisfying the h~2-parallelogram condition when the lowest-order Raviart-Thomas space is employed in the mixed covolume method.The authors prove O(h~2) accuracy between the approximate velocity or pressure and a suitable projection of the real velocity or pressure in the L~2 norm.Numerical experiments illustrating the theoretical results are provided.     

GLOBAL SUPERCONVERGENCE ANALYSIS OF WILSON ELEMENT FOR SOBOLEV AND VISCOELASTICITY TYPE EQUATIONS

In this paper the Wilson nonconforming finite element is employed to solve Sobolev and viscoelasticity type equations. By means of post-processing technique, global superconvergence estimates are obtained for quasi-uniform rectangular meshes. Finally, an error correction scheme is presented.