Global stability of multiple solutions in plane sudden-expansion flow

2012 ◽  
Vol 702 ◽  
pp. 378-402 ◽  
Author(s):  
Daniel Lanzerstorfer ◽  
Hendrik C. Kuhlmann

AbstractThe two-dimensional, incompressible flow in a plane sudden expansion is investigated numerically for a systematic variation of the geometry, covering expansion ratios (steps-to-outlet heights) from $0. 25$ to $0. 95$. By means of a three-dimensional linear stability analysis global temporal modes are scrutinized. In a symmetric expansion the primary bifurcation is stationary and two-dimensional, breaking the mirror symmetry with respect to the mid-plane. The secondary asymmetric flow experiences a secondary instability to different three-dimensional modes, depending on the expansion ratio. For a moderately asymmetric expansion only one of the two secondary flows (the connected branch) is realized at low Reynolds numbers. Since the perturbed secondary flow does not deviate much from the symmetric secondary flow, both secondary stability boundaries are very close to each other. For very small and very large expansion ratios an asymptotic behaviour is found for suitably scaled critical Reynolds numbers and wavenumbers. Representative instabilities are analysed in detail using an a posteriori energy transfer analysis to reveal the physical nature of the instabilities. Depending on the geometry, pure centrifugal and elliptical amplification processes are identified. We also find that the basic flow can become unstable due to the effects of flow deceleration, streamline convergence and high shear stresses, respectively.

Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


2005 ◽  
Vol 128 (4) ◽  
pp. 671-679 ◽  
Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations failed to capture completely the total expansion effect on the flow, which is influenced by both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When two-dimensional simulations were performed using the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations were compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, thus sustaining laminar flow symmetry to higher Reynolds numbers. Last, and most important, when the logarithm of the critical Reynolds number was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


2011 ◽  
Vol 693 ◽  
pp. 1-27 ◽  
Author(s):  
Daniel Lanzerstorfer ◽  
Hendrik C. Kuhlmann

AbstractThe two-dimensional, incompressible flow over a backward-facing step is considered for a systematic variation of the geometry covering expansion ratios (step to outlet height) from 0.25 to 0.975. A global temporal linear stability analysis shows that the basic flow becomes unstable to different three-dimensional modes depending on the expansion ratio. All critical modes are essentially confined to the region behind the step extending downstream up to the reattachment point of the separated eddy. An energy-transfer analysis is applied to understand the physical nature of the instabilities. If scaled appropriately, the critical Reynolds number approaches a finite asymptotic value for very large step heights. In that case centrifugal forces destabilize the flow with respect to an oscillatory critical mode. For moderately large expansion ratios an elliptical instability mechanism is identified. If the step height is further decreased the critical mode changes from oscillatory to stationary. In addition to the elliptical mechanism, the strong shear in the layer emanating from the sharp corner of the step supports the amplification process of the critical mode. For very small step heights the basic state becomes unstable due to the lift-up mechanism, which feeds back on itself via the recirculating eddy behind the step, resulting in a steady critical mode comprising pronounced slow and fast streaks.


1994 ◽  
Vol 261 ◽  
pp. 169-198 ◽  
Author(s):  
F. J. Wang ◽  
G. A. Domoto

The hydrodynamic instability of a viscous incompressible flow with a free surface is studied both numerically and experimentally. While the free-surface flow is basically two-dimensional at low Reynolds numbers, a three-dimensional secondary flow pattern similar to the Taylor vorticies between two concentric cylinders appears at higher rotational speeds. The secondary flow has periodic velocity components in the axial direction and is characterized by a distinct spatially periodic variation in surface height similar to a standing wave. A numerical method, using boundary-fitted coordinates and multigrid methods to solve the Navier–Stokes equations in primitive variables, is developed to treat two-dimensional free-surface flows. A similar numerical technique is applied to the linearized three-dimensional perturbation equations to treat the onset of secondary flows. Experimental measurements have been obtained using light sheet techniques to visualize the secondary flow near the free surface. Photographs of streak lines were taken and compared to the numerical calculations. It has been shown that the solution of the linearized equations contains most of the important features of the nonlinear secondary flows at Reynolds number higher than the critical value. The experimental results also show that the numerical method predicts well the onset of instability in terms of the critical wavenumber and Reynolds number.


2013 ◽  
Vol 727 ◽  
pp. 533-581 ◽  
Author(s):  
D. P. Wall ◽  
M. Nagata

AbstractThree-dimensional exact, finite-amplitude solutions are presented for the problem of channel flow subject to a system rotation about a spanwise axis. The solutions are of travelling wave form, and may bifurcate as tertiary flows from the two-dimensional streamwise-independent secondary flow, or as secondary flows directly from the basic flow. For the tertiary flows, we consider solutions of spanwise superharmonic and subharmonic type. We distinguish flows on the basis of symmetry, originating eigenmode and major solution branch, and thus identify 15 distinct flows: 5 superharmonic tertiary, 5 subharmonic tertiary and 5 secondary flows. The tertiary flows all feature a single layer of vortical structures in the spanwise–wall-normal plane, the secondary flows feature single-, double-, triple- or quadruple-layer flow structures in this plane. All flows feature low-speed streamwise-orientated streaks in the streamwise velocity component and/or pulses of low-speed streamwise velocity. The streaks may be sinusoidal or varicose. Sinusoidal streaks are flanked by staggered streamwise vortices, varicose streaks and pulses are flanked by aligned vortices. A comparison with previous simulation and experimental studies finds that the simplest three-dimensional flows observed previously correspond to superharmonic tertiary flows bifurcating from the upper branch of the secondary flow. The mean absolute vorticity of the present flows is also considered. A flattening of the profile of this vorticity is observed in the central region of the channel for two-dimensional secondary and many of the three-dimensional flows, with two-step profiles also observed. This phenomenon is attributed to mixing of the vorticity across zones of the channel in which streamwise vortex structures exist, and is demonstrated by a two-dimensional model. The phenomenon appears to be distinct to that observed in fully turbulent rotating channel flows.


2002 ◽  
Vol 473 ◽  
pp. 167-190 ◽  
Author(s):  
DWIGHT BARKLEY ◽  
M. GABRIELA M. GOMES ◽  
RONALD D. HENDERSON

Results are reported from a three-dimensional computational stability analysis of flow over a backward-facing step with an expansion ratio (outlet to inlet height) of 2 at Reynolds numbers between 450 and 1050. The analysis shows that the first absolute linear instability of the steady two-dimensional flow is a steady three-dimensional bifurcation at a critical Reynolds number of 748. The critical eigenmode is localized to the primary separation bubble and has a flat roll structure with a spanwise wavelength of 6.9 step heights. The system is further shown to be absolutely stable to two-dimensional perturbations up to a Reynolds number of 1500. Stability spectra and visualizations of the global modes of the system are presented for representative Reynolds numbers.


1996 ◽  
Vol 118 (2) ◽  
pp. 347-352 ◽  
Author(s):  
R. G. Dominy ◽  
D. A. Kirkham

Interturbine diffusers provide continuity between HP and LP turbines while diffusing the flow upstream of the LP turbine. Increasing the mean turbine diameter offers the potential advantage of reducing the flow factor in the following stages, leading to increased efficiency. The flows associated with these interturbine diffusers differ from those in simple annular diffusers both as a consequence of their high-curvature S-shaped geometry and of the presence of wakes created by the upstream turbine. It is shown that even the simplest two-dimensional wakes result in significantly modified flows through such ducts. These introduce strong secondary flows demonstrating that fully three-dimensional, viscous analysis methods are essential for correct performance modeling.


2000 ◽  
Vol 411 ◽  
pp. 1-38 ◽  
Author(s):  
C. ROSS ETHIER ◽  
SUJATA PRAKASH ◽  
DAVID A. STEINMAN ◽  
RICHARD L. LEASK ◽  
GREGORY G. COUCH ◽  
...  

Numerical and experimental techniques were used to study the physics of flow separation for steady internal flow in a 45° junction geometry, such as that observed between two pipes or between the downstream end of a bypass graft and an artery. The three-dimensional Navier–Stokes equations were solved using a validated finite element code, and complementary experiments were performed using the photochromic dye tracer technique. Inlet Reynolds numbers in the range 250 to 1650 were considered. An adaptive mesh refinement approach was adopted to ensure grid-independent solutions. Good agreement was observed between the numerical results and the experimentally measured velocity fields; however, the wall shear stress agreement was less satisfactory. Just distal to the ‘toe’ of the junction, axial flow separation was observed for all Reynolds numbers greater than 250. Further downstream (approximately 1.3 diameters from the toe), the axial flow again separated for Re [ges ] 450. The location and structure of axial flow separation in this geometry is controlled by secondary flows, which at sufficiently high Re create free stagnation points on the model symmetry plane. In fact, separation in this flow is best explained by a secondary flow boundary layer collision model, analogous to that proposed for flow in the entry region of a curved tube. Novel features of this flow include axial flow separation at modest Re (as compared to flow in a curved tube, where separation occurs only at much higher Re), and the existence and interaction of two distinct three-dimensional separation zones.


Author(s):  
A. Perdichizzi ◽  
V. Dossena

This paper describes the results of an experimental investigation of the three-dimensional flow downstream of a linear turbine cascade at off-design conditions. The tests have been carried out for five incidence angles from −60 to +35 degrees, and for three pitch-chord ratios: s/c = 0.58,0.73,0.87. Data include blade pressure distributions, oil flow visualizations, and pressure probe measurements. The secondary flow field has been obtained by traversing a miniature five hole probe in a plane located at 50% of an axial chord downstream of the trailing edge. The distributions of local energy loss coefficients, together with vorticity and secondary velocity plots show in detail how much the secondary flow field is modified both by incidence and cascade solidity variations. The level of secondary vorticity and the intensity of the crossflow at the endwall have been found to be strictly related to the blade loading occurring in the blade entrance region. Heavy changes occur in the spanwise distributions of the pitch averaged loss and of the deviation angle, when incidence or pitch-chord ratio is varied.


Author(s):  
S. Friedrichs ◽  
H. P. Hodson ◽  
W. N. Dawes

The endwall film-cooling cooling configuration investigated by Friedrichs et al. (1996, 1997) had in principle sufficient cooling flow for the endwall, but in practice, the redistribution of this coolant by secondary flows left large endwall areas uncooled. This paper describes the attempt to improve upon this datum cooling configuration by redistributing the available coolant to provide a better coolant coverage on the endwall surface, whilst keeping the associated aerodynamic losses small. The design of the new, improved cooling configuration was based on the understanding of endwall film-cooling described by Friedrichs et al. (1996, 1997). Computational fluid dynamics were used to predict the basic flow and pressure field without coolant ejection. Using this as a basis, the above described understanding was used to place cooling holes so that they would provide the necessary cooling coverage at minimal aerodynamic penalty. The simple analytical modelling developed in Friedrichs et al. (1997) was then used to check that the coolant consumption and the increase in aerodynamic loss lay within the limits of the design goal. The improved cooling configuration was tested experimentally in a large scale, low speed linear cascade. An analysis of the results shows that the redesign of the cooling configuration has been successful in achieving an improved coolant coverage with lower aerodynamic losses, whilst using the same amount of coolant as in the datum cooling configuration. The improved cooling configuration has reconfirmed conclusions from Friedrichs et al. (1996, 1997); firstly, coolant ejection downstream of the three-dimensional separation lines on the endwall does not change the secondary flow structures; secondly, placement of holes in regions of high static pressure helps reduce the aerodynamic penalties of platform coolant ejection; finally, taking account of secondary flow can improve the design of endwall film-cooling configurations.


Sign in / Sign up

Export Citation Format

Share Document