Environmental Fluid Mechanics at UC Berkeley



Home

 

Research

 

Teaching

 

Publications

 

Students

 

Awards

 

Department website

 

Fotini Katopodes Chow

 

Research interests

My current research interests are in performing large-eddy simulations (LES) of the atmospheric boundary layer, with a focus on the development and testing of new turbulence models and improved boundary conditions for flow over complex terrain (such as mountainous and urban areas).

Current projects include studies of boundary layer processes in steep valleys, turbulence in urban boundary layer flow, source inversion for contaminant dispersion in urban areas, and development of coupled models for improving the representation of land-atmosphere interactions, among others. The Environmental Fluid Mechanics at Berkeley page describes further research in our group.

See the environmental engineering alumni newsletter or the College of Engineering Lab Notes for a non-technical overview of some aspects of my work.

 

Selected projects (see this page too)

Flow over complex terrain - Riviera Valley

Results from simulations of flow in the Riviera Valley at 150 m horizontal resolution. Contours of along-valley flow in meters per second (red: up-valley, into the page, blue: down-valley, out of the page). Vectors show flow in the cross-valley direction (not all vectors are shown).

Increases in available computational power now allow high-resolution simulations of flow over complex terrain, but the appropriate numerical and physical parameters required by such simulations are not generally known. The influence of parameterizations such as those used for turbulence, soil moisture, solar radiation, surface roughness, the configuration of initial conditions, lateral boundary conditions, and the choice of numerical grids is highly situation dependent. This project investigates the steps necessary to achieve accurate large- eddy simulations of flow in highly complex terrain. Specifically, we examine the flow and temperature fields in the Riviera Valley, located in the Alps in southern Switzerland. The simulation results are verified through comparisons to surface and radiosonde observations in the Riviera Valley, obtained during the Mesoscale Alpine Programme (MAP) Riviera Project. We evaluate the model sensitivity to changes in parameterizations such as those listed above. The boundary layer processes in our simulated valley are also studied, including comparisons to aircraft flight data, descriptions of along-valley wind transitions and secondary cross-valley circulations, and a heat budget analysis. See Chow et al. 2005 and Weigel et al. 2005 for details.

 

Source inversion for urban contaminant dispersion

The ability to determine the source of a contaminant plume in urban environments is crucial for emergency response applications. Locating the source based on downwind concentration measurements, however, is complicated by the presence of buildings which can divert flow in unexpected directions. High-resolution flow simulations are now possible for predicting plume evolution in complex urban geometries, where contaminant dispersion is affected by the flow around individual buildings. Using stochastic sampling algorithms and Bayesian inference together with a high-resolution CFD model, we have developed methods to reconstruct an atmospheric release event to determine the plume source and release rate in urban environments based on point measurements of concentration. Event reconstruction algorithms are applied first for flow around a prototype isolated building (a cube), and then using observations and flow conditions from Oklahoma City during the Joint URBAN 2003 field campaign. Stochastic sampling methods (Markov Chain Monte Carlo) are used to extract likely source term parameters, taking into consideration measurement and forward model errors. See Chow et al. 2006 for details.

 

Explicit filtering and turbulence modeling for LES

Vertical profiles of normalized mean horizontal wind speed (<U>) for neutral boundary layer simulations using an eddy-viscosity closure (Smagorinsky) and the dynamic reconstruction model (DRM) with a near-wall stress model. The traditional (Smagorinsky) closure model overestimates the flow by 10% at heights z/H ~ 0.1.

The equations for large-eddy simulation (LES) are obtained by applying a low-pass filter to the Navier–Stokes equations. This filtering operator divides the flow into so-called resolved and subfilter-scale (SFS) motions. When the equations are solved on a discrete grid, a discretization operator is applied to the equations as well, which further divides the turbulent flow field; the subfilter scales are divided into resolved SFS and unresolved SFS regions. The RSFS contribution can be theoretically reconstructed (e.g. using series expansions), and the SGS stress must be modeled. In this investigation, we study the influence of numerical errors on LES of turbulent channel flow, as well as the influence of the filtering approach and the reconstruction level on the turbulence models. The results demonstrate that improvements can be obtained for a given resolution and code by using explicit filtering and the dynamic reconstruction model (DRM). A similar approach is taken for LES of atmospheric boundary layer flow, again with significant improvement in comparisons with similarity theory (log velocity profile). In the atmospheric boundary layer, surface roughness becomes important, and its contribution must be included in the total stress distribution. These ideas are also extended for simulations of flow over an isolated hill in Scotland. See Gullbrand and Chow 2003, Chow et al. 2005, and Chow and Street 2004 for details.