UFR 1-06 Best Practice Advice: Difference between revisions
m (Dave.Ellacott moved page SilverP:UFR 1-06 Best Practice Advice to UFR 1-06 Best Practice Advice) |
|||
(9 intermediate revisions by the same user not shown) | |||
Line 5: | Line 5: | ||
}} | }} | ||
=Axisymmetric buoyant far-field plume | = Axisymmetric buoyant far-field plume = | ||
Underlying Flow Regime 1-06 | |||
= Best Practice Advice = | = Best Practice Advice = | ||
Line 54: | Line 52: | ||
Hypothesis (GGDH) instead of the Standard Gradient Diffusion Hypothesis | Hypothesis (GGDH) instead of the Standard Gradient Diffusion Hypothesis | ||
(SGDH). More refined models have been suggested which could further | (SGDH). More refined models have been suggested which could further | ||
improve model predictions, see for example Hossain & Rodi.</li> | improve model predictions, see for example Hossain & Rodi [[UFR_1-06_References|[8]]].</li> | ||
<li>Do not neglect the buoyancy source term in the | <li>Do not neglect the buoyancy source term in the | ||
''&epsilon''-equation as this can lead to problems with | ''ε''-equation as this can lead to problems with | ||
numerical stability.</li> | numerical stability.</li> | ||
Line 75: | Line 73: | ||
jet. However, the standard ''k – ε'' model predicts | jet. However, the standard ''k – ε'' model predicts | ||
the spreading rate for round jets to be 15% ''higher '' than for | the spreading rate for round jets to be 15% ''higher '' than for | ||
the plane jets.</li> | the plane jets [[UFR_1-06_References|[67]]].</li> | ||
</ul> | </ul> | ||
=== Application Uncertainties === | === Application Uncertainties === | ||
In Van Maele & Merci’s calculations, they substituted | |||
the streamwise normal stress in the GGDH expression for the production | |||
due to buoyancy, ''G,'' with the turbulent kinetic energy. They | |||
justified this ad-hoc modification on the basis that that the | |||
streamwise normal stress was underpredicted using a linear ''k – ε'' | |||
model and therefore using ''k'' was more | |||
appropriate. However, this modification may not be appropriate in more | |||
complex flows where the flow direction is not aligned to one of the | |||
coordinate axes. The effect on Van Maele & Merci’s | |||
results would have been to increase the buoyancy production term | |||
''G''. One could therefore expect to see a slightly less | |||
significant difference between SGDH and GGDH approaches without this | |||
ad-hoc modification. | |||
Van Maele & Merci also did not report any tests that had been performed | |||
to assess the influence of the entrainment boundaries on the flow | |||
predictions. Ideally, simulations should have been performed using a | |||
smaller or larger domain to demonstrate that the presence of the | |||
entrainment boundaries had no effect on the solution. | |||
=== Recommendations for further work === | === Recommendations for further work === | ||
It is recommended that further experiments and a systematic | |||
re-evaluation of available data be undertaken to establish with | |||
confidence the self-similar behaviour of axisymmetric buoyant plumes, | |||
similar to the exercise undertaken 30 years ago by Chen & Rodi [[UFR_1-06_References|[1]]] | |||
Both List [[UFR_1-06_References|[6]]] | |||
and Dai ''et al.'' [[UFR_1-06_References|[10]]] | |||
have called into question whether the measurements by | |||
George ''et al.'' [[UFR_1-06_References|[3]]] were carried out sufficiently far | |||
from the source, in the region where self-similar behaviour exists. | |||
This issue was, however, dismissed by | |||
Shabbir & George [[UFR_1-06_References|[11]]][[UFR_1-06_References|[34]]]. | |||
Line 91: | Line 123: | ||
{{ACContribs | {{ACContribs | ||
|authors=Simon Gant | |authors=Simon Gant | ||
|organisation= | |organisation=UK Health & Safety Laboratory | ||
}} | }} | ||
© copyright ERCOFTAC 2010 | © copyright ERCOFTAC 2010 |
Latest revision as of 19:15, 11 February 2017
Axisymmetric buoyant far-field plume
Underlying Flow Regime 1-06
Best Practice Advice
Best Practice Advice for the UFR
Key Physics
The key physics to be captured in this UFR is the self-similar behaviour of a spreading axisymmetric buoyant plume.
Numerical Modelling Issues
- The flow can be treated as axisymmetric.
- For a grid-independent resolution of the George et al. [3] plume, at least (40 × 100) grid nodes should be used in the (radial × axial) directions. At least 10 nodes should be used radially to resolve the plume source.
- Discretization schemes should be at least second-order accurate.
- For further advice on boundary conditions, see Test Case
Physical Modelling
- Use of the standard k – ε model with or without the common Standard Gradient Diffusion Hypothesis (SGDH) for the production term due to buoyancy, G, will probably result in overprediction of the centreline mean parameters, and underprediction of the spreading rate of the plume.
- To obtain more accurate plume predictions with a k – ε model, use the Generalized Gradient Diffusion Hypothesis (GGDH) instead of the Standard Gradient Diffusion Hypothesis (SGDH). More refined models have been suggested which could further improve model predictions, see for example Hossain & Rodi [8].
- Do not neglect the buoyancy source term in the ε-equation as this can lead to problems with numerical stability.
- If you are considering only the far-field region of a buoyant plume, where density differences are small, the Boussinesq approximation can be used. If, however, your flow domain includes the region nearer the source of buoyancy where density differences are appreciable, avoid using the Boussinesq approximation.
- For cases where buoyancy is not as strong as in a plume, in the limit of a non-buoyant axisymmetric jet, be aware of the limitations of the standard k – ε model. The spreading rate of a non-buoyant round jet is 15% lower than for a two-dimensional, plane jet. However, the standard k – ε model predicts the spreading rate for round jets to be 15% higher than for the plane jets [67].
Application Uncertainties
In Van Maele & Merci’s calculations, they substituted the streamwise normal stress in the GGDH expression for the production due to buoyancy, G, with the turbulent kinetic energy. They justified this ad-hoc modification on the basis that that the streamwise normal stress was underpredicted using a linear k – ε model and therefore using k was more appropriate. However, this modification may not be appropriate in more complex flows where the flow direction is not aligned to one of the coordinate axes. The effect on Van Maele & Merci’s results would have been to increase the buoyancy production term G. One could therefore expect to see a slightly less significant difference between SGDH and GGDH approaches without this ad-hoc modification.
Van Maele & Merci also did not report any tests that had been performed
to assess the influence of the entrainment boundaries on the flow
predictions. Ideally, simulations should have been performed using a
smaller or larger domain to demonstrate that the presence of the
entrainment boundaries had no effect on the solution.
Recommendations for further work
It is recommended that further experiments and a systematic re-evaluation of available data be undertaken to establish with confidence the self-similar behaviour of axisymmetric buoyant plumes, similar to the exercise undertaken 30 years ago by Chen & Rodi [1] Both List [6] and Dai et al. [10] have called into question whether the measurements by George et al. [3] were carried out sufficiently far from the source, in the region where self-similar behaviour exists. This issue was, however, dismissed by Shabbir & George [11][34].
Contributed by: Simon Gant — UK Health & Safety Laboratory
© copyright ERCOFTAC 2010