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.     

Transient Free Convection Fluid Flow in Domains Partially Filled with Porous Media

The problem of transient free convection in domains partly filled with porous substrates is investigated analytically using Laplace transformation technique. Four configurations are considered which are subject to an isothermal heating boundary condition. The Brinkman-extended Darcy model is adopted to describe the hydrodynamics behavior of the porous domain.     

Heat Transfer Characteristics of Reciprocating Flows in Channels Partially Filled with Porous Medium

This paper presents the effect of introducing a porous medium on the flow regime and heat transfer of a two-dimensional channel through which the flow is reciprocating. The channel is discretely heated from above and is insulated in the bottom which can simulate a cooling mechanism for compact circuit boards. In this ideal geometry, a fully developed reciprocating flow is established via oscillating pressure gradient. In side boundaries, velocity and temperature are assumed to be periodic. A certain volume of this channel is occupied by a porous medium which is shown to be an effecting tool for augmentation of heat transfer. At first, Momentum equations of the domain are solved analytically (Brinkman-extended Darcy model is used for porous region) and then the energy equation is solved numerically using alternating direction implicit (ADI) method. Finally a case study is investigated for a high-porous and high-conductive medium (Aluminum alloy T-6201) and the enhancing effect and optimization criteria are discussed in the result section.     

Numerical Simulation of Turbulent Flow and Heat Transfer in a Three-Dimensional Channel Coupled with Flow Through Porous Structures

This study investigates numerically the turbulent flow and heat transfer characteristics of a T-junction mixing, where a porous media flow is vertically discharged in a 3D fully developed channel flow. The fluid equations for the porous medium are solved in a pore structure level using an Speziale, Sarkar and Gatski turbulence model and validated with open literature data. Overall, two types of porous structures, consisted of square pores, are investigated over a wide range of Reynolds numbers: an in-line and a staggered pore structure arrangement. The flow patterns, including the reattachment length in the channel, the velocity field inside the porous medium as well as the fluctuation velocity at the interface, are found to be strongly affected by the velocity ratio between the transversely interacting flow streams. In addition, the heat transfer examination of the flow domain reveals that the temperature distribution in the porous structure is more uniform for the staggered array. The local heat transfer distributions inside the porous structure are also studied, and the general heat transfer rates are correlated in terms of area-averaged Nusselt number accounting for the effects of Reynolds number, velocity ratio as well as the geometrical arrangement of the porous structures.     

Analytical Investigation of the Effect of Viscous Dissipation on Couette Flow in a Channel Partially Filled with a Porous Medium

An analytical study of viscous dissipation effect on the fully developed forced convection Couette flow through a parallel plate channel partially filled with porous medium is presented. A uniform heat flux is imposed at the moving plate while the fixed plate is insulated. In the fluid-only region the flow field is governed by Navier–Stokes equation while the Brinkman-extended Darcy law relationship is considered in the fully saturated porous medium. The interface conditions are formulated with an empirical constant β due to the stress jump boundary condition. Fluid properties are assumed to be constant and the longitudinal heat conduction is neglected. A closed-form solution for the velocity and temperature distributions and also the Nusselt number in the channel are obtained and the viscous dissipation effect on these profiles is briefly investigated.     

Analytical and Numerical Study on Thermally Developing Forced Convective Flow in a Channel Filled with a Highly Porous Medium Under Local Thermal Non-Equilibrium

Transport in Porous Media - An analytical and numerical study was made on thermally developing forced convective flow in a channel filled with a fluid-saturated porous medium, subject to constant...     

An Analysis of Thermal Radiation in Porous Media Under Local Thermal Non-equilibrium

Transport in Porous Media - The present work investigates the thermal radiation transport inside porous media under local thermal non-equilibrium conditions. Two different geometrical situations, a...     

Darcy Flow in a Wavy Channel Filled with a Porous Medium

Flow in channels bounded by wavy or corrugated walls is of interest in both technological and geological contexts. This paper presents an analytical solution for the steady Darcy flow of an incompressible fluid through a homogeneous, isotropic porous medium filling a channel bounded by symmetric wavy walls. This packed channel may represent an idealized packed fracture, a situation which is of interest as a potential pathway for the leakage of carbon dioxide from a geological sequestration site. The channel walls change from parallel planes, to small amplitude sine waves, to large amplitude nonsinusoidal waves as certain parameters are increased. The direction of gravity is arbitrary. A plot of piezometric head against distance in the direction of mean flow changes from a straight line for parallel planes to a series of steeply sloping sections in the reaches of small aperture alternating with nearly constant sections in the large aperture bulges. Expressions are given for the stream function, specific discharge, piezometric head, and pressure.     

Physical Clogging of Uniformly Graded Porous Media Under Constant Flow rates

This study investigated the physical clogging of uniformly graded porous media under constant flow rates using natural porous media and suspensions. The porous media selected for this experimental study was a fine-to-medium sandy soil fractioned into thirteen uniformly graded beds: seven unisize beds and six uniform beds. The physical clogging of the beds was studied using two types of silt suspensions as along with two suspension concentrations and three water discharges. It was found that the permeability reduction due to physical clogging \([(K_\mathrm{i} - K_\mathrm{t})/K_\mathrm{i}]\) increased with decreasing \({D}_{15}/{d}_{85}\) ratios until a critical value of \({D}_{15}/{d}_{85}\), after which a surface mat of suspension was formed on the porous media. It was also found that the value of reduced permeability at any time (at any number of pore volumes of injected suspension-laden water), \(K_\mathrm{t}\), is directly proportional to square of \({D}_{15}\) and inversely proportional to \({C}_{\mathrm{u}}\) of the porous media and \({d}_{85}\) of suspensions. The effects of suspension type and flow rates on physical clogging seemed to depend on the size of the pores in the porous media.     

A Local Thermal Non-Equilibrium Analysis of Fully Developed Forced Convective Flow in a Tube Filled with a Porous Medium

A local thermal non-equilibrium model has been considered for the case of thermally fully developed flow within a constant heat flux tube filled with a porous medium. Exact temperature profiles for the fluid and solid phases are found after combining the two individual energy equations and then transforming them into a single ordinary differential equation with respect to the temperature difference between the solid phase and the wall subject to constant heat flux. The exact solutions for the case of metal-foam and air combination reveal that the local thermal equilibrium assumption may fail for the case of constant heat flux wall. The Nusselt number is presented as a function of the Peclet number, which shows a significant increase due to both high stagnant thermal conductivity and thermal dispersion resulting from the presence of the metal-foam.     

Conjugate Phase Change Heat Transfer in an Inclined Compound Cavity Partially Filled with a Porous Medium: A Deformed Mesh Approach

In this paper, the melting process of a PCM inside an inclined compound enclosure partially filled with a porous medium is theoretically addressed using a novel deformed mesh method. The sub-domain area of the compound enclosure is made of a porous layer and clear region. The right wall of the enclosure is adjacent to the clear region and is subject to a constant temperature of T_c. The left wall, which is connected to the porous layer, is thick and thermally conductive. The thick wall is partially subject to the hot temperature of T_h. The remaining borders of the enclosure are well insulated. The governing equations for flow and heat transfer, including the phase change effects and conjugate heat transfer at the thick wall, are introduced and transformed into a non-dimensional form. A deformed grid method is utilized to track the phase change front in the solid and liquid regions. The melting front movement is controlled by the Stefan condition. The finite element method, along with Arbitrary Eulerian–Lagrangian (ALE) moving grid technique, is employed to solve the non-dimensional governing equations. The modeling approach and the accuracy of the utilized numerical approach are verified by comparison of the results with several experimental and numerical studies, available in the literature. The effect of conjugate wall thickness, inclination angle, and the porous layer thickness on the phase change heat transfer of PCM is investigated. The outcomes show that the rates of melting and heat transfer are enhanced as the thickness of the porous layer increases. The melting rate is the highest when the inclination angle of the enclosure is 45°. An increase in the wall thickness improves the melting rate.

     

Flow of a Weakly Conducting Fluid in a Channel Filled with a Porous Medium

We investigate the fully developed flow in a fluid-saturated porous medium channel with an electrically conducting fluid under the action of a parallel Lorentz force. The Lorentz force varies exponentially in the vertical direction due to low fluid electrical conductivity and the special arrangement of the magnetic and electric fields at the lower plate. Exact analytical solutions are derived for fluid velocity and the results are presented in figures. All these flows are new and are presented for the first time in the literature.     

Role of Thermal Diffusion on Heat and Mass Transfer in Porous Media

In this paper, thermal diffusion phenomena in a porous cavity are investigated. The Brinkman model, coupled with the energy and the mass balance equations was solved numerically using a finite element techniques. A two-component system was included in the model. Different models were investigated to demonstrate the importance of the Soret effect with the presence of gravity vector. We do not take into consideration the pressure effect in the thermal diffusion. Even with such simplification to the problem, results reveal that the thermal diffusion is important and drives a strong convection. A series of convection cells are observed and steady-state solutions are obtained. Asymmetric solutions are obtained for various cases of dual-porosity porous media. Variations in the gravity vector indicated that the convection patterns, as well as the role of Soret coefficient, are profoundly impacted. Finally, the importance of including thermal diffusion in petroleum reservoir simulation is discussed.     

Thermal Instability in a Horizontal Porous Channel with Horizontal Through Flow and Symmetric Wall Heat Fluxes

The linear thermoconvective instability of the basic parallel flow in a plane and horizontal porous channel is investigated. The boundary walls are assumed to be impermeable and subject to symmetric and uniform heat fluxes. The wall heat fluxes produce either a net heating or a net cooling of the fluid saturated porous medium. A horizontal mass flow rate is externally impressed leading to a stationary basic state with a temperature gradient inclined to the vertical. A region of possibly unstable thermal stratification exists either in the lower half-channel (boundary heating), or in the upper half-channel (boundary cooling). The convective instability of the basic flow is governed by the Rayleigh number and by the Péclet number. In the case of boundary heating, the thermal instability arises when the Rayleigh number exceeds its critical value, that depends on the Péclet number. The change of the critical Rayleigh number as a function of the Péclet number is determined numerically for arbitrary normal modes oblique to the basic flow direction. The most dangerous modes are the longitudinal rolls, with a wave vector perpendicular to the basic velocity. There exists a minimum value of the Péclet number, 19.1971, below which no linear instability is detected.     

Buoyant Flow with Viscous Heating in a Vertical Circular Duct Filled with a Porous Medium

Combined, forced, and free flow in a vertical circular duct filled with a porous medium is investigated according to the Darcy–Boussinesq model. The effect of viscous dissipation is taken into account. It is shown that a thermal boundary condition compatible with fully developed and axisymmetric flow is either a linearly varying wall temperature in the axial direction or, only in the case of uniform velocity profile, an axial linear-exponential wall temperature change. The case of a linearly varying wall temperature corresponds to a uniform wall heat flux and includes the uniform wall temperature as a special case. A general analytical solution procedure is performed, by expressing the seepage velocity profile as a power series with respect to the radial coordinate. It is shown that, for a fixed thermal boundary condition, i.e., for a prescribed slope of the wall temperature, and for a given flow rate, there exist two solutions of the governing balance equations provided that the flow rate is lower than a maximum value. When the maximum value is reached, the dual solutions coincide. When the flow rate is higher than its maximum, no axisymmetric solutions exist. E. Magyari is on leave from the Institute of Building Technology, ETH—Zürich.     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

Fully Developed Mixed Convection Flow Between Inclined Parallel Plates Filled with a Porous Medium

A fully developed mixed convection flow between inclined parallel flat plates filled with a porous medium is considered through which there is a constant flow rate and with heat being supplied to the fluid by the same uniform heat flux on each plate. The equations governing this flow are made non-dimensional and are seen to depend on two dimensionless parameters, a mixed convection parameter λ and the Péclet number Pe, as well as the inclination γ of the plates to the horizontal. The velocity and temperature profiles are obtained in terms of λ, Pe and γ when the channel is inclined in an upwards direction as well as for horizontal channels. The limiting cases of small and large λ and small Pe are considered with boundary-layer structures being seen to develop on the plates for large values of λ.