Stokes方程的压力梯度局部投影间断有限元法

本文对定常的Stokes方程提出了一种新的间断有限元法,通过将通常的间断Galerkin有限元法与压力梯度局部投影相结合,建立了一个稳定的间断有限元格式,对速度和压力的任意分片多项式空间P_l(K),P_m(K)的间断有限元逼近证明了解的存在唯一性,给出了关于速度和压力的L~2范数的最优误差估计.     

可压缩Navier-Stokes方程的压力梯度局部投影间断有限元法

将压力梯度投影与间断有限元法相结合,对可压缩线性化N-S方程提出了一种稳定化间断有限元格式．证明了此格式在速度和压力有限元空间无需满足B-B型条件的情况下,解的存在性和唯一性,以及相应的误差估计．     

拟等级网格下非线性延迟微分方程间断有限元法

本文利用间断有限元法求解非线性延迟微分方程,在拟等级网格下.给出非线性延迟微分方程间断有限元解的整体收敛阶和局部超收敛阶,数值实验验证了理论结果的正确性.     

非定常Navier-Stokes方程的稳定化特征有限元法

1引言特征线有限元法是求解对流扩散问题的有效方法。在处理对流占优问题时,表现出了很好的稳定性[8]。对于求解Navier-Stokes方程,文[9]建立了特征有限元格式,并进行了详细分析,但得到的收敛阶O(h~m △t (h~(m 1)/△t))只是拟丰满的。文[10]对此作了非线性稳定性的进一步分析,给出了关于速度和压力的最优误差估计。但目前所有的特征有限元法都要求有限元空间满足inf-sup条件,这就排除了工程实际应用计算方便的低阶有     

四段式有限元法分析热耦合的流体-固体相互作用问题

给出了一种流(体)-热-结构综合的分析方法,固体中的热传导耦合了粘性流体中的热对流,因而在固体中产生热应力.应用四段式有限元法和流线逆风Petrov-Galerkin法分析热粘性流动,应用Galerkin法分析固体中的热传导和热应力.应用二阶半隐式Crank-Nicolson格式对时间积分,提高了非线性方程线性化后的计算效率.为了简化所有有限元公式,采用3节点的三角形单元,对所有的变量:流体的速度分量、压力、温度和固体的位移,使用同阶次的插值函数.这样做的主要优点是,使流体-固体介面处的热传导连接成一体.数个测试问题的结果表明,这种有限元法是有效的,且能加深对流(体)-热-结构相互作用现象的理解.     

解Stokes方程的一种无散稳定二次有限体积法格式

在三角形网格上构造了一种求解Stokes方程的Lagrange二次有限体积法格式.取连续的二次有限元空间与间断的线性有限元空间分别作为Stokes方程的速度项与压力项的试探空间,从而保证了离散方程的速度解在宏元三角形单元上满足局部质量守恒性,且有限元空间对自然满足所谓的inf-sup条件.采用特殊的有限体积法映射与对偶剖分,求解Stokes方程的Lagrange二次有限体积法格式等价于相对应的有限元法格式,因此确保了有限体积法格式的无条件(无需约束三角形网格的几何形状)稳定性和关于速度项的最优阶H1范数的误差估计.最后,数值实验展示了理论结果的正确性以及有限体积法的数值模拟在计算流体力学中的有效性.     

对流扩散方程的混合时间间断时空有限元方法

构造并分析二阶对流扩散方程的混合时间间断时空有限元格式．利用混合有限元方法将二阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散低阶方程．证明数值解的稳定性、存在唯一性和收敛性．最后通过数值结果验证该算法的有效性和可行性．     

一类四阶抛物型积分-微分方程的混合间断时空有限元法

构造四阶抛物型积分-微分方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明离散解的稳定性,存在唯一性和收敛性.     

带广义边界条件的四阶抛物型方程的混合间断时空有限元法

构造具有广义边界条件的四阶线性抛物型方程的混合间断时空有限元格式,利用混合有限元方法将高阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散方程,证明了离散解的存在唯一性,稳定性和收敛性,并给出数值算例验证了方法的有效性.     

手性障碍电磁散射的PML有限元计算

该文建立了手性障碍电磁散射问题的二维模型, 给出问题的有限元分析, 并利用结合PML(perfectly matched layers)技术的有限元法进行数值模拟.     

Numerical analysis of two partitioned methods for uncoupling evolutionary MHD flows

Magnetohydrodynamics (MHD) studies the dynamics of electrically conducting fluids, involving Navier–Stokes (NSE) equations in fluid dynamics and Maxwell equations in eletromagnetism. The physical processes of fluid flows and electricity and magnetism are quite different and numerical simulations of each subprocess can require different meshes, time steps, and methods. In most terrestrial applications, MHD flows occur at low‐magnetic Reynold numbers. We introduce two partitioned methods to solve evolutionary MHD equations in such cases. The methods we study allow us at each time step to call NSE and Maxwell codes separately, each possibly optimized for the subproblem's respective physics. Complete error analysis and computational tests supporting the theory are given.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1083–1102, 2014     

Automated FEM discretizations for the Stokes equation

Current FEM software projects have made significant advances in various automated modeling techniques. We present some of the mathematical abstractions employed by these projects that allow a user to switch between finite elements, linear solvers, mesh refinement and geometry, and weak forms with very few modifications to the code. To evaluate the modularity provided by one of these abstractions, namely switching finite elements, we provide a numerical study based upon the many different discretizations of the Stokes equations. AMS subject classification (2000)  74S05, 65Y99, 35Q30     

A posteriori error analysis for locally conservative mixed methods

In this work we present a theoretical analysis for a residual-type error estimator for locally conservative mixed methods. This estimator was first introduced by Braess and Verfürth for the Raviart-Thomas mixed finite element method working in mesh-dependent norms. We improve and extend their results to cover any locally conservative mixed method under minimal assumptions, in particular, avoiding the saturation assumption made by Braess and Verfürth. Our analysis also takes into account discontinuous coefficients with possibly large jumps across interelement boundaries. The main results are applied to the nonconforming finite element method and the interior penalty discontinuous Galerkin method as well as the mixed finite element method.

     

On the relationship between finite volume and finite element methods applied to the Stokes equations

We investigate the relationship between finite volume and finite element approximations for the lower‐order elements, both conforming and nonconforming for the Stokes equations. These elements include conforming, linear velocity‐constant pressure on triangles, conforming bilinear velocity‐constant pressure on rectangles and their macro‐element versions, and nonconforming linear velocity‐constant pressure on triangles and nonconforming rotated bilinear velocity‐constant pressure on rectangles. By applying the relationship between the two methods, we obtain the convergence finite volume solutions for the Stokes equations. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 440–453, 2001.     

Weak Galerkin finite element methods for Parabolic equations

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal‐order error estimates in both H¹ and L² norms are established. Numerical tests are performed and reported. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A note on the design of hp-version interior penalty discontinuous Galerkin finite element methods for degenerate problems

^** Email: emmanuil.georgoulis{at}mcs.le.ac.uk^*** Email: al{at}maths.strath.ac.uk We consider a variant of the hp-version interior penalty discontinuousGalerkin finite element method (IP-DGFEM) for second-order problemsof degenerate type. We do not assume uniform ellipticity ofthe diffusion tensor. Moreover, diffusion tensors of arbitraryform are covered in the theory presented. A new, refined recipefor the choice of the discontinuity-penalization parameter (thatis present in the formulation of the IP-DGFEM) is given. Makinguse of the recently introduced augmented Sobolev space framework,we prove an hp-optimal error bound in the energy norm and anh-optimal and slightly p-suboptimal (by only half an order ofp) bound in the L² norm (the latter, for the symmetric versionof the IP-DGFEM), provided that the solution belongs to an augmentedSobolev space.     

A PARALLEL ITERATIVE DOMAIN DECOMPOSITION ALGORITHM FOR ELLIPTIC PROBLEMS

1.IntroductionNolloverlappillgdomaindecolllpositionnletllodshavereceivedalotofattentionlenlsilllldallowefficielltparallelisnl.F'Orarecentdevelopmelltofthesemethods,werefertot…     

Convergence analysis of an adaptive edge element method for Maxwell's equations

An efficient and reliable a posteriori error estimate is derived for solving three-dimensional static Maxwell's equations by using the edge elements of first family. Based on the a posteriori error estimates, an adaptive finite element method is constructed and its convergence is established. Compared with the existing results, an important advantage of the new theory lies in its feature that the usual marking of elements based on the oscillation is not needed in our adaptive algorithm, while the linear convergence of the algorithm can be still demonstrated in terms of the reduction of the energy-norm error and the oscillation. Numerical examples are provided which demonstrate the effectiveness and robustness of the adaptive methods.     

IMPLICIT-EXPLICIT MULTISTEP FINITE ELEMENT-MIXED FINITE ELEMENT METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device consists of a Poisson equa-tion for the electric potential and of two nonlinear parabolic equations for the electrondensity and hole density.The electric potential equation is discretized by a mixed finiteelement method. The electron and hole density equations are treated by implicit-explicitmultistep finite element methods. The schemes are very efficient. The optimal order errorestimates both in time and space are derived.     

A IPN×IPN Spectral Element Projection Method for the Unsteady Incompressible Navier-Stokes Equations

In this paper, we present a PN×PN spectral element method and a detailed comparison with existing methods for the unsteady incompressible Navier-Stokes equations. The main purpose of this work consists of: (i) detailed comparison and discussion of some recent developments of the temporal discretizations in the frame of spectral element approaches in space; (ii) construction of a stable PN×PN method together with a PN→PN-2 post-filtering. The link of different methods will be clarified. The key feature of our method lies in that only one grid is needed for both velocity and pressure variables, which differs from most well-known solvers for the Navier-Stokes equations. Although not yet proven by rigorous theoretical analysis, the stability and accuracy of this one-grid spectral method are demonstrated by a series of numerical experiments.