Full Content is available to subscribers

Subscribe/Learn More  >

Parallel Dynamic Large Eddy Simulation of Turbulent Flow Around MIRA Model

[+] Author Affiliations
H. G. Choi, S. W. Kang, J. Y. Yoo

Seoul National University, Seoul, Korea

Paper No. FEDSM2002-31100, pp. 1007-1012; 6 pages
  • ASME 2002 Joint U.S.-European Fluids Engineering Division Conference
  • Volume 1: Fora, Parts A and B
  • Montreal, Quebec, Canada, July 14–18, 2002
  • Conference Sponsors: Fluids Engineering Division
  • ISBN: 0-7918-3615-0 | eISBN: 0-7918-3600-2
  • Copyright © 2002 by ASME


For the large scale computation of turbulent flows around an arbitrarily shaped body, a parallel LES (large eddy simulation) code has been recently developed in which domain decomposition method is adopted. METIS and MPI (message passing interface) libraries are used for domain partitioning and data communication between processors, respectively. For unsteady computation of the incompressible Navier-Stokes equation, 4-step splitting finite element algorithm [1] is adopted and Smagorinsky or dynamic LES model can be chosen for the modeling of small eddies in turbulent flows. For the outlet (open) boundary condition, a Dirichlet boundary condition for the pressure is proposed. For the validation and performance-estimation of the parallel code, a three-dimensional laminar flow generated by natural convection inside a cube has been solved. We have confirmed that our code gives accurate results compared with previous studies. Regarding the speed-up of the code, the present parallel code with parallel block-Jacobi preconditioner is about 50 times faster than the corresponding serial code with 64 processors when approximately one million grid points are used. Most of the CPU time is consumed in solving elliptic type pressure equation. For the validation of LES models, turbulent channel flows are simulated at Re = 180, which is based on the channel half height and friction velocity using 51 × 71 × 71 grid system. It has been shown that our results agree well with the well-known results by Kim et al. [2] with less grid points than used by them in terms of time-averaged velocity field and velocity fluctuation. Lastly, we have solved the turbulent flow around MIRA (Motor Industry Research Association) model at Re = 1.6 × 106 which is based on the model height and inlet free stream velocity. Both Smagorinsky and dynamic models are tested, comparing estimated drag coefficients and pressure distribution along the model surface with the existing experimental data [3]. With the help of the parallel code developed in this study, we are able to obtain a unsteady solution of the turbulent flow field around a vehicle discretized by approximately three million grid points within two weeks when 32 IBM-SP2-processors are used. The calculated drag coefficient agrees better with the experimental result [3] than those using two equation turbulence models [4].

Copyright © 2002 by ASME



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In