多孔介质中可压缩可混溶驱动问题的特征—有限元方法

在有界区域上多孔介质中可压缩可混溶的油、水两相驱动问题是由非线性偏微分方程组的初、边值问题所决定.Douglas和Roberts曾提出其数学模型并研究了半离散化方法.本文对压力方程采用有限元和混合元两种方法.对饱和度方程采用特征—有限元方法.此方法的截断误差较标准有限元小的多,随之饱和度的计算更加精确.且可用     

可压缩可混溶驱动问题的有限元配置法

多孔介质中可压缩可混溶驱动问题是非线性抛物系统,压力方程和饱和度方程用有限元配置方法来求解,证明了配置解的存在唯一性,最后得到了最优阶的误差估计．     

多孔介质中可压缩可混溶驱动问题的有限体积元法

有界区域上多孔介质中可压缩可混溶驱动问题由两个非线性抛物型方程耦合而成：压力方程和饱和度方程均是抛物型方程．运用有限体积元法对两个方程进行数值分析,给出了全离散有限体积元格式,并通过详细的理论分析,得到了近似解与原问题真解的最优H^1模误差估计。     

求解多孔介质中可压缩混溶驱动问题的一类满足极大值原理的数值方法

多孔介质中可压缩混溶驱动问题是用非线性抛物型方程组来描述的．用Ｐｏｔｅｍｐａ格式求其数值解．证明了构造的求解方法满足极大值原理,从而可以保证饱和度的数值解在［０,１］范围内这一物理特性,同时还得到了解的收敛性．     

二维可压缩可混溶流体驱动问题的特征差分法

二维可压缩可混溶流体驱动问题的特征差分法顾海明（青岛化工学院计算机系）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＤＩＦＦＥＲＥＮＣＥＭＥＴＨＯＤＳＦＯＲＴＷＯ－ＰＨＡＳＥＦＬＯＷＩＮＰＯＲＯＵＳＭＥＤＩＡ￥ＧｕＨａｔ－ｉｎｉｎｇ（ＱｉｎｇｄａｏＩｎｓｔｉｔｕｔｅｏ...     

多孔介质中不可压缩流体的可混溶驱动问题的配置法

有了以上准备之后,我们给出有限元剖分,以h_c,h_p分别表示饱和度方程和压力方程的空间剖分步长,Δt_c,Δt_p分别表示饱和度方程和压力方程的时间剖分步长,且使:     

多孔介质中两相可压缩混溶流体驱动问题的交替方向法

多孔介质中二相可压缩流体驱动问题可以由一组非线性方程组来描述,它包含一个压力方程和一个浓度方程.文中对压力方程提出了稳定化校正格式,对浓度方程提出了二阶迎风交替方向差分格式,并结合双线性插值,给出了最优的L~2估计.结果表明该格式有一定的理论意义和实际应用价值.     

多孔介质中两相不可压缩可混溶驱动问题的动态混合元方法及其特征修正格式

许多依赖时间的问题涉及到局部化现象,如突出的前沿位置、激波、边界层等, 其位置随时间而变动.多孔介质中两相不可压缩可混溶驱动问题是一典型的、有代表性 的\"局部化现象\"问题,其数学模型为耦合非线性偏微分方程组的初边值问题.为减轻数 值解在局部前沿位置的数值振荡,提高解的精确性,本文给出了该问题的动态混合元格 式和沿特征线修正的动态混合元格式,证明了其收敛性,并给出了误差估计.     

多孔介质中可压缩可混溶驱动问题的共轭梯度迭代法的L2模误差估计

2006年3月 高等学校计算数学学报 1数学模型 多孔介质中可压缩可混溶驱动问题的模型是两个非线性抛物型方程:压力方程和饱 和度方程.Douglass和Roberts曾提出其数学模型并研究了半离散化方法[“一”}.袁益让对 此模型研究了特征一有限元方法[s]和差分法10]. 本人对可压缩可混溶驱动问题的模型曾研究了共扼梯度迭代解与原问题真解的最优 阶H‘模误差估计阁.其中饱和度方程的弥散项为一甲·(D(劝甲c)，而本文讨论的是D(司 情况下的尸模误差估计.就护模而言，对此模型目前尚未有人讨论过.从本文可看到， 由于饱和度方程中含有拭c)鬓这一项，…     

两相渗流驱动问题的体积有限元数值计算和分析

本文在一般的三角形剖分上对两相渗流驱动提出了全离散体积有限元 ,并分析了带有弥散项时格式的收敛性 ,得到H1 模的最优估计 .     

含弥散核废料污染问题特征有限元方法

研究可压缩核废料污染问题,考虑包含分子扩散和弥散的一般情形．从问题的原始数学模型出发,构造了特征线有限元全离散格式．对通常使用的构造误差方程的方法进行了改进,证明了格式的最佳L∞（J;H1（Ω））收敛性．     

油水驱动半正定问题迎风混合元格式及其数值分析

研究了不可压缩油水两相渗透流驱动问题.在扩散矩阵仅是半正定的假设条件下,提出了迎风混合元方法.混合元方法近似压力方程,饱和度方程的对流项用Godunov迎风格式来处理,扩散项则用推广的混合元来逼进,并推导出格式的误差估计.此种格式的优越性表现在两个方面:首先是饱和度方程的扩散矩阵仅是半正定的;二是摒弃了特征格式所限制的周期性条件,更适用于实际问题.     

Mixed Method for Compressible Miscible Displacement with Dispersion in Porous Media

Compressible miscible displacement of one fluid by another in porous media is modelled by a nonlinear parabolic system. A finite element procedure is introduced to approximate the concentration of one fluid and the pressure of the mixture. The concentration is treated by a Galerkin method while the pressure is treated by a parabolic mixed finite element method. The effect of dispersion, which is neglected in [1], is considered. Optimal order estimates in L2 are derived for the errors in the approximate solutions.     

A Hybrid WENO Method with Modified Ghost Fluid Method for Compressible Two-Medium Flow Problems

In this paper, we develop a simplified hybrid weighted essentially non-oscillatory (WENO) method combined with the modified ghost fluid method (MGFM) [31] to simulate the compressible two-medium flow problems. The MGFM can turn the two-medium flow problems into two single-medium cases by defining the ghost fluids state in terms of the predicted the interface state, which makes the material interface “invisible”. For the single medium flow case, we adapt between the linear upwind scheme and the WENO scheme automatically by identifying the regions of the extreme points for the reconstruction polynomial as same as the hybrid WENO scheme [55]. Instead of calculating their exact locations, we only need to know the regions of the extreme points based on the zero point existence theorem, which is simpler for implementation and saves computation time. Meanwhile, it still keeps the robustness and has high efficiency. Extensive numerical results for both one and two dimensional two-medium flow problems are performed to demonstrate the good performances of the proposed method.     

Hybrid Finite Element Method for Two-phase Miscible Displacement in Porous Media

ＨｙｂｒｉｄＦｉｎｉｔｅＥｌｅｍｅｎｔＭｅｔｈｏｄｆｏｒＴｗｏ┐ｐｈａｓｅＭｉｓｃｉｂｌｅＤｉｓｐｌａｃｅｍｅｎｔｉｎＰｏｒｏｕｓＭｅｄｉａ＊）ＬｉａｎｇＤｏｎｇ（梁栋）ＣｈｅｎｇＡｉｊｉｅ（程爱杰）（ＤｅｐａｒｔｍｅｎｔｏｆＭａｔｈｅｍａｔｉｃｓ，...     

A Steady Weak Solution of the Equations of Motion of a Viscous Incompressible Fluid through Porous Media in a Domain with a Non-Compact Boundary

We assume that Ω is a domain in ℝ² or in ℝ³ with a non-compact boundary, representing a generally inhomogeneous and anisotropic porous medium. We prove the weak solvability of the boundary-value problem, describing the steady motion of a viscous incompressible fluid in Ω. We impose no restriction on sizes of the velocity fluxes through unbounded components of the boundary of Ω. The proof is based on the construction of appropriate Galerkin approximations and study of their convergence. In Sect. 4, we provide several examples of concrete forms of Ω and prescribed velocity profiles on ∂Ω, when our main theorem can be applied.     

A Method for Numerical Simulation of Haline Convective Flows in Porous Media as Applied to Geology

Computational Mathematics and Mathematical Physics - A numerical code for simulating haline convective flows in porous media based on the finite difference method on a staggered nonuniform grid is...     

基于二次插值的非饱和土壤水流问题的特征差分方法及数值模拟

本文针对一类非饱和土壤水流问题,提出了基于二次插值的特征差分格式,得到了严谨的L2模误差估计.并作了数值试验,指明方法的有效性.