Quantifying turbulence model uncertainty in Reynolds-averaged Navier–Stokes simulations of a pin-fin array. Part 2: Scalar transport. (15th September 2020)

• Record Type: Journal Article
• Authors: Hao, Zengrong; Gorlé, Catherine
• Is Part Of: Computers & fluids. Volume 209(2020)
• Publication Date: 2020-09-15
• Subjects: Turbulence model form uncertainties -- Turbulent scalar flux model -- Epistemic uncertainty quantification (EUQ) -- Reynolds stress realizability; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2020.104642
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: Uncertainty in LEV assumption for Reynolds stress is propagated to scalar transport. Scalar flux models involving turbulence anisotropy exhibit good bounding behaviors. Bounding behaviors depend on both momentum transport and relevant stress components. Abstract: The accuracy of Reynolds-averaged Navier-Stokes results for turbulent scalar transport are affected by epistemic uncertainty in the Reynolds stress model in two ways: by altering the mean velocity field that advects the scalar, and by altering the inputs required for scalar flux models. We investigate these effects by propagating uncertainty in the Reynolds stress model to the prediction of scalar quantities in simulations of a pin-fin heat exchanger. The Reynolds stress model uncertainty is quantified by perturbing the anisotropy of the predicted stress tensor towards the three limiting realizable states. This uncertainty is then propagated to the scalar turbulence transport via the conservation law for the mean scalar and the turbulent scalar flux model. We consider three different scalar flux models that explicitly take the Reynolds stresses as input, and evaluate the performance of the models by verifying if high-fidelity large eddy simulation results are encompassed by the model predictions. The results indicate that the predicted uncertainty depends on both the general level of momentum transport and the most-relevant stress component, which is affected by the anisotropy perturbations. The predictions provide a similar bounding behavior for the overall temperature field as for the momentum field, but fail in bounding the local heat transfer rates in several locations. The bounding behaviors are further analyzed in terms of the predicted uncertainty in the flux magnitude, direction, and divergence to identify opportunities for further improvement of the proposed methods.