多孔介质中卡森流体的分形分析

研究了非牛顿流体中的卡森流体在多孔介质中的流动特性.基于服从分形分布的弯曲毛细管束模型,运用分形几何理论推导出了该流体在多孔介质中流动的流量、流速、启动压力梯度和有效渗透率的分形解析解.模型中的每一个参数都有明确的物理意义,它将卡森流体在多孔介质中的流动特性与多孔介质的微结构参数有机联系起来.文中给出了卡森流体的流速、启动压力梯度和有效渗透率随着各影响因素的变化趋势,并进行了讨论.所得分形模型可以更深刻地理解卡森流体在多孔介质中流动的内在物理机理. 关键词： 多孔介质 卡森流体 分形     

多孔介质模拟用于气膜冷气流动研究的尝试

本文采用多孔介质模拟与带真实气膜孔流动模拟方法,对平板气膜孔模型在低速条件下进行对比研究,以确定适用于涡轮冷却叶片气膜冷气与主流流动影响研究的多孔介质模拟方法。结果表明,采用流场参数进行多孔介质关键参数迭代计算,可以简化多孔介质方法在气膜冷气流动模拟中的应用。多孔介质模拟方法可以在较小的网格量下,获得较为详细的流场信息,对气膜冷气的流量、冷气射流速度、掺混平均的速度剖面、掺混损失、充分掺混区的平均气膜冷却效率、绝热壁温等进行较为准确的模拟。     

基于Biot-喷射流统一模型Maxwell流体饱和孔隙介质中的弹性波

推广Biot-Tsiklauri声学模型的同时借鉴Dvorkin和Nur的工作,建立了具有任意孔径分布并顾及喷射流动机制的非牛顿流体饱和孔隙介质声学模型,研究了非牛顿流体(Maxwell流体)饱和孔隙介质中的弹性波的衰减和频散特性.着重讨论充孔隙Maxwell流体的非牛顿流效应对弹性波的频散和衰减的影响.研究表明,饱和流体的非牛顿流效应和喷射流动机制均是引起弹性波波频散和衰减的重要因素.依据非牛顿流体(Maxwell流体)饱和各向同性孔隙介质的Biot-喷射流声学模型,喷射流动只影响纵波的频散和衰减,而饱和流体的非牛顿流效应不仅影响纵波,而且还影响横波的频散和衰减.     

多相非牛顿流体驱替过程的格子Boltzmann模拟

将非牛顿流体模型引入多相格子Boltzmann方法（LBM）,研究多相非牛顿流体的驱替过程.首先以牛顿流体驱替过程作对照,比较非牛顿流体的驱替特性.然后基于简单结构模型,分析不同界面张力下,非牛顿流体的驱替能力和驱替过程中被驱替相的形态变化规律,揭示界面张力的影响机理.结果表明：界面张力在驱替过程中起阻碍作用,其它参数相同情况下,界面张力越低,流体的驱替能力越强.最后,结合四参数随机生长（QSGS）结构生成方法,模拟不同压差、不同界面张力以及不同润湿性条件下多孔复杂结构内的非牛顿多相驱替过程,分析这些参数对过程的影响.     

随机四参数法生成多孔介质及渗流模拟

本文基于多孔介质内流动的广泛应用背景,模拟了多孔介质渗流.采用随机四参数生长法生成了多孔介质模型,并对模型中的微小固相进行了过滤处理,将随机四参数生长法生成的多孔介质结合格子玻尔兹曼方法D2Q9模型模拟了压力驱动的多孔介质中的渗流.通过将模拟结果与著名的CK理论公式的计算结果进行比较,发现在孔隙率高时误差较小(cp=0...     

多孔介质模型在管壳式换热器数值模拟中的应用

本文详细讨论了多孔介质模型在管壳式换热器数值模拟中的应用,开发了一套能自动生成多孔介质特性参数的通用程序。该程序主要基于三维交错网格及SIMPLE算法,然后运用该模型,采用改进的κ-ε模型和壁面函数法,对换热器壳侧的湍流流动进行了数值模拟。计算结果与换热器冷态实验结果符合良好,表明该模型和计算方法是切实可行的。     

生物过滤器中非均匀性流动的数值研究

生物过滤技术因其具有有效性、低成本和环境友好等优点引起了人们的广泛关注.该技术主要通过生物过滤器去除含有H2S等废气的有毒有害气体.运用格子Boltzmann方法对三种生物过滤器模型中多孔介质的非均匀性流动进行了数值模拟.数值模拟结果表明,多孔介质的性质和进口流动条件对临界Rayleigh数有显著影响,临界Rayleigh数随着多孔介质的孔隙度和Darcy数的增大而逐渐变小,并随着进口Reynolds数的增大而逐渐变大.所得结果可望为生物过滤器的优化设计提供一个合理的理论依据.     

基于几何重建的气液两相流孔隙尺度模拟

本文基于二维micro CT扫描图像进行数值重构得到三维多孔介质几何模型,采用Shan-Chen多相格子Boltzmann方法(LBM)建立适用不规则几何边界区域内气液两相流动数值模型,在此基础上对复杂多孔介质内的两相流动过程进行孔隙尺度模拟。结果表明,气-液两相流动受限于多孔介质的复杂孔隙结构,气液表面张力、壁面润湿特性以及进出口压差对两相流动特性将产生明显影响。     

含内热源多孔介质–自由流通道内非达西强制对流

基于多孔介质局部非热平衡模型,对考虑内热源条件下的多孔介质–自由流耦合通道内非达西对流换热特性进行研究。多孔介质区内流体运动方程采用Darcy-Brinkman-Forchheimer模型,利用有限差分法获得通道内各区域流体运动速度、流固相温度分布及努塞尔数,并进一步分析了相关参数对流体流动传热的影响。结果表明：在文中研究参数下,惯性参数F_f对通道内各区域流体速度及温度分布的影响仅在达西数Da大于10^-3时需要考虑;增加F_f或降低固相内热源W_s绝对值会使流固两相温差减小,且改变固相内热源换热方向会使多孔介质区内流固相发生温度分岔现象;固相内热源对通道换热效果影响较大且更为复杂,不同惯性参数F_f下,考虑W_s时可能会使Nu出现奇异点。     

多孔介质强化池内沸腾换热特性

为探讨多孔介质对强化池内沸腾传热的机理,在38×3mm的容器内上,应用多相流VOF模型,耦合多孔介质流动与传热模型,通过添加用户自定义程序实现了对光滑和多孔加热壁面的池内沸腾传热过程进行数值模拟。数值模拟结果与实验结果相符合,所建立的数学模型为多孔介质的沸腾传热研究提供了依据。     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

Elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media

This paper studies the elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media with the nonzero boundary slip velocity for pore size distribution. The coefficient bF_m(ω) that measures the deviation from Poiseuille flow friction in such media is presented. Based on this coefficient, we investigate the properties of elastic waves by calculating their phase velocities and attenuation coefficients as functions of frequency and the behaviour of the dynamic permeability. The study shows that the pore size distribution removes oscillations in all physical quantities in the non-Newtonian regime. Consideration of the nonzero boundary slip effect in non-Newtonian (Maxwell) fluid-saturated porous media results in (a) an overall increase of the dynamic permeability, (b) an increase of phase velocities of fast Biot waves and shear waves except in the low frequency domain and an overall increase of phase velocity of slow Biot waves and (c) an overall increase of the attenuation of three Biot waves in the intermediate frequency domain except in the deeply non-Newtonian regime. The study also shows that the attenuation coefficient of slow Biot waves is small in the deeply non-Newtonian regime at higher frequency, which encourages us to detect slow Biot waves in oil-saturated porous rock.     

确定多孔介质流动参数的格子Boltzmann方法

本文采用D2Q9不可压缩Lattice-Boltzmann模型模拟了多孔介质中的流场,该模型从介观层次上描述了多孔介质孔隙中的流动特性,验证了Darcy定律,给出了一种计算渗透率的简便方法,本文还运用图像处理方法来构造多孔介质,该方法构造的二维多孔介质能较好地体现实际多孔介质的几何特征。     

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效，用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。     

Elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media

Abstract

This paper studies the elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media with the nonzero boundary slip velocity for pore size distribution. The coefficient bF _m (ω) that measures the deviation from Poiseuille flow friction in such media is presented. Based on this coefficient, we investigate the properties of elastic waves by calculating their phase velocities and attenuation coefficients as functions of frequency and the behaviour of the dynamic permeability. The study shows that the pore size distribution removes oscillations in all physical quantities in the non-Newtonian regime. Consideration of the nonzero boundary slip effect in non-Newtonian (Maxwell) fluid-saturated porous media results in (a) an overall increase of the dynamic permeability, (b) an increase of phase velocities of fast Biot waves and shear waves except in the low frequency domain and an overall increase of phase velocity of slow Biot waves and (c) an overall increase of the attenuation of three Biot waves in the intermediate frequency domain except in the deeply non-Newtonian regime. The study also shows that the attenuation coefficient of slow Biot waves is small in the deeply non-Newtonian regime at higher frequency, which encourages us to detect slow Biot waves in oil-saturated porous rock.     

Analysis of the flow of non-Newtonian visco-elastic fluids in fractal reservoir with the fractional derivative

In both the oil reservoir engineering and seepage flow mechanics, heavy oil with relaxation property shows non-Newtonian rheological characteristics. The relationship between shear rate g& and shear stress t is nonlinear. Because of the relaxation phenomena of heavy oil flow in porous media, the equation of motion can be written as[1] 2,rrvpqkppqtrrtll秏骣+=-+琪抖桫 (1) where lv and lp are velocity relaxation and pressure retardation times. For most porous media, the above motion equation (1)…     

变形双重介质分形油藏非达西流动分析

考虑与实际生产相符的介质的双重特性和分形特征,并考虑介质的变形,引入双重分形介质渗透率模数,建立应力敏感地层双重分形介质系统的流动方程.采用Douglas-Jones预估-校正法获得无限大地层定产量生产时变形双重分形介质模型的数值解和无限大地层定压生产时分形介质双孔模型的数值解,作出了典型的压力曲线图版.     

On the analytic solution of the steady flow of a fourth grade fluid

The steady flow of a fourth grade fluid is a problem belonging to non-Newtonian fluid mechanics and deserves to be more widely studied than it has been to date. In the non-linear regime the literature is scarce. We develop a formulation suitable for solution of hydrodynamic equation containing non-linear rheological effects of fourth grade fluids. The homotopy analysis method (HAM) is used to investigate the flow of a fourth grade fluid past a porous plate. Explicit analytic solution is given. The non-linear effects on the velocity distribution is shown and discussed. Comparison of the present analysis is also made with the existing results in the literature.     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

Influence of system size and solvent flow on the distribution of wormlike micelles in a contraction-expansion geometry

Viscoelastic wormlike micelles are formed by surfactants assembling into elongated cylindrical structures. These structures respond to flow by aligning, breaking and reforming. Their response to the complex flow fields encountered in porous media is particularly rich. Here we use a realistic mesoscopic Brownian Dynamics model to investigate the flow of a viscoelastic surfactant (VES) fluid through individual pores idealized as a step expansion-contraction of size around one micron. In a previous study, we assumed the flow field to be Newtonian. Here we extend the work to include the non-Newtonian flow field previously obtained by experiment. The size of the simulations is also increased so that the pore is much larger than the radius of gyration of the micelles. For the non-Newtonian flow field at the higher flow rates in relatively large pores, the density of the micelles becomes markedly non-uniform. In this case, we find that the density in the large, slowly moving entry corner regions is substantially increased.