Search results

The purpose of this paper is to investigate thermally and hydrodynamically fully developed convection in a duct of rectangular cross-section containing a porous medium and fluid layer.

The Darcy–Brinkman–Forchheimer flow model is adopted. A finite difference method of second-order accuracy with the Southwell-over-relaxation method is deployed to solve the non-dimensional momentum and energy conservation equations under physically robust boundary conditions.

It is found that the presence of porous structure and different immiscible fluids exert a significant impact on controlling the flow. Graphical results for the influence of the governing parameters i.e. Grashof number, Darcy number, porous media inertia parameter, Brinkman number and ratios of viscosities, thermal expansion and thermal conductivity parameters on the velocity and temperature fields are presented. The volumetric flow rate, skin friction and rate of heat transfer at the left and right walls of the duct are also provided in tabular form. The numerical solutions obtained are validated with the published study and excellent agreement is attained.

To the author’s best knowledge this study original in developing the numerical code using FORTRAN to assess the fluid properties for immiscible fluids. The study is relevant to geothermal energy systems, thermal insulation systems, resin flow modeling for liquid composite molding processes and hybrid solar collectors.

The purpose of this paper is to provide a numerical study of conjugate heat transfer by mixed convection and conduction in a lid-driven enclosure with thick vertical porous layer. The effect of the relevant parameters: Richardson number (Ri=0.1, 1, 10) and thermal conductivity ratio (Rk=0.1, 1, 10, 100) are investigated.

The studied system is a two dimensional lid-driven enclosure with thick vertical porous layer. The left vertical wall of the enclosure is allowed to move in its own plane at a constant velocity. The enclosure is heated from the right vertical wall isothermally. The left and the right vertical walls are isothermal but temperature of the outside of the right vertical wall is higher than that of the left vertical wall. Horizontal walls are insulated. The governing equations are solved by finite volume method and the SIMPLE algorithm.

From the finding results, it is observed that: for the two studied cases, heat transfer rate along the hot wall is a decreasing function of thermal conductivity ratio irrespective of Richardson numbers contrary to the heat transfer rate along the fluid-porous layer interface which is an increasing function of thermal conductivity ratio. At forced convection dominant regime, the difference between heat transfer rate for upward and downward moving wall is insensitive to the thermal conductivity ratio. For downward moving wall, average Nusselt number is higher than that of upward moving wall.

Some applications: building applications, furnace design, nuclear reactors, air solar collectors.

From the bibliographic work and the authors’ knowledge, the conjugate mixed convection in lid-driven partially porous enclosures has not yet been investigated which motivates the present work that represent a continuation of the preceding investigations.

This paper aims to provide a limited, but selective bibliography on modelling heat and mass transfer in composite fluid‐porous domains.

Since the pioneer study by Beavers and Joseph, the problem of interface continuity and/or jump conditions at a fluid‐porous interface has been of interest to the fluid mechanics and heat and mass transfer community. The paper is concerned both with numerical simulations of heat and fluid flow in such systems, and with the linear stability problems.

The one‐ and two‐domain formulations are equivalent. Using the Darcy‐Brinkman extension instead of the Darcy model reduces the number of ad hoc parameters in this configuration.

The problem of double diffusive convection has still to be solved and analyzed.

The discussion on the interface conditions is of great relevance to many industrial and practical situations.

The important question of the macroscopic formulation of the problem is tackled in the paper.

The purpose of this paper is to implement the numerical analysis based on finite volume method to compare the effects of stress-jump (SJ) and stress-continuity (SC) conditions on flow structure around and through a porous circular cylinder.

In this study, a steady flow of a viscous, incompressible fluid around and through a porous circular cylinder of diameter “D,” using Darcy-Brinkman-Forchheimer’s equation in the porous region, is discussed. The SJ condition proposed by Ochoa-Tapia and Whitaker is applied at the porous-fluid interface and compared with the traditional interfacial condition based on the SC condition in fluid and porous media. 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 < 40 and 10-7 < Da < 10-2, respectively, and the porosities are e=0.45, 0.7 and 0.95.

Results show that the SJ condition leads to a much smaller boundary layer within porous medium near the interface as compared to the SC condition. Two interfacial conditions yield similar results with decrease in porosity.

There is no published research in the literature about the effects of important parameters, such as Porosity and Darcy numbers on different fluid-porous interface conditions for a porous cylinder and comparison the effects of SJ and SC conditions on flow structure around and through a porous circular cylinder.

This paper aims to consider numerical analysis of laminar double-diffusive natural convection inside a non-homogeneous closed medium composed of a saturated porous matrix and a clear binary fluid under spatial sinusoidal heating/cooling on one side wall and uniform salting.

The domain of interest is a partially square porous enclosure with sinusoidal wall heating and cooling. The fluid flow, heat and mass transfer dimensionless governing equations associated with the corresponding boundary conditions are discretized using the finite volume method. The resulting algebraic equations are solved by an in-house FORTRAN code and the SIMPLE algorithm to handle the non-linear character of conservation equations. The validity of the in-house FORTRAN code is checked by comparing the current results with previously published experimental and numerical works. The effect of the porous layer thickness, the spatial frequency of heating and cooling, the Darcy number, the Rayleigh number and the porous to fluid thermal conductivity ratio is analyzed.

The results demonstrate that for high values of the spatial frequency of heating and cooling (f = 7), temperature contours show periodic variations with positive and negative values providing higher temperature gradient near the thermally active wall. In this case, the temperature variation is mainly in the porous layer, while the temperature of the clear fluid region is practically the same as that imposed on the left vertical wall. This aspect can have a beneficial impact on thermal insulation. Besides, the porous to fluid thermal conductivity ratio, Rk, has practically no effect on Shhot wall, contrary to Nuinterface where a strong increase is observed as Rk is increased from 0.1 to 100, and much heat transfer from the hot wall to the clear fluid via the porous media is obtained.

The findings are useful for devices working on double-diffusive natural convection inside non-homogenous cavities.

The authors believe that the presented results are original and have not been published elsewhere.

To study the combined effects of permeability and couplestresses on the performance characteristics of a secant‐shaped porous slider bearing lubricated with Stokes couplestress fluid.

The modified Reynolds type equation governing the fluid film pressure is derived on the basis of Stokes microcontinuum theory of couplestress fluids by using Beavers‐Joshep slip boundary conditions at the fluid‐porous interface. The modified Reynolds equation is solved analytically and closed form expressions are obtained for the fluid film pressure, load carrying capacity, frictional force and centre of pressure.

The bearing characteristics are computed for various values of the couple stress parameter, slip parameter and the permeability parameters. It is found that, the effect of couple stress is to increase the load carrying capacity and to decrease the coefficient of friction as compared to the Newtonian case. However, the effect of permeability parameter is to decrease the load carrying capacity. It is also found that, the effect of couple stresses is to shift the center of pressure towards the outlet edge.

The end effects are neglected in the analysis and these can be included in the study by considering full three‐dimensional problem.

Reduction in the load carrying capacity due to the presence of porous facing can be compensated by the use of lubricants containing additives of proper size. As a result the bearing performance can be improved.

This paper provides closed form expressions for the bearing characteristics and are analyzed with respect to non‐dimensional parameters, viz. couplestress parameter, slip parameter and permeability parameter. This paper offers help to design engineers to design efficient bearing systems.

The purpose of this study is a numerical analysis of transient natural convection in a square partially porous cavity with a heat-generating and heat-conducting element using the local thermal non-equilibrium model under the effect of cooling from the vertical walls. It should be noted that this research deals with a development of passive cooling system for the electronic devices.

The domain of interest is a square cavity with a porous layer and a heat-generating element. The vertical walls of the cavity are kept at constant cooling temperature, while the horizontal walls are adiabatic. The heat-generating solid element is located on the bottom wall. A porous layer is placed under the clear fluid layer. The governing equations, formulated in dimensionless stream function, vorticity and temperature variables with corresponding initial and boundary conditions, are solved using implicit finite difference schemes of the second order accuracy. The governing parameters are the Darcy number, viscosity variation parameter, porous layer height and dimensionless time. The effects of varying these parameters on the average total Nusselt number along the heat source surface, the average temperature of the heater, the fluid flow rate inside the cavity and on the streamlines and isotherms are analyzed.

The results show that in the case of local thermal non-equilibrium the total average Nusselt number is an increasing function of the interphase heat transfer coefficient and the porous layer thickness, while the average heat source temperature decreases with the Darcy number and viscosity variation parameter.

An efficient numerical technique has been developed to solve this problem. The originality of this work is to analyze unsteady natural convection within a partially porous cavity using the local thermal non-equilibrium model in the presence of a local heat-generating solid element. The results would benefit scientists and engineers to become familiar with the analysis of convective heat transfer in enclosures with local heat-generating heaters and porous layers, and the way to predict the heat transfer rate in advanced technical systems, in industrial sectors including transportation, power generation, chemical sectors and electronics.

The purpose of this paper is to report on the flow past a porous square cylinder, implementing the stress jump treatments for the porous‐fluid interface.

The numerical method was developed for flows involving an interface between a homogenous fluid and a porous medium. It is based on the finite volume method with body‐fitted and multi‐block grids. The Brinkman‐Forcheimmer extended model was used to govern the flow in the porous medium region. At its interface, a shear stress jump that includes the inertial effect was imposed, together with a continuity of normal stress.

The present model is validated by comparing with those for the flow around a solid circular cylinder. Results for flow around porous square cylinder are presented with flow configurations for different Darcy number, 10⁻² to 10⁻⁵, porosity from 0.4 to 0.8, and Reynolds number 20 to 250. The flow develops from steady to unsteady periodic vortex shedding state. It was found that the stress jump interface condition can cause flow instability. The first coefficient β has a more noticeable effect whereas the second coefficient β₁ has very small effect, even for Re=200. The effects of the porosity, Darcy number, and Reynolds number on lift and drag coefficients, and the length of circulation zone or shedding period are studied.

The present study implements the numerical method based on finite volume method with a collocated variable arrangement to treat the stress jump condition.

In this work, a new two-phase version of the finite element-based Artificial Compressibility (AC) Characteristic-Based Split (CBS) algorithm is developed and applied for the first time to heat and mass transfer phenomena in porous media with associated phase change. The purpose of this study is to provide an alternative for the theoretical analysis and numerical simulation of multiphase transport phenomena in porous media. Traditionally, the more complex Separate Flow Model was used in which the vapour and liquid phases were considered as distinct fluids and mathematically described by the conservation laws for each phase separately, resulting in a large number of governing equations.

Even though the adopted mathematical model presents analogies with the conventional multicomponent mixture flow model, it is characterized by a considerable reduction in the number of the differential equations for the primary variables. The fixed-grid numerical formulation can be applied to the resolution of general problems that may simultaneously include a superheated vapour region, a two-phase zone and a sub-cooled liquid region in a single physical domain with irregular and moving phase interfaces in between. The local thermal non-equilibrium model is introduced to consider the heat exchange between fluid and solid within the porous matrix.

The numerical model is verified considering the transport phenomena in a homogenous and isotropic porous medium in which water is injected from one side and heated from the other side, where it leaves the computational domain in a superheated vapour state. Dominant forces are represented by capillary interactions and two-phase heat conduction. The obtained results have been compared with the numerical data available in the scientific literature.

The present algorithm provides a powerful routine tool for the numerical modelling of complex two-phase transport processes in porous media.

For the first time, the stabilized AC-CBS scheme is applied to the resolution of compressible viscous flow transport in porous materials with associated phase change. A properly stabilized matrix inversion-free procedure employs an adaptive local time step that allows acceleration of the solution process even in the presence of large source terms and low diffusion coefficients values (near the phase change point).

The purpose of this paper is to focus on the numerical analysis of transient free convection heat transfer in partially porous cylindrical domains. The authors analyze the dependence of velocity and temperature fields on the geometry, by analyzing transient flow behavior for different values of cavity aspect ratio and radii ratio; both inner and outer radius are assumed variable in order to not change the difference r_o-r_i. Moreover, several Darcy numbers have been considered.

A dual time-stepping procedure based on the transient artificial compressibility version of the characteristic-based split algorithm has been adopted in order to solve the transient equations of the generalized model for heat and fluid flow through porous media. The present model has been validated against experimental data available in the scientific literature for two different problems, steady-state free convection in a porous annulus and transient natural convection in a porous cylinder, showing an excellent agreement.

For vertically divided half porous cavities, with Rayleigh numbers equal to 3.4×106 for the 4:1 cavity and 3.4×105 for the 8:1 cavity, the numerical results show that transient oscillations tend to disappear in presence of cylindrical geometry, differently from what happens for rectangular one. The magnitude of this phenomenon increases with radii ratio; the porous layer also affects the stability of velocity and temperature fields, as oscillations tend to decrease in presence of a porous matrix with lower value of the Darcy number.

A proper analysis of partially porous annular cavities is fundamental for the correct estimation of Nusselt numbers, as the formulas provided for rectangular domains are not able to describe these problems.

The proposed model represents a useful tool for the study of transient natural convection problems in porous and partially porous cylindrical and annular cavities, typical of many engineering applications. Moreover, a fully explicit scheme reduces the computational costs and ensures flexibility.

This is the first time that a fully explicit finite element scheme is employed for the solution of transient natural convection in partially porous tall annular cavities.