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

2014 ◽  
Vol 136 (7) ◽  
Author(s):  
Abhro Pal ◽  
Kameswararao Anupindi ◽  
Yann Delorme ◽  
Niranjan Ghaisas ◽  
Dinesh A. Shetty ◽  
...  

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.

2005 ◽  
Vol 73 (1) ◽  
pp. 79-87
Author(s):  
Shari J. Kimmel-Klotzkin ◽  
Fadi P. Deek

Numerical simulations of turbulent convection under the influence of rotation will help understand mixing in oceanic flows. Though direct numerical simulations (DNS) can accurately model rotating convective flows, this method is limited to small scale and low speed flows. A large eddy simulation (LES) with the Smagorinsky subgrid scale model is used to compute the time evolution of a rotating convection flow generated by a buoyancy source of finite size at a relatively high Rayleigh number. Large eddy simulations with eddy viscosity models have been used successfully for other rotating convective flows, so the Smagorinsky model is a reasonable starting point. These results demonstrate that a LES can be used to model larger scale rotating flows, and the resulting flow structure is in good agreement with DNS and experimental results. These results also demonstrate that the qualitative behavior of vorticies which form under the source depend on the geometry of the flow. For source diameters that are small compared to the size of the domain, the vortices propagate away from the source. On the other hand, if the ratio of source diameter to domain size is relatively large, the vortices are constrained beneath the source. Though the results are qualitatively similar to a direct numerical simulation (DNS) and other LES, in this simulation the flow remains laminar much longer than the DNS predicts. This particular flow is complicated by the turbulence transition between the convective plume and the quiescent ambient fluid, and an eddy viscosity model is inadequate to accurately model this type of flow. In addition, the Smagorinsky model is not consistent in a noninertial reference frame. Thus the Smagorinsky model is not the optimal choice for this type of flow. In particular, the estimation model has demonstrated better results for other types of rotating flows and is the recommended subgrid scale model for future work.


Author(s):  
Gorazd Medic ◽  
Om Sharma

Flow over three low-pressure turbine airfoils presented in [1] is analyzed for a range of Reynolds numbers (30,000 to 150,000) by means of large-eddy simulation. Baseline computational grid for these 2D linear cascade configurations consisted of 35 millions cells, and additional finer grids of 70 millions cells were used for grid sensitivity studies. For these low Reynolds number flows, this represents a quasi-DNS resolution which minimizes the role of the subgrid-scale model — however, WALE subgrid-scale model [7] was still employed. The configurations were analyzed for low free-stream turbulence intensity, as well as for 4% turbulence intensity at free-stream. Laminar separation exists on the suction side, and, depending on the Reynolds number, the flow at the outer edge of the separation either transitions, and the separation closes before the trailing edge, or not. Detailed comparisons to measurements are presented for computed surface pressure and total pressure losses over the range of Reynolds numbers for all three airfoils; these show that LES analyses are able to capture the main trends across all three geometries.


Sign in / Sign up

Export Citation Format

Share Document