Liquid Flow in Partially Saturated Porous Media

We consider a picture for the filtration of a liquid in a partiallysaturated porous medium, leading to a two-phase one-dimensionalfree boundary problem of the following type: The liquid pressuresatisfies an elliptic equation in the saturated region and anon-linear parabolic equation in the unsaturated region, whilepressure and velocity are continuous across the interface. This scheme reduces to the study of the non-linear parabolicfree boundary problem in the unsaturated phase with cauchy dataprescribed on the free boundary, for such a problem existence,uniqueness and continuous dependence theorems are proved.     

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

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

饱和多孔介质中骨架的应变局部化萌生条件

应用饱和多孔介质控制方程和Liapunov稳定理论,导出了固相应力和有效应力描述的多孔介质骨架应变局部化的萌生条件．不同应力形式表达的多孔介质基体的控制方程,相应的应变局部化萌生条件的表达形式也不尽相同,其原因源于骨架本构中固液两相之间相互作用的不同描述．应用得出的Terzaghi有效应力描述的应变局部化萌生条件,可以理论解释多孔介质中固、液两相不同相对运动出现的破坏方式,如管涌、滑坡和泥石流．应用简单算例说明了应变局部化条件的具体实施方法．     

多孔介质中两相不可压缩不易混溶渗流问题的特征配置法

多孔介质中两相不可压缩不易混溶渗流问题是非线性偏微分方程的耦合系统,其中压力方程是椭圆的用配置法逼近,而饱和度方程是对流占优的抛物方程,用特征配置法来逼近,并且证明了数值解的存在唯一性,最后得到了最优的误差估计.     

饱和多孔介质的超粘弹性本构理论研究

为了建立能考虑固体材料、多孔固体与流体可逆和不可逆变形的饱和多孔介质超粘弹性理论,以多孔固相为参考构型,以有效应力、材料真实应力和流相真实孔压作为状态变量,结合混合物均匀化响应原理获得各项均符合热力学功共轭特征的饱和多孔介质能量平衡方程,根据非平衡热力学熵分解理论求得熵流和熵产.结果表明,超弹塑性理论是该理论的一个特例;多孔固体的总变形可分为固相间隙和材料变形两部分,间隙应变与Terzaghi有效应力构成功共轭对,材料应变与材料真实应力构成功共轭对.饱和多孔介质的自由能可分为固相和流相两部分.当固相间隙和材料变形解耦时,固相所含的自由能又可分为间隙和材料两部分.证明了Skempton有效应力不是饱和多孔介质的基本应力状态变量.     

多层非线性渗流耦合系统的特征分数步差分方法

对多层非线性渗流耦合系统提出适合并行计算的特征分数步差分格式, 利用变分形式、能量方法、粗细网格配套、分片双二次插值、差分算子乘积交换性、高阶差分算子的分解、先验估计的理论和技巧, 得到收敛性的最佳阶的l₂误差估计. 该方法已成功的应用到多层油资源评估的生产实际中.     

含液体非弹性多孔介质中波传播过程的失稳与逸散性

在基于Biot理论的饱和-非饱和多孔介质的动力-渗流模型中计及流固惯性耦合效应。对单轴应变的一维情况讨论了饱和和非饱和多孔介质中波传播过程的驻值失稳和逸散性,分析了流固粘性耦合,流固惯性耦合,流固混合体的压缩性,孔隙水饱和度,及固体骨架在高应变速率下材料粘弹塑性本构行为等因素的影响。该工作将对克服饱和与非饱和多孔介质在强动荷载下波传播过程的数值求解困难提供理论上的依据和启示。     

Theory and Application of Characteristic Finite Element Domain Decomposition Procedures for Coupled System of Dynamics of Fluids in Porous Media

For a coupled system of multiplayer dynamics of fluids in porous media,the characteristic finiteelement domain decomposition procedures applicable to parallel arithmetic are put forward.Techniques suchas calculus of variations,domain decomposition,characteristic method,negative norm estimate,energy methodand the theory of prior estimates are adopted.Optimal order estimates in L~2 norm are derived for the error inthe approximate solution.     

On the Description of Partially Saturated Soils by the Theory of Porous Media

In this contribution, a multi‐phase soil model based on the Theory of Porous Media (TPM) is presented. The model is fully coupled in the following constitutive phases: An elasto‐plastic or elasto‐viscoplastic solid skeleton, a materially incompressible pore‐liquid (water) and a materially compressible pore‐gas (air). The interaction of the solid skeleton and the pore‐fluids is specified by a capillary pressure‐saturation relation, whereas the mobilities of the fluid phases in the pore‐space of the solid skeleton are described by the so‐called relative permeabilities. Finally, a gravity governed initial‐boundary‐value problem solved by the FE method is presented.     

Wave Attenuation in Effective Models Describing Porous and Fractured Media Saturated by a Fluid

Wave attenuation is introduced in the effective model of media that consists of alternating elastic and fluid layers. This attenuation is due to the friction on the boundaries between elastic and fluid layers and is described by additional terms in equations of the effective model. An investigation of these equations allows one to derive expressions of the attenuation coefficients for every body wave propagating along the layers. Bibliography: 9 titles.__________Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 297, 2003, pp. 216–229.     

非线性多孔介质流动模型问题的分歧

本文讨论一类多孔介质流动模型非线性边值问题．对这类问题以及广泛的一类算子方程，我们得到了渐近歧点的存在性     

含有溶质的流体在两层多孔介质中的渗流问题

本文讨论含有溶质的流体在两层多孔介质中的渗流问题，即（θ（x，U）t=（K（x，U）Ux-K（x，U））x，（x，t）∈GT，（θ（x，U）V（x，t）t=（DθVx）x-（V（KUx-K））x，（x，t）∈GT，U（x，0）=U0（x），V（x，0）=V0（x），0≤x≤2，U（0，t）-h0（t），U（2，t）=h2（t），0≤t≤T，V（0，t）=g0（t），V（2，t）=g2（t），0≤t≤T。其中θ（x,U)=θ1(x,U)，当（x,t)∈D1=｛0≤x≤1,0≤t≤T｝；θ（x,U)=θ2(x,U)当（x,t）∈D2＋1｛1＜x≤2，0≤t≤T｝。K（x,U)=K1(x,U)当（x,t)∈D1;K(x,U)=K2(x,U)，当(x,t)∈D2。θi，Ki分别是Di上的介质含水率及水力传导率，V是溶质的浓度，此外还要求U，V，K（x,U)(Ux-1)及DθVx V(KUx-K)在x=1连续。     

Convection Problems in Anisotropic Porous Media with Nonhomogeneous Porosity and Thermal Diffusivity

Convection problem in anisotropic and inhomogeneous porous media has been analysed. In particular, the effect of variable permeability and thermal diffusivity with respect to the vertical direction, has been studied. Linear and nonlinear stability analysis of the conduction solution have been performed.     

A Parallel Implementation of the Algebraic Multigrid Method for Solving Problems in Dynamics of Viscous Incompressible Fluid

An algorithm for improving the scalability of the multigrid method used for solving the system of difference equations obtained by the finite volume discretization of the Navier–Stokes equations on unstructured grids with an arbitrary cell topology is proposed. It is based on the cascade assembly of the global level; the cascade procedure gradually decreases the number of processors involved in the computations. Specific features of the proposed approach are described, and the results of solving benchmark problems in the dynamics of viscous incompressible fluid are discussed; the scalability and efficiency of the proposed method are estimated. The advantages of using the global level in the parallel implementation of the multigrid method which sometimes makes it possible to speed up the computations by several fold.     

三维非线性对流扩散问题的数值方法在渗流力学的应用

对三维非线性对流扩散问题提出一类适合并行计算的二阶迎风分数步差分格式,采用分数步技术,将三维问题化为连续解3个一维问题计算．利用变分形式、能量方法、差分算子乘积交换性、高阶差分算子的分解、微分方程先验估计的理论和技巧,得到收敛性的最佳阶的误差估计．该方法已成功的应用油资源运移聚集渗流力学数值模拟计算、海水入侵预测和防治的工程实践中．     

不可压饱和多孔弹性梁的大挠度非线性数学模型

在孔隙流体仅存在沿梁轴线方向扩散的假定下,建立了微观不可压饱和多孔弹性梁大挠度问题的非线性数学模型．利用Galerkin截断法,研究了固定端不可渗透、自由端可渗透的饱和多孔弹性悬臂梁在自由端突加集中载荷作用下的非线性弯曲,得到了梁骨架的挠度、弯矩以及孔隙流体压力等效力偶等的时间响应和沿轴线的分布．比较了大挠度非线性和小挠度线性理论的结果,揭示了两者间的差异．研究发现大挠度理论的结果小于相应的小挠度理论结果,并且,大挠度理论的结果趋于其稳态值的时间小于相应的小挠度理论结果趋于其稳态值的时间．     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.