LES Simulation of Backflow Vortex Structure at the Inlet of an Inducer

[+] Author and Article Information
Nobuhiro Yamanishi

 Japan Aerospace Exploration Agency, Ibaraki 305 8505, Japanyamanishi.nobuhiro@jaxa.jp

Shinji Fukao

 Mitsubishi Heavy Industries, Hyogo 676 8686, Japanshinji̱fukao@mhi.co.jp

Xiangyu Qiao

 ANSYS, Kanagawa 220 6216, Japanxiangyu.qiao@ansys.com

Chisachi Kato

 The University of Tokyo, Tokyo 153 8505, Japanckato@iis.u-tokyo.ac.jp

Yoshinobu Tsujimoto

 Osaka University, Osaka 560 8531, Japantujimoto@me.es.osaka-u.ac.jp

J. Fluids Eng 129(5), 587-594 (Sep 28, 2006) (8 pages) doi:10.1115/1.2717613 History: Received August 25, 2005; Revised September 28, 2006

Turbopump inducers often have swirling backflow under a wide range of flow rates because they are designed with a certain angle of attack even at the design point in order to attain high cavitation performance. When the flow rate is decreased, the backflow region extends upstream and may cause various problems by interacting with upstream elements. It is also known that the backflow vortex structure occurs in the shear layer between the main flow and the swirling backflow. Experimental studies on the backflow from an inducer have given us insight into the characteristics of backflow vortex structure, but the limited information has not lead to the complete understanding of the phenomena. Numerical studies based on Reynolds-averaged Navier-Stokes (RANS) computations usually deteriorate when the flow field of interest involves large-scale separations, as shown by a previous study by Tsujimoto (2005). On the other hand, the numerical approach using the Large Eddy Simulation (LES) technique has the potential to predict unsteady flows and/or flow fields that include regions of large-scale separation much more accurately than RANS computations does in general. The present paper describes the application of the LES code developed by one of the authors (Kato) to further understand the backflow vortex structure at the inlet of an inducer. First, the internal flow of the inducer was simulated, as a way to evaluate the validity of the proposed method, under a wide range of inlet flow coefficients. The static pressure peformance and the length of the backflow region was compared with measured values, and good agreement was obtained. Second, using the validated LES code, the fundamental characteristics of the backflow vortex was investigated in detail. It was found that the backflow vortices are formed in a circumferentially twisted manner at the boundary between the swirling backflow and the straight inlet flow. Also, the backflow vortices rotate in the same direction as the inducer, but with half of the circumferential flow velocity in the backflow region. Another finding was that the backflow region expands toward the center of the flow field and the number of vortices decrease, as the flow coefficient decreases. To the best of our knowledge, this is the first computation of the backflow at the inducer inlet to achieve quantitative agreement with measured results, and give new findings to the complicated three-dimensional structure of the backflow, which was very limited under experimental studies.

Copyright © 2007 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

Schematic of backflow vortex

Grahic Jump Location
Figure 2

Inducer geometry

Grahic Jump Location
Figure 3

Computational region

Grahic Jump Location
Figure 4

Computational mesh

Grahic Jump Location
Figure 5

Static pressure performance

Grahic Jump Location
Figure 6

Axial location of upstream edge of backflow

Grahic Jump Location
Figure 7

Pressure contours and velocity vector in the r-θ cross section at z∕Dt=−0.50 (ϕ=0.06, after 55 revolutions)

Grahic Jump Location
Figure 8

Backflow vortex at the shear layer between main flow and swirling backflow (ϕ=0.06, after 55 revolutions)

Grahic Jump Location
Figure 9

Circumferential velocity and pressure distribution at various axial locations (ϕ=0.06, after 55 revolutions)

Grahic Jump Location
Figure 10

Axial velocity distribution at r∕Dt=0.50 (tip) (ϕ=0.06, after 55 revolutions)

Grahic Jump Location
Figure 11

Pressure distribution at r∕Dt=0.50 (tip) (ϕ=0.06, after 55 revolutions)

Grahic Jump Location
Figure 12

Rotation process of backflow vortex (ϕ=0.06, from 54 to 55 revolutions)

Grahic Jump Location
Figure 13

Backflow vortex (ϕ=0.05, after 55 revolutions)

Grahic Jump Location
Figure 14

Backflow vortex (ϕ=0.07, after 45 revolutions)

Grahic Jump Location
Figure 15

Normalized propagation velocity and maximum circumferential velocity (z∕Dt=−0.18, except for ϕ=0.078: z∕Dt=−0.10)

Grahic Jump Location
Figure 16

Number of vortices

Grahic Jump Location
Figure 17

Profile of vortex filaments projected to the meridional plane




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