Comparison of dns of compressible and incompressible turbulent droplet-laden heated channel flow with phase transition

In this paper a turbulent channel flow with dispersed droplets is examined. The dispersed phase is allowed to have phase transition, which leads to heat and mass transfer between the phases, and correspondingly modulates turbulent flow properties. As a point of reference we examine the flow of water droplets in air, containing also the vapor of water. The key element of this study concerns the treatment of the carrier phase as either a compressible or an incompressible fluid. We compare simulation results obtained with a pseudo-spectral discretization for the incompressible flow to those obtained with a finite volume approach for the compressible flow. The compressible formulation is not tailored for low Mach flow and we need to resort to a Mach number that is artificially high for simulation feasibility. We discuss differences in fluid flow, heat- and mass transfer and dispersed droplet properties. The main conclusion is that both formulations give a good general correspondence. Flow properties such as velocity fields agree very closely, while heat transfer as characterized by the Nusselt number differs by around 25%. Droplet sizes are shown to be slightly larger, particularly in the center of the channel, in case the compressible formulation is chosen. A low-Mach compressible formulation is required for a fully quantitative comparison

The effect of phase transitions on the droplet size distribution in homogeneous isotropic turbulence

We investigate the dynamics of an ensemble of discrete aerosol droplets in a homogeneous, isotropic turbulent flow. Our focus is on the stationary distribution of droplet sizes that develops as a result of evaporation and condensation effects. For this purpose we simulate turbulence in a domain with periodic boundary conditions using pseudo-spectral discretization. We solve in addition equations for the temperature and for a scalar field, which represents the background humidity against which the size of the droplets evolves. We apply large-scale forcing of the velocity field to reach a statistically steady state. The droplets are transported by the turbulent field while exchanging heat and mass with the evolving temperature and humidity fields. In this Euler-Lagrange framework, we assume the droplets volume fraction to be sufficiently low to allow one-way coupling of the droplets and turbulence dynamics. The motion of the droplets is time-accurately tracked. The Stokes drag force is included in the equation of motion of the individual droplets. The responsiveness of the droplets to small turbulent scales is directly related to the size of the individual spherical droplets. We perform direct numerical simulation to ultimately obtain the probability density function of the evolving radius of the droplets at different points in time with characteristic heat and mass transfer parameters. We determine the gradual convergence of the distribution function to its statistically stationary state for forced homogeneous, isotropic turbulence

Improvement of heat- and mass transfer modeling for single iron particles combustion using resolved simulations

In this work, we use a boundary layer resolved model to improve a Lagrangian point particle model to simulate the combustion of single iron particles. By resolving the full boundary layer, mass and heat transfer are accurately modeled, including Stefan flow. Therefore, the model is suitable to improve point particle models. This work focuses on the first stage of iron combustion, which lasts up to the maximum temperature. Temperature- and composition-dependent properties are used and phase transitions from solid to liquid and liquid to gas are taken into account. The Nusselt and Sherwood correlations are investigated in conditions typical for iron particle combustion. It is found that the 1/2-film temperature is the best film rule to use to model heat- and mass transfer for iron particle combustion. The boundary layer resolved model is used to validate the point particle models. Then, the model is systematically elaborated by including a temperature-dependent particle density, slip velocity and Stefan flow. The individual and combined effect of these phenomena on the burn duration are investigated. Including all these effects decreases the time to maximum temperature by around 25%. Furthermore, it is shown that if one neglects physical phenomena like slip and Stefan flow, but uses the 1/3-film rule instead of the 1/2-film rule, errors cancel and still reasonable agreement is obtained with experiments