Hex-PORTMAN: Heat Flux Splitting in Porous Materials for Thermal Management

Lead Research Organisation: Queen's University Belfast
Department Name: Sch Mechanical and Aerospace Engineering

Abstract

Thermal management plays a vital role in determining the efficiency, safety and reliability of technological development in a plethora of industries including aerospace, automotive, computing and renewable energy sectors. The developments in these industries have culminated in a considerable surge in the power densities, which goes hand-in-hand with the increase of generated heat flux and subsequent undesirable temperature rise in system components. Porous materials (i.e. solids, which are permeated by a network of pores) have been demonstrated to be competitive microfluidic materials for effective cooling in high heat flux applications because of their fluid permeability and high surface area, which augments the heat transfer from hot surfaces to the cooling fluid passing through the porous media.

Past studies have theoretically investigated the flow and thermal characteristics of the porous media systems for thermal management using the volume-averaged approach, which is a popular low-cost engineering approach for studying transport in porous media. However, after more than a decade of research, this problem has still not been resolved. This is primarily because the splitting mechanism of the external heat flux between the solid and fluid phases in the porous media is unknown and determination of the thermal boundary condition for volume-averaged solvers remains a scientific challenge. This ambitious project will, for the first time, address this fundamental problem of flow and heat transfer in porous media systems through a comprehensive series of experimental and modelling studies.

This project will benefit from partnership with world-renowned scientists: Prof Kambiz Vafai (KV)-University of California Riverside, Dr Mahdi Azarpeyvand (MA)-University of Bristol and Prof Kamel Hooman (KH)-University of Queensland, with the involvement of one PDRA and four PhD students. KV is a world-leading scientist in the field of transport in porous media and will bring his key knowledge in understanding the heat flux splitting in the porous media. MA and KH will support the project for experimental measurements of the velocity field in the system. This project is also of direct relevance to industry with the involvement of UK-based companies (Glen Dimplex, B9 Energy, and BL Refrigeration) who will be deploying the fully validated volume-averaged solver developed in the project for the purpose of thermal management using porous materials in application to electronics cooling, energy storage and solar photovoltaic systems, respectively.

Planned Impact

Engagement with industry during and after this project will be the key to delivering economic impact and eventually getting products to market that will generate profit within the UK. The first project beneficiaries are the UK manufacturing sectors working in the field of thermal management in electronics for ground and space applications. In this sector, implementation of our low-cost modelling tool will allow heat exchange devices such as metal foam heat sinks to be designed more reliably. It has the potential to increase the cooling capacity of the electronic devices by 15% compared to the conventional finned heat sinks, which will increase the life-span of the devices and lower costs for end-users. This is clearly recognised by our industrial project partner, Glen Dimplex who will use the developed solver for their future design of thermal management devices for their electronics appliances. Additionally, the project has attracted the interests of UK-based manufactures in Clean Energy: (1) B9 Energy has recognised the project's benefit in manufacturing a novel liquid piston gas compressor (LPGC) for application to compressed energy storage system. As demonstrated in the letter by B9, deployment of the project's solver allows for an optimum design of a LPGC technology, which can have about 18% higher efficiency than the conventional LPGC system; (2) BL Refrigeration shows high interests in utilising the project's outcomes for an optimum design of cooling systems for Solar Photovoltaic (PV) panels. Such a cooling technology has the potential to increase the electricity production of PV panels by 10% over its 25 years life for a typical domestic household in the UK climate.

Societal impact will be realised by the application of the project's outcomes to large-scale applications in energy storage and solar energy. A recent study showed that by 2030, UK can meet emission reduction goals and electricity need predominantly using solar and energy storage along wind that can provide more than 60% of the total UK's electricity. Solar alone is estimated to create about 50,000 jobs per year, contributing £25.5 billion in GVA to the UK economy and putting £425 million back into consumers pockets. Thermal management using cost-effective materials such as metal foams have been specifically identified as technology drivers in solar PV and energy storage systems. This project will enhance our fundamental understanding of the thermal characteristics of flow in porous materials, which will contribute to this growing industry by improving the efficiency of the solar PVs and compressed air energy storage systems using metal foams. Additionally, the environmental impact of the project will be realised through supporting the UK's target in reducing GHG emissions by 80% by 2050, since instead of burning fossil fuels, solar PVs can be deployed for electricity production and compressed air energy storage for storing the surplus electricity for use on demand.

The additional benefits of the project to the general public will be through contributing to the training and education of UK engineers with the necessary skills in the area of experimental and computational fluid dynamics, which is an area of priority for EPSRC. The PI will directly benefit by expanding his research ambition, developing a strong international network, and strengthen his industrial and academic connections for future grant applications, which will lead to the future training of many more students and PDRAs. Pupils and students inspired by the research will seek careers in engineering. Clean and renewable energies are naturally inspiring to many students and the wider public. If successful, we envisage that the research findings will be incorporated into the undergraduate course at Queen's University, where each year, six final year students will benefit by spreading the project findings via interacting with the PI's current activities in supervising final year projects.
 
Description At the first stage of this study, a pore-scale large eddy simulation investigation has been conducted to better understand the flow features of a channel partially filled with a porous block. For this purpose, an scalable open-source solver and utility based on OpenFOAM have been developed that is capable to simulate passive scalar convection and diffusion for pore-scale turbulence modelling (RANS and LES) of porous structures.

Numerical results illustrate that some portion of the fluid entering the porous blocks leave the porous structure to the clear region (positive leakage) through the porous-fluid interface. Moreover, it is argued that the boundary layer at the porous-fluid interface is continuously disrupted by the positive leakage across the interface. This is due to the negative vertical pressure gradient in the porous region, causing the flow to move up into the clear region through the interface. For the first time, it is explained that positive leakage has a noticeable effect on the coherent structures, turbulent intensity, and turbulent kinetic energy above the porous-fluid interface. Moreover, the results show that although the fluid flow within the porous block experiences severe gradient, the turbulent velocity fluctuations are limited due to the small-scale pore size there.

Moreover, the complexity of flow physics and numerical methods have opened several important research questions especially about the:
Grid resolution assessment in pore-scale large eddy simulation (LES) of porous media.
Mesh generation of complex geometry for porous structures.
Development of a numerical method for increasing the speed up of the computational expensive LES approach.
Origination of complex turbulent flow structures within and over a porous block.

Therefore, in the project, we are in close collaboration with different leading scientists from the University of Bristol (UoB-UK), University of Queensland (UQ-Australia), and University of California Riverside (UCR-US) in the field of transport in porous media. For code development and high-performance computing based on OpenFOAM, we also benefit from the technical advice of other scientists from Chalmers University of Technology (Sweden) and University of Windsor (Canada) that has expanded the research networks greatly.
Exploitation Route • Development of open-source code based on OpenFOAM that is scalable on HPC clusters
• Presenting the achievements in seminars and workshops. For example, the preliminary achievements of this project are going to be presented in a talk at the University of Manchester in May 2021.
• Being Present at different (face-to-face and virtual) conferences. For example, we have presented our pore-scale methodology for porous media modelling in the virtual "International Symposium on Numerical Methods in Heat and Mass Transfer 2020, ISNMHMT2020 11-13 Dec, 2020, Ningbo, China" and, two conference papers have been submitted for UKHTC2021(https://www.ukhtc2021.org/) and ATE-HEFAT 2021(https://hefat2021.org/).
• Submitting journal papers. For instance, a journal paper about the characteristic flow features of a channel partially filled with a porous block is now under preparation to be submitted to share our new findings with other researchers in this field
Sectors Aerospace, Defence and Marine,Agriculture, Food and Drink,Electronics,Energy,Environment

 
Description 'This grant has been transferred to the University of Manchester, please see EP/T012242/2
Sector Energy
Impact Types Economic

 
Title Turbulent kinetic energy post-process utility for LES in OpenFOAM 
Description Description: The developed utility in OpenFOAM source code calculates and writes some mean fields from an LES calculation. The following fields are created: - resLES: resolutness of an LES - TKEMean: mean turbulent kinetic energy - TKEMeanProd: production term of the mean turbulent kinetic energy - turbDiffusionMean: turbulent diffusion term of the mean turbulent kinetic energy - SGSDiffusionMean: subgrid-scale diffusion term of the mean turbulent kinetic energy - viscDiffusionMean: viscous diffusion term of the mean turbulent kinetic energy Fields UMean, UPrime2Mean, turbDiffMean, SGSDiffMean, and kMean must exist. 
Type Of Material Data analysis technique 
Year Produced 2021 
Provided To Others? No  
Impact This utility can be used for the post-processing of LES results. Especially, the influence of each term in the budget of turbulent kinetic energy (TKE) can be extracted from the LES results by this utility. It can also be used for validation of any sub-grid scale LES model based on the balance of different terms in the TKE transport equation. 
 
Title Passive scalar solver in OpenFOAM based on the RANS and LES turbulence models 
Description Description: Transient solver for incompressible, turbulent flow of Newtonian fluids. Heading Solver details: The solver uses the PIMPLE (merged PISO-SIMPLE) algorithm. Sub-models include: - turbulence modelling, i.e. laminar, RAS or LES; - run-time selectable MRF and finite volume options, e.g. explicit porosity. Solver capability: 1. This solver can also be used for Post-processing of turbulence results. 2. It is the solver plus two essential utilities to calculate budgets of the mean turbulent kinetic energy equation, turbulent heat flux, and wall heat flux for the incompressible solver. 3. Comments have been added to the main solver which should make it easy to understand the terms of TKE budgets and heat flux. 4. Both solver and utilities ("mjPostProcessTKEBudgestsLES" and "mjPostProcessMeanWallHeatFluxIncompressible") have been compiled on OF v2012 and v2006. 5-1- The implementation requires a well-established averaged Velocity field (UMean and TMean) to exist in the time folder. The process of averaging of TKE budgets happens within the main solver and UMean and TMean are NOT updated within the code. So a reasonable averaged field must exist prior to beginning of simulation using "mjPassiveScalarPimpleFoam". To provide this averaged field, one can perform simulation for several flow-through time using the native "mjPassiveScalarPimpleFoam" solver. 5-2-Three of TKE budgets are calculated afterward the simulation through the "mjPostProcessTKEBudgestsLES" utility provided. This utility requires averaged fields, including: -UMean -kMean -UPrime2Mean -turbDiffMean -SGSDiffMean If any of the above fields would not exist, the utility will be terminated. 5-3-turbulent heat flux (UPrime*TPrime) should be time-averaged. 5-4-wallHeatFlux are calculated afterward the simulation through the "mjPostProcessMeanWallHeatFluxIncompressible" utility provided. This utility requires averaged fields, including: -TMean -alphatMean NOTE: The solver is now under validation procedure. 
Type Of Technology New/Improved Technique/Technology 
Year Produced 2021 
Impact The solver has been developed in OpenFOAM code that is very general in the CFD community and can be extended by other researchers. It also supports different turbulence models that make it very helpful for different research fields. 
 
Description A NUMERICAL STUDY OF ACTIVE AND PASSIVE HEAT TRANSFER ENHANCEMENT IN LATENT HEAT THERMAL ENERGY STORAGE DEVICES 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Study participants or study members
Results and Impact In this talk, I am going to present my latest findings on the application of porous media in "Thermal Energy Storage" design.

Latent Heat Thermal Energy Storage (LHTES) devices implement phase-change materials (PCM) to store and release the thermal energy from the latent heat of fusion of a material. PCMs have high energy densities and can be effective at storing a large amount of thermal energy. This is especially useful when used in conjunction with renewable energy sources to balance energy production with demand. One of the disadvantages of PCMs, however, is their poor thermal conductivity, inhibiting the LHTES device's melting and solidification performance. This numerical study aims to research both active and passive heat transfer methods to improve the melting performance of a PCM, making them more effective when used in thermal energy storage.
Year(s) Of Engagement Activity 2022
URL https://www.ukhtc2021.org/
 
Description Flow Interaction Between Porous and Non-porous region in a Channel Partially Filled with a Porous Block: Pore-scale LES Study, Proceedings of the 7th World Congress on Momentum, Heat and Mass Transfer (MHMT'22), 07 - 09 April 2022 (Virtual Conference) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact In this talk, I want to discuss the flow behaviour inside and outside of a porous block that is interested in industry for flow control and heat transfer enhancement.
Abstract of the talk:
The present study investigates fluid flow interaction between porous and non-porous regions in a channel partially filled with a porous block. For this purpose, a detailed pore-scale large eddy simulation is utilized. Flow visualization shows that some portion of the fluid entering the porous blocks is pushed upwards to the porous-fluid interface and leaves the porous region; this phenomenon is called flow leakage. Spectral analysis of vertical velocity and correlation coefficients confirm the flow leakage. Below the porous interface, the magnitude of correlation coefficients exposes a strong positive correlation between vertical velocity fluctuations that reveals the upward tendency of flow in the porous region. This trend is also observed across the porous interface which confirms momentum transfer through the porous interface. Moreover, spectral analysis of vertical velocity reveals that the dominant frequencies within the porous region exist in the non-porous region where the flow leakage is pronounced. This observation shows strong momentum transfer between the porous and non-porous regions due to the flow leakage.
Year(s) Of Engagement Activity 2022
URL https://mhmtcongress.com/
 
Description International Symposium on Numerical Methods in Heat and Mass Transfer 2020, ISNMHMT2020 11-13 Dec, 2020, Ningbo, China 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact Turbulent flow in porous media is ubiquitous in nature and its understanding is central for unravelling the physics underlying the natural phenomena in food drying, soil evaporation as well as fostering technological development in engineering applications including packed bed reactors, energy storage and electronics cooling. In these natural and engineering applications, a step change in the fundamental understanding of flow in composite porous-fluid systems, which consists of a fluid-saturated porous medium and a flow passing over it is crucial for characterisation and diagnostic analysis of such complex problems.
In this talk, the flow behaviour for a channel partially filled with a porous block was explained to the audience of the symposium. Moreover, I explained different methodologies for modelling porous media namely: volume averaging method and pore-scale simulation. the main objectives of the talk were as follows:
1. The application of pore-scale CFD simulation of porous medium
2. Difference between flow behaviour of a channel partially filled with a porous block and a full porous channel.
The CFD results were very interesting for the audience since the pore-scale flow pattern within the porous structure had been analysed directly.
Year(s) Of Engagement Activity 2020
 
Description Large Eddy Simulation of Turbulent Flow Over Fluid-Saturated Porous Media, 17th UK Heat Transfer Conference, UKHTC2021, (Virtual Conference) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Study participants or study members
Results and Impact At this conference, I am going to talk about the following concepts that are related to heat transfer enhancement using porous materials.
Abstract of the talk:
In this paper, the heat and fluid flow behaviour for a channel partially filled with two different porous blocks and a solid block is investigated numerically using a detailed pore-scale large eddy simulation (LES). Results show that a significant portion of the fluid entering the porous blocks leaks from the porous region to the non-porous region through the porous-fluid interface. Detailed flow patterns visualization inside the porous blocks indicate that the flow leakage leads to the creation of counter-rotating vortex pairs of fluid flow within and above the porous blocks that result in the formation of organized hairpin structures. Investigation of coherent structures show that the legs and the head of hairpins are generated separately by different mechanisms; the legs are produced by flow leakage while the heads are created by vorticity rollup at the shear layer above the porous interface. Finally, results show that the heat transfer mechanisms over the porous media are governed by the formation of organized hairpin structures.
Year(s) Of Engagement Activity 2022
URL https://www.ukhtc2021.org/
 
Description Pore Scale Large Eddy Simulation of Flow in a Channel Partially Filled with Porous Block, 15th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics (ATE-HEFAT 2021), 25-28 July 2021 (Best paper award, Virtual Conference) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Professional Practitioners
Results and Impact In this talk, I presented my latest finding on the application of porous material in both natural and engineering applications.
Abstract of the talk:
Present work examines the turbulent feature of a flow passing over a fluid-saturated porous media in a channel using a detailed pore-scale large eddy simulation approach. Two different porous configurations (i.e., a packed bed of sphere and a porous block formed from rectangular cross-section ligament) are analysed and the results are compared against those obtained for the case when a solid block with the same size of the porous block is inserted in the channel. Results show that some portion of the fluid entering the porous blocks leave the porous structure to the clear region (positive leakage) through the porous-fluid interface. This is due to the negative vertical pressure gradient inside the porous region. The flow leakage leads to the formation of counter-rotating vortex pairs of the flow within and above the porous region, altering the coherent structures of the flow above the interface. Results indicate that although the fluid flow within the porous is highly vortical and experiences intense gradient, the turbulence fluctuations are limited because of the existence of porous ligaments and small-scale pores size. In comparison to the solid block case, the magnitude of the turbulent kinetic energy is decreased drastically above the porous block, since the reverse flow on the porous-fluid interface is deteriorated by the positive leakage as compared with the solid block.
Year(s) Of Engagement Activity 2021
URL https://hefat2021.org/
 
Description Pore-scale large-eddy simulations of turbulent flow in a composite porous-fluid system, UK Fluids Conference, 8-10 September 2021 (Virtual Conference) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach National
Primary Audience Postgraduate students
Results and Impact At this conference, I presented my latest findings of the flow behaviour in and over a porous block for researchers in this field of study. The conference was held virtually.
Abstract of the talk:
In the present study, the flow features for a channel partially filled with two different porous blocks and a solid block are investigated numerically using a detailed pore-scale large eddy simulation (LES) approach. The effect of the porous block on the coherent structures and large reverse flow patterns along the top surface of the blocks are investigated and compared with the solid block. Moreover, the characteristics of the wake flow and the flow leakage from the porous region to the clear flow through the porous-fluid interface are examined in detail.
Pore-scale LES results show that some portion of the fluid entering the porous blocks is pushed upwards toward the porous-fluid interface due to the negative vertical pressure gradient inside the porous block and leaves the porous structures to the clear region (positive leakage). This tendency of flow leads to the formation of counter-rotating vortex pairs (CRVPs) of the fluid flow within and above the porous block and changes the coherent structures above the porous-fluid interface. It is also argued that the magnitudes of turbulent kinetic energy and turbulent velocity fluctuations are decreased drastically over the porous block as compared with the solid block since the reverse flow on the porous-fluid interface is deteriorated by the positive leakage.
Year(s) Of Engagement Activity 2021
URL https://www.southampton.ac.uk/~assets/doc/comms%20and%20marketing/2021_EngineeR_AbstractBooklet.pdf