22–26 Jul 2024
CICG - GENEVA, Switzerland
Europe/Zurich timezone

Modelling heat ingress in LH2 containment systems – benchmarking improved thermal resistance network models against the finite element method

25 Jul 2024, 14:00
2h
Poster area

Poster area

Poster Presentation (120m) ICEC 12: Thermal properties and numerical studies Thu-Po-3.4

Speaker

Sindre Stenen Blakseth (SINTEF Energy Research / Norwegian University of Science and Technology)

Description

Hydrogen is projected to play an important role in the decarbonization of industries that are not readily electrifiable, such as aviation, shipping, and high-energy process industries [1]. To this end, effective solutions for hydrogen transport and storage are needed at a scale far exceeding current capacity. As the hydrogen value chain expands, a wide range of computational modelling tools are required to ensure that the deployed technology is satisfactory in terms of thermo-dynamic efficiency, economic profitability, and safety. In the present work, we discuss and compare two modelling tools designed to predict heat ingress and boil-off generation in liquid hydrogen containment systems.

In the first tool, the finite element method (FEM) is used to solve the steady heat equation on a tetrahedral mesh that accurately represents the geometry of the containment system. This tool, which we have implemented using the open-source NGSolve library [2], provides high predictive accuracy, but it is also computationally intensive. These characteristics imply that the tool is most suited for detailed assessments of the containment system in isolation, as computational cost is then not a critical concern.

For more comprehensive studies, where several value chain components are considered simultaneously, the computational cost of modelling each component must be kept at a minimum. Then, an alternative to the FEM-based tool is desirable. This motivates the second modelling tool we consider: a state-based thermal network model that exploits the analogy between electrical and thermal conduction [3]. In this model, the various components of the containment system are represented by nodes connected by thermal resistances that are computed, e.g., using exact solutions of the heat equation for simple geometries, such as spherical shells and cylinders. With this approach, some local information about the temperature field within the system is lost, but high-level performance indicators such as total heat ingress can still be predicted with reasonable accuracy.

In this work, we present novel augmentations of the network model that significantly improve the model's predictive accuracy without any adverse increase in model complexity. One such augmentation is a semi-analytic estimate of the temperature field surrounding the cold spot at the intersection between the support structure and the outer tank. We also discuss how to handle geometrically complex support structure components such as a connection between a support skirt and the inner tank. The FEM tool is used to provide empirical error bounds for the network model across a range of material choices and support structure geometries. Based on these bounds, we establish that the augmented network model is sufficiently accurate for use in more comprehensive studies of liquid hydrogen carrier ships, including optimization of ship design and reliquefaction capacity. Results from a case study on 40 000 m3 capacity LH2 tanks – conducted in the LH2 Pioneer project (Research Council of Norway grant 320233) – will be presented in order to demonstrate the efficacy of the modelling tools discussed.

References
[1] IPCC, 2022: Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [P.R. Shukla, J. Skea, R. Slade, A. Al Khourdajie, R. van Diemen, D. McCollum, M. Pathak, S. Some, P. Vyas, R. Fradera, M. Belkacemi, A. Hasija, G. Lisboa, S. Luz, J. Malley, (eds.)]. Cambridge University Press, Cambridge, UK and New York, NY, USA. doi: 10.1017/9781009157926
[2] J. Schoeberl, C++11 implementation of finite elements in NGSolve, Vienna University of Technology - Institute for Analysis and Scientific Computing, (2014).
[3] A. Majumdar, T. Steadman, J. Maroney, J. Sass, J. Fesmire, Numerical modeling of propellant boil-off in a cryogenic storage tank, in: AIP Conference Proceedings, Vol. 985, American Institute of Physics, (2008), pp. 1507–1514

Submitters Country Norway

Primary authors

Sindre Stenen Blakseth (SINTEF Energy Research / Norwegian University of Science and Technology) Ailo Aasen (SINTEF Energy Research) Magnus Aa. Gjennestad (SINTEF Energy Research)

Co-author

André Massing (Norwegian University of Science and Technology)

Presentation materials