Search results

A combination of highly conductive porous media and nanofluids is an efficient way for improving thermal performance of relevant applications. For precisely predicting the flow and thermal transport of nanofluids in porous media, the purpose of this paper is to explore the inter-phase coupling numerical methods.

Based on the lattice Boltzmann (LB) method, this study combines the convective flow, non-equilibrium thermal transport and phase interactions of nanofluids in porous matrix and proposes a new multi-phase LB model. The micro-scale momentum and heat interactions are especially analyzed for nanoparticles, base fluid and solid matrix. A set of three-phase LB equations for the flow/thermal coupling of base fluid, nanoparticles and solid matrix is established.

Distributions of nanoparticles, velocities for nanoparticles and the base fluid, temperatures for three phases and interaction forces are analyzed in detail. Influences of parameters on the nanofluid convection in the porous matrix are examined. Thermal resistance of nanofluid convective transport in porous structures are comprehensively discussed with the models of multi-phases. Results show that the Rayleigh number and the Darcy number have significant influences on the convective characteristics. The result with the three-phase model is mildly larger than that with the local thermal non-equilibrium model.

This paper first creates the multi-phase theoretical model for the complex coupling process of nanofluids in porous structures, which is useful for researchers and technicians in fields of thermal science and computational fluid dynamics.

This paper describes the numerical solutions of type‐IV shock‐on‐shock interactions in hypersonic thermochemical nonequilibrium air flows around blunt bodies. The Navier‐Stokes equation solver for a chemically reacting and vibrationally relaxing gas mixture was applied to the present problem, where the concepts of the Advection Upstream Splitting Method (AUSM) and the Lower‐Upper Symmetric Gauss‐Seidel (LU‐SGS) method were basically employed along with the two‐temperature thermochemical model of Park. The aerodynamic heating with or without the shock‐on‐shock interaction to a sphere and circular cylinders are simulated for a hypersonic nonequilibrium flow. The numerical results show that typical type‐IV shock‐on‐shock interaction pattern with a supersonic jet structure is also formed in a high‐enthalpy thermochemical nonequilibrium flow, and that the contribution of convective heat flux in the translational/rotational mode to the total heat flux is dominant. Furthermore, the inherent unsteadiness of nonequilibrium type‐IV shock‐on‐shock interaction flowfield is discussed briefly.

The purpose of this study is to numerically analyze the mixed convection and entropy generation in an annulus with a rotating heated inner cylinder for single-wall carbon nanotube (SWCNT)–water nanofluid flow using local thermal nonequilibrium (LTNE) model. An examination of the system behavior is presented considering the heat-generating solid phase inside the porous layer partly filled at the inner surface of the outer cylinder.

The discretized governing equations for nanofluid and porous layer by means of the finite volume method are solved by using the SIMPLE algorithm.

It is found that the buoyancy force and rotational effect have an important impact on the change of the strength of streamlines and isotherms for nanofluid flow. The minimum average Nusselt number on the inner cylinder is obtained at Ra$_E$ = 10$^4$, and the minimum total entropy generation is found at Re = 400 for given parameters. The entropy generation minimization is determined in case of different nanoparticle volume fractions. It is observed that at the same external Rayleigh numbers, the LTNE condition obtained with internal heat generation is very different from that without heat generation.

To the best of the authors’ knowledge, there is no previous paper presenting mixed convection and entropy generation of SWCNT–water nanofluid in a porous annulus under LTNE condition. The addition of nanoparticles to based fluid leads to a decrease in the value of minimum total entropy generation. Thus, using nanofluid has a significant role in the thermal design and optimization of heat transfer applications.

This study aims to investigate lattice Boltzmann (LB) simulation of the fluid flow and heat transfer characteristics of a heated porous elliptic cylinder in uniform flow based on the two-domain scheme. In the present research, the effect of axis ratio (1 ≤ AR ≤ 2), Reynolds number (5 ≤ Re ≤ 40) and Darcy number (10⁻⁴ ≤ Da ≤ 10⁻²) are studied.

To perform the LB simulation based on the two-domain scheme, the nonequilibrium extrapolation method is modified to model the heat transfer interfacial conditions required at the curved interface.

The results show that the axis ratio as well as Reynolds and Darcy numbers significantly affect the fluid flow and heat transfer characteristics of the porous elliptic cylinder. It is shown that for AR > 1, the phenomenon of detached recirculating zone occurs at much higher Darcy numbers compared with the case of the porous circular cylinder (AR = 1). The results show that the location of maximum temperature within the cylinder moves downstream when the Reynolds number, Darcy number and axis ratio increase. It is also concluded that the average Nusselt number of a porous elliptic cylinder is always lower than that of a porous circular cylinder.

The LB simulation of forced convection from a porous cylinder in uniform flow with a curved interface based on the two-domain scheme has not been studied yet.

This paper aims to establish a lattice Boltzmann (LB) method for solid-liquid phase transition (SLPT) from the pore scale to the representative elementary volume (REV) scale. By applying this method, detailed information about heat transfer and phase change processes within the pores can be obtained, while also enabling the calculation of larger-scale SLPT problems, such as shell-and-tube phase change heat storage systems.

Three-dimensional (3D) pore-scale enthalpy-based LB model is developed. The computational input parameters at the REV scale are derived from calculations at the pore scale, ensuring consistency between the two scales. The approaches to reconstruct the 3D porous structure and determine the REV of metal foam were discussed. The implementation of conjugate heat transfer between the solid matrix and the solid−liquid phase change material (SLPCM) for the proposed model is developed. A simple REV-scale LB model under the local thermal nonequilibrium condition is presented. The method of bridging the gap between the pore-scale and REV-scale enthalpy-based LB models by the REV is given.

This coupled method facilitates detailed simulations of flow, heat transfer and phase change within pores. The approach holds promise for multiscale calculations in latent heat storage devices with porous structures. The SLPT of the heat sinks for electronic device thermal control was simulated as a case, demonstrating the efficiency of the present models in designing and optimizing SLPT devices.

A coupled pore-scale and REV-scale LB method as a numerical tool for investigating phase change in porous materials was developed. This innovative approach allows for the capture of details within pores while addressing computations over a large domain. The LB method for simulating SLPT from the pore scale to the REV scale was given. The proposed method addresses the conjugate heat transfer between the SLPCM and the solid matrix in the enthalpy-based LB model.

The onset of convection in a ferrofluid-saturated porous layer has been investigated using a local thermal nonequilibrium (LTNE) model by allowing the solid phase to transfer heat via a Cattaneo heat flux theory while the fluid phase to transfer heat via usual Fourier heat-transfer law. The flow in the porous medium is governed by modified Brinkman-extended Darcy model. The instability of the system is discussed exactly for stress-free boundaries, while for rigid-ferromagnetic/paramagnetic boundaries the results are obtained numerically using the Galerkin method. The presence of Cattaneo effect introduces oscillatory convection as the preferred mode of instability contrary to the occurrence of instability via stationary convection found in its absence. Besides, oscillatory ferroconvection is perceived when the solid thermal relaxation time parameter exceeds a threshold value and increase in its value is to hasten the oscillatory onset. The effect of different boundary conditions on the instability of the system is noted to be qualitatively same. The paper aims to discuss these issues.

The investigators would follow the procedure of Straughan (2013) to obtain the expression for Rayleigh number. The Brinkman-extended Darcy model is used to describe the flow in a porous medium. The investigators have used a Galerkin method to obtain the numerical results for rigid-ferromagnetic/paramagnetic boundaries, while the instability of the system is discussed exactly for stress-free boundaries.

The Cattaneo–LTNE porous ferroconvection has been analyzed for different velocity and magnetic boundary conditions. The Brinkman-extended Darcy model is used to describe the flow in a porous medium. The effect of different types of velocity and magnetic boundary conditions on the instability of the system has been highlighted. The instability of the system is discussed exactly for stress-free boundaries, while for rigid-ferromagnetic/paramagnetic boundaries the results are obtained numerically using the Galerkin method.

The novelty of the present paper is to combine LTNE and second sound effects in solids on thermal instability of a ferrofluid-saturated porous layer by retaining the usual Fourier heat-transfer law in the ferrofluid. The Brinkman-extended Darcy model is used to describe the flow in a porous medium. The effect of different types of velocity and magnetic boundary conditions on the instability of the system is discussed.

Partially‐decoupled upwind‐based total‐variation‐diminishing (TVD) finite‐difference schemes for the solution of the conservation laws governing two‐dimensional non‐equilibrium vibrationally relaxing and chemically reacting flows of thermally‐perfect gaseous mixtures are presented. In these methods, a novel partially‐decoupled flux‐difference splitting approach is adopted. The fluid conservation laws and species concentration and vibrational energy equations are decoupled by means of a frozen flow approximation. The resulting partially‐decoupled gas‐dynamic and thermodynamic subsystems are then solved alternately in a lagged manner within a time marching procedure, thereby providing explicit coupling between the two equation sets. Both time‐split semi‐implicit and factored implicit flux‐limited TVD upwind schemes are described. The semi‐implicit formulation is more appropriate for unsteady applications whereas the factored implicit form is useful for obtaining steady‐state solutions. Extensions of Roe's approximate Riemann solvers, giving the eigenvalues and eigenvectors of the fully coupled systems, are used to evaluate the numerical flux functions. Additional modifications to the Riemann solutions are also described which ensure that the approximate solutions are not aphysical. The proposed partially‐decoupled methods are shown to have several computational advantages over chemistry‐split and fully coupled techniques. Furthermore, numerical results for single, complex, and double Mach reflection flows, as well as corner‐expansion and blunt‐body flows, using a five‐species four‐temperature model for air demonstrate the capabilities of the methods.

The Fay and Riddell (F–R) formula is an empirical equation for estimating the stagnation-point heat flux on noncatalytic and fully catalytic surfaces, based on an assumption of equilibrium. Because of its simplicity, the F–R has been used extensively for reentry flight design as well as ground test facility applications. This study aims to investigate the uncertainties of the F-R formula by considering velocity gradient, chemical species at the boundary layer edge, and the thermochemical nonequilibrium (NEQ) behind the shock layer under various hypersonic NEQ flow environments.

The stagnation-point heat flux calculated with the F–R formula was evaluated by comparison with thermochemical NEQ calculations and existing flight experimental values.

The comparisons showed that the F–R underestimated the noncatalytic heat flux, because of the chemical composition at the surface. However, for fully catalytic heat flux, the F–R results were similar to values of surface heat flux from thermochemical NEQ calculations, because the F–R formula overestimates the diffusive heat flux. When compared with the surface heat flux results obtained from flight experimental data, the F–R overestimated the fully catalytic heat flux. The error was 50% at most.

The results provided guidelines for the F–R calculations under hypersonic flight conditions and for determining the approximate error range for noncatalytic and fully catalytic surfaces.

The purpose of this work is the study of the steady free convection in a square differentially heated cavity filled by a Brinkman bidisperse porous medium. An appropriate mathematical model considering the Brinkman, momentum and energy interphase terms is proposed. The dependence of the stream functions, isotherms and of the Nusselt numbers on the governing parameters is analysed.

The both phases of flow and heat transfer are solved numerically using a modified finite difference technique. The algebraic system obtained after discretization is solved using the SOR method. The results are found to be in a significant agreement with the ones presented by the literature for a Darcy bidisperse porous medium and a Brinkman monodisperse porous medium.

The effects of the governing parameters on stream functions, isotherms and Nusselt numbers are discussed. It has been found that in the case of the Brinkman bidisperse model, the Nusselt numbers decrease compared to the Darcy model, and this behaviour is significant in comparison to the Brinkman monodisperse case.

A mathematical model for the free convection inside a cavity filled by a non-Darcy bidisperse porous medium, based on the Brinkman equation, is used. The effect of Darcy number, Rayleigh number, modified inter-phase heat transfer parameter, modified thermal conductivity ratio and the inertial parameters is studied. To the best of the authors’ knowledge, this problem has not been studied before, and the results are new and original.

Efficient thermal management of central processing unit (CPU) cooling systems is vital in the context of advancing information technology and the demand for enhanced data processing speeds. This study aims to explore the thermal performance of a CPU cooling setup using a cylindrical porous metal foam heat sink.

Nanofluid flow through the metal foam is simulated using the Darcy–Brinkman–Forschheimer equation, accounting for magnetic field effects. The temperature distribution is modeled through the local thermal equilibrium equation, considering viscous dissipation. The problem’s governing partial differential equations are solved using the similarity method. The CPU’s hot surface serves as a solid wall, with nanofluid entering the heat sink as an impinging jet. Verification of the numerical results involves comparison with existing research, demonstrating strong agreement across numerical, analytical and experimental findings. Ansys Fluent® software is used to assess temperature, velocity and streamlines, yielding satisfactory results from an engineering standpoint.

Investigating critical parameters such as Darcy number (10⁻⁴ ≤ Da_D ≤ 10⁻²), aspect ratio (0.5 ≤ H/D ≤ 1.5), Reynolds number (5 ≤ Re_D,bf ≤ 3500), Eckert number (0 ≤ EC_bf ≤ 0.1) , porosity (0.85 ≤ ε ≤ 0.95), Hartmann number (0 ≤ Ha_D,bf ≤ 300) and the volume fraction of nanofluid (0 ≤ φ ≤ 0.1) reveals their impact on fluid flow and heat sink performance. Notably, Nusselt number will reduce 45%, rise 19.2%, decrease 14.1%, and decrease 0.15% for Reynolds numbers of 600, with rising porosity from 0.85 to 0.95, Darcy numbers from 10⁻⁴ to 10⁻², Eckert numbers from 0 to 0.1, and Hartman numbers from 0 to 300.

Despite notable progress in studying thermal management in CPU cooling systems using porous media and nanofluids, there are still significant gaps in the existing literature. First, few studies have considered the Darcy–Brinkman–Forchheimer equation, which accounts for non-Darcy effects and the flow and geometric interactions between coolant and porous medium. The influence of viscous dissipation on heat transfer in this specific geometry has also been largely overlooked. Additionally, while nanofluids and impinging jets have demonstrated potential in enhancing thermal performance, their utilization within porous media remains underexplored. Furthermore, the unique thermal and structural characteristics of porous media, along with the incorporation of a magnetic field, have not been fully investigated in this particular configuration. Consequently, this study aims to address these literature gaps and introduce novel advancements in analytical modeling, non-Darcy flow, viscous dissipation, nanofluid utilization, impinging jets, porous media characteristics and the impact of a magnetic field. These contributions hold promising prospects for improving CPU cooling system thermal management and have broader implications across various applications in the field.