Please use this identifier to cite or link to this item: http://bura.brunel.ac.uk/handle/2438/32634
Title: Effective thermal conductivity prediction model for vacuum insulation cores
Authors: Kaushik, D
Nunna, A
Singh, H
Keywords: thermal conductivity;vacuum insulation;powder material;finite element;MATLAB
Issue Date: 7-Jan-2026
Publisher: Elsevier
Citation: Kaushik, D., Nunna, A. and Singh, H. (2026) 'Effective thermal conductivity prediction model for vacuum insulation cores', Results in Engineering, 29, 109035, pp. 1 - 10. doi: 10.1016/j.rineng.2026.109035.
Abstract: Accurate prediction of thermal conductivity of porous granular materials enables the identification and rapid optimisation of new composite core materials for Vacuum Insulation Panels (VIPs). To date, no computer model has reported the use of multi-sized particles with a combined Finite Element Analysis (FEA) and Discrete Element Method (DEM) approach to predict thermal conductivity of VIP core. To fill this knowledge gap, we propose a FEA and DEM based thermal conductivity prediction model for powdery composites, particularly suited for VIPs. The geometry of the model was formed using random packing of multisized spherical particles using DEM and a MATLAB PDE-based thermal model solver was used to obtain a solution for the generated geometry. The model can predict effective thermal conductivity whilst accounting for the VIP-specific fundamental heat exchange phenomena. The results from the model are validated against those obtained from accompanying thermal conductivity measurements performed using the Transient Hot Wire method, with results falling within the error range for temperatures <491.15 K. Effective thermal conductivity of a perlite and Silicon Carbide (SiC) core at 0.1 mbar over a temperature range of 303 K to 803 K, with the proportion of perlite varying from 100 % to 50 % (by weight), as predicted by the model, is presented. The thermal conductivity of the 50 % perlite-50 % SiC composite had the lowest rate of increase of thermal conductivity of 46.3 %, with the value increasing from 12.1 mW m-1 K-1 at 303 K to 17.7 mW m-1 K-1 at 803 K.
Description: Data availability: Data will be made available on request.
Supplementary materials are available online at: https://www.sciencedirect.com/science/article/pii/S2590123026000782#sec0012 .
URI: https://bura.brunel.ac.uk/handle/2438/32634
DOI: https://doi.org/10.1016/j.rineng.2026.109035
Other Identifiers: ORCiD: Harjit Singh https://orcid.org/0000-0003-3448-1175
Article number: 109035
Appears in Collections:Dept of Mechanical and Aerospace Engineering Research Papers

Files in This Item:
File Description SizeFormat 
FullText.pdfCopyright © 2026 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license ( https://creativecommons.org/licenses/by/4.0/ )3.49 MBAdobe PDFView/Open


This item is licensed under a Creative Commons License Creative Commons