Optimisation of crust freezing in meat processing via Computational Fluid Dynamics

TL;DR

This paper presents a computationally efficient CFD model for continuous impingement freezing of meat products, enabling detailed analysis of freezing behavior and optimal process parameters with minimal computational resources.

Contribution

A new simplified CFD solver integrated into OpenFOAM for modeling phase change in continuous meat freezing processes, capturing complex freezing dynamics efficiently.

Findings

Maximum freezing achieved at jet distance H=7.2D

Optimal conveyor speed corresponds to Pe_s=8244

Model captures nonlinear freezing behavior due to phase change

Abstract

In this work a numerical model for two-dimensional axisymmetric continuous freezing by impingement of processed meat or similar products in food industry moving along a conveyor belt is presented. The model represents a more computationally efficient alternative to solve conjugate heat transfer between a fluid and a solid, accompanied by phase change in some constituents of the solid phase. In the model presented here it is assumed that the solid can be represented as an homogeneous medium, with its thermophysical properties depending on the temperature. The impingement freezing model is conceived to be valid for highly processed vegetarians products or meat such as sausages, mince or ham freezing. Furthermore, this approach is much simpler in terms of computational cost whilst it still captures the complexity of continuous freezing under industrial setting. The methodology is implemented as a new solver in the widely used open-source Computational Fluid Dynamics (CFD) library OpenFOAM. Overall, highly non-linear freezing behaviour was found due to the phase change inside the solid and the associated heat of fusion. We studied the effect of high fluid Reynolds numbers as well as investigating the optimal distance between the jet and the solid surface for different speeds of the conveyor. We found that the maximum freezing is obtained positioning the jet at a distance H = 7.2D (where D is the diameter of the impinging jet) and setting the speed of the conveyor such that the P\'eclet number of the solid is $\text{Pe}_{\text{s}}=8244$. The methodology developed allows to obtain detailed insight on the freezing process for various impingement configurations at a minimum computational cost using a freely available open-source tool.