0
Research Papers

Large Eddy Simulation of Transitional Flow in an Idealized Stenotic Blood Vessel: Evaluation of Subgrid Scale Models

[+] Author and Article Information
Abhro Pal, Yann Delorme, Niranjan Ghaisas, Dinesh A. Shetty, Steven H. Frankel

School of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907

Kameswararao Anupindi

School of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907
e-mail: kamesh.a@gmail.com

1Corresponding author.

Manuscript received January 28, 2013; final manuscript received April 4, 2014; accepted manuscript posted May 7, 2014; published online May 22, 2014. Editor: Victor H. Barocas.

J Biomech Eng 136(7), 071009 (May 22, 2014) (8 pages) Paper No: BIO-13-1048; doi: 10.1115/1.4027610 History: Received January 28, 2013; Revised April 04, 2014; Accepted May 07, 2014

In the present study, we performed large eddy simulation (LES) of axisymmetric, and 75% stenosed, eccentric arterial models with steady inflow conditions at a Reynolds number of 1000. The results obtained are compared with the direct numerical simulation (DNS) data (Varghese et al., 2007, “Direct Numerical Simulation of Stenotic Flows. Part 1. Steady Flow,” J. Fluid Mech., 582, pp. 253–280). An inhouse code (WenoHemo) employing high-order numerical methods for spatial and temporal terms, along with a 2nd order accurate ghost point immersed boundary method (IBM) (Mark, and Vanwachem, 2008, “Derivation and Validation of a Novel Implicit Second-Order Accurate Immersed Boundary Method,” J. Comput. Phys., 227(13), pp. 6660–6680) for enforcing boundary conditions on curved geometries is used for simulations. Three subgrid scale (SGS) models, namely, the classical Smagorinsky model (Smagorinsky, 1963, “General Circulation Experiments With the Primitive Equations,” Mon. Weather Rev., 91(10), pp. 99–164), recently developed Vreman model (Vreman, 2004, “An Eddy-Viscosity Subgrid-Scale Model for Turbulent Shear Flow: Algebraic Theory and Applications,” Phys. Fluids, 16(10), pp. 3670–3681), and the Sigma model (Nicoud et al., 2011, “Using Singular Values to Build a Subgrid-Scale Model for Large Eddy Simulations,” Phys. Fluids, 23(8), 085106) are evaluated in the present study. Evaluation of SGS models suggests that the classical constant coefficient Smagorinsky model gives best agreement with the DNS data, whereas the Vreman and Sigma models predict an early transition to turbulence in the poststenotic region. Supplementary simulations are performed using Open source field operation and manipulation (OpenFOAM) (“OpenFOAM,” http://www.openfoam.org/) solver and the results are inline with those obtained with WenoHemo.

FIGURES IN THIS ARTICLE
<>
Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

Grid 3 used for WenoHemo simulations

Grahic Jump Location
Fig. 2

Comparison of mean axial velocity profiles for various grid sizes considered. Dotted line represents grid 1, Solid line represents grid 2, dashed-dotted line represents grid 3. DNS data is represented by solid dots. The scale represents 5 velocity units (u/Um) for each axial unit (x/D).

Grahic Jump Location
Fig. 3

(a) Normalized mean axial velocity profiles (u¯/Um) at indicated locations along the axial direction. Solid lines correspond to WenoHemo simulation, solid dots indicate DNS [9]. Scale represents 1 axial unit (x/D) equals 5 velocity units (u¯/Um). (b) Contours of normalized vorticity magnitude (|ω→|D/Um) on XY-plane at Z = 0.

Grahic Jump Location
Fig. 4

Contours of normalized vorticity magnitude (|ω→|D/Um) on the XZ-plane, at Y = 0 for DNS and LES simulations as indicated

Grahic Jump Location
Fig. 5

Normalized mean axial velocity (u¯/Um) profiles at indicated locations. Lines show present simulations using WenoHemo and dots indicate the DNS result.

Grahic Jump Location
Fig. 9

Instantaneous coherent structures, identified by the λ2 criterion defined by Ref. [30], colored by instantaneous normalized vorticity magnitude |ω→|D/Um. The inset shows a close-up view of the structures.

Grahic Jump Location
Fig. 10

Turbulent energy spectra E11 for the indicated SGS models and S represents nondimensional frequency

Grahic Jump Location
Fig. 11

Mesh of approximately 1.0 × 106 cells used in OpenFOAM simulations

Grahic Jump Location
Fig. 12

Normalized mean axial velocity (u¯/Um) profiles at the indicated locations for the SGS models considered. Lines denote result from present simulation using OpenFOAM and dots indicate the DNS result.

Grahic Jump Location
Fig. 13

TKE profiles at the indicated locations for the SGS models considered. Lines indicate the present simulation result obtained using OpenFOAM and dots indicate the DNS result.

Grahic Jump Location
Fig. 6

RMS of velocity fluctuations at indicated axial locations. Solid lines indicate LES with Smagorinsky Model using WenoHemo, solid dots indicate DNS [9].

Grahic Jump Location
Fig. 7

Turbulent kinetic energy profiles at indicated axial locations. Solid line indicates LES using WenoHemo, solid dots indicate DNS [9].

Grahic Jump Location
Fig. 8

Contour plots of time averaged SGS activity 〈Asgs〉 (top), and vorticity magnitude |ω→|D/Um (bottom) predicted by WenoHemo for eccentric model

Tables

Errata

Discussions

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 Journal Articles
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