Onset of Buoyancy-driven Instability in Porous Media Melted from Below

When a nonhomogeneous solid is melting from below, convection may be induced in a thermally–unstable melt layer. In this study, the onset of buoyancy-driven convection during time-dependent melting is investigated by using similarly transformed disturbance equations. The critical Darcy–Rayleigh numbers based on the melt-layer thickness, Ra _H,c, are found numerically for various conditions. For small superheats, the present predictions show that Ra _H,c is located between 27.1 and 4π ² and it approaches the well-known results of the original Horton–Rogers–Lapwood problem. However, for high superheats, it is dependent on the phase change rate λ and the relation of Ra _H,c λ = 25.89 is shown.     

Onset of Buoyancy-Driven Convection in Porous Media Saturated with Cold Water Cooled from Above

When porous media saturated with initially stagnant cold water around the density maximum temperature are cooled from above, convection may be induced in an unstable lower layer. In this study, the onset of buoyancy-driven convection during time-dependent cooling is investigated using the propagation theory, which transforms disturbance equations similarly, and also considering the density inversion effect. The critical Darcy–Rayleigh number Ra _D,c is found as a function of the dimensionless density maximum temperature θ _max. For Ra _D > Ra _D,c the dimensionless critical time τ _c to mark the onset of instability is presented as a function of Ra _D and θ _max. These critical conditions are compared with previous theoretical results.     

Analytical Solution for Forced Convection in a Semi-Circular Channel Filled with a Porous Medium

Fully developed laminar forced convection inside a semi-circular channel filled with a Brinkman-Darcy porous medium is studied. Analytical solutions for flow and constant flux heat transfer are found using a mixture of Cartesian and cylindrical coordinates. The problem depends on a parameter s, which is proportional to the inverse square of the Darcy number. Velocity boundary layers exist when s is large. Both friction factor-Reynolds number product and Nusselt number are determined. Closed form expressions for the clear fluid () limit are found. Rare analytical solutions not only describe fundamental channel flows, but also serve as a check for more complicated numerical solutions.     

Mixed Convection Boundary-Layer Flow in a Porous Medium Filled with Water Close to its Maximum Density

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water close to its maximum density is considered for uniform wall temperature and outer flow. The problem can be reduced to similarity form and the resulting equations are examined in terms of a mixed convection parameter λ and a parameter δ which measures the difference between the ambient temperature and the temperature at the maximum density. Both assisting (λ > 0) and opposing flows (λ < 0) are considered. A value δ₀ is found for which there are dual solutions for a range λ_c < λ < 0 of λ (the value of λ_c dependent on δ) and single solutions for all λ ≥ 0. Another value of δ₁ of δ, with δ₁ > δ₀, is found for which there are dual solutions for a range 0 < λ < λ_c of positive values of λ, with solutions for all λ≤ 0. There is also a range δ₀ <  δ < δ₁ where there are solutions only for a finite range of λ, with critical points at both positive and negative values of λ, thus putting a finite limit on the range of existence of solutions.     

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 λ.     

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.     

Reply to Comments on “Non-Darcian Forced Convection Flow of Viscous Dissipating Fluid over a Flat Plate Embedded in a Porous Medium”

In this note, we reply to the comments by Rees and Magyari (2007) on our article (Aydin and Kaya 2007). They mainly stated that the thermal boundary conditions we defined at the edge of the boundary layer were incompatible with the energy equation. This is questionable and therefore we will discuss it below. They disclosed that our results were in error. However, this is quite misleading. Scientifically, they cannot reach such a conclusion without comparing our results with what they thought to be correct. In fact, this misleading and unproven statement will be shown not to be correct in the following.     

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two-dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left-hand vertical wall has temperature T _h and the right-hand vertical wall is maintained at temperature T _c (T _h > T _c) and the horizontal walls are insulated. At time t > 0, the left-hand vertical wall temperature is suddenly raised to which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.     

Vortex Instability of the Asymptotic Dissipation Profile in a Porous Medium

In this note we consider the thermoconvective stability of the recently-discovered asymptotic dissipation profile (ADP). The ADP is a uniform thickness, parallel-flow boundary layer which is induced by a cold surface in a warm saturated porous medium in the presence of viscous dissipation. We have considered destabilisation in the form of stream-wise vortex disturbances. The critical wavenumber and Rayleigh number for the onset of convection have been determined for all angles of the cooled surface between the horizontal and the vertical for which the ADP exists. The paper closes with a presentation of some strongly nonlinear computations of steady vortices.     

Free Convection from a Vertical Plate in a Porous Medium Subjected to a Porous Medium Subjected to a Sudden Change in Surface Heat Flux

An analysis is made of the transient free convection from a vertical flat plate which is embedded in a fluid-saturated porous medium. It is assumed that for time a steady state temperature or velocity has been obtained in the boundary-layer which occurs due to a uniform flux dissipation rate . Then at time the heat flux on the plate is suddenly changed to and maintained at this value for 0$$ " align="middle" border="0"> . An analytical solution has been obtained for the temperature/velocity field for small times in which the transport effects are confined within an inner layer adjacent to the plate. These effects cause a new steady boundary layer. A numerical solution of the full boundary-layer equations is then obtained for the whole transient from to the steady state, firstly by means of a step-by-step method and then by a matching technique. The transition between the two distinct solution methods is always observed to occur very near to the turning point of the plate surface temperature, a time at which the fluid temperature is close to its steady state profile. The solution obtained using the step-by-step method shows excellent agreement with the small time analytical solution. Results are presented to illustrate the occurrence of transients from both small and large increases and decreases in the levels of existing energy inputs.     

Numerical Simulations of Flows in a Heterogeneous Porous Medium

Numerical simulations are employed to investigate the fluid flow and pressure loss in a heterogeneous block within a composite porous medium. The mean permeability of the heterogeneous block is seen to affect the overall effective flux significantly. The heterogeneous parameters of the permeability field, such as the correlation length and variance, affect it quite differently. Because of the channelling effects, the effective flux depends strongly on the realization of the permeability for larger correlation length. Under a specific permeability field, higher effective flux results from smaller variances. The influences of the inertial factor are found to be insignificant within the range of practical interests.     

Mixed Convection Boundary-Layer Flow Near the Stagnation Point on a Vertical Surface in a Porous Medium: Brinkman Model with Slip

The steady boundary-layer flow near the stagnation point on an impermeable vertical surface with slip that is embedded in a fluid-saturated porous medium is investigated. Using appropriate similarity variables, the governing system of partial differential equations is transformed into a system of ordinary differential equations. This system is then solved numerically. The features of the flow and the heat transfer characteristics for different values of the governing parameters, namely, the Darcy–Brinkman, Γ, mixed convection, λ, and slip, γ, parameters, are analysed and discussed in detail for the cases of assisting and opposing flows. It is found that dual solutions exist for assisting flows, as well as those usually reported in the literature for opposing flows. A stability analysis of the steady flow solutions encountered for different values of the mixed convection parameter λ is performed using a linear temporal stability analysis. This analysis reveals that for γ  =  0 (slip absent) and Γ  =  1 the lower solution branch is unstable while the upper solution branch is stable.     

Effects of Temperature-Dependent Viscosity on Forced Convection Inside a Porous Medium

Considering the exponential viscosity–temperature relation, effect of temperature-dependent viscosity on forced convection of a liquid through a porous medium, bounded by isoflux parallel plates, is investigated numerically based on the general model of momentum transfer. Local effects of viscosity variation on the distribution of velocity and temperature are analyzed. Moreover, global aspects of the problem are investigated where corrections are proposed for total pressure drop and the fully developed Nusselt number, in the form of out/in viscosity ratio. Results are obtained over a wide range of permeabilities from clear (of solid material) fluid to very low permeability, where for constant properties one expects a nearly slug flow.     

Fourier Representation of Intermittent Flow in a Porous Medium

The field measurements and numerical results for intermittent flow regime in a sandy soil show that the time distributions of the soil water flux q(z, t), and the soil water content θ(z, t)at various depths are periodic in nature, where t is time and z is the depth (i.e., at the surface z = 0 and at depths z = − 5, − 10, − 15 cm, etc). The period of q(z, t) and θ(z, t) variations are generally determined by the sum of the duration of pulse and the duration between the initiation of two consecutive pulses of water at the soil surface. Fourier series models have been given for q(z, t) and θ(z, t) variations. The predicted Fourier results for these variations have been compared with the experimentally verified numerical results—designated as observed values. The results show that the amplitudes of these variations were damped exponentially with depth, and the phase shift increased linearly with depth.     

The Effect of Strong Heterogeneity on the Onset of Convection in a Porous Medium

The effect of strong heterogeneity on the onset of convection induced by a vertical density gradient in a saturated porous medium governed by Darcy’s law is discussed. The general case, where there is heterogeneity in both the vertical and horizontal directions, and where there is heterogeneity in permeability, thermal conductivity, and applied temperature gradient, is considered. A computer package has been developed to implement an algorithm giving a criterion for instability. This package is applied to the case of a cube partitioned into octants and to the cases where the permeability and thermal conductivity vary continuously across a cube in either a linear or a quadratic manner.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Free Convection in Spherical Annular Sectors Filled with a Porous Medium

Calculated free convection flows and heat transfer are presented for concentric spherical annular sectors, filled with a porous medium. Two isothermal walls and an adiabatic radial wall at the sector angle define the sectors. The governing equations (in the stream function and temperature formulation) are solved numerically using ADI (alternative direction implicit) finite-difference method. Over the range of geometric parameters examined, the obtained results for spherical annuli and low Rayleigh number Ra. As Ra increases, multicellular flows develop for small values of the aspect ratio parameter . In addition, analytical solutions of the governing equations were obtained for small values of Ra (1) and it was shown that these solutions agree well with those obtained numerically. Significant differences in the local heat transfer rates on the inner and outer walls of the spherical annuli were observed from these solutions.     

Modelling the Effect of Surface Tension on Immiscible Displacement in a Saturated Porous Medium

Comments are made about the model employed by Chen and Vafai for forced convection in a porous medium channel, with a surface tension effect at the moving interface between two fluids, when one fluid is displaced by the other. A simple situation is analysed, and the circumstances under which surface tension effects are important in this case are clarified.     

Attenuated Wave Field in Fluid-Saturated Porous Medium with Excitations of Multiple Sources

This research addresses the investigation of an elastic wave field in a homogeneous and isotropic porous medium which is fully saturated by a Newtonian viscous fluid. A new methodology is developed for describing the wave field in the medium excited by multiple energy sources. To quantify the relative displacements between the fluid and solid of the medium, the governing equations of the elastic wave propagation are derived in the form of displacements specially. The velocities and attenuation of the waves are considered as functions of viscosity and frequency. Making use of the Hankel function and the moving-coordinate method, a model of the wave motion with multiple cylindrical wave sources is built. Making use of the model established in this research, the relative displacement between the fluid and the solid can be quantified, and the wave field in the porous media can then be determined with the given energy sources. Numerical simulations of cylindrical waves from multiple energy sources propagating in the porous medium saturated by viscous fluid are performed for demonstrating the practicability of the model developed.     

Interaction of Transverse Heterogeneity and Thermal Development of Forced Convection in a Porous Medium

The classical Graetz methodology is applied to investigate the thermal development of forced convection in a parallel plate channel filled by a saturated porous medium whose permeability and thermal conductivity vary in the transverse direction. It was found that there is a significant interaction between heterogeneity and thermal development.