两相可混溶驱动问题中近似压力方程的一种新方法

1　介　　绍ΩR2为凸多边形区域,Ω上的两相可混溶驱动问题可由以下微分方程系统来描述a)-.[a(x,c)(p-r(c)]=.u=q,b)φ(x)ct+u.c-(Dc)=(c-c)q=g(c),(1.1)其中a(x,c)=-k(x)μ(c),k(x)为介质的渗透率,μ(c)为流体的粘度,p为流体的压力,φ(x)为介质的孔隙度,c为一相流体的体积浓度,q为外部流体的体积流速,且满足相容性条件∫Ωqdx=0.D是2×2阶矩阵,D=φ(x)[dmI+|u|(dlE(u)+dtE⊥(u))],E(u)=(uiuj/|u|2)2×2,dm为分子扩散系数,dl,dt分别为横向、纵向弥散系数.系统的边界条件、初始条件:n为边界单位外法向a)u.n=0,(x,t)∈Ω×Jb)2i,j=1Dij(…     

裂缝—孔隙介质中地下水污染问题的混合元方法

1 引言 地下水污染问题是一类重要的环境问题.由于地质结构往往是裂缝-孔隙的双重介质,给这类实际问题的数值模拟带来极大困难.本文从描述裂缝-孔隙双重介质的均匀化模型出发,研究问题的数值方法,并对其进行理论分析,对解决地下水污染问题具有重要的 理论和实用意义. 设Ω为R~2中具有光滑边界的有界区域,令J=(0,T].由达西定律和质量守恒原理有以下模型: (1.1) (1.2) (1.3)初边值条件为 (1.4) (1.5)其中p为流体压力,u为Darcy速度,c,c′分别为污染物的浓度和介质表面的吸附浓度,s_1(x),s_2(x)分别为流动水和非流动水中相对贮水率,D为扩散矩阵,为孔隙度,I为单位阵,dm为分子扩散系数,d_l,d_t为纵向和横向弥散系数,a(x)为交换系数,q为源汇项,c~*为源汇项处的浓度值,在q<0处,c~*=c,在q>0处,c~*为源汇项处污染物的已知浓度.γ为的外法向量.相容性条件为|q(x,t)dx=0,x∈Ω,t∈J.     

不可压混流驱动问题的最小二乘混合元方法

1引言考虑多孔介质中两相不可压缩可混溶渗流驱动问题,它是由一组非线性耦合的椭园型压力方程和抛物型浓度方程组成：dVV。—一山人V什）gVV却）一q,VEn,（．1）＆,,。＿．、。。—一。x）＿＋u·grade－dlv（D（u）grade）一（1－c）q－,xEn,tEJ,（1．2）＆”－－’”””‘”－”””——－’——,、—’一其中a（）一a（x,c）一是（x）／卢（c）,J一［0,Ti,DcyR‘为水平油藏区域．方程式（1．l）一（1．2）中各物理量的意义如下：广为流体压力,c为流体的浓度,u为流体的Darer速度,叶为源汇项,／一—。x（q,O）,…     

扩散方程的随机Dirichlet问题

令D。表示d+1维欧氏空间R。d的有界子集.旨在用概率方法利用时空布朗运动探讨D。上如下扩散方程的随机D irich let问题:12Δu(x。(t))+q(x。(t))u(x。(t))=tu(x。(t)),x。(t)∈D。(*)其中q是给定的定义在D。上的有界Ho。lder连续函数.本文解决了上述扩散方程(*)的随机Dirichlet问题的解在S3内存在性及唯一性问题.     

具混合边界条件的趋化-流体耦合模型解的整体存在性

本文考虑如下趋化-流体耦合模型的混合边值问题:{nt + u· ▽n =△nm-x▽· (n▽c),ct+u·▽c=/△c-cn,ut +u· ▽u =△u-▽π +n▽ψ,▽·u=0.主要研究其在空间有界二维区域上解的整体存在性及一致有界性问题.首先,证明慢扩散情形(m＞1)混合非齐次边值问题的一致有界弱解的整体存在...     

非线性双曲型方程的混合元方法

1　引　　言考虑下述非线性双曲型方程的混合问题:c(x,u)utt-.(a(x,u)u)=f(x,u,t),　　x∈Ω,t∈J,(1.1)u(x,0)=u0(x),　　x∈Ω,(1.2)ut(x,0)=u1(x),　　x∈Ω,(1.3)u(x,t)=-g(x,t),　　(x,t)∈Ω×J,(1.4)其中ΩR2是一具有Lipschitz边界Ω的有界区域,J=[0,T],0相似文献   

一类反应扩散方程组定性分析

§1引言 以往处理反应扩散方程组: (x=(x_1,X_2,…,x_n)∈Ω,u=(u_1,U_2,…,u_n))的工作,大都是考虑D为对角矩阵的情形。对于这种情形,只要反应项f(x,t,u)有较好的性质,可以用与处理单个方程类似的办法处理。在f(x,t,u)为拟单调函数的假定下,鲍家馼(C.V.Pao)用上下解法考虑了这类方程组的定解问题,见[1]、[2]、[3]。     

多孔介质两相驱动问题的特征差分和Schwarz区域分裂并行算法

本文考虑多孔介质两相驱动问题的数值解法。用混合元方法求解压力方程，可同时得到速度和压力的近似；对浓度方程，给出了两类特征差分与Ｓｃｈｗａｒｚ型区域分裂引结合的数值格式，以减小对流项产生的数值弥散，减小所处理问题的规模和实行并行计算。     

关于“问题征解”的解答

《数学学习》1990年第一期上刊登有“问题征解”一则:若u在有界闭区域D上连续,在D内部有一、二阶偏导数,满足且u在D边界上之值为0.试证u在D上只恒为0.     

海仑数组的探求

如果三角形的三边长为整数且面积亦为整数，则称之为海仑三角形．海仑三角形的三边长所构成的数组（a，b，c）称之为海会数组．本文对海会数组进行新的探索．假定D＞0，D不是平方数，c是非0整数．设x=u，y=V是不定方程x~2-Dy~2=c的一个解，那么就称u＋v是它的一个解．其中当u≥0，v≥0时，最小的一个叫做基本解．再设x＋y是Pell方程X~2-Dy~2=1的任意一个解，则容易验证（u十v）（X十y）（=ux＋uyD＋（ux＋ut）也是x~2-Dy~2=c的解．设三角形三边长分别为a，从一a十…，C，其中p为奇数（可正可负）．则其面积为由于这个关于C’的M次方程…     

多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.     

Numerical analysis of incompressible wormhole propagation with mass-preserving characteristic mixed finite element procedure

In this article, we construct a new combined characteristic mixed finite element procedure to simulate the incompressible wormhole propagation. In this procedure, we use the classical mixed finite element method to solve the pressure equation and a modified mass-preserving characteristic finite element method for the solute transport equation, and solve the porosity function straightly by the given concentration. This combined method not only keeps mass balance globally but also preserves maximum principle for the porosity. We considered the corresponding convergence and derive the optimal L²-norm error estimate. Finally, we present some numerical examples to confirm theoretical analysis.

     

Characteristic finite element analysis of pattern formation dynamical model in polymerizing actin flocks

To illustrate the characteristic states of actin wave formation numerically, such as, spots, spirals and traveling waves, we propose a finite element method for simulating pattern formation in polymerizing actin flocks. Here we use the traditional Galerkin finite element method to solve the average filament polarization equation, and a mass-conservative characteristic finite element method to solve the density of actin equation, which deals well with the convection-dominated problem and keeps the mass balance globally. We study the convergence of this method and give the corresponding error estimate under some regularity assumptions of the solution. Finally some numerical results show the pattern formation in polymerizing actin flocks.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

Monte Carlo method via a numerical algorithm to solve a parabolic problem

This paper is intended to provide a numerical algorithm consisted of the combined use of the finite difference method and Monte Carlo method to solve a one-dimensional parabolic partial differential equation. The numerical algorithm is based on the discretize governing equations by finite difference method. Due to the application of the finite difference method, a large sparse system of linear algebraic equations is obtained. An approach of Monte Carlo method is employed to solve the linear system. Numerical tests are performed in order to show the efficiency and accuracy of the present work.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

Combined use of the finite element and finite superelement methods

Results of the theoretical and numerical studies of an algorithm based on the combined use of the finite element and finite superelement methods are presented. Estimates of the errors for one of the variants of the method applied to solving the Laplace equation are obtained. The method can be used to solve a problem concerning the skin layer appearing due to high velocities.     

Solution of Burgers' equation for large Reynolds number using finite elements with moving nodes

Burgers' equation often arises in the mathematical modelling used to solve problems in fluid dynamics involving turbulence. Numerical difficulties arise in the solution for the case of large Reynolds number. To obtain high accuracy, finite element methods are important. The aim of this paper is to summarize relevant past work and to use a moving node finite element method to obtain a solution of Burgers' equation under certain prescribed conditions. The results for high Reynolds number are compared with accurate results obtained by other authors.     

抛物方程时间连续有限元全离散格式的超收敛性

利用张量积分解和时间方向单元正交分解,证明了线性抛物型方程的时间连续全离散有限元在单元节点和内部的特征点的超收敛性.并用连续有限元计算了非线性Schrodinger方程,验证了能量的守恒性.计算结果与理论相吻合.     

Decoupled characteristic stabilized finite element method for time‐dependent Navier–Stokes/Darcy model

In this article, we propose and analyze a new decoupled characteristic stabilized finite element method for the time‐dependent Navier–Stokes/Darcy model. The key idea lies in combining the characteristic method with the stabilized finite element method to solve the decoupled model by using the lowest‐order conforming finite element space. In this method, the original model is divided into two parts: one is the nonstationary Navier–Stokes equation, and the other one is the Darcy equation. To deal with the difficulty caused by the trilinear term with nonzero boundary condition, we use the characteristic method. Furthermore, as the lowest‐order finite element pair do not satisfy LBB (Ladyzhen‐Skaya‐Brezzi‐Babuska) condition, we adopt the stabilized technique to overcome this flaw. The stability of the numerical method is first proved, and the optimal error estimates are established. Finally, extensive numerical results are provided to justify the theoretical analysis.