Effects of Time-Periodic Thermal Boundary Conditions and Internal Heating on Heat Transport in a Porous Medium

The effects of time-periodic boundary temperatures and internal heating on Nusselt number in the Bénard–Darcy convective problem has been considered. The amplitudes of temperature modulation at the lower and upper surfaces are considered to be very small. By performing a weakly non-linear stability analysis, the Nusselt number is obtained in terms of the amplitude of convection, which is governed by the non-autonomous Ginzburg–Landau equation, derived for the stationary mode of convection. The effects of internal Rayleigh number, amplitude and frequency of modulation, thermo-mechanical anisotropies, and Vadasz number on heat transport have been analyzed and depicted graphically. Increasing values of internal Rayleigh number results in the enhancement of heat transport in the system. Further, the study establishes that the heat transport can be controlled effectively by a mechanism that is external to the system.     

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.     

Effect of Discrete Heating on Natural Convection in a Rectangular Porous Enclosure

The main objective of this article is to study the effect of discrete heating on free convection heat transfer in a rectangular porous enclosure containing a heat-generating substance. The left wall of the enclosure has two discrete heat sources and the right wall is isothermally cooled at a lower temperature. The top and bottom walls, and the unheated portions of the left wall are adiabatic. The vorticity–stream function formulation of the governing equations is numerically solved using an implicit finite difference method. The effects of aspect ratio, Darcy number, heat source length, and modified Rayleigh number on the flow and heat transfer are analyzed. The numerical results reveal that the rate of heat transfer increases as the modified Rayleigh number and the Darcy number increases, but decreases on increasing the aspect ratio. The average heat transfer rate is found to be higher at the bottom heater than at the top heater in almost all considered parameter cases except for ε = 0.5. Also, the maximum temperature takes place generally at the top heater except for the case ε = 0.5, where the maximum temperature is found at the bottom heater. Further, the numerical results reveal that the maximum temperature decreases with the modified Rayleigh number and increases with the aspect ratio.     

Convective Flows in a Melt During Crystal Growth under Periodic Temperature Field Conditions

The results of an analysis of thermo-gravitational convection in Czochralski single-crystal growth from a melt under periodic temperature field conditions are presented. The numerical modeling is based on solving the unsteady 3-D Navier-Stokes and heat conduction equations in the Boussinesq approximation. It is shown that using different heating regimes for the crystal medium provides additional opportunities for controlling heat and mass transfer.     

Convective Flows on Reactive Surfaces in Porous Media

A model for the convective flow in a fluidsaturated porous medium containing a reactive component is considered. This component undergoes an exothermic reaction (modelled by a first order mechanism) on an impermeable bounding surface, the resulting heat released driving the convective flow. Large Rayleigh number flow near a stagnation point is treated in detail by first considering the steady states. Multiple solution branches and critical points arising from a hysteresis bifurcation are identified. The form that these solution branches take depends on whether or not the effects of reactant consumption are included. An initialvalue problem is then discussed. This shows that both the lower (slow reaction) and upper (fast reaction) solution branches are stable (and the ultimate state of the system). When the parameter values are such that there is no steady state, the solution develops a finitetime singularity, the nature of which is analysed.     

Staggered Finite Volume Modeling of Transport Phenomena in Porous Materials with Convective Boundary Conditions

The aim of this article is to present a one-dimensional finite volume modeling of mass transport phenomena in partially saturated open porous media. The finite volume model is based on a staggered solution strategy which permits to solve the equations sequentially. This appropriate partitioning reduces thus the size of the system to be solved at each time step. Therefore, two iterative algorithms are proposed to solve the discretized problem. Moreover, an original treatment of the convective boundary condition, based on a suitable linearization to derive the algebraic form of this condition, is also presented. The computational model is then specifically used for studying the drying process of a cementitious material. For this purpose, the two proposed algorithms are compared in terms of computational efficiency. Thereafter, the convective boundary condition is analyzed with regard to the drying kinetics. Finally, a one-dimensional drying experimental test is simulated.     

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.     

Multi-layered Porous Foam Effects on Heat Transfer and Entropy Generation of Nanofluid Mixed Convection Inside a Two-Sided Lid-Driven Enclosure with Internal Heating

Transport in Porous Media - Mixed convection of Cu-water nanofluid inside a two-sided lid-driven enclosure with an internal heater, filled with multi-layered porous foams is studied numerically and...     

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.     

Turbulence Characteristics in Flows Over Solid and Porous Square Ribs Mounted on Porous Walls

PIV measurements have been performed for turbulent flows in a rib-mounted channel whose bottom wall is made of a porous layer. The ratio of the rib and channel heights is fixed at 0.5. The effects of the wall and rib permeability are investigated focusing on the separating and reattaching flows at the bulk Reynolds number of 10³???10⁴. Three kinds of foamed ceramics are employed as the porous media. They have the same porosity of 0.8 but each permeability is different from the others. Its normalized values by the rib height are 0.89 × 10^???4, 1.47 × 10^???4 and 3.87 × 10^???4. Two kinds of square cylinder ribs: an impermeable smooth solid rib or a permeable porous rib which is made of the same porous medium as that for the bottom wall are used. The obtained turbulent velocity fields of the solid rib flows indicate that the turbulent intensity behind the rib becomes weak and the recirculation bubble in the clear channel tends to vanish as the the wall permeability increases. In the porous rib flow, the recirculation and the reattachment point shift downstream and turbulence becomes weaker due to the bleeding flow through the rib. In the higher permeability cases, the recirculation bubble hardly exists due to the flows through not only the bottom wall but also the porous rib. From the measurements, it is suggested that in the solid rib flows, a reverse flow region exists inside the porous wall whereas in porous rib flows, such reverse flow does not exist at higher permeability.     

Some Pitfalls in the Modelling of Convective Flows in Porous Media

Some modelling deficiencies in various recent papers, on convective flows in porous media, are pointed out and discussed. These deficiencies are related to the Forchheimer coefficient, mixed double-diffusive convection, magnetohydrodynamic mixed convection, convection in a rarefied gas, and geophysical phenomena.     

Onset of Convection with Internal Heating in a Weakly Heterogeneous Porous Medium

Linear stability analysis is applied to the onset of convection due to internal heating in a porous medium with weak vertical and horizontal heterogeneity. It is found that the effect of horizontal heterogeneity of each of permeability and thermal conductivity is slightly destabilizing. Increase of permeability in the upward direction is destabilizing and increase in the downward direction is stabilizing, and the reverse is true for increase of conductivity.     

Mixed Convection Boundary-Layer Flow on a Vertical Surface in a Porous Medium with a Constant Convective Boundary Condition

The mixed convection boundary-layer flow on a vertical surface heated convectively is considered when a constant surface heat transfer parameter is assumed. The problem is seen to be chararterized by a mixed convection parameter $\gamma $ γ . The flow and heat transfer near the leading edge correspond to forced convection solution and numerical solutions are obtained to determine how the solution then develops. The solution at large distances is obtained and this identifies a critical value $\gamma _c$ γ c of the parameter $\gamma $ γ . For $\gamma > \gamma _c$ γ > γ c a solution at large distances is possible and this is approached in the numerical integrations. For $\gamma <\gamma _c$ γ < γ c the numerical solution breaks down at a finite distance along the surface with a singularity, the nature of which is discussed.     

Onset of Convection with Internal Heating in a Porous Medium Saturated by a Nanofluid

Linear stability analysis was applied to the onset of convection due to internal heating in a porous medium saturated by a nanofluid. A model in which the effects of thermophoresis and Brownian motion are taken into account is employed. We utilized more realistic boundary conditions than in the previous work on this subject; now the nanofluid particle fraction is allowed to adapt to the temperature profile induced by the internal heating, subject to the requirement that there is zero perturbation flux across a boundary. The results show that the presence of the nanofluid particles leads to increased instability of the system. We identified two combinations of dimensionless parameters that are the major controllers of convection instability in the layer.     

Numerical Simulation of Forced Convective Heat Transfer Past a Square Diamond-Shaped Porous Cylinder

Fluid flow and heat transfer around and through a porous cylinder is an important issue in engineering applications. In this paper a numerical study is carried out for simulating the fluid flow and forced convection heat transfer around and through a square diamond-shaped porous cylinder. The flow is two-dimensional, steady, and laminar. Conservation laws of mass, momentum, and heat transport equations are applied in the clear region and Darcy–Brinkman–Forchheimer model for simulating the flow in the porous medium has been used. Equations with the relevant boundary conditions are numerically solved using a finite volume approach. In this study, Reynolds and Darcy numbers are varied within the ranges of $1<Re<45$ and $10^{-6}<Da<10^{- 2}$ , respectively. The porosity $(\varepsilon )$ is 0.5. This paper presents the effect of Reynolds and Darcy numbers on the flow structure and heat transfer characteristics. Finally, these parameters are compared among solid and porous cylinder. It was found that the drag coefficient decreases and flow separation from the cylinder is delayed with increasing Darcy number. Also the size of the thermal plume decreases by decreasing Darcy number.     

Convection in a Lid-Driven Heat-Generating Porous Cavity with Alternative Thermal Boundary Conditions

Mixed convection flow in a two-sided lid-driven cavity filled with heat-generating porous medium is numerically investigated. The top and bottom walls are moving in opposite directions at different temperatures, while the side vertical walls are considered adiabatic. The governing equations are solved using the finite-volume method with the SIMPLE algorithm. The numerical procedure adopted in this study yields a consistent performance over a wide range of parameters that were 10⁻⁴ ≤ Da ≤ 10⁻¹ and 0 ≤ Ra _I ≤ 10⁴. The effects of the parameters involved on the heat transfer characteristics are studied in detail. It is found that the variation of the average Nusselt number is non-linear for increasing values of the Darcy number with uniform or non-uniform heating condition.     

Natural Convection Boundary-Layer Flow in a Porous Medium with Temperature-Dependent Boundary Conditions

The natural convection boundary-layer flow on a surface embedded in a fluid-saturated porous medium is discussed in the case when the wall heat flux is related to the wall temperature through a power-law variation. The flow within the porous medium is assumed to be described by Darcy’s law and the Boussinesq approximation is assumed for the density variations. Two cases are discussed, (i) stagnation-point flow and (ii) flow along a vertical surface. The possible steady states are considered first with the governing partial equations reduced to ordinary differential equations by similarity transformations and these latter equations further transformed to previously studied free-convection problems. This identifies values of the exponent N in the power-law wall temperature variation, N = 3/2 for stagnation-point flows and 3/2 ≤ N ≤ 3 for the vertical surface, where similarity solutions do not exist. Time development for stagnation-point flows is seen to depend on N, for N <  3/2 the steady state is approached at large times, for N ≥ 3/2 a singularity develops at finite time leading to thermal runaway. Numerical solutions for the vertical surface, where the temperature-dependent boundary condition becomes more significant as the solution develops, show that, for N < 3/2, the corresponding similarity solution is approached, whereas for N >  3/2 the solution breaks down at a finite distance along the surface.     

Heat Transfer in Porous Microchannels with Second-Order Slipping Boundary Conditions

Transport in Porous Media - The paper presents results of a study of forced convection in a vertical flat and circular microchannels incorporating porous medium under slip boundary conditions of...     

Non-Darcy Mixed Convection in a Vertical Porous Channel with Boundary Conditions of Third Kind

Combined free and forced convection flow in a parallel plate vertical channel filled with porous matrix is analyzed in the fully developed region with boundary conditions of third kind. The flow is modeled using the Brinkman?CForchheimer-extended Darcy equations. The plates exchange heat with an external fluid. Both conditions of equal and different reference temperatures of the external fluid are considered. Governing equations are solved numerically by shooting technique that uses classical explicit Runge?CKutta scheme and Newton?CRaphson method as a correction scheme and analytically using perturbation series method for Darcy model. The velocity field, the temperature field and Nusselt numbers are obtained for governing parameters such as porous parameter, inertia term and perturbation parameter for equal and unequal Biot numbers and are displayed graphically. The dimensionless mean velocity and bulk temperature are also determined. It is found that the numerical solutions agree for small values of the perturbation parameter in the absence of the inertial forces.