Nonsimilar Boundary Layer Analysis of Free Convection Heat Transfer Over a Vertical Cylinder in Bidisperse Porous Media

This work presents a boundary layer analysis for the free convection heat transfer from a vertical cylinder in bidisperse porous media with constant wall temperature. A boundary layer analysis and the two-velocity two-temperature formulation are used to derive the nonsimilar governing equations. The transformed governing equations are solved by the cubic spline collocation method to yield computationally efficient numerical solutions. The effects of inter-phase heat transfer parameter, modified thermal conductivity ratio, and permeability ratio on the heat transfer and flow characteristics are studied. Results show that an increase in the modified thermal conductivity ratio and the permeability ratio can effectively enhance the free convection heat transfer of the vertical cylinder in a bidisperse porous medium. Moreover, the thermal nonequilibrium effects are strong for low values of the inter-phase heat transfer parameter.     

多孔介质平板通道强迫对流中热局部非平衡时的热应力

根据微观不可压饱和多孔介质热-力-流相互作用的一般理论,在固相骨架小变形的假定下,考虑固相和流相相互作用的粘性耗散,研究了多孔介质平板通道强迫对流热局部非平衡的热应力问题.建立了问题的热-力数学模型,根据饱和多孔介质的平衡方程,在固相骨架只存在横向位移的假定下,求解了固相骨架的位移和相应的热应力,数值考察了各种物性参数对热应力分布的影响,讨论了热局部平衡模型的适用性.     

GURTIN-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF A NON-LOCAL THERMAL EQUILIBRIUM SATURATED POROUS MEDIUM

Based on the porous media theory and by taking into account the effects of the pore fluid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fluid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fluid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small deformation of the solid phase, small velocity of the fluid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fluid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles,especially Hu-Washizu type variational principles, for the initial boundary value problems of dynamic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fluid-saturated elastic porous media, which have been discussed previously.     

Thermo-Osmosis Effect in Saturated Porous Medium

In this paper, the thermo-poroelasticity theory is used to investigate the quasi-static response of temperatures, pore pressure, stress, displacement, and fluid flux around a cylindrical borehole subjected to impact thermal and mechanical loadings in an infinite saturated poroelastic medium. It has been reported in literatures that coupled flow known as thermo-osmosis by which flux is driven by temperature gradient, can significantly change the fluid flux in clay, argillaceous and many other porous materials whose permeability coefficients are very small. This study presents a mathematical model to investigate the coupled effect of thermo-osmosis in saturated porous medium. The energy balance equations presented here fulfill local thermal non-equilibrium condition (LTNE) which is different from the local thermal equilibrium transfer theory, accounting for that temperatures of solid and fluid phases are not the same and governed by different heat transfer equations. Analytical solutions of temperatures, pore pressure, stress, displacement, and fluid flux are obtained in Laplace transform space. Numerical results for a typical clay are used to investigate the effect of thermo-osmosis. The effects of LTNE on temperatures, pore pressure, and stress are also studied in this paper.     

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

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

TEMPERATURE PROFILES OF LOCAL THERMAL NONEQUILIBRIUM FOR THERMAL DEVELOPING FORCED CONVECTION IN POROUS MEDIUM PARALLEL PLATE CHANNEL

Based on the two-energy equation model, taking into account viscous dissipation due to the interaction between solid skeleton and pore fluid flow, temperature expressions of the solid skeleton and pore fluid flow are obtained analytically for the thermally developing forced convection in a saturated porous medium parallel plate channel, with walls being at constant temperature. It is proved that the temperatures of the two phases for the local thermal nonequilibrium approach to the temperature derived from the one-energy equation model for the local thermal equilibrium when the heat exchange coefficient goes to infinite. The temperature profiles are shown in figures for different dimensionless parameters and the effects of the parameters on the local thermal nonequilibrium are revealed by parameter study.     

Pore-Scale Simulation of Interphase Multicomponent Mass Transfer for Subsurface Flow

A simulation framework is proposed to simulate multicomponent multiphase flow in porous media at the pore scale. It solves equations for the species concentrations in the framework of the volume-of-fluid approach including thermodynamics equilibrium at the fluid/fluid interface. Particular attention is paid to the derivation of the boundary condition for the concentration at the solid walls. The method is validated by comparison with analytical solutions of simple setups. Then, the approach is used to investigate and upscale mass transfer across interfaces in different configurations, including the drainage of water in a tube by a gas carrying a contaminant, mass transfer in thin films, and mass transfer in complex porous structures under various flow conditions.     

Modeling and Simulation of Local Thermal Non-equilibrium Effects in Porous Media with Small Thermal Conductivity

The two-equation model in porous media can describe the local thermal non-equilibrium (LTNE) effects between fluid and solid at REV scale, with the temperature differences in a solid particle neglected. A multi-scale model has been proposed in this study. In the model, the temperature differences in a solid particle are considered by the coupling of the fluid energy equation at REV scale with the heat conduction equation of a solid particle at pore scale. The experiments were conducted to verify the model and numerical strategy. The multi-scale model is more suitable than the two-equation model to predict the LTNE effects in porous media with small thermal conductivity. The effects of particle diameter, mass flow rate, and solid material on the LTNE effects have been investigated numerically when cryogenic nitrogen flows through the porous bed with small thermal conductivity. The results indicate that the temperature difference between solid center and fluid has the same trend at different particle diameters and mass flow rates, while the time to reach the local thermal equilibrium is affected by solid diameter dramatically. The results also show that the temperature difference between solid center and surface is much greater than that between solid surface and fluid. The values of \( \rho {\text{c}} \) for different materials have important influence on the time to reach the local thermal equilibrium between solid and fluid.     

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.     

Pore Network Analysis of Resistivity Index for Water-Wet Porous Media

Pore network analysis is used to investigate the effects of microscopic parameters of the pore structure such as pore geometry, pore-size distribution, pore space topology and fractal roughness porosity on resistivity index curves of strongly water-wet porous media. The pore structure is represented by a three-dimensional network of lamellar capillary tubes with fractal roughness features along their pore-walls. Oil-water drainage (conventional porous plate method) is simulated with a bond percolation-and-fractal roughness model without trapping of wetting fluid. The resistivity index, saturation exponent and capillary pressure are expressed as approximate functions of the pore network parameters by adopting some simplifying assumptions and using effective medium approximation, universal scaling laws of percolation theory and fractal geometry. Some new phenomenological models of resistivity index curves of porous media are derived. Finally, the eventual changes of resistivity index caused by the permanent entrapment of wetting fluid in the pore network are also studied.Resistivity index and saturation exponent are decreasing functions of the degree of correlation between pore volume and pore size as well as the width of the pore size distribution, whereas they are independent on the mean pore size. At low water saturations, the saturation exponent decreases or increases for pore systems of low or high fractal roughness porosity respectively, and obtains finite values only when the wetting fluid is not trapped in the pore network. The dependence of saturation exponent on water saturation weakens for strong correlation between pore volume and pore size, high network connectivity, medium pore-wall roughness porosity and medium width of the pore size distribution. The resistivity index can be described succesfully by generalized 3-parameter power functions of water saturation where the parameter values are related closely with the geometrical, topological and fractal properties of the pore structure.     

Validation of the Local Thermal Equilibrium Assumption in Forced Convection of Non-Newtonian Fluids through Porous Channels

In this paper, we assess the validity of the local thermal equilibrium assumption in the non-Newtonian forced convection flow through channels filled with porous media. For this purpose, the problem is solved numerically using local thermal non-equilibrium and non-Darcian models. Numerical solutions obtained over broad ranges of representative dimensionless parameters are utilized to map conditions at which the local thermal equilibrium assumption can or cannot be employed. The circumstances of a higher modified Peclet number, a lower modified Biot number, a lower fluid-to-solid thermal conductivity ratio, a lower power-law fluid index, and a lower microscopic and macroscopic frictional flow resistance coefficients, are identified as unfavorable circumstances for the local thermal equilibrium (LTE) condition to hold. Quantitative LTE validity maps that reflect the proportional effect of each parameter as related to others are presented.     

Analysis of hydrodynamic and thermal dispersion in porous media by means of a local approach

A pore scale analysis is implemented in this numerical study to investigate the behavior of microscopic inertia and thermal dispersion in a porous medium with a periodic structure. The macroscopic characteristics of the transport phenomena are evaluated with an averaging technique of the controlling variables at a pore scale level in an elementary cell of the porous structure. The Darcy–Forchheimer model describes the fluid motion through the porous medium while the continuity and Navier–Stokes equations are applied within the unit cell. An average energy equation is employed for the thermal part of the porous medium. The macroscopic pressure loss is computed in order to evaluate the dominant microscopic inertial effects. Local fluctuations of velocity and temperature at the pore scale are instrumental in the quantification of the thermal dispersion through the total effective thermal diffusivity. The numerical results demonstrate that microscopic inertia contributes significantly to the magnitude of the macroscopic pressure loss, in some instances with as much as 70%. Depending on the nature of the porous medium, the thermal dispersion may have a marked bearing on the heat transfer, particularly in the streamwise direction for a highly conducting fluid and certain values of the Peclet number.     

A pore scale study on turbulent combustion in porous media

This paper presents pore scale simulation of turbulent combustion of air/methane mixture in porous media to investigate the effects of multidimensionality and turbulence on the flame within the pores of porous media. In order to investigate combustion in the pores of porous medium, a simple but often used porous medium consisting of a staggered arrangement of square cylinders is considered in the present study. Results of turbulent kinetic energy, turbulent viscosity ratio, temperature, flame speed, convective heat transfer and thermal conductivity are presented and compared for laminar and turbulent simulations. It is shown that the turbulent kinetic energy increases from the inlet of burner, because of turbulence created by the solid matrix with a sudden jump or reduction at the flame front due to increase in temperature and velocity. Also, the pore scale simulation revealed that the laminarization of flow occurs after flame front in the combustion zone and turbulence effects are important mainly in the preheat zone. It is shown that turbulence enhances the diffusion processes in the preheat zone, but it is not enough to affect the maximum flame speed, temperature distribution and convective heat transfer in the porous burner. The dimensionless parameters associated with the Borghi–Peters diagram of turbulent combustion have been analyzed for the case of combustion in porous media and it is found that the combustion in the porous burner considered in the present study concerns the range of well stirred reactor very close to the laminar flame region.     

Review on Scale Dependent Characterization of the Microstructure of Porous Media

The paper discusses local porosity theory and its relation with other geometric characterization methods for porous media such as correlation functions and contact distributions. Special emphasis is placed on the charcterization of geometric observables through Hadwigers theorem in stochastic geometry. The four basic Minkowski functionals are introduced into local porosity theory, and for the first time a relationship is established between the Euler characteristic and the local percolation probabilities. Local porosity distributions and local percolation probabilities provide a scale dependent characterization of the microstructure of porous media that can be used in an effective medium approach to predict transport.     

Evaluation of the Darcy's law performance for two-fluid flow hydrodynamics in a particle debris bed using a lattice-Boltzmann model

In the present paper, multiphase flow dynamics in a porous medium are analyzed by employing the lattice-Boltzmann modeling approach. A two-dimensional formulation of a lattice-Boltzmann model, using a D2Q9 scheme, is used. Results of the FlowLab code simulation for single phase flow in porous media and for two-phase flow in a channel are compared with analytical solutions. Excellent agreement is obtained. Additionally, fluid-fluid interaction and fluid-solid interaction (wettability) are modeled and examined. Calculations are performed to simulate two-fluid dynamics in porous media, in a wide range of physical parameters of porous media and flowing fluids. It is shown that the model is capable of determining the minimum body force needed for the nonwetting fluid to percolate through the porous medium. Dependence of the force on the pore size, and geometry, as well as on the saturation of the nonwetting fluid is predicted by the model. In these simulations, the results obtained for the relative permeability coefficients indicate the validity of the reciprocity for the two coupling terms in the modified Darcy's law equations. Implication of the simulation results on two-fluid flow hydrodynamics in a decay-heated particle debris bed is discussed. Received on 1 December 1999     

Conjugate Natural Convection in a Porous Enclosure with Non-Uniform Heat Generation

Effects of a conductive wall on natural convection in a square porous enclosure having internal heating at a rate proportional to a power of temperature difference is studied numerically in this article. The horizontal heating is considered, where the vertical walls heated isothermally at different temperatures while the horizontal walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and finite difference method is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (0 ???Ra ???1000), the internal heating and the local exponent parameters (0 ????? ???5), (1 ????? ???3), the wall to porous thermal conductivity ratio (0.44 ???Kr ???9.9) and the ratio of wall thickness to its width (0.02 ???D ???0.5). The results are presented to show the effect of these parameters on the fluid flow and heat transfer characteristics. It is found a strong internal heating can generate significant maximum fluid temperature more than the conductive solid wall. Increasing value thermal conductivity ratio and/or decreasing the thickness of solid wall can increase the maximum fluid temperature. It is also found that at very low Rayleigh number, the heat transfer across the porous enclosure remain stable for any values of the thermal conductivity ratio.     

Mixed Convection Boundary Layer Flow Past a Horizontal Circular Cylinder Embedded in a Bidisperse Porous Medium

Adopting a two-temperature and two-velocity model, appropriate to a bidisperse porous medium (BDPM) proposed by Nield and Kuznetsov (2008), the classical steady, mixed convection boundary layer flow about a horizontal, isothermal circular cylinder embedded in a porous medium has been theoretically studied in this article. It is shown that the boundary layer analysis leads to expressions for the flow and heat transfer characteristics in terms of an inter-phase momentum parameter, a thermal diffusivity ratio, a thermal conductivity ratio, a permeability ratio, a modified thermal capacity ratio, and a buoyancy or mixed convection parameter. The transformed partial differential equations governing the flow and heat transfer in the f-phase (the macro-pores) and the p-phase (the remainder of the structure) are solved numerically using a very efficient implicit finite-difference technique known as Keller-box method. A good agreement is observed between the present results and those known from the open literature in the special case of a traditional Darcy formulation (monodisperse system).     

Non-Darcian Effects on Buoyancy-Induced Heat Transfer in a Partially Divided Square Enclosure with Internal Heat Generation

The present investigation addresses non-Darcian effects on the buoyancy-induced heat transfer in a partially divided square enclosure with internal heat generation. The generalized model of the momentum equation, which is also known as the Forchheimer–Brinkman extended Darcy model, which takes into account boundary and inertia effects, was used in representing the fluid motion inside the porous layer. The local thermal equilibrium condition was assumed to be valid for the range of the thermophysical parameters considered in the present investigation. The transport equations were solved using the finite element formulation based on the Galerkin method of weighted residuals. The validity of the numerical code used was ascertained by comparing our results with previously published results. Results were obtained in terms of streamlines, isotherms, and Nusselt number for various geometrical parameters specifying the height and width of the partition. In addition, the effects of external and internal Rayleigh numbers and Darcy number were highlighted in the proposed study.