AbstractThis work presents a multi‐phase approach to model freezing processes within a coupled finite element formulation using the Theory of Porous Media (TPM), see also [1]. The focus of the study is to enhance the numerical stability for the volume fraction ice as is it highly sensitive to phase discontinuities. To address this issue, a discontinuous approximation is employed for the volume fraction ice while the remaining degrees of freedom are still continuously approximated. As a result, the corresponding degrees of freedom are independently solved for each element, as the element nodes are not interconnected with other elements. Academic examples are provided within a standard Galerkin procedure to highlight the optimization effect, comparing the results obtained with a continuous and discontinuous approximation for the volume fraction ice. Additionally, the implementation of the discontinuous approximation allows for the application of static condensation, resulting in reduced computation time.
Read full abstract