Implementation, demonstration and validation of a user-defined wall-function for direct precipitation fouling in ANSYS Fluent

TL;DR

This paper presents a new user-defined wall-function for simulating precipitation fouling in CFD, implemented in ANSYS Fluent, validated against experimental and detailed CFD data for heat exchangers.

Contribution

It introduces a computationally efficient 1D boundary-layer model as a UDF in ANSYS Fluent for industry-scale fouling simulations.

Findings

Good agreement with experimental deposition rates

Model reduces computational cost

Effective for industry-scale heat exchanger fouling

Abstract

In a previous paper (Johnsen et al., 2015) and presentation (Johnsen et al., 2016), we developed and demonstrated a generic modelling framework for the modelling of direct precipitation fouling from multi-component fluid mixtures that become super-saturated at the wall. The modelling concept involves the 1-dimensional transport of the fluid species through the turbulent boundary layer close to the wall. The governing equations include the Reynolds-averaged (RANS) advection-diffusion equations for each fluid species, and the axial momentum and energy equations for the fluid mixture. The driving force for the diffusive transport is the local gradient in the species' chemical potential. Adsorption mechanisms are not modelled per se, but the time-scale of adsorption is reflected in the choice of Dirichlet boundary conditions for the depositing species, at the fluid-solid interface. In this paper, the modelling framework is implemented as a user-defined function (UDF) for the CFD software ANSYS Fluent, to act as a wall boundary condition for mass-transfer to the wall. The subgrid, 1-dimensional formulation of the model reduces the computational cost associated with resolving the fine length-scales at which the boundary-layer mass transfer is determined, and allows for efficient modelling of industry-scale heat exchangers suffering from fouling. The current paper describes the modelling framework, and demonstrates and validates its applicability in a simplified 2D heat exchanger geometry (experimental and detailed CFD modelling data by P\"a\"akk\"onen et al. (2012, 2016)). By tuning the diffusivity, only, good agreement with the experimental data and the detailed CFD model was obtained, in terms of area-averaged deposition rates.