The multidimensional positive definite advection transport algorithm: nonoscillatory option
Journal of Computational Physics
Modeling a no-slip flow boundary with an external force field
Journal of Computational Physics
A flexible inner-outer preconditioned GMRES algorithm
SIAM Journal on Scientific Computing
MPDATA: a finite-difference solver for geophysical flows
Journal of Computational Physics
Large-eddy simulations of convective boundary layers using nonoscillatory differencing
Physica D - Special issue originating from the 18th Annual International Conference of the Center for Nonlinear Studies, Los Alamos, NM, May 11&mdash ;15, 1998
An all-scale anelastic model for geophysical flows: dynamic grid deformation
Journal of Computational Physics
Journal of Computational Physics
MPDATA: An edge-based unstructured-grid formulation
Journal of Computational Physics
Iterated upwind schemes for gas dynamics
Journal of Computational Physics
DNS of buoyancy-dominated turbulent flows on a bluff body using the immersed boundary method
Journal of Computational Physics
On numerical realizability of thermal convection
Journal of Computational Physics
Prediction of wall-pressure fluctuation in turbulent flows with an immersed boundary method
Journal of Computational Physics
Pores resolving simulation of Darcy flows
Journal of Computational Physics
EULAG, a computational model for multiscale flows: An MHD extension
Journal of Computational Physics
An unstructured-mesh atmospheric model for nonhydrostatic dynamics
Journal of Computational Physics
Hi-index | 31.48 |
We perform large-eddy simulations (LES) of the flow past a scale model of a complex building. Calculations are accomplished using two different methods to represent the edifice. The first method employs the standard Gal-Chen and Somerville terrain-following coordinate transformation, common in mesoscale atmospheric simulations. The second method uses an immersed boundary approach, in which fictitious body forces in the equations of motion are used to represent the building by attenuating the flow to stagnation within a time comparable to the time step of the model. Both methods are implemented in the same hydrodynamical code (EULAG) using the same nonoscillatory forward-in-time (NFT) incompressible flow solver based on the multidimensional positive definite advection transport algorithms (MPDATA). The two solution methods are compared to wind tunnel data collected for neutral stratification. Profiles of the first- and second-order moments at various locations around the model building show good agreement with the wind tunnel data. Although both methods appear to be viable tools for LES of urban flows, the immersed boundary approach is computationally more efficient. The results of these simulations demonstrate that, contrary to popular opinion, continuous mappings such as the Gal-Chen and Somerville transformation are not inherently limited to gentle slopes. Calculations for a strongly stratified case are also presented to point out the substantial differences from the neutral boundary layer flows.