Search

The computational cost of large eddy simulation (LES) increases rapidly with the Reynolds number when applied to attached boundary layers. This problem can be avoided by use of a Reynolds-averaged Navier-Stokes (RANS) model in the inner part of the boundary layer, which reduces the computational cost drastically. Such hybrid LES/RANS methods yield accurate results in general, but suffer from an artificial buffer layer and a shift in the velocity profile around the modeling interface. This velocity shift can be removed by use of additional forcing, but the results are very sensitive to the forcing amplitude. The present paper proposes a feedback algorithm which efficiently finds the appropriate amplitude and thus yields accurate flow statistics. The feedback algorithm is relatively robust, both in that it is insensitive to the values of the parameters involved and that it yields accurate results with different forcing fields and for different Reynolds numbers. It is argued that the feedback algorithm is consistent with the underlying assumptions of hybrid LES/RANS and that it does not introduce additional empiricism into the method. [ABSTRACT FROM AUTHOR]

The release of chemical, biological, radiological, or nuclear (CBRN) agents by terrorists or rogue states in a North American city (densely populated urban centre) and the subsequent exposure, deposition and contamination are emerging threats in an uncertain world. The modeling of the transport, dispersion, deposition and fate of a CBRN agent released in an urban environment is an extremely complex problem that encompasses potentially multiple space and time scales. The availability of high-fidelity, time-dependent models for the prediction of a CBRN agent's movement and fate in a complex urban environment can provide the strongest technical and scientific foundation for support of Canada's more broadly based effort at advancing counter-terrorism planning and operational capabilities.The objective of this paper is to report the progress of developing and validating an integrated, state-of-the-art, high-fidelity multi-scale, multi-physics modeling system for the accurate and efficient prediction of urban flow and dispersion of CBRN (and other toxic) materials discharged into these flows. Development of this proposed multi-scale modeling system will provide the real-time modeling and simulation tool required to predict injuries, casualties and contamination and to make relevant decisions (based on the strongest technical and scientific foundations) in order to minimize the consequences of a CBRN incident in a populated centre. [ABSTRACT FROM AUTHOR]

A numerical simulation has been performed of the disturbed flow through and over a two-dimensional array of rectangular buildings immersed in a neutrally stratified deep rough-walled turbulent boundary-layer flow. The model used for the simulation was the steady-state Reynolds-averaged Navier–Stokes equations with linear and non-linear eddy viscosity formulations for the Reynolds stresses. The eddy viscosity was determined using a high-Reynolds number form of the k–ϵ turbulence-closure model with the boundary conditions at the wall obtained with a standard wall-function approach. The resulting system of partial differential equations was solved using the SIMPLE algorithm in conjunction with a non-orthogonal, colocated, cell-centered, finite volume procedure. The predictive capabilities of the high-resolution computational fluid dynamics (CFD) simulations of urban flow are validated against a very detailed and comprehensive wind tunnel data set. Vertical profiles of the mean streamwise velocity and the turbulence kinetic energy are presented and compared to those measured in the wind tunnel simulation.It is found that the performance of all the turbulence models investigated is generally good—most of the qualitative features in the disturbed turbulent flow field through and over the building array are correctly reproduced. The quantitative agreement is also fairly good (especially for the mean velocity field). Overall, the non-linear k–ϵ model gave the best performance among four different turbulence closure models examined. The turbulence energy levels within the street canyons and in the exit region downstream of the last building were underestimated by all four turbulence closure models. This appears to contradict the ‘stagnation point anomaly’ associated with the standard k–ϵ model which is a result of the excessive turbulence energy production due to normal straining. A possible explanation for this is the inability of the present models to account properly for the effects of secondary strains on the turbulence and/or for the effects of large-scale flapping of the strong shear layer at the canopy top.The results of the high-resolution CFD simulations have been used to diagnose values of the drag coefficient to be used in a distributed drag force representation of the obstacles in the array. Comparisons of the measured spatially-averaged time-mean mean velocity and turbulence kinetic energy in the array with predictions of the disturbed flow using the distributed drag force approach have been made. [Copyright &y& Elsevier]

Abstract: The present study is focused on large eddy simulations (LES) that use a statistical (RANS) turbulence model near solid walls, and on the artificial buffer layer that is formed at the interface between these two modeling regions. Additional forcing is used to trigger resolved motions in the LES region more quickly, and leads to improved results in several ways. The study investigates the artificial buffer layer and how it changes with the use of forcing in an in-depth manner, with the purpose of increased understanding of the increasingly popular hybrid LES/RANS group of methods. The artificial buffer layer is shown to extend from below the modeling interface to well above it, in fact up to 20% of the boundary layer thickness for the cases studied here. The artificial buffer layer is found to be similar to the true buffer layer in many aspects, including a high correlation between the streamwise and wall normal velocity components in the ‘superstreaks’. This indicates that while the superstreaks are highly anisotropic and have unphysical length scales, they still contribute to the resolved shear stress. The forcing does not remove the artificial buffer layer, but it does reduce its extent and increases the resolved shear stress. This increase is mainly associated with increased fluctuations of the wall normal velocity. A simple, low-dimensional forcing model is proposed and tested, with favorable results. The model is simple to implement and easily generalized to more complex geometries. [Copyright &y& Elsevier]