UFR 1-07 Best Practice Advice: Difference between revisions
Jump to navigation
Jump to search
| Line 46: | Line 46: | ||
* 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<ref name="ftn13">S. Tieszen, Private Communication, March 2010.</ref> 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. | * 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<ref name="ftn13">S. Tieszen, Private Communication, March 2010.</ref> 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 === | === Physical Modelling === | ||
| Line 53: | Line 54: | ||
* The baroclinic torque is non-zero and therefore should not be neglected. | * The baroclinic torque is non-zero and therefore should not be neglected. | ||
* The Boussinesq approximation, where flow properties are assumed not to vary as a function of temperature or composition, and where buoyancy is only included as an additional body-force term in the momentum equations, should not be used. The Boussinesq approximation is only appropriate for modelling small density differences, equivalent to a temperature difference less than around 15°C in air [[UFR_1-07_References#80|[80]]]. | |||
=== Application Uncertainties === | === Application Uncertainties === | ||
Revision as of 10:49, 13 July 2010
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[1] 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.
- The Boussinesq approximation, where flow properties are assumed not to vary as a function of temperature or composition, and where buoyancy is only included as an additional body-force term in the momentum equations, should not be used. The Boussinesq approximation is only appropriate for modelling small density differences, equivalent to a temperature difference less than around 15°C in air [80].
Application Uncertainties
Recommendations for further work
Footnotes
- ↑ S. Tieszen, Private Communication, March 2010.
Contributed by: Simon Gant — UK Health & Safety Laboratory
© copyright ERCOFTAC 2010