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

在有界区域上多孔介质中可压缩可混溶的油、水两相驱动问题是由非线性偏微分方程组的初、边值问题所决定.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 High-Order Kernel-Free Boundary Integral Method for Incompressible Flow Equations in Two Space Dimensions

This paper presents a fourth-order kernel-free boundary integral method for the time-dependent, incompressible Stokes and Navier-Stokes equations defined on irregular bounded domains. By the stream function-vorticity formulation, the incompressible flow equations are interpreted as vorticity evolution equations. Time discretization methods for the evolution equations lead to a modified Helmholtz equation for the vorticity, or alternatively, a modified biharmonic equation for the stream function with two clamped boundary conditions. The resulting fourth-order elliptic boundary value problem is solved by a fourth-order kernel-free boundary integral method, with which integrals in the reformulated boundary integral equation are evaluated by solving corresponding equivalent interface problems, regardless of the exact expression of the involved Green's function. To solve the unsteady Stokes equations, a four-stage composite backward differential formula of the same order accuracy is employed for time integration. For the Navier-Stokes equations, a three-stage third-order semi-implicit Runge-Kutta method is utilized to guarantee the global numerical solution has at least third-order convergence rate. Numerical results for the unsteady Stokes equations and the Navier-Stokes equations are presented to validate efficiency and accuracy of the proposed method.     

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.     

Mixed Laguerre–Legendre Spectral Method for Incompressible Flow in an Infinite Strip

This paper concerns the mixed Laguerre–Legendre spectral approximation and its application to numerical simulation of incompressible flow in an infinite strip. Some approximation results in weighted Sobolev spaces are given. A Laguerre–Legendre spectral scheme for the stream function form of Navier–Stokes equations is constructed. The stability and the convergence of the proposed scheme are proved. The numerical experiments show the high accuracy of this method. The main techniques used in this paper are also applicable to other nonlinear partial differential equations in an infinite strip.     

A Novel Localized Meshless Collocation Method for Numerical Simulation of Flume Dynamic Characteristics北大核心CSCD

局部边界节点法是一种基于非奇异半解析基函数和移动最小二乘原理的新型无网格配点技术,该方法把每个节点处的未知变量表示为该点对应的局部子域内节点处物理量的线性组合,该文基于局部边界节点法对数值波浪水槽进行了研究。首先,通过基准算例确定了Laplace算子非奇异半解析基函数的合理形状参数值。进一步,基于合理的参数选取,用较少的离散节点即可成功模拟波浪传播行为,将得到的数值结果与其他文献数值结果比较,可以发现局部边界节点法用更少的局部点即可得到较好的数值结果。最后,以保护近海岸建筑物为目标,模拟了水下防波堤对波浪传播的影响。结果表明,当波浪与梯形防波堤发生作用后,波峰变得比较陡峭,而波谷变得相对比较平坦,为近海岸防波堤的相关研究和设计提供了数值参考。     

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

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