A Smooth Embedding Domain Method Based On the Penalty Approach

It is well known that convergence of the fictitious domain formulation with boundary Lagrange multipliers is slow due to the lower global regularity of its solution. This article presents a smoothed variant of this approach which is based on a formulation in the form of a state constraint optimal control problem. The convergence rate is increased as seen from a model example.     

对流扩散方程的有限体积-有限元方法的误差估计

本文结合有限体积方法和有限元方法处理非线性对流扩散问题,非线性对流项利用有限体积方法处理,扩散项利用有限元方法离散,并给近似解的误差估计。     

On the Convergence of Finite Element Method for Second Order Elliptic Interface Problems

The purpose of this paper is to study the convergence of finite element approximation to the exact solution of general self-adjoint elliptic equations with discontinuous coefficients. Due to low global regularity of the solution, it is difficult to achieve optimal order of convergence with classical finite element methods [Numer. Math. 1998; 79:175–202]. In this paper, an isoparametric type of discretization is used to prove optimal order error estimates in L ² and H ¹ norms when the global regularity of the solution is low. The interface is assumed to be of arbitrary shape and is smooth for our purpose. Further, for the purpose of numerical computations, we discuss the effect of numerical quadrature on finite element solution, and the related optimal order estimates are also established.     

On Fictitious Domain Formulations for Maxwell''s Equations

We consider fictitious domain-Lagrange multiplier formulations for variational problems in the space H(curl: Ω{\bf)} derived from Maxwell's equations. Boundary conditions and the divergence constraint are imposed weakly by using Lagrange multipliers. Both the time dependent and time harmonic formulations of the Maxwell's equations are considered, and we derive well-posed formulations for both cases. The variational problem that arises can be discretized by functions that do not satisfy an a-priori divergence constraint.     

Analysis on a Finite Volume Element Method for Stokes Problems

Finite volume element method for the Stokes problem is considered. We use a conforming piecewise linear function on a fine grid for velocity and piecewise constant element on a coarse grid for pressure. For general triangulation we prove the equivalence of the finite volume element method and a saddle-point problem, the inf-sup condition and the uniqueness of the approximation solution. We also give the optimal order H^1 norm error estimate. For two widely used dual meshes we give the L^2 norm error estimates, which is optimal in one case and quasi-optimal in another ease. Finally we give a numerical example.     

Maximum Norm Estimates for Finite Volume Element Method for Non-selfadjoint and Indefinite Elliptic Problems

In this paper, we establish the maximum norm estimates of the solutions of the finite volume element method (FVE) based on the P1 conforming element for the non-selfadjoint and indefinite elliptic problems.     

非自共轭和不定椭圆问题的有限体积元方法的一致收敛性

In this paper, we prove the existence, uniqueness and uniform convergence of the solution of finite volume element method based on the P1 conforming element for non-selfadjoint and indefinite elliptic problems under minimal elliptic regularity assumption.     

Stokes问题的区域分解法

本文主要讨论了Ｓｔｏｋｅｓ问题的非重迭型两仓区域性情形的区域分解算法，首先讨论了连续情形，然后将区域分解算法应用到Ｓｔｏｋｅｓ问题的非协调离散情形。     

Shape Optimization and Fictitious Domain Approach for Solving Free Boundary Problems of Bernoulli Type

This contribution deals with an efficient method for the numerical realization of the exterior and interior Bernoulli free boundary problems. It is based on a shape optimization approach. The state problems are solved by a fictitious domain solver using boundary Lagrange multipliers.     

Uniform Convergence Analysis for Singularly Perturbed Elliptic Problems with Parabolic Layers

In this paper,using Lin's integral identity technique,we prove the optimal uniform convergence θ(N_x~(-2)In~2N_x N_y~(-2)In~2N_y) in the L~2-norm for singularly per- turbed problems with parabolic layers.The error estimate is achieved by bilinear fi- nite elements on a Shishkin type mesh.Here N_x and N_y are the number of elements in the x- and y-directions,respectively.Numerical results are provided supporting our theoretical analysis.     

正交各向异性弹性力学平面问题的样条虚边界元法

采用域外奇点技术并根据问题的边界条件,建立了正交各向异性弹性力学平面问题的非奇异虚边界积分方程,然后采用性态优越的B样条函数去逼近未知虚荷载函数,并采用性能稳定的最小二乘边界子段法去消除边界余量,据此获得积分方程的数值解.数值算例表明:该方法具有相当高的精度和良好的数值稳定性,且计算工作量少.文中引言部分还对域外奇点法的发展作了系统的评述.     

无穷凹角区域各向异性问题的自然边界元法

本文研究无穷凹角区域上一类各向异性问题的自然边界元法.利用自然边界归化原理,获得该问题的Poisson积分公式和自然积分方程,给出了自然积分方程的数值方法,以及逼近解的收敛性和误差估计,最后给出了数值例子,以示方法的可行性和有效性.     

Shape Optimization for Semi-Linear Elliptic Equations Based on an Embedding Domain Method

We study a class of shape optimization problems for semi-linear elliptic equations with Dirichlet boundary conditions in smooth domains in ℝ². A part of the boundary of the domain is variable as the graph of a smooth function. The problem is equivalently reformulated on a fixed domain. Continuity of the solution to the state equation with respect to domain variations is shown. This is used to obtain differentiability in the general case, and moreover a useful formula for the gradient of the cost functional in the case where the principal part of the differential operator is the Laplacian. Online publication 23 January 2004.     

Nonoverlapping Domain Decomposition Method with Mixed Element for Elliptic Problems

1.IntroductionDomaindecompositionasanewmethodofcomputationalmathematics,waJsdevel-opedsincethedevelopmentofparallelcomputersandmultiprocessorsupercomputers-Usingdomaindecompositionwecandecreasethescaleoftheproblemandimplementthesub-problemsonparallelcomputer.Fromatechnicalpointofviewmostofdo-maindecompositionmethodsconsideredsofarhavebeendealingwithfiniteelementmethods.In[1,2]ZhangandHuanghavegivenakindofnonoverlappingdomaindecompositionprocedurewithpiecewiselinearfiniteelementapproximation.…     

Cascadic Multigrid for Finite Volume Methods for Elliptic Problems

In this paper, some effective cascadic multigrid methods are proposed for solving the large scale symmetric or nonsymmetric algebraic systems arising from the finite volume methods for second order elliptic problems. It is shown that these algorithms are optimal in both accuracy and computational complexity. Numerical experiments are reported to support our theory.     

自适应有限元方法在凹角域线性椭圆方程的应用(续)

作为序列文章自适应有限元方法在凹角域线性椭圆方程的应用的第三篇,在本文我们将给出并详细论证一个重要结论即 |?(u(x)-U(x))|≤Ch(x)|x|^β-2,|x|≥C′h且进一步分析说明在本序列文章的第一部分和地二部分得出方法都是以此为基础作出的。     

有界对称域的特征

We give a necessary and sufficient condition for a domain to be biholomorphic to a bounded symmetric domain.     

Fictitious domain method with penalty for an incompressible fluid

Fictitious domain method shows great advantages when handling problems with complex and constantly varying domains. In this article, we propose an algorithm which extends the fictitious domain method by introducing penalties. Test results with the numerical examples of backward facing step problem and the flow around steady and dynamic cylinder problem show that the algorithm we propose is highly efficient for solving incompressible fluid problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

抛物型方程在变动区域中的Galarkin有限元法及其对Stefan问题的应用

文中叙述了具变系数抛物型方程在变动区域中的有限元半离散格式和全离散θ格式,给出了误差界和稳定性结果,并应用于求解Stefan问题。