A stochastic SGS model with application to turbulent channel flow with a passive scalar

Author(s):  
Linus Marstorp ◽  
Geert Brethouwer ◽  
Arne V. Johansson
2013 ◽  
Vol 721 ◽  
pp. 541-577 ◽  
Author(s):  
Amin Rasam ◽  
Geert Brethouwer ◽  
Arne V. Johansson

AbstractIn Marstorpet al. (J. Fluid Mech., vol. 639, 2009, pp. 403–432), an explicit algebraic subgrid stress model (EASSM) for large-eddy simulation (LES) was proposed, which was shown to considerably improve LES predictions of rotating and non-rotating turbulent channel flow. In this paper, we extend that work and present a new explicit algebraic subgrid scalar flux model (EASSFM) for LES, based on the modelled transport equation of the subgrid-scale (SGS) scalar flux. The new model is derived using the same kind of methodology that leads to the explicit algebraic scalar flux model of Wikströmet al. (Phys. Fluids, vol. 12, 2000, pp. 688–702). The algebraic form is based on a weak equilibrium assumption and leads to a model that depends on the resolved strain-rate and rotation-rate tensors, the resolved scalar-gradient vector and, importantly, the SGS stress tensor. An accurate prediction of the SGS scalar flux is consequently strongly dependent on an accurate description of the SGS stresses. The new EASSFM is therefore primarily used in connection with the EASSM, since this model can accurately predict SGS stresses. The resulting SGS scalar flux is not necessarily aligned with the resolved scalar gradient, and the inherent dependence on the resolved rotation-rate tensor makes the model suitable for LES of rotating flow applications. The new EASSFM (together with the EASSM) is validated for the case of passive scalar transport in a fully developed turbulent channel flow with and without system rotation. LES results with the new model show good agreement with direct numerical simulation data for both cases. The new model predictions are also compared to those of the dynamic eddy diffusivity model (DEDM) and improvements are observed in the prediction of the resolved and SGS scalar quantities. In the non-rotating case, the model performance is studied at all relevant resolutions, showing that its predictions of the Nusselt number are much less dependent on the grid resolution and are more accurate. In channel flow with wall-normal rotation, where all the SGS stresses and fluxes are non-zero, the new model shows significant improvements over the DEDM predictions of the resolved and SGS quantities.


Author(s):  
Kyoungyoun Kim ◽  
Radhakrishna Sureshkumar

A direct numerical simulation (DNS) of viscoelastic turbulent channel flow with the FENE-P model was carried out to investigate turbulent heat transfer mechanism of polymer drag-reduced flows. The configuration was a fully-developed turbulent channel flow with uniform heat flux imposed on both walls. The temperature was considered as a passive scalar. The Reynolds number based on the friction velocity (uτ) and channel half height (δ) is 125 and Prandtl number is 5. Consistently with the previous experimental observations, the present DNS results show that the heat-transfer coefficient was reduced at a rate faster than the accompanying drag reduction rate. Statistical quantities such as root-mean-square temperature fluctuations and turbulent heat fluxes were obtained and compared with those of a Newtonian fluid flow. Budget terms of the turbulent heat fluxes were also presented.


1992 ◽  
Vol 114 (3) ◽  
pp. 598-606 ◽  
Author(s):  
N. Kasagi ◽  
Y. Tomita ◽  
A. Kuroda

A direct numerical simulation (DNS) of the fully developed thermal field in a two-dimensional turbulent channel flow of air was carried out. The isoflux condition was imposed on the two walls so that the local mean temperature increased linearly in the streamwise direction. With any buoyancy effect neglected, temperature was considered as a passive scalar. The computation was executed on 1,589,248 grid points by using a spectral method. The statistics obtained were root-mean-square temperature fluctuations, turbulent heat fluxes, turbulent Prandtl number, and dissipation time scales. They agreed fairly well with existing experimental and numerical simulation data. Each term in the budget equations of temperature variance, its dissipation rate, and turbulent heat fluxes was also calculated. It was found that the temperature fluctuation θ′ was closely correlated with the streamwise velocity fluctuation u′, particularly in the near-wall region. Hence, the distribution of budget terms for the streamwise and wall-normal heat fluxes, u′θ′ and v′θ′, were very similar to those for the two Reynolds stress components, u′u′ and u′v′.


AIChE Journal ◽  
2005 ◽  
Vol 51 (7) ◽  
pp. 1938-1950 ◽  
Author(s):  
V. K. Gupta ◽  
R. Sureshkumar ◽  
B. Khomami

Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner

Abstract In this paper a new method is introduced to model the transport of entropy waves and equivalence ratio fluctuations in turbulent flows. The model is based on the Navier-Stokes equations and includes a transport equation for a passive scalar, which may stand for entropy or equivalence ratio fluctuations. The equations are linearized around the mean turbulent fields, which serve as the input to the model in addition to a turbulent eddy viscosity, which accounts for turbulent diffusion of the perturbations. Based on these inputs, the framework is able to predict the linear response of the flow velocity and passive scalar to harmonic perturbations that are imposed at the boundaries of the computational domain. These in this study are fluctuations in the passive scalar and/or velocities at the inlet of a channel flow. The code is first validated against analytic results, showing very good agreement. Then the method is applied to predict the convection, mean flow dispersion and turbulent mixing of passive scalar fluctuations in a turbulent channel flow, which has been studied in previous work with Direct Numerical Simulations (DNS). Results show that our code reproduces the dynamics of coherent passive scalar transport in the DNS with very high accuracy and low numerical costs, when the DNS mean flow and Reynolds stresses are provided. Furthermore, we demonstrate that turbulent mixing has a significant effect on the transport of the passive scalar fluctuations. Finally, we apply the method to explain experimental observations of transport of equivalence ratio fluctuations in the mixing duct of a model burner.


Sign in / Sign up

Export Citation Format

Share Document