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

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

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

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

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

藏东南波密地区岩石风化速率及其影响因素分析

根据解剖学发现的人体小腿骨间膜血管和胶原纤维有序排列的特点,建立三维组织间隙液渗 流模型,控制方程用Brinkman方程和连续性方程,使用Fluent软件进行数值模拟. 计算结 果显示,组织液在组织间隙中流动的方向总体是与平行毛细血管的方向相一致的;组织间隙 多孔率会影响速度的分布,当多孔率增加时,空间的速度趋于均匀,从而使得最大速度变小, 正常生理范围内多孔率的改变对于组织液流动的影响较小. 但当多孔率很小时(组织纤维化), 会大大影响组织液流动的均匀性;组织间隙平行胶原的存在,会使得组织间隙速度在空间上 的分布趋于均匀;另外,间隙流场随血压的增加而增加,随组织压的增加而降低,这与临床 和经络研究的发现相一致.     

基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究

认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键. 微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率. 本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响. 结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.      

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

底部热加载下两层多孔介质热流耦合现象研究

底部热加载条件下非线性流动传热问题的求解一直是多孔介质研究领域的难题。采用实验测试及数值模拟的方法,对底部热加载方式下两层多孔介质内热流耦合对流传热解的特性进行了研究。研究结果表明:在小瑞利数工况下,两层多孔介质内骨架内材料比例对温度场分布和非线性特性有重要影响;确定了非线性分叉、震荡解的特征值,并计算出了不同材料比例下的分叉及震荡所对应的临界瑞利数;通过实验验证了底部热加载时两层多孔介质内温度随时间呈非稳态震荡变化。结论可为实现热流分层、局部削弱或强化传热提供参考。     

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

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

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究

认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键.微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率.本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响.结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.     

Effect of an Inserted Porous Layer Located at a Wall of a Parallel Plate Channel on Forced Convection Heat Transfer

A theoretical study is performed on heat and fluid flow in partially porous medium filled parallel plate channel. A uniform symmetrical heat flux is imposed onto the boundaries of the channel partially filled with porous medium. The dimensional forms of the governing equations are solved numerically for different permeability and effective thermal conductivity ratios. Then, the governing equations are made dimensionless and solved analytically. The results of two approaches are compared and an excellent agreement is observed, indicating correctness of the both solutions. An overall Nusselt number is defined based on overall thermal conductivity and difference between the average temperature of walls and mean temperature to compare heat transfer in different channels with different porous layer thickness, Darcy number, and thermal conductivity ratio. Moreover, individual Nusselt numbers for upper and lower walls are also defined and obtained. The obtained results show that the maximum overall Nusselt number is achieved for thermal conductivity ratio of 1. At specific values of Darcy number and thermal conductivity ratio, individual Nusselt numbers approach to infinity since the value of wall temperatures approaches to mean temperature.     

Fully Developed Forced Convection in a Parallel Plate Channel with a Centered Porous Layer

In this study, fully developed heat and fluid flow in a parallel plate channel partially filled with porous layer is analyzed both analytically and numerically. The porous layer is located at the center of the channel and uniform heat flux is applied at the walls. The heat and fluid flow equations for clear fluid and porous regions are separately solved. Continues shear stress and heat flux conditions at the interface are used to determine the interface velocity and temperature. The velocity and temperature profiles in the channel for different values of Darcy number, thermal conductivity ratio, and porous layer thickness are plotted and discussed. The values of Nusselt number and friction factor of a fully clear fluid channel (Nu _cl = 4.12 and fRe _cl = 24) are used to define heat transfer increment ratio (e_th = Nu_p/Nu_cl)({\varepsilon _{\rm th} =Nu_{\rm p}/Nu_{\rm cl})} and pressure drop increment ratio (e_p = fRe_p/fRe_cl )({\varepsilon_{\rm p} =fRe_{\rm p}/fRe_{\rm cl} )} and observe the effects of an inserted porous layer on the increase of heat transfer and pressure drop. The heat transfer and pressure drop increment ratios are used to define an overall performance (e = e_th/e_p)({\varepsilon = \varepsilon_{\rm th}/\varepsilon_{\rm p})} to evaluate overall benefits of an inserted porous layer in a parallel plate channel. The obtained results showed that for a partially porous filled channel, the value of e{\varepsilon} is highly influenced from Darcy number, but it is not affected from thermal conductivity ratio (k _r) when k _r > 2. For a fully porous material filled channel, the value of e{\varepsilon} is considerably affected from thermal conductivity ratio as the porous medium is in contact with the channel walls.     

Thermal Non-equilibrium Natural Convection in a Square Enclosure with Heat-Generating Porous Layer on Inner Walls

Differentially heated enclosure with heat-generating porous layer on inner walls is studied computationally for non-Darcy flow and thermal non-equilibrium models. In this study, this problem is investigated for different internal and external Rayleigh numbers, Darcy numbers, porosity-scaled thermal conductivity ratio, solid-/fluid-scaled heat transfer coefficient and dimensionless thickness of the porous layer. The results indicate that the dimensionless thickness of the porous layer has an important effect on the heat transfer in the enclosure. It was found that the thermal non-equilibrium model is needed for small values of the porosity-scaled thermal conductivity ratio and the solid-/fluid-scaled heat transfer coefficient. It is shown that the convection of heat due to internal heat generation is increased in the enclosure when the ratio of internal Rayleigh number to external Rayleigh number is larger.     

Analysis of fluid flow and heat transfer in a channel with intermittent heated porous blocks

A numerical study of forced convection enhancement in a channel intermittently heated is presented in this work. The use of porous blocks mounted on the heated parts of the channel to improve thermal performance is investigated. In order to account for the inertia, drag and boundary effects, the Brinkman-Forchheimer-extended Darcy model is adopted for the flow inside the porous regions. The effects of several parameters such as Darcy number, the block dimensions, the number of blocks and the thermal conductivity ratio are documented. The results show that the blocks may alter substantially the flow pattern depending on the permeability of the porous medium, and may improve the heat transfer and reduce the wall temperature under certain circumstances. Received on 9 June 1997     

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.     

Numerical Study of Heat Transfer in a Rectangular Channel with Porous Covering Obstacles

A numerical study is performed to analyze steady laminar forced convection in a channel in which discrete heat sources covered with porous material are placed on the bottom wall. Hydrodynamic and heat transfer results are reported. The flow in the porous medium is modeled using the Darcy–Brinkman–Forchheimer model. A computer program based on control volume method with appropriate averaging for diffusion coefficient is developed to solve the coupling between solid, fluid, and porous region. The effects of parameters such as Reynolds number, Prandtl number, inertia coefficient, and thermal conductivity ratio are considered. The results reveal that the porous cover with high thermal conductivity enhances the heat transfer from the solid blocks significantly and decreases the maximum temperature on the heated solid blocks. The mean Nusselt number increases with increase of Reynolds number and Prandtl number, and decrease of inertia coefficient. The pressure drop along the channel increases rapidly with the increase of Reynolds number.     

Forced convection from an isolated heat source in a channel with porous medium

A numerical study is made of flow and heat transfer characteristics of forced convection in a channel that is partially filled with a porous medium. The flow geometry models convective cooling process in a printed circuit board system with a porous insert.The channel walls are assumed to be adiabatic. Comprehensive numerical solutions are acquired to the governing Navier-Stokes equations, using the Brinkman-Forchheimer-extended Darcy model for the regions of porous media. Details of flow and thermal fields are examined over ranges of the principal parameters; i.e., the Reynolds number Re, the Darcy number Da (≡K/H²), the thickness of the porous substrate S, and the ratio of thermal conductivities R_k (≡k_eff/k). Two types of the location of the porous block are considered. The maximum temperature at the heat source and the associated pressure drop are presented for varying Re, Da, S, and R_k. The results illustrate that as S increases or Da decreases, the fluid flow rate increases. Also, as R_k increases for fixed Da, heat transfer rates are augmented. Explicit influences of Re on the flow and heat transport characteristics are also scrutinized. Assessment is made of the utility of using a porous insert by cross comparing the gain in heat transport against the increase in pressure drop.     

Generalized Plain Couette Flow and Heat Transfer in a Composite Channel

An analytical study of fluid flow and heat transfer in a composite channel is presented. The channel walls are maintained at different constant temperatures in such a way that the temperatures do not allow for free convection. The upper plate is considered to be moving and the lower plate is fixed. The flow is modeled using Darcy–Lapwood–Brinkman equation. The viscous and Darcy dissipation terms are included in the energy equation. By applying suitable matching and boundary conditions, an exact solution has been obtained for the velocity and temperature distributions in the two regions of the composite channel. The effects of various parameters such as the porous medium parameter, viscosity ratio, height ratio, conductivity ratio, Eckert number, and Prandtl number on the velocity and temperature fields are presented graphically and discussed.     

The Use of Porous Fins for Heat Transfer Augmentation in Parallel-Plate Channels

In this study, a steady, fully developed laminar forced convection heat augmentation via porous fins in isothermal parallel-plate duct is numerically investigated. High-thermal conductivity porous fins are attached to the inner walls of two parallel-plate channels to enhance the heat transfer characteristics of the flow under consideration. The Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous fins. This study reports the effect of several operating parameters on the flow hydrodynamics and thermal characteristics. This study demonstrates, mainly, the effects of porous fin thickness, Darcy number, thermal conductivity ratio, Reynolds number, and microscopic inertial coefficient on the thermal performance of the present flow. It is found that the highest Nusselt number is achieved at fully filled porous duct which requires the highest pumping pressure. The results show that using porous fins requires less pumping pressure with comparable high heat augmentation weight against fully filled porous duct. It is found that higher Nusselt numbers are achieved by increasing the microscopic inertial coefficient (A), the Reynolds number (Re), and the thermal conductivity of the porous substrate k ₂. The results show that heat transfer can be enhanced (1) with the use of high thermal conductivity fins, (2) by decreasing the Darcy number, and (3) by increasing microscopic inertial coefficient.     

Mixed Convection Heat Transfer in the Annulus Between Two Concentric Vertical Cylinders Using Porous Layers

A numerical investigation of the steady-state, laminar, axi-symmetric, mixed convection heat transfer in the annulus between two concentric vertical cylinders using porous inserts is carried out. The inner cylinder is subjected to constant heat flux and the outer cylinder is insulated. A finite volume code is used to numerically solve the sets of governing equations. The Darcy–Brinkman–Forchheimer model along with Boussinesq approximation is used to solve the flow in the porous region. The Navier–Stokes equation is used to describe the flow in the clear flow region. The dependence of the average Nusselt number on several flow and geometric parameters is investigated. These include: convective parameter, λ, Darcy number, Da, thermal conductivity ratio, K _r, and porous-insert thickness to gap ratio (H/D). It is found that, in general, the heat transfer enhances by the presence of porous layers of high thermal conductivity ratios. It is also found that there is a critical thermal conductivity ratio on which if the values of Kr are higher than the critical value the average Nusselt number starts to decrease. Also, it found that at low thermal conductivity ratio (K _r ≈ 1) and for all values of λ the porous material acts as thermal insulation.