Calculation of Laminar Separated Flow in Symmetric Two-Dimensional Diffusers

1995 ◽  
Vol 117 (4) ◽  
pp. 612-616 ◽  
Author(s):  
Yeng-Yung Tsui ◽  
Chia-Kang Wang

This study is concerned with numerical analysis of laminar separated flow in symmetric, two-dimensional, straight-walled diffusers. With Reynolds numbers Re = 56 and 114 and expansion ratios ER = 3 and 4, totally, there are four cases considered. At the low Reynolds number and the low expansion ratio the flow in the diffuser is nearly symmetric to the center line, irrespective of the diffusion angle. As Reynolds number or expansion ratio increases, a large recirculation region forms at one side wall and a small one at the other side. For the case with Re = 114 and ER = 4 the small recirculating flow disappears at small diffusion angles and a third recirculating flow appear in the same side of the small main recirculation region for large diffusion angles. The pressure recovery reaches its peak value somewhere downstream of the reattachment point of the large recirculating flow. The effectiveness of the diffuser deteriorates as the diffusion angle increases, apart from that at Re = 56 the effectiveness increases from θ = 15 to 30 deg. Symmetric flow solutions can be obtained by incorporating a symmetric relaxation method. The pressure recovery is higher for the symmetric flow than that for the asymmetric flow owing to the weaker recirculating strength in the former.

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.


1988 ◽  
Vol 110 (3) ◽  
pp. 289-296 ◽  
Author(s):  
F. Durst ◽  
J. C. F. Pereira

This paper presents results of numerical studies of the impulsively starting backward-facing step flow with the step being mounted in a plane, two-dimensional duct. Results are presented for Reynolds numbers of Re = 10; 368 and 648 and for the last two Reynolds numbers comparisons are given between experimental and numerical results obtained for the final steady state flow conditions. In the computational scheme, the convective terms in the momentum equations are approximated by a 13-point quadratic upstream weighted finite-difference scheme and a fully implicit first order forward differencing scheme is used to discretize the temporal derivatives. The computations show that for the higher Reynolds numbers, the flow starts to separate on the lower and upper corners of the step yielding two disconnected recirculating flow regions for some time after the flow has been impulsively started. As time progresses, these two separated flow regions connect up and a single recirculating flow region emerges. This separated flow region stays attached to the step, grows in size and approaches, for the time t → ∞, the dimensions measured and predicted for the separation region for steady laminar backward-facing flow. For the Reynolds number Re = 10 the separation starts at the bottom of the backward-facing step and the separation region enlarges with time until the steady state flow pattern is reached. At the channel wall opposite to the step and for Reynolds number Re = 368, a separated flow region is observed and it is shown to occur for some finite time period of the developing, impulsively started backward-facing step flow. Its dimensions change with time and reduce to zero before the steady state flow pattern is reached. For the higher Reynolds number Re = 648, the secondary separated flow region opposite to the wall is also present and it is shown to remain present for t → ∞. Two kinds of the inlet conditions were considered; the inlet mean flow was assumed to be constant in a first study and was assumed to increase with time in a second one. The predicted flow field for t → ∞ turned out to be identical for both cases. They were also identical to the flow field predicted for steady, backward-facing step flow using the same numerical grid as for the time-dependent predictions.


2011 ◽  
Vol 681 ◽  
pp. 411-433 ◽  
Author(s):  
HEMANT K. CHAURASIA ◽  
MARK C. THOMPSON

A detailed numerical study of the separating and reattaching flow over a square leading-edge plate is presented, examining the instability modes governing transition from two- to three-dimensional flow. Under the influence of background noise, experiments show that the transition scenario typically is incompletely described by either global stability analysis or the transient growth of dominant optimal perturbation modes. Instead two-dimensional transition effectively can be triggered by the convective Kelvin–Helmholtz (KH) shear-layer instability; although it may be possible that this could be described alternatively in terms of higher-order optimal perturbation modes. At least in some experiments, observed transition occurs by either: (i) KH vortices shedding downstream directly and then almost immediately undergoing three-dimensional transition or (ii) at higher Reynolds numbers, larger vortical structures are shed that are also three-dimensionally unstable. These two paths lead to distinctly different three-dimensional arrangements of vortical flow structures. This paper focuses on the mechanisms underlying these three-dimensional transitions. Floquet analysis of weakly periodically forced flow, mimicking the observed two-dimensional quasi-periodic base flow, indicates that the two-dimensional vortex rollers shed from the recirculation region become globally three-dimensionally unstable at a Reynolds number of approximately 380. This transition Reynolds number and the predicted wavelength and flow symmetries match well with those of the experiments. The instability appears to be elliptical in nature with the perturbation field mainly restricted to the cores of the shed rollers and showing the spatial vorticity distribution expected for that instability type. Indeed an estimate of the theoretical predicted wavelength is also a good match to the prediction from Floquet analysis and theoretical estimates indicate the growth rate is positive. Fully three-dimensional simulations are also undertaken to explore the nonlinear development of the three-dimensional instability. These show the development of the characteristic upright hairpins observed in the experimental dye visualisations. The three-dimensional instability that manifests at lower Reynolds numbers is shown to be consistent with an elliptic instability of the KH shear-layer vortices in both symmetry and spanwise wavelength.


Author(s):  
Aya Kito ◽  
Kazuaki Sugawara ◽  
Hiroyuki Yoshikawa ◽  
Terukazu Ota

The direct numerical simulation methodology was employed to analyze the unsteady features of a three-dimensional separated flow and heat transfer around a downward step in a rectangular channel, and to clarify systematically the channel expansion ratio effects upon them. Numerical calculations were carried out using the finite difference method. The Reynolds number Re based on the mean velocity at inlet and the step height was varied from 300 to 1000. The channel expansion ratio ER is 1.5, 2.0 and 3.0 under a step aspect ratio of 36.0. It is found that the flow is steady upto Re = 500 but becomes sensibly unsteady at Re = 700 for all the three expansion ratios. In the case of ER = 2.0, the separated shear layer is most unstable. In the case of ER = 1.5, the longitudinal vortices formed near the side walls of channel are strongest. Nusselt number reaches its maximum in the reattachment flow region and also in the neighborhood of the side wall, and their locations depend greatly upon ER and Re.


Author(s):  
Majid Nabavi ◽  
Luc Mongeau

In this study, two-dimensional laminar incompressible and turbulent compressible flow through the planar diffuser (gradual expansion) for different divergence half angles of the diffuser (θ), and different Reynolds numbers (Re) was numerically studied. The effects of θ on the critical Reynolds number at which the onset of asymmetric flow is observed, were investigated. In the laminar flow regime, it was observed that for every values of θ, there is a critical Re beyond which the flow is asymmetric. However, in the turbulent flow regime, for θ ≥ 20°, even at low Reynolds number the flow is asymmetric. Only for θ ≤ 10°, symmetric flow was observed below a critical Re.


1997 ◽  
Vol 351 ◽  
pp. 167-199 ◽  
Author(s):  
S. BALACHANDAR ◽  
R. MITTAL ◽  
F. M. NAJJAR

The properties of the time- and span-averaged mean wake recirculation region are investigated in separated flows over several different two-dimensional bluff bodies. Ten different cases are considered and they divide into two groups: cylindrical geometries of circular, elliptic and square cross-sections and the normal plate. A wide Reynolds number range from 250 to 140000 is considered, but in all the cases the attached portion of the boundary layer remains laminar until separation. The lower Reynolds number data are from direct numerical simulations, while the data at the higher Reynolds number are obtained from large-eddy simulation and the experimental work of Cantwell & Coles (1983), Krothapalli (1996, personal communication), Leder (1991) and Lyn et al. (1995). Unlike supersonic and subsonic separations with a splitter plate in the wake, in all the cases considered here there is strong interaction between the shear layers resulting in Kármán vortex shedding. The impact of this fundamental difference on the distribution of Reynolds stress components and pressure in relation to the mean wake recirculation region (wake bubble) is considered. It is observed that in all cases the contribution from Reynolds normal stress to the force balance of the wake bubble is significant. In fact, in the cylinder geometries this contribution can outweigh the net force from the shear stress, so that the net pressure force tends to push the bubble away from the body. In contrast, in the case of normal plate, owing to the longer wake, the net contribution from shear stress outweighs that from the normal stress. At higher Reynolds numbers, separation of the Reynolds stress components into incoherent contributions provides more insight. The behaviour of the coherent contribution, arising from the dominant vortex shedding, is similar to that at lower Reynolds numbers. The incoherent contribution to Reynolds stress, arising from small-scale activity, is compared with that of a canonical free shear layer. Based on these observations a simple extension of the wake model (Sychev 1982; Roshko 1993a, b) is proposed.


1986 ◽  
Vol 108 (1) ◽  
pp. 93-97 ◽  
Author(s):  
R. H. Buckholz

Questions concerning the functional role of spanwise wing corrugation in living systems are experimentally investigated. Attention was initially directed to this problem by observation of the irregular shape of many insect wings as well as other studies indicating higher lift on these wings. First, a flow visualization scheme was used to observe and photograph streamlines around two different wing sections. One of these, a sheet metal model with geometry matching that of a butterfly wing, was studied at a chord Reynolds number of 1500 and at a Reynolds number of 80 based on corrugation depth. A steady-state recirculation region near the model leading edge was found, and the separated flow region above this recirculation zone formed a laminar reattachment to the model. A second thicker wing was corrugated on the upper surface. Closed streamlines inside these upper surface corrugations were photographed at Reynolds numbers of 8000 and 3800 based on chord length, and 200 and 90 based on corrugation depth. Reductions in pressures on the corrugated upper wing surface relative to a smooth upper wing surface were then measured.


1967 ◽  
Vol 27 (1) ◽  
pp. 97-109 ◽  
Author(s):  
D. B. Spalding

A power-law relation is derived between the Stanton number and the Reynolds number, expressing the law of heat transfer for a wall adjacent to a region of turbulent separated flow. The derivation is based on Prandtl's (1945) proposal for the laws of dissipation, diffusion and generation of turbulent kinetic energy. The constants appearing in these laws are determined by reference to experimental data for the hydrodynamic properties of the constant-stress and the linear-stress layers.The agreement between the resulting predictions and the experimental data of other workers is sufficiently good to suggest that the actual mechanism of heat transfer from separated flows has much in common with that which is postulated. Closer agreement can be expected only after the present one-dimensional analysis has been superseded by a two-dimensional one.


2010 ◽  
Vol 648 ◽  
pp. 225-256 ◽  
Author(s):  
B. E. STEWART ◽  
M. C. THOMPSON ◽  
T. LEWEKE ◽  
K. HOURIGAN

A study investigating the flow around a cylinder rolling or sliding on a wall has been undertaken in two and three dimensions. The cylinder motion is specified from a set of five discrete rotation rates, ranging from prograde through to retrograde rolling. A Reynolds number range of 20–500 is considered. The effects of the nearby wall and the imposed body motion on the wake structure and dominant wake transitions have been determined. Prograde rolling is shown to destabilize the wake flow, while retrograde rotation delays the onset of unsteady flow to Reynolds numbers well above those observed for a cylinder in an unbounded flow.Two-dimensional simulations show the presence of two recirculation zones in the steady wake, the lengths of which increase approximately linearly with the Reynolds number. Values of the lift and drag coefficient are also reported for the steady flow regime. Results from a linear stability analysis show that the wake initially undergoes a regular bifurcation from a steady two-dimensional flow to a steady three-dimensional wake for all rotation rates. The critical Reynolds number Rec of transition and the spanwise wavelength of the dominant mode are shown to be highly dependent on, but smoothly varying with, the rotation rate of the cylinder. Varying the rotation from prograde to retrograde rolling acts to increase the value of Rec and decrease the preferred wavelength. The structure of the fully evolved wake mode is then established through three-dimensional simulations. In fact it is found that at Reynolds numbers only marginally (~5%) above critical, the three-dimensional simulations indicate that the saturated state becomes time dependent, although at least initially, this does not result in a significant change to the mode structure. It is only at higher Reynolds numbers that the wake undergoes a transition to vortex shedding.An analysis of the three-dimensional transition indicates that it is unlikely to be due to a centrifugal instability despite the superficial similarity to the flow over a backward-facing step, for which the transition mechanism has been speculated to be centrifugal. However, the attached elongated recirculation region and distribution of the spanwise perturbation vorticity field, and the similarity of these features with those of the flow through a partially blocked channel, suggest the possibility that the mechanism is elliptic in nature. Some analysis which supports this conjecture is undertaken.


2005 ◽  
Vol 109 (1097) ◽  
pp. 325-335
Author(s):  
M. Ali ◽  
T. Fujiwara

Abstract A numerical study on mixing of air and hydrogen is performed by solving two-dimensional full Navier-Stokes equations. The main stream is air of Mach 5 entering through the configured inlet of the combustor and gaseous hydrogen is injected from the configured jet on the side wall. Supersonic mixing and diffusion mechanisms of a transverse hydrogen jet in two-dimensional finite air streams have been analyzed and discussed. The computed results are compared with the experimental data and show good agreement. For an otherwise fixed combustor geometry, the air inlet width and injection angle are varied to study the physics of mixing and flow field characteristics. On the effect of inlet width variation, two competing phenomena have been observed: (i) upstream of injector the strength of recirculation is higher for wider inlet and consequently the mixing increases, and (ii) downstream, the diffusion of hydrogen decreases with the increase of inlet width and eventually mixing decreases. As a result, in far downstream the mixing efficiency increases up to certain inlet width and then decreases for further increment of inlet width. For the variation of injection angle results show that upstream of injector the mixing is dominated by recirculation and downstream the mixing is dominated by mass concentration of hydrogen. Upstream recirculation is dominant for injecting angle 60° and 90°. Incorporating the various effects, perpendicular injection shows the maximum mixing efficiency and its large upstream recirculation region has a good flame holding capability.


Sign in / Sign up

Export Citation Format

Share Document