Search results

The purpose of this study is to validate whether the local thermal equilibrium for unsteady state is an appropriate assumption for the porous media with closed pores. It also compares the transient temperatures between the pore scale and volume averaged approaches to prove that the volume averaged method is an appropriate technique for the heat transfer in closed-cell porous media. The interfacial heat transfer coefficient for the closed-cell porous media is also discussed in details.

The governing equations for the pore scale and continuum domains are given. They are solved numerically for the pore scale and volume-averaged domains. The results are compared and discussion was done. The performed discussions and explanations are supported with figure and graphics.

A local thermal non-equilibrium exits for the closed-cell porous media in which voids are filled with water during the unsteady heat transfer process. Local thermal non-equilibrium condition exists in the cells under high temperature gradient and it disappears when the heat transfer process becomes steady-state. Although a local thermal equilibrium exists in the porous media in which the voids are filled with air, a finite value for heat transfer coefficient is found. The thermal diffusivity of air and solid phase are close to each other and hence a local thermal equilibrium exists.

The study is done only for the closed-cell porous media and for Rayleigh number till 10⁵. Two common working fluids as water and air are considered.

There are many applications of porous media with closed pores particularly in the industry, such as the closed-cell metal foam or the closed cells in porous materials such as foods and plastic-based insulation material. The obtained results are important for transient heat transfer in closed-cell porous materials.

The obtained results are important from the transient application of heat transfer in the closed-cell material existing in nature and industry.

The authors’ literature survey shows that it is the first time the closed-cell porous media is discussed from local thermal non-equilibrium point of view and it is proved that the local thermal non-equilibrium can exist in the closed-cell porous media. Hence, two equations as solid and fluid equations should be used for unsteady heat transfer in a closed-cell porous medium.

The purpose of this paper is to develop an empiricism free, first principle‐based model to simulate fluid flow and heat transfer through porous media.

Conventional approaches to the problem are reviewed. A multi‐scale approach that makes use of the sample simulations at the individual pore levels is employed. The effect of porous structures on the global fluid flow is accounted for via local volume averaged governing equations, while the closure terms are accounted for via averaging flow characteristics around the pores.

The performance of the model has been tested for an isothermal flow case. Good agreement with experimental data were achieved. Both the permeability and Ergun coefficient are shown to be flow properties as opposed to the empirical approach which typically results in constant values of these parameters independent of the flow conditions. Hence, the present multi‐scale approach is more versatile and can account for the possible changes in flow characteristics.

Further validation including non‐isothermal cases is necessary. Current scope of the model is limited to incompressible flows. The methodology can accommodate extension to compressible flows.

This paper proposes a method that eliminates the dependence of the numerical porous media simulations on empirical data. Although the model increases the fidelity of the simulations, it is still computationally affordable due to the use of a multi‐scale methodology.

Regularities imposed by an external wall on the random distribution of particles in a packed bed lead to anisotropic wall effects of the bed. The resulting deviations of near—wall porosity from the average bed—porosity markedly affect the average velocity profile, not only near the wall but also some distance into the bed itself. The effect of wall channelling due to such variations in porosity near external boundaries is predicted by means of two different numerical solution methods for a unified model for granular porous media. The results are shown to compare favourably to experimental and numerical results reported in the literature.

The purpose of this paper is to focus on modeling buoyancy driven viscous flow and heat transfer through saturated packed pebble‐beds via a set of homogeneous volume‐averaged conservation equations in which local thermal disequilibrium is accounted for.

The local thermal disequilibrium accounted for refers to the solid and liquid phases differing in temperature in a volume‐averaged sense, which is modeled by describing each phase with its own governing equation. The partial differential equations are discretized and solved via a vertex‐centered edge‐based dual‐mesh finite volume algorithm. A compact stencil is used for viscous terms, as this offers improved accuracy compared to the standard finite volume formulation. A locally preconditioned artificial compressibility solution strategy is employed to deal with pressure incompressibility, whilst stabilisation is achieved via a scalar‐valued artificial dissipation scheme.

The developed technology is demonstrated via the solution of natural convective flow inside a heated porous axisymmetric cavity. Predicted results were in general within 10 per cent of experimental measurements.

This is the first instance in which both artificial compressibility and artificial dissipation is employed to model flow through saturated porous materials.

The purpose of this study is to investigate the applicability of volume average which is extensively used for analyzing the heat and fluid flow (both for single-phase and solid/liquid-phase change) in a closed cell porous medium numerically.

Heat conduction equations for the solid frame and fluid (or phase change material) are solved for pore scale and volume average approaches. The study mainly focuses on the effect of porosity and the number of porous media unit cell on the agreement between the results of the pore scale and volume average approaches.

It is observed for the lowest porosity values such as 0.3 and the number of porous media unit cell as 4 in heat transfer direction, the results between two approaches may be questionable for the single-phase fluid. By increasing the number of porous media unit cell in heat transfer direction, the agreement between two approaches becomes better. In general, for high porosity values (such as 0.9) the agreement between the results of two approaches is in the acceptable range both for single-phase and solid/liquid-phase change. Two charts on the applicability of volume average method for single-phase and solid/liquid-phase change are presented.

The authors’ literature survey shows that it is the first time the applicability of volume average which is extensively used for analyzing the heat and fluid flow in a closed cell porous medium is investigated numerically.

The purpose of this paper is to introduce a novel computational method to evaluate damage accumulation in a solder joint of an electronic package, when exposed to operating temperature environment. A procedure to implement the method is suggested, and a discussion of the method and its possible applications is provided in the paper.

Methodologically, interpolated response surfaces based on specially designed finite element (FE) simulation runs, are employed to compute a damage metric at regular time intervals of an operating temperature profile. The developed method has been evaluated on a finite-element model of a lead-free PBGA256 package, and accumulated creep strain energy density has been chosen as damage metric.

The method has proven to be two orders of magnitude more computationally efficient compared to FE simulation. A general agreement within 3 percent has been found between the results predicted with the new method, and FE simulations when tested on a number of temperature profiles from an avionic application. The solder joint temperature ranges between +25 and +75°C.

The method can be implemented as part of reliability assessment of electronic packages in the design phase.

The method enables increased accuracy in thermal fatigue life prediction of solder joints. Combined with other failure mechanisms, it may contribute to the accuracy of reliability assessment of electronic packages.

The purpose of this study is to determine interfacial convective heat transfer coefficient numerically, for a porous media consisting of square blocks in inline arrangement under mixed convection heat transfer.

The continuity, momentum and energy equations are solved in dimensionless form for a representative elementary volume of porous media, numerically. The velocity and temperature fields for different values of porosity, Ri and Re numbers are obtained. The study is performed for the range of Ri number from 0.01 to 10, Re number from 100 to 500 and porosity value from 0.51 to 0.96. Based on the obtained results, the value of the interfacial convective heat transfer coefficient is calculated by using volume average method.

It was found that at low porosities (such as 0.51), the interfacial Nusselt number does not considerably change with Ri and Re numbers. However, for porous media with high Ri number and porosity (such as 10 and 0.51, respectively), secondary flows occur in the middle of the channel between rods improving heat transfer between solid and fluid, considerably. It is shown that the available correlations of interfacial heat transfer coefficient suggested for forced convection can be used for mixed convection for the porous media with low porosity (such as 0.51) or for the flow with low Ri number (such as 0.01).

To the best of the authors’ knowledge, there is no study on determination of interfacial convective heat transfer coefficient for mixed convection in porous media in literature. The present study might be the first study providing an accurate idea on the range of this important parameter, which will be useful particularly for researchers who study on mixed convection heat transfer in porous media, macroscopically.

This paper aims to tackle the problem of some ambiguity of the momentum equation formulation in the commonly used macroscopic models of two‐phase solid/liquid region, developing during alloy solidification. These different appearances of the momentum equation are compared and the issue is addressed of how the choice of the particular form affects velocity and temperature fields.

Attention is focused on the ensemble averaging method, which, owing to its stochastic nature, is a new promising tool for setting up the macroscopic transport equations in highly inhomogeneous multiphase micro‐ and macro‐structures, with morphology continuously changing in time when the solidification proceeds. The basic assumptions of the two other continuum models, i.e. based on the classical mixture theory and on the volume‐averaging technique, are also unveiled. These three different forms of the momentum equation are then compared analytically and their impact on calculated velocity and temperature distribution in the mushy zone is studied for the selected test problem of binary alloy solidification driven by diffusion and thermal natural convection in a square mould.

It is found that a chosen appearance of the momentum equation mildly affects temporal velocity/temperature, and shapes of the phase interface at longer times of the solidification.

This mainly results from small variations of the liquid fraction across the mushy zone and from a low solidification rate, and it may change drastically when anisotropic properties of the mushy zone, solutal convection, different phase densities and cooling conditions are considered. Therefore, further comprehensive study is needed.

The paper addresses how the different focus of the momentum equation for liquid flow is compared.

The purpose of this paper is to simulate two macrosegregation benchmarks with a newly developed stabilized finite element algorithm based on a semi-implicit pressure correction scheme.

A streamline-upwind/Petrov–Galerkin (SUPG) stabilized finite element algorithm is developed for the coupled conservation equations of mass, momentum, energy and species. A semi-implicit pressure correction method combined with SUPG stabilization technique is proposed to solve the convection flow during solidification. An analytically derived enthalpy method is adopted to solve the energy conservation equation. The nonlinearities of the energy and species equations are tackled by Newton–Raphson method. Two macrosegregation benchmarks considering the solidification of an Al-4.5 per cent Cu alloy and a Sn-10 per cent Pb alloy are simulated.

A very good agreement is achieved by comparison with the classical finite volume method and a novel meshless method for the Al-4.5 per cent Cu alloy solidification benchmark. Moreover, a unique reference numerical solution has been obtained. Besides, it is demonstrated that the stabilized finite element algorithm can capture the flow instability and channel segregation during solidification of the Sn-10 per cent Pb alloy.

A semi-implicit pressure correction method combined with SUPG stabilization technique is adopted to develop robust stabilized finite element algorithm for the macrosegregation model. A new enthalpy formulation for heat transfer problems with phase change is derived analytically.

The purpose of this study is to analyze heat transfer for solid–liquid phase change in two inclined cavities assisted with open cell and closed cell porous structures for enhancement of heat transfer and compare them.

The heat transfer analysis is done numerically. The set of conservation equations for mass, momentum and energy for phase change material (PCM) and conduction heat transfer equation for metal frame are solved. Furthermore, temperature and solid–liquid fraction distributions for a cavity filled only with PCM are also obtained for comparison. The porosity is 0.9 for both porous structures. Rayleigh number and inclination angle change from 1 to 10⁸, and from −90° to 90°, respectively.

The present study reveals that the use of closed cell structures not only can make phase change faster than open cell structure (except for Ra = 10⁸ and = 90°) but also provide more stable process. The use of a closed cell porous structure in a cavity with PCM can reduce melting period up to 55% more than a cavity with an open cell porous structure. The rate of this additional enhancement depends on Rayleigh number and inclination angle.

To the best of the authors’ knowledge, this is the first time that the comparison between closed cell and open cell porous structures for heat transfer enhancement in a solid/liquid phase change process is reported. Authors believe that the present study will lead more attentions on the use of closed cell porous structures.