带滑动边界条件的Stokes方程的边界元逼近

前言 带滑动边界条件的Stokes方程,在诸如具有自由表面或具有大攻角的流体模型中起着重要作用。在电镀或容器壁可与流体起化学反应等流动问题中,经典的Stokes问题的不滑动边界条件不再成立,而滑动边界条件才是适当的物理模型。关于这一类实际问题,已有一些数值结果,但仅有很少的工作是就一般问题进行系统的分析。在文献[9]中,R.Verfuth就带滑动边界条件的定常Navier-Stokes方程给出了一种混合有限     

解Stokes问题的P1非协调四边形元的稳定化方法

本文研究了用(~P)_1-Q_0元(其中(~P)_1表示P_1非协调四边形元)解Stokes问题的非协调混合有限元稳定化逼近方法.(~P)_1-Q_0元不满足LBB条件(见[7,14] ),因而其不能直接用来求解Stokes问题.受[3] 的启发,我们提出了一种用(~P)_1-Q_0元解Stokes问题的稳定化方法,证明了这种方法的稳定性和离散问题解的存在唯一性,得到了最优误差估计.文章最后给出的数值算例验证了我们的理论结果.     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

非定常Stokes方程的降阶稳定化CN有限体积元外推模型

利用稳定化的Crank-Nicolson(CN)有限体积元方法和特征投影分解方法,建立非定常Stokes方程的一种自由度很少、精度足够高的降阶稳定化CN有限体积元外推模型,并给出这种降阶稳定化CN有限体积元外推模型解的误差估计和算法的实现.最后用数值例子说明数值结果与理论结果相吻合,并阐明这种降阶稳定化CN有限体积元外推模型的优越性.     

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

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

求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

Stokes问题的非协调区域分解法

1 引言 区域分解法是近年来新崛起的偏微分方程的数值解法.由于区域分解法比通常的数值解法有其独特的优点,再加上并行机的迅速发展,故国内外大批数值分析学家竞相投入这一研究行列.现在这方面的研究工作已有许多[1,2,8,9,10,12,13,15,16,18,19,20],并且在实际计算方面也开始理论化,系统化.但到目前为止,对非协调元的非重迭型区域分解法的研究还甚少,尤其是用来解决Stokes问题.对此,我们提出本文,目的在于利用在实际计算中经常采用的一类非协调元(C-R元),配合区域分解法这一新思想来处理Stokes问题,并得到了所给定算法的几何收敛性结果. 对Stokes问题连续情形下的区域分解算法的讨论已在文[8]中有所涉及,这里不再叙述.本文主要讨论Stokes问题非协调离散情形下的区域分解算法.     

非定常Stokes方程的稳定化全离散有限体积元格式

本文研究非定常Stokes方程的有限体积元方法,给出一种基于两个局部高斯积分的稳定化全离散格式,并给其有限体积元解的误差分析.     

关于Stokes和线性Navier-Stokes方程的广义维数分裂迭代方法

本文研究了鞍点问题的迭代法.在Benzi等人提出的维数分裂(DS)迭代方法的基础上,提出了具有三个参数的广义维数分裂(GDS)迭代法,该方法包含了DS迭代法,理论分析表明该方法是无条件收敛的.通过对有限差分法和有限元法离散的Stokes问题及有限元法离散的Oseen问题的数值结果表明,本文所给方法是有效的.     

Stokes波高阶谐波系数递推的数值计算*

本文用W.H.Hui提出的方法,在半物理平面内重新表述了Stokes波的数学模型和边界条件,提出了两种更有效的数值计算方法来获得Stokes波高阶谐波系数,并可递推至无穷.通过小参数转换,重新得到了Cokelet(1977)的波速和半波高的摄动展开式.     

Singularities of the quad curl problem

We consider the quad curl problem in smooth and non smooth domains of the space. We first give an augmented variational formulation equivalent to the one from [25] if the datum is divergence free. We describe the singularities of the variational space which correspond to the ones of the Maxwell system with perfectly conducting boundary conditions. The edge and corner singularities of the solution of the corresponding boundary value problem with smooth data are also characterized. We finally obtain some regularity results of the variational solution.     

THE ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE COMPLEX AMPLITUDE IN A COUPLED BAY-RIVER SYSTEM

We consider the numerical approximations of the complex amplitude in a coupled bayriver system in this work. One half-circumference is introduced as the artificial boundary in the open sea, and one segment is introduced as the artificial boundary in the river if the river is semi-infinite. On the artificial boundary a sequence of high-order artificial boundary conditions are proposed. Then the original problem is solved in a finite computational domain, which is equivalent to a variational problem. The numerical approximations for the original problem are obtained by solving the variational probiem with the finite element method. The numerical examples show that the artificial boundary conditions given in this work are very effective.     

THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the up-stream artificial boundary and two vertical sides are introduced as the downstream arti-ficial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.     

Signorini接触问题的边界变分不等式

本以Signorini接触问题为背景，讨论了变分不等式与边值问题的等价性，利用Green公式，基本解和基本解法向导数的性质，将区域型变分不等式化成等价的边界型变分不等式，并证明了边界变分不等式解的存在唯一性，为使用边界元方法数值求解提供理论依据。     

周期裂纹削弱的无限长板条的应力分析

作出了周期裂纹削弱的无限长板条的应力分析．假设这些裂纹均在水平位置,又板条承受 y方向的拉伸力p．此时边值问题归结为一个复杂混合边值问题．发现,对此问题言,特征展开变分原理方法 ( eigenfunction expansion variational method,简称为EEVM)是非常有效的．研究了裂纹端的应力强度因子和T-应力．从拉伸力作用下的弹性变形考虑,开裂板条可等价于一不开裂的正交异性板条．还分析了等价正交异性板条的弹性性质．最后给出了算例和数值结果．     

平面弹性力学问题等价的间接变量边界积分方程

用变分法确立平面弹性力学外边值问题的确切形式.在此基础之上,导出外边值问题的等价的间接变量边界积分方程.一些实例表明传统的惯用的直接变量边界积分方程与原边值问题是不等价的.     

Alternating iteration and elliptic variational inequalities

Summary Free boundary value problems, too complicated for formulation as a variational inequality, are broken up into two problems on overlapping regions. On one region the problem is treated as an ordinary boundary value problem; on the second region, the free boundary part of the problem is reduced to a variational inequality. By solving the two problems successively it is shown that under certain conditions the successive solutions converge to a single function that gives a solution of the original problem. Application to a filtration problem is given.     

Asymptotic analysis of the linear Boltzmann equation

The stiff boundary value problem for the time-independent linear Boltzmann equation is considered in the variational formulation in terms of the coercive symmetric bilinear form. It is shown that it is asymptotically equivalent to the related variational diffusion problem when the stiffness parameter tends to zero. Higher order perturbations are also discussed.     

ARTIFICIAL BOUNDARY METHOD FOR THE THREE-DIMENSIONAL EXTERIOR PROBLEM OF ELASTICITY

The exact boundary condition on a spherical artificial boundary is derived for thethree-dimensional exterior problem of linear elasticity in this paper. After this bound-ary condition is imposed on the artificial boundary, a reduced problem only defined in abounded domain is obtained. A series of approximate problems with increasing accuracycan be derived if one truncates the series term in the variational formulation, which isequivalent to the reduced problem. An error estimate is presented to show how the errordepends on the finite element discretization and the accuracy of the approximate problem.In the end, a numerical example is given to demonstrate the performance of the proposedmethod.     

电动力学电磁场边值问题的广义变分原理

给出了线性各项异性电磁场边值问题的广义虚功原理表达式,运用钱伟长教授提出的方法建立了该问题的广义变分原理,可直接反映该问题的全部特征,即4个Maxwell方程、2个场强-位势方程、2个本构方程和8个边界条件．继而导出了一族有先决条件的广义变分原理．作为例证,导出了两个退化形式的广义变分原理,和已知的广义变分原理等价．此外还导出了两个修正的广义变分原理,可为该问题提供杂交有限元模型．建立的各广义变分原理可为电磁场边值问题的有限元应用提供更为完善的理论基础．