Numerical description of axisymmetric blue whirls over liquid-fuel pools

This numerical investigation is focused on determining the structures of blue whirls, recently found to occur in laboratory investigations of fire whirls when the circulation becomes sufficiently large to produce a vortex breakdown that drastically shortens the fire whirl and correspondingly reduces residence times, so that the yellow flames turn blue. The computations address axisymmetric configurations for round pools of liquid fuels flush with and at the center of a larger solid horizontal disc, at the outer edge of which vanes of adjustable angles cause the entrained air to enter with a controllable azimuthal component of velocity. The nondimensionlized conservation equations employed include realistic Lewis numbers with temperature-dependent transport coefficients and a one-step chemical-kinetic approximation that correctly reproduces laminar burning velocities. Buoyancy and radiant energy transport from the flames to the liquid surface are both taken into account, the latter being found to be essential for the blue whirl. Along with the vaporization-equilibrium and energy-conservation boundary conditions at the fuel surface, inflow boundary conditions are provided by a recently developed solution for the boundary-layer flow over the solid disc, while zero-gradient outflow conditions are applied above the whirl. Controlling nondimensional parameters, besides Reynolds, Damköhler, and Froude numbers, are a ratio of radiant to convective energy flux and a ratio of azimuthal to inward radial flow velocity in the boundary layer at the edge of the disc. The computed conditions for the onset of the blue whirl, as well as the computed structure of the whirl itself, bear close resemblance to what was found experimentally.