Large-eddy simulation of turbulent flow using the finite element method

The equations of motion describing turbulent flows (in both the low and high Reynolds-number regimes) are well established. However, present day computers cannot meet the enormous computational requirement for numerically solving the governing equations for common engineering flows in the high Reynolds number turbulent regime. The characteristics that make turbulent, high Reynolds number flows difficult to simulate is the extreme range of time and space scales of motion. Most current engineering calculations are performed using semi-empirical equations, developed in terms of the flow mean (average) properties. These turbulence{open_quote} models{close_quote} (semi-empirical/analytical approximations) do not explicitly account for the eddy structures and thus, the temporal and spatial flow fluctuations are not resolved. In these averaging approaches, it is necessary to approximate all the turbulent structures using semi-empirical relations, and as a result, the turbulence models must be tailored for specific flow conditions and geometries with parameters obtained (usually) from physical experiments. The motivation for this research is the development of a finite element turbulence modeling approach which will ultimately be used to predict the wind flow around buildings. Accurate turbulence models of building flow are needed to predict the dispersion of airborne pollutants. The building flow turbulence models used today are not capable of predicting the three-dimensional separating and reattaching flows without the manipulation of many empirical parameters. These empirical parameters must be set by experimental data and they may vary unpredictably with building geometry, building orientation, and upstream flow conditions