Comparison of Turbulence Modeling Strategies for Indoor Flows

Turbulence modeling techniques are compared for the simulation of low speed indoor air flow in a simple room. The effect of inlet turbulence intensity on the flow field is investigated using the constant coefficient large eddy simulation (LES) model with uniform mean inlet conditions at several levels of inlet turbulence intensities. The results show significant differences between the simulations with laminar inflow conditions and those in which turbulence was introduced at the inlet. For simulations with turbulent inlet conditions, it is noticed that the jet transitions to a state of fully developed turbulence wherein the dynamics of the flow become nearly insensitive to any further increase in the level of inlet turbulence. For laminar flow conditions, it is seen that the jet slowly spreads and mixes with the quiescent room air. As a result, the jet reaches a fully developed turbulent state further away from the inlet relative to the simulations with inlet turbulence. The effect of using experimental inlet profiles is also investigated. It is seen that, close to the inlet, the flow is sensitive to the inflow details, whereas further away from the inlet, these effects become less pronounced. The results from the constant coefficient and the dynamic LES models are compared. The most noticeable differences in the flow occur at the locations where the subgrid-scale's contribution to the turbulent kinetic energy is highest. Finally, the results from the dynamic LES and the k-∈ models are compared. It is found that there are significant differences between the two models for the zero inlet turbulence limit where the flow is most probably transitional in nature and turbulence has not yet reached a fully developed state. It is seen that in the laminar inflow case the k-∈ model predicts a fully turbulent jet very close to the inlet and thus fails to capture the slow development of the jet found in LES. Accordingly, the k-∈ model results are nearly insensitive to the level of inlet turbulence especially far from the origin of the flow. It is also seen that for cases with nonzero inlet turbulence level, the k-∈ model predicts the general features of the mean flow reasonably well; however, the k-∈ model overpredicts the jet spreading rate and the turbulent kinetic energy close to the inlet. Furthermore, the k-∈ model under predicts the turbulence level near the corner of the ceiling as it fails to capture the complicated mean velocity and turbulent kinetic energy, most likely because of the highly intermittent flow pattern found there in LES.

[1]  K. G. Budden.  The basic equations , 1985 .

[2]  P. Moin,et al.  Simulation of spatially evolving turbulence and the applicability of Taylor's hypothesis in compressible flow , 1992 .

[3]  Zhao Zhang,et al.  Evaluation of Various Turbulence Models in Predicting Airflow and 1 Turbulence in Enclosed Environments by CFD : Part-1 : 2 Summary of Prevalent Turbulence Models 3 4 , 2007 .

[4]  Jianbo Jiang.  Experimental and Numerical Study of Air Flows in a Full Scale Room , 2009 .

[5]  Yi Jiang,et al.  Using Large Eddy Simulation to Study Airflows in and around Buildings , 2003 .

[6]  M. Germano,et al.  Turbulence: the filtering approach , 1992, Journal of Fluid Mechanics.

[7]  I. Celik,et al.  Random Flow Generation Technique for Large Eddy Simulations and Particle-Dynamics Modeling , 2001 .

[8]  Qingyan Chen.  COMPARISON OF DIFFERENT k-ε MODELS FOR INDOOR AIR FLOW COMPUTATIONS , 1995 .

[9]  J. Ferziger,et al.  A new non-eddy viscosity subgrid-scale model and its application to channel flow , 1995 .

[10]  P. Nielsen,et al.  The Selection of Turbulence Models for Prediction of Room Airflow , 1998 .

[11]  C. Meneveau,et al.  A Lagrangian dynamic subgrid-scale model of turbulence , 1994, Journal of Fluid Mechanics.

[12]  Peter V. Nielsen,et al.  Modifications of the v2-f Model for Computing the Flow in a 3D Wall Jet , 2003 .

[13]  P. Moin,et al.  A dynamic localization model for large-eddy simulation of turbulent flows , 1995, Journal of Fluid Mechanics.

[14]  Caskey,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS I . THE BASIC EXPERIMENT , 1962 .

[15]  Rainer Friedrich,et al.  ANALYSING TURBULENT BACKWARD-FACING STEP FLOW WITH THE LOWPASS-FILTERED NAVIER-STOKES EQUATIONS , 1990 .

[16]  J. Deardorff.  Three-dimensional numerical study of the height and mean structure of a heated planetary boundary layer , 1974 .

[17]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[18]  D. Marr,et al.  Velocity measurements in the breathing zone of a moving thermal manikin within the indoor environment , 2007 .

[19]  B. Launder,et al.  Lectures in mathematical models of turbulence , 1972 .

[20]  H. Michel,et al.  Spectral dynamics in the B800 band of LH2 from Rhodospirillum molischianum: a single-molecule study , 2004 .

[21]  Qingyan Chen,et al.  Large eddy simulation of indoor airflow with a filtered dynamic subgrid scale model , 2000 .

[22]  P. Durbin.  SEPARATED FLOW COMPUTATIONS WITH THE K-E-V2 MODEL , 1995 .

[23]  S. Orszag,et al.  Renormalization group analysis of turbulence. I. Basic theory , 1986 .

[24]  Wei Chen,et al.  Large Eddy Simulation of Natural and Mixed Convection Airflow Indoors with Two Simple Filtered Dynamic Subgrid Scale Models , 2000 .

[25]  P. Moin,et al.  Direct numerical simulation of turbulent flow over a backward-facing step , 1997, Journal of Fluid Mechanics.

[26]  Zhiqiang John Zhai,et al.  Evaluation of Various Turbulence Models in Predicting Airflow and Turbulence in Enclosed Environments by CFD: Part 2—Comparison with Experimental Data from Literature , 2007 .

[27]  Peter V. Nielsen,et al.  Large Eddy Simulations of the Flow in Three-Dimensional Ventilated Room , 1996 .

[28]  U. Schumann.  Subgrid Scale Model for Finite Difference Simulations of Turbulent Flows in Plane Channels and Annuli , 1975 .

[29]  B. A. Kader.  Temperature and concentration profiles in fully turbulent boundary layers , 1981 .

[30]  M. Wolfshtein.  The velocity and temperature distribution in one-dimensional flow with turbulence augmentation and pressure gradient , 1969 .

[31]  D. Wilcox.  Turbulence modeling for CFD , 1993 .

[32]  Jean-Pierre Bertoglio,et al.  Assessment of the vortex method for Large Eddy Simulation inlet conditions , 2006 .

[33]  S. Orszag,et al.  Renormalization group analysis of turbulence. I. Basic theory , 1986, Physical review letters.

[34]  S. Pope.  Ten questions concerning the large-eddy simulation of turbulent flows , 2004 .

[35]  Qingyan Chen,et al.  A Procedure for Verification, Validation, and Reporting of Indoor Environment CFD Analyses , 2002 .

[36]  S. Pope.  Turbulent Flows: FUNDAMENTALS , 2000 .

[37]  Peter V. Nielsen,et al.  The Velocity Characteristics of Ventilated Rooms , 1978 .

[38]  B. Launder,et al.  Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc , 1974 .

[39]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[40]  Qingyan Chen,et al.  Significant questions in predicting room air motion , 1992 .

[41]  M. Kopera.  Direct numerical simulation of turbulent flow over a backward-facing step , 2011 .

[42]  Mglc Marcel Loomans,et al.  The measurement and simulation of indoor air flow , 1998 .

[43]  K. Lilly.  The representation of small-scale turbulence in numerical simulation experiments , 1966 .

[44]  John D. Spengler,et al.  Indoor Air Quality Handbook , 2000 .