0
Research Papers: Techniques and Procedures

A Numerical Simulation Algorithm of the Inviscid Dynamics of Axisymmetric Swirling Flows in a Pipe

[+] Author and Article Information
J. Granata, L. Xu, Z. Rusak

Department of Mechanical,
Aerospace, and Nuclear Engineering,
Rensselaer Polytechnic Institute,
Troy, NY 12180-3590

S. Wang

Department of Mathematics,
University of Auckland,
38 Princes Street,
Auckland 1142, New Zealand

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received September 22, 2015; final manuscript received March 16, 2016; published online June 6, 2016. Assoc. Editor: Praveen Ramaprabhu.

J. Fluids Eng 138(9), 091402 (Jun 06, 2016) (14 pages) Paper No: FE-15-1682; doi: 10.1115/1.4033321 History: Received September 22, 2015; Revised March 16, 2016

Current simulations of swirling flows in pipes are limited to relatively low Reynolds number flows (Re < 6000); however, the characteristic Reynolds number is much higher (Re > 20,000) in most of engineering applications. To address this difficulty, this paper presents a numerical simulation algorithm of the dynamics of incompressible, inviscid-limit, axisymmetric swirling flows in a pipe, including the vortex breakdown process. It is based on an explicit, first-order difference scheme in time and an upwind, second-order difference scheme in space for the time integration of the circulation and azimuthal vorticity. A second-order Poisson equation solver for the spatial integration of the stream function in terms of azimuthal vorticity is used. In addition, when reversed flow zones appear, an averaging step of properties is applied at designated time steps. This adds slight artificial viscosity to the algorithm and prevents growth of localized high-frequency numerical noise inside the breakdown zone that is related to the expected singularity that must appear in any flow simulation based on the Euler equations. Mesh refinement studies show agreement of computations for various mesh sizes. Computed examples of flow dynamics demonstrate agreement with linear and nonlinear stability theories of vortex flows in a finite-length pipe. Agreement is also found with theoretically predicted steady axisymmetric breakdown states in a pipe as flow evolves to a time-asymptotic state. These findings indicate that the present algorithm provides an accurate prediction of the inviscid-limit, axisymmetric breakdown process. Also, the numerical results support the theoretical predictions and shed light on vortex dynamics at high Re.

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

References

Leibovich, S. , 1984, “ Vortex Stability and Breakdown: Survey and Extension,” AIAA J., 22(9), pp. 1192–1206. [CrossRef]
Sarpkaya, T. , 1995, “ Turbulent Vortex Breakdown,” Phys. Fluids, 7(10), pp. 2301–2303. [CrossRef]
Novak, F. , and Sarpkaya, T. , 2000, “ Turbulent Vortex Breakdown at High Reynolds Numbers,” AIAA J., 38(5), pp. 825–834. [CrossRef]
Faler, J. H. , and Leibovich, S. , 1977, “ Disrupted States of Vortex Flow and Vortex Breakdown,” Phys. Fluids, 20(9), pp. 1385–1400. [CrossRef]
Grabowski, W. J. , and Berger, S. A. , 1976, “ Solutions of the Navier–Stokes Equations for Vortex Breakdown,” J. Fluid Mech., 75(3), pp. 525–544. [CrossRef]
Spall, R. E. , Gatski, T. B. , and Grosch, C. E. , 1987, “ A Criterion for Vortex Breakdown,” Phys. Fluids, 30(11), pp. 3434–3440. [CrossRef]
Spall, R. E. , Gatski, T. B. , and Ash, R. L. , 1990, “ The Structure and Dynamics of Bubble-Type Vortex Breakdown,” Proc. R. Soc. London, Ser. A, 429(1877), pp. 613–637. [CrossRef]
Darmofal, D. L. , 1996, “ Comparisons of Experimental and Numerical Results for Axisymmetric Vortex Breakdown in Pipes,” Comput. Fluids, 25(4), pp. 353–371. [CrossRef]
Beran, P. S. , and Culick, F. E. C. , 1992, “ The Role of Non-Uniqueness in the Development of Vortex Breakdown in Tubes,” J. Fluid Mech., 242, pp. 491–527. [CrossRef]
Beran, P. S. , 1994, “ The Time-Asymptotic Behavior of Vortex Breakdown in Tubes,” Comput. Fluids, 23(7), pp. 913–937. [CrossRef]
Lopez, J. M. , 1994, “ On the Bifurcation Structure of Axisymmetric Vortex Breakdown in a Constricted Pipe,” Phys. Fluids, 6(11), pp. 3683–3693. [CrossRef]
Leibovich, S. , and Kribus, A. , 1990, “ Large Amplitude Wave Trains and Solitary Waves in Vortices,” J. Fluid Mech., 216, pp. 459–504. [CrossRef]
Vyazmina, E. , Nichols, J. W. , Chomaz, J. M. , and Schmid, P. J. , 2009, “ The Bifurcation Structure of Viscous Steady Axisymmetric Vortex Breakdown With Open Lateral Boundaries,” Phys. Fluids, 21(7), p. 074107. [CrossRef]
Meliga, P. , and Gallaire, F. , 2011, “ Control of Axisymmetric Vortex Breakdown in a Constricted Pipe: Nonlinear Steady States and Weakly Nonlinear Asymptotic Expansions,” Phys. Fluids, 23(8), p. 084102. [CrossRef]
Snyder, D. O. , and Spall, R. E. , 2000, “ Numerical Simulation of Bubble-Type Vortex Breakdown Within a Tube-and-Vane Apparatus,” Phys. Fluids, 12(3), pp. 603–608. [CrossRef]
Wang, S. , and Rusak, Z. , 1996, “ On the Stability of an Axisymmetric Rotating Flow in a Pipe,” Phys. Fluids, 8(4), pp. 1007–1016. [CrossRef]
Wang, S. , and Rusak, Z. , 1997, “ The Dynamics of a Swirling Flow in a Pipe and Transition to Axisymmetric Vortex Breakdown,” J. Fluid Mech., 340, pp. 177–223. [CrossRef]
Liang, H. , and Maxworthy, T. , 2005, “ An Experimental Investigation of Swirling Jets,” J. Fluid Mech, 525, pp. 115–159. [CrossRef]
Benjamin, T. B. , 1962, “ Theory of the Vortex Breakdown Phenomenon,” J. Fluid Mech., 14(04), pp. 593–629. [CrossRef]
Gallaire, F. , Chomaz, J.-M. , and Huerre, P. , 2004, “ Closed-Loop Control of Vortex Breakdown: A Model Study,” J. Fluid Mech., 511, pp. 67–93. [CrossRef]
Rusak, Z. , Wang, S. , Xu, L. , and Taylor, S. , 2012, “ On the Global Nonlinear Stability of a Near-Critical Swirling Flow in a Long Finite-Length Pipe and the Path to Vortex Breakdown,” J. Fluid Mech., 712, pp. 295–326. [CrossRef]
Wang, S. , and Rusak, Z. , 1997, “ The Effect of Slight Viscosity on a Near-Critical Swirling Flow in a Pipe,” Phys. Fluids, 9(7), pp. 1914–1927. [CrossRef]
Keller, J. J. , Egli, W. , and Exley, W. , 1985, “ Force- and Loss-Free Transitions Between Flow States,” Z. Angew. Math. Phys., 36(6), pp. 854–889. [CrossRef]
Rusak, Z. , 1996, “ Axisymmetric Swirling Flow Around a Vortex Breakdown Point,” J. Fluid Mech., 323, pp. 79–105. [CrossRef]
Mattner, T. W. , Joubert, P. N. , and Chong, M. S. , 2002, “ Vortical Flow, Part 1: Flow Through a Constant Diameter Pipe,” J. Fluid Mech., 463, pp. 259–291. [CrossRef]
Umeh, C. O. U. , Rusak, Z. , Gutmark, E. , Villalva, R. , and Cha, D. J. , 2010, “ Experimental and Computational Study of Nonreacting Vortex Breakdown in a Swirl-Stabilized Combustor,” AIAA J., 48(11), pp. 2576–2585. [CrossRef]
Rusak, Z. , and Lamb, D. , 1999, “ Prediction of Vortex Breakdown in Leading-Edge Vortices Above Slender Delta Wings,” J. Aircr., 36(4), pp. 659–667. [CrossRef]
Batchelor, G. K. , 1956, “ On Steady Laminar Flow With Closed Streamlines at Large Reynolds Number,” J. Fluid Mech., 1(02), pp. 177–190. [CrossRef]
Squire, H. B. , 1956, “ Rotating Fluids,” Surveys in Mechanics, G. K. Batchelor and R. M. Davies, Cambridge University Press, Cambridge, UK, pp. 139–161.
Long, R. R. , 1953, “ Steady Motion Around a Symmetrical Obstacle Moving Along the Axis of a Rotating Liquid,” J. Meteorol., 10(3), pp. 197–203. [CrossRef]
Rusak, Z. , 1998, “ The Interaction of Near-Critical Swirling Flows in a Pipe With Inlet Azimuthal Vorticity Perturbations,” Phys. Fluids, 10(7), pp. 1672–1684. [CrossRef]
Gallaire, F. , and Chomaz, J.-M. , 2004, “ The Role of Boundary Conditions in a Simple Model of Incipient Vortex Breakdown,” Phys. Fluids, 16(2), pp. 274–286. [CrossRef]
Leclaire, B. , Sipp, D. , and Jacquin, L. , 2007, “ Near-Critical Swirling Flow in a Contracting Duct: The Case of Plug Axial Flow With Solid-Body Rotation,” Phys. Fluids, 19(9), p. 091701. [CrossRef]
Leclaire, B. , and Sipp, D. , 2010, “ A Sensitivity Study of Vortex Breakdown Onset to Upstream Boundary Conditions,” J. Fluid Mech., 645, pp. 81–119. [CrossRef]
Buntine, J. D. , and Saffman, P. G. , 1995, “ Inviscid Swirling Flows and Vortex Breakdown,” Proc. R. Soc. London, Ser. A, 449(1935), pp. 139–153. [CrossRef]
Xu, L. , 2012, “ Vortex Flow Stability, Dynamics and Feedback Stabilization,” Ph.D. dissertation, Rensselaer Polytechnic Institute, Troy, NY.
Granata, J. , 2014, “ An Active Feedback Flow Control Theory of the Vortex Breakdown Process,” Ph.D. dissertation, Rensselaer Polytechnic Institute, Troy, NY.
Rusak, Z. , Granata, J. , and Wang, S. , 2015, “ An Active Feedback Flow Control Theory of the Axisymmetric Vortex Breakdown Process,” J. Fluid Mech., 774, pp. 488–528. [CrossRef]
Webster, B. E. , Shephard, M. S. , Rusak, Z. , and Flaherty, J. E. , 1994, “ Automated Adaptive Time-Discontinuous Finite-Element Method for Unsteady Compressible Airfoil Aerodynamics,” AIAA J., 32(4), pp. 748–757. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

Time history of Amin (solid line), Aout (dashed line), and Amax (dotted line) computed from the meshes 100 × 100, 150 × 150, 200 × 200, and 300 × 300; for ω = 0.95 and δ = −4

Grahic Jump Location
Fig. 3

A snapshot of the flow evolution for ω = 0.95 and δ = −4 at (a) t = 50, (b) t ∼ 90, and at (c) time-asymptotic state, for meshes 100 × 100, 150 × 150, 200 × 200, and 300 × 300. The results from all meshes are essentially the same for all time.

Grahic Jump Location
Fig. 4

Comparison of the profile of stream function at the outlet from four meshes 100 × 100, 150 × 150, 200 × 200, 300 × 300 and the theoretical prediction of Wang and Rusak [17] (heavy line); for ω = 0.95

Grahic Jump Location
Fig. 5

Evolution of σ* as a function of t* for meshes 100 × 100 and 200 × 200. Also shown for reference are the linear stability growth rate according to Wang and Rusak [16] (dashed–dotted line) and the growth rate according to the weakly nonlinear model of Rusak et al. [21] (dashed line); for ω = 1.1 and δ = −4.

Grahic Jump Location
Fig. 6

Time history of Amin (solid line), Aout (dashed line), and Amax (dotted line) computed from the meshes 100 × 100 and 200 × 200; for ω = 1.1 and δ = −4

Grahic Jump Location
Fig. 7

Snapshots of the flow evolution for ω = 1.1 and δ = −4 at (a) t = 15 and at (b) time-asymptotic state, for meshes 100 × 100 and 200 × 200. The results from all meshes are essentially the same for all time.

Grahic Jump Location
Fig. 8

Comparison of the profile of stream function at the outlet from two meshes 100 × 100 (circles), 200 × 200 (triangles) and the theoretical prediction of Wang and Rusak [17] (solid line); for ω = 1.1

Grahic Jump Location
Fig. 9

Stream function contour-line snapshots of flow evolution from a perturbed columnar state to a breakdown state when ω = 0.8 and δ = −15 at (a) t = 0, (b) t = 66, (c) t = 85, (d) t = 105, (e) t = 134, (f) t = 155 (time-asymptotic state)

Grahic Jump Location
Fig. 10

Stream function contour-lines of time-asymptotic breakdown states at (a) ω = 0.8, (b) ω = 0.78, (c) ω = 0.76, and (d) ω = 0.74

Grahic Jump Location
Fig. 11

Stream function contour-line snapshots of flow evolution from a breakdown state at ω = 0.74 to a columnar state at ω = 0.735 at (a) t = 0, (b) t = 80, (c) t = 180, and at (d) time-asymptotic state

Grahic Jump Location
Fig. 1

Evolution of σ* as a function of t* for meshes 100 × 100, 150 × 150, 200 × 200, and 300 × 300. Also shown for reference are the linear stability growth rate according to Wang and Rusak [16] (dashed–dotted line) and the growth rate according to the weakly nonlinear model of Rusak et al. [21] (dashed line); for ω = 0.95 and δ = −4.

Grahic Jump Location
Fig. 12

Stream function contour-line snapshots of flow evolution from a near-columnar state to a breakdown state when ω = 1.2 and δ = −4 at (a) t = 0, (b) t = 10, (c) t = 25, (d) t = 40, (e) t = 60 and at (f) time-asymptotic state

Grahic Jump Location
Fig. 13

Comparison between the numerical results and predictions according to the theory of vortex breakdown

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