渗流力学的近况和展望

渗流力学研究流体在多孔介质的运动规律.渗流是指流体在多孔介质内的流动;孔隙介质、裂缝-孔隙介质以及各种类型的毛细管体系等均属多孔介质.渗流力学是流体力学与多孔介质理论和表面物理化学等学科交叉渗透产生的一个独立的学科领域,是多种工程技术的理论基础.     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

基于侵入混沌多项式法的随机多孔介质内顺磁性流体热磁对流不确定度量化

目前流体流动与传热问题的研究大都基于确定性工况条件,而现实流体流动与传热问题中存在着大量不确定性因素,计算流体力学的不确定性量化提供了一种理解流体物性、边界条件与初始条件等不确定性因素对模拟结果影响的能力.为揭示随机多孔介质内顺磁性流体热磁对流的传播规律与演化特征,本文发展了一种基于侵入式多项式混沌展开法的热磁对流不确...     

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型,介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展.通过线性稳定性理论,分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, DarcyBrinkman-Oldroyd以及Darcy-Brinkman-Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响.利用弱非线性分析方法,揭示失稳临界点附近热对流流动的分叉情况,以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式.采用数值模拟方法,研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的,而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞,最后发展为混沌状态.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

多孔介质-自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质-自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究.在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响.结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差.当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反.较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.     

基于孔隙弹性的水力压裂耦合方法与压裂状态参数数值研究

本文系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂FEM的数值实现方法。基本架构为：(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式; (3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式; (4)提出全流量加载的耦合分析方案。最后作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性。     

基于孔隙弹性耦合的水力压裂数值模拟研究

系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂的FEM数值实现方法.基本架构为:(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式;(3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式;(4)提出全流量加载的耦合分析方案.作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性.     

基于孔隙弹性的水力压裂耦合方法与数值模拟研究

本文系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂的FEM数值实现方法。基本架构为:(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式;(3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式;(4)提出全流量加载的耦合分析方案。最后作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性。     

Influence of heat transfer and fluid flow on crack growth in multilayered porous/dense materials using XFEM: Application to Solid Oxide Fuel Cell like material design

An advanced numerical model is developed to investigate the influence of heat transfer and fluid flow on crack propagation in multi-layered porous materials. The fluid flow, governed by the Navier–Stokes and Darcy’s law, is discretized with the nonconforming Crouzeix–Raviart (CR) finite element method. A combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods is used to solve the advection–diffusion heat transfer equation in the flow channel and in the fluid phase within the porous material. The crack is assumed to affect only the heat diffusion within the porous layer, therefore a time splitting technique is used to solve the heat transfer in the fluid and the solid phases separately. Thus, within the porous material, the crack induces a discontinuity of the temperature at the crack surfaces and a singularity of the flux at the crack tip. Conduction in the solid phase is solved using the eXtended Finite Element Method (XFEM) to better handle the discontinuities and singularities caused by the cracks. The XFEM is also used to solve the thermo-mechanical problem and to track the crack propagation. The multi-physics model is implemented then validated for the transient regime, this necessitated a post processing treatment in which, the stress intensity factors (SIF) are computed for each time step. The SIFs are then used in the crack propagation criterion and the crack orientation angle. The methodology seems to be robust accurate and the computational cost is reduced thanks to the XFEM.     

三维椭球形缺陷启裂载荷随倾角和压缩度变化的规律

本文基于扩展有限元法模拟不同尺寸和倾角三维椭球形缺陷的启裂行为,随后引入压缩度描述缺陷扁平程度,最后讨论缺陷倾角和压缩度对其启裂载荷的影响。结果表明：当缺陷形状呈扁椭球形时,随倾角(缺陷倾向与力方向的夹角)增加,启裂载荷先减小后增加。当缺陷形状呈近球形时,启裂载荷随倾角增加变化轻微。当缺陷倾角为0°时,其启裂载荷随压缩度增加先减小后增加。当倾角为15°,30°,45°,60°和75°时,启裂载荷随压缩度增加而减小。

     

板壳式换热器波纹倾角对换热和阻力性能影响

研究了波纹倾角（α=45°,60°,75°）对板壳式换热器单流道内单相流动与换热过程的影响,分析了波纹流道内的速度分布、湍动能分布、压力分布和温度分布;基于范宁摩擦阻力因子f和Nu与Re的关系,提出了板壳式换热器不同波纹倾角下换热特性和流动阻力特性的预测关联式;通过计算不同波纹倾角下的综合性能评价指标（performance evaluation criteria,PEC）和面积质量因子（j/f）综合评价了流动与换热性能。结果表明：波纹倾角是影响圆形板间流体流动形态的因素之一,随着波纹倾角增大,会出现“十字交叉流”向“曲折流”转变;所提出的关联式能够很好地预测板壳式换热器内阻力与换热性能,阻力特性偏差在±14%范围内,换热特性偏差在±7%范围内;α为45° 时,j/f较大,流道内阻力相对小;α为60° 时,PEC较大,流道内换热性能相对强。

     

Nonlinear convection in a non-Darcy porous medium

In this paper, the natural convection in a non-Darcy porous medium is studied using a temperature-concentration-dependent density relation. The effect of the two parameters responsible for the nonlinear convection is analyzed for different values of the inertial parameter, dispersion parameters, Rayleigh number, Lewis number, Soret number, and Dufour number. In the aiding buoyancy, the tangential velocity increases steeply with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter when the inertial effect is zero. However, when the inertial effect is non-zero, the effect of the nonlinear temperature parameter and the nonlinear concentration parameter on the tangential velocity is marginal. The concentration distribution varies appreciably and spreads in different ranges for different values of the double dispersion parameters, the inertial effect parameter, and also for the parameters which control the nonlinear temperature and the nonlinear concentration. Heat and mass transfer varies extensively with an increase in the nonlinear temperature parameter and the nonlinear concentration parameter depending on Dacry and non-Darcy porous media. The variation in heat and mass transfer when all the effects, i.e., the inertial effect, double dispersion ef- fects, and Soret and Dufour effects, are simultaneously zero and non-zero. The combined effects of the nonlinear temperature parameter, the nonlinear concentration parameter and buoyancy are analyzed. The effect of the nonlinear temperature parameter and the nonlinear concentration parameter and also the cross diffusion effects on heat and mass transfer are observed to be more in Darcy porous media compared with those in non- Darcy porous media. In the opposing buoyancy, the effect of the temperature parameter is to increase the heat and mass transfer rate, whereas that of the concentration parameter is to decrease.     

应变率对含裂隙红砂岩裂纹扩展模式及破碎特征的影响

以含不同倾角预制裂纹的长方形板状红砂岩为研究对象,采用分离式霍普金森压杆沿试样宽度方向施加冲击荷载,使用高速摄像机记录裂纹扩展过程,获得试样的裂纹路径特征以及动态压缩强度和动态弹性模量,利用筛分统计法分析试样碎块分布特征,结合分形理论定量描述试样破碎程度及特点,探讨中应变率条件下含裂隙试样裂纹扩展模式与动态力学性质和破碎程度的相互关系。研究结果表明,应变率较高时试样会更多地出现远场裂纹和离层裂纹,并且相比相关低应变率实验结果,中应变率范围内试样破坏模式及裂纹分布情况随应变率的变化规律是不同的。随着应变率的提高,试样大体上从1条拉伸裂纹的临界破坏演变成X形剪切裂纹为主的复杂裂隙网络,并且不同角度预制裂隙对于这种裂纹扩展模式的演变有重要影响。在预制裂纹倾角一定的情况下,岩样动态压缩强度和动态弹性模量表现出明显的应变率效应,不同角度预制裂纹对于试样的应变率敏感性有显著影响。随裂纹倾角的增大,试样的动态强度、动态弹性模量和分形维数表现出的变化趋势具有一定的相似性,大体呈现先减小后增大的趋势,裂纹倾角为45°的试样的动态压缩强度、动态弹性模量和分形维数均为最小。随应变率的升高,不同预制裂纹倾角的试样碎块分布更加分散,应变率越高,预制裂纹倾角对于岩石冲击破碎程度、分形维数的影响越显著。

     

Steady and Unsteady Heat Transfer in a Channel Partially Filled with Porous Media Under Thermal Non-Equilibrium Condition

Steady and pulsatile flow and heat transfer in a channel lined with two porous layers subject to constant wall heat flux under local thermal non-equilibrium (LTNE) condition is numerically investigated. To do this, a physical boundary condition in the interface of porous media and clear region of the channel is derived. The objective of this work is, first, to assess the effects of local solid-to-fluid heat transfer (a criterion indicating on departure from local thermal equilibrium (LTE) condition), solid-to-fluid thermal conductivity ratio and porous layer thickness on convective heat transfer in steady condition inside a channel partially filled with porous media; second, to examine the impact of pulsatile flow on heat transfer in the same channel. The effects of LTNE condition and thermal conductivity ratio in pulsatile flow are also briefly discussed. It is observed that Nusselt number inside the channel increases when the problem is tending to LTE condition. Therefore, careless consideration of LTE may lead to overestimation of heat transfer. Solid-to-fluid thermal conductivity ratio is also shown to enhance heat transfer in constant porous media thickness. It is also revealed that an increase in the amplitude of pulsation may result in enhancement of Nusselt number, while Nusselt number has a minimum in a certain frequency for each value of amplitude.     

多层平行裂隙型多孔介质通道内流体流动传热特性

基于Brinkman-extended Darcy模型和局部热平衡模型,对多层平行裂隙型多孔介质通道内的流动传热特性进行研究.获得了多层平行裂隙型多孔介质通道内各区域的速度场、温度场、摩擦系数及努塞尔数解析解,并分析了裂隙层数、达西数、空心率、有效热导率之比等对通道内流动传热特性的影响.结果表明:达西数较小时,通道多孔介质层内会出现不随高度变化的达西速度,此达西速度会随裂隙层数的增加而增大,但却不受各裂隙层下多孔介质层位置变化的影响.增加裂隙层数会减弱空心率对压降的影响,会使通道内流体压降升高,但升高程度会逐渐降低.增大热导率之比或减小空心率会使多裂隙通道内出现阶梯式温度分布,而在较小热导率之比或较大空心率时多裂隙情况下的温度分布曲线会趋于一致.此外,当热导率之比较小时,多层裂隙通道内的传热效果在任何空心率下都要优于单裂隙情况,当热导率之比较大时,存在临界空心率使各裂隙层数通道内的传热效果相同,且多裂隙通道内继续增加裂隙层数对传热强度影响不大.     

爆破中双线型聚能药包最佳成缝角度

为探究聚能张开角对双线型聚能药包结构炸药有效利用率和聚能效应的影响,通过瞬时爆轰假说理论对有效聚能炸药边界方程进行推导,分析不同聚能张开角聚能装药结构炸药的有效利用率;通过水泥砂浆物理模型试验,研究不同聚能张开角预裂孔成缝规律;采用LS-DYNA数值模拟软件,建立不同聚能张开角数值模型,揭示不同聚能张开角的双线型聚能结构药包射流的侵彻过程。研究结果表明聚能张开角为75°时,炸药产生聚能效应的有效利用率最大;聚能结构药包聚能槽张开角为75°时,预裂孔成缝效果明显优于聚能槽张开角为60°的聚能结构药包,沿聚能槽方向应力集中效应和侵彻深度最佳,炮孔壁上岩石单元最先达到应力峰值。针对聚能张开角为75°的双线型聚能结构药包开展了不同岩性预裂爆破现场试验,板岩和白云岩两种不同岩性,在孔距增大20%的条件下,双线型聚能预裂爆破效果优于常规预裂爆破。     

Experimental Investigations of Natural Convection in Wire-Woven Bulk Kagome

Natural convection in an open and highly anisotropic cellular material, wire-woven bulk Kagome (WBK) is experimentally characterized. A series of pressure drop and heat transfer experiments are conducted. In particular, effect of inclination angle on the heat transfer rate for porous-cored heat sinks (WBK and metallic foams), as well as smooth plate was experimentally investigated. There exists an optimal inclination angle for smooth plate and porous-cored heat sinks in association with the maximum heat transfer rate (Nusselt number). An advance in optimal inclinational angle was observed in porous media during the full inclination angle range from \(0^{\circ }\) (horizontal) to \(90^{\circ }\) (vertical) due to the enhanced lateral conduction and a further advance was found in WBK specimen. A strong anisotropic heat transfer existed in terms of different orientations of WBK specimen: O-B orientation with a bigger cross-sectional surface area density (blockage ratio) reveals a higher heat transfer rate than O-A orientation. Further, in comparison with isotropic metallic foams with a given porosity, the WBK which is positioned at vertical orientation can dissipate more heat than the foams due to the higher permeability resulting from less pressure drop in the WBK formed by an assembly of cylindrical wires.     

Numerical and experimental investigation of convective drying in unsaturated porous media with bound water

The heat and mass transfer in an unsaturated wet cylindrical porous bed packed with quartz particles was investigated theoretically for relatively low convective drying rates. Local thermodynamic equilibrium was assumed in the mathematical model describing the multi-phase flow in the unsaturated porous media using the energy and mass conservation equations to describe the heat and mass transfer during the drying. The drying model included convection and capillary transport of the free water, diffusion of bound water, and convection and diffusion of the gas. The numerical results indicated that the drying process could be divided into three periods, the temperature rise period, the constant drying rate period and the decreasing drying rate period. The numerical results agreed well with the experimental data verifying that the mathematical model can evaluate the drying performance of porous media for low drying rates. The effects of drying conditions such as the ambient temperature, the relative humidity, and the velocity of the drying air, on the drying process were evaluated by numerical solution.