banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Transient growth analysis of flow through a sudden expansion in a circular pipe
Rent this article for


Image of FIG. 1.
FIG. 1.

Sketch illustrating the evolution of a perturbation through an expanding pipe. Small inlet perturbations are amplified in the region of the separated axisymmetric shear layer, but eventually decay downstream. Hence, even though the flow is linearly stable, it supports very strong transient growth of perturbations.

Image of FIG. 2.
FIG. 2.

Geometry of the expanding pipe. The computational flow domain is illustrated with the cylindrical coordinate system and the inlet and outlet lengths indicated (not to scale).

Image of FIG. 3.
FIG. 3.

A spectral-element mesh used in this study. For the case illustrated there are 563 elements, an inflow length of , and an outflow length of . When required, meshes with more elements and outflow lengths up to have been used. The two-dimensional mesh (a) used for linear analysis is extended to three dimensions (b) for nonlinear analysis (DNS) where Fourier expansions are used in azimuth.

Image of FIG. 4.
FIG. 4.

Contour plot of the base-flow streamfunction at , showing the separation and reattachment of the flow, and the recirculation region behind the expansion. Contours are drawn at intervals of 0.125 in the core of the flow and at intervals of 0.02 in the recirculation region.

Image of FIG. 5.
FIG. 5.

Relationship between downstream reattachment point and Reynolds number for base flows with a fully developed inlet profile up to . Points are the computed value of and the solid line shows the best-fit proportionality given by . The dotted lines indicate reattachment lengths for the base flows at (corresponding to Fig. 4) and .

Image of FIG. 6.
FIG. 6.

(a) Energy growth envelopes at for azimuthal mode numbers (as indicated) up to . (b) Enlargement of (a) for small . Curves for through are in decreasing monotonic order on the right-hand side.

Image of FIG. 7.
FIG. 7.

Contours of optimal transient energy growth as a function of time horizon and Reynolds number for (a) azimuthal mode number and (b) azimuthal mode number .

Image of FIG. 8.
FIG. 8.

Energy growth under linear evolution at , for three different initial conditions corresponding to optimal perturbations at , 70, and 110. Circles denote the optimal growth envelope. Each linear evolution curve touches the envelope at its respective value, as indicated by a filled circle.

Image of FIG. 9.
FIG. 9.

Evolution of optimal initial disturbance in the mode visualized through contours/isosurfaces of azimuthal velocity at from (bottom) in time intervals of four units in the spatial range of . The panel labeled shows the evolved disturbance at its maximum growth: here the spatial range is and the isosurface levels are two orders of magnitude larger than in the other panels.

Image of FIG. 10.
FIG. 10.

Evolution of optimal initial disturbance in the mode visualized through isosurfaces of radial velocity at from (bottom) in time intervals of two units in the spatial range of .

Image of FIG. 11.
FIG. 11.

Physical interpretation of the maximal disturbance at . Shown is a linear superposition of the base flow with the optimal disturbance at the time of maximum growth, . The disturbance has a relative energy magnitude of 5% compared with the base flow. The visualization is a semitransparent isosurface of azimuthal vorticity at a level highlighting the separated shear layer.

Image of FIG. 12.
FIG. 12.

Optimal energy growth , as a function of Re for each of the first four azimuthal modes. Above , increases exponentially with Reynolds number at a rate of 0.45 orders of magnitude for each increase in 100 in Re. Inset shows detail at small Re. Curves intersect at the critical Reynolds number of each mode .

Image of FIG. 13.
FIG. 13.

Location of the centroid of the disturbance at its maximum growth (points connect by solid lines), compared to the location of the reattachment point (dashed line). The disturbance reaches its maximum approximately five diameters upstream of reattachment point.

Image of FIG. 14.
FIG. 14.

Leading four growth envelopes for the optimal and suboptimal perturbations for the mode at . The circles on the dotted line at show the first eight leading eigenvalues at the point of the peak optimal growth.

Image of FIG. 15.
FIG. 15.

Isosurfaces of azimuthal velocity for the four leading perturbations at , evolved to time . (a) and (b) are the optimal and first suboptimal modes. (c) and (d) are the next pair of suboptimal modes (see Fig. 14).

Image of FIG. 16.
FIG. 16.

Modal energy in a noisy inflow DNS of the expanding-pipe flow at (a) , (b) , and (c) . Modal energies visible above the noise floor of are labeled. Initially the flow is the steady axisymmetric base flow. The axisymmetric energy is off the scale of the figure.

Image of FIG. 17.
FIG. 17.

Visualization of noise-driven flow and optimal perturbation at (a) and (b) in the range of . The upper half of each pipe shows the optimal perturbation at the respective Reynolds number, evolved to the point of maximum growth. The lower half shows a snapshot of the noise-driven flow for comparison.

Image of FIG. 18.
FIG. 18.

Isosurfaces of streamwise (top), radial (middle) and azimuthal (bottom) velocity components for the noise-driven flow showing the downstream disturbance induced by the stochastic forcing at the same time instant as Fig. 17(b).

Image of FIG. 19.
FIG. 19.

Energy of noise-driven flow through the line for the (solid lines) and (dotted lines). [(a)–(c)] , (d) , and (e) at 1/4 the noise level of (d). Vertical lines indicate , the centroid of the optimal linear perturbation at the corresponding value of Re.

Image of FIG. 20.
FIG. 20.

Radial profiles of velocity component standard deviations comparing the noise-driven simulation (left) with the linear analysis (right). The three velocity components each normalized to their peak value are shown: streamwise (solid line), radial velocity (dashed line), and azimuthal velocity (dotted line).

Image of FIG. 21.
FIG. 21.

Dependence of maximum transient energy growth on bulk Reynolds numbers for three separated flows: two-dimensional disturbances in a backward-facing step flow (Ref. 23); steady stenotic flow (Ref. 9); and the present 1–2 axisymmetric expansion, compared with Hagen–Poiseuille flow (Refs. 44 and 47).


Generic image for table
Table I.

Effect of domain length on representative quantities at . The reattachment point, leading eigenvalue, and optimal growth for are given for different inflow and outflow lengths. The polynomial order is .

Generic image for table
Table II.

Streamwise position of base flow reattachment point (in units of step height, ) at , as a function of spectral-element polynomial order . The outflow length is .

Generic image for table
Table III.

Leading eigenvalues from a linear stability analysis of flows in the 1–2 axisymmetric expansion for Reynolds numbers indicated. All values correspond to azimuthal mode number . All are real and negative and hence the axisymmetric base flow is linearly stable up to at least .

Generic image for table
Table IV.

Table of critical Reynolds numbers , maximum growth values , and corresponding time horizons for maximum growth, (for and as indicated) for each of the first four azimuthal modes, .

Generic image for table
Table V.

Characteristics of the optimal perturbations at the time of maximum growth for modes at Reynolds numbers indicated. Along with the centroid location , of the evolved perturbation energy, as well as the local axial wavelength and temporal frequency St of the perturbation.

Generic image for table
Table VI.

Statistics for DNS with stochastic inflow forcing at . Here we define the centroid based on the distribution of turbulent kinetic energy in the outflow section of the pipe.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Transient growth analysis of flow through a sudden expansion in a circular pipe