UFR 1-07 Best Practice Advice

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Unsteady Near-Field Plumes

Underlying Flow Regime 1-07

Best Practice Advice

Best Practice Advice for the UFR

Key Physics

The key physics of this UFR is the transient, unsteady behaviour in the near-field of a turbulent buoyant helium-air plume. The flow features two key instabilities. Firstly, the Rayleigh-Taylor instability related to the presence of dense fluid above less-dense fluid, which gives rise to fingers or spikes of dense fluid separated by rising bubbles of lighter fluid. Secondly, the Kelvin-Helmholtz instability related to the shear-layer interface between the rising plume and the ambient fluid, which produces roll-up vortex sheets on the boundary between the two layers of fluid travelling at different velocities. The flow is very challenging to predict using CFD, due to the sharp density gradients at the plume exit which produce flow conditions where small scales of turbulent motion feed into the larger scales.

Numerical Modelling

• For LES, the flow cannot be treated as two-dimensional or axisymmetric. Full three-dimensional time-dependent simulations must be performed.
• For simulation of the selected UFR test case, open boundaries should be used on all sides of the flow domain except for the floor. Constant pressure boundaries may be used, although if a fully-compressible code is used, care will need to be taken to ensure that the boundaries are non-reflective.
• For simulation of the selected UFR test case, the domain should extend at least 4 metres radially and vertically to minimize any effects of the open boundaries on the development of the plume. Ideally, tests should be performed to ensure that the location of the open boundaries has no significant effect on the results.
• The finest mesh should be used given the available computing resources. The results discussed above suggest that a mesh of around 4 million nodes should give good agreement with the experiments in terms of mean flow quantities, but may still be insufficient for good predictions of fluctuations or RMS values. Tieszen noted that at least 75 cells across the base diameter of the plume are necessary to avoid significant differences in the vertical centreline velocity compared to the measured values. Ideally, a grid-dependence study should be undertaken to investigate the magnitude of these effects.

Physical Modelling

• Either the fully-compressible or the low-Mach-number form of the Navier-Stokes equations can be used. The fully-compressible N-S equations require more careful treatment to avoid acoustic waves reflecting back into the domain from open boundaries. Furthermore, they will require a very short time-step, based on the speed of sound instead of the local flow speed, unless special treatments are used. For details of a fully-compressible N-S treatment, see DesJardin et al. [1].

• The baroclinic torque is non-zero and therefore should not be neglected.

Footnotes

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Contributed by: Simon Gant — UK Health & Safety Laboratory

© copyright ERCOFTAC 2010