Bifurcation analysis of two-dimensional laminar flow through sudden expansion channel

In this work, bifurcation characteristics of unsteady, viscous, Newtonian laminar flow in two-dimensional sudden expansion and sudden contraction-expansion channels have been studied for different values of expansion ratio. The governing equations have been solved using finite volume method and FLUENT software has been employed to visualize the simulation results. Three different mesh studies have been performed to calculate critical Reynolds number (Recr) for different types of bifurcation phenomena. It is found that Recr decreases with the increase in expansion ratio (ER).

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.


1997 ◽  
Vol 330 ◽  
pp. 169-188 ◽  
Author(s):  
N. ALLEBORN ◽  
K. NANDAKUMAR ◽  
H. RASZILLIER ◽  
F. DURST

Two-dimensional laminar flow of an incompressible viscous fluid through a channel with a sudden expansion is investigated. A continuation method is applied to study the bifurcation structure of the discretized governing equations. The stability of the different solution branches is determined by an Arnoldi-based iterative method for calculating the most unstable eigenmodes of the linearized equations for the perturbation quantities. The bifurcation picture is extended by computing additional solution branches and bifurcation points. The behaviour of the critical eigenvalues in the neighbourhood of these bifurcation points is studied. Limiting cases for the geometrical and flow parameters are considered and numerical results are compared with analytical solutions for these cases.


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.


2018 ◽  
Vol 7 (4.35) ◽  
pp. 148 ◽  
Author(s):  
Nur Irmawati Om ◽  
Rozli Zulkifli ◽  
P. Gunnasegaran

The influence of utilizing different nanofluids types on the liquid cold plate (LCP) is numerically investigated. The thermal and fluid flow performance of LCP is examined by using pure ethylene glycol (EG), Al2O3-EG and CuO-EG. The volume fraction of the nanoparticle for both nanofluid is 2%. The finite volume method (FVM) has been used to solved 3-D steady state, laminar flow and heat transfer governing equations. The presented results indicate that Al2O3-EG able to provide the lowest surface temperature of the heater block followed by CuO-EG and EG, respectively. It is also found that the pressure drop and friction factor are higher for Al2O3-EG and CuO-EG compared to the pure EG.


A class of two-dimensional channels, with walls whose radius of curvature is uniformly large relative to local channel width, is described, and the velocity field of laminar flow through these channels is obtained as a power series in the small curvature parameter. The leading term is the Jeffery-Hamel solution considered in part I, and it is shown here how the higher-order terms are found. Terms of the third approximation have been computed. The theory is applied to two examples, for one of which experimental results are available and confirm the theoretical values with fair accuracy.


2019 ◽  
Vol 24 (2) ◽  
pp. 387-410
Author(s):  
Md. Shahjada Tarafder ◽  
M. Al Mursaline

Abstract This paper deals with the numerical simulation of a turbulent flow around two-dimensional bodies by the finite volume method with non-orthogonal body-fitted grid. The governing equations are expressed in Cartesian velocity components and solution is carried out using the SIMPLE algorithm for collocated arrangement of scalar and vector variables. Turbulence is modeled by the k- ε turbulence model and wall functions are used to bridge the solution variables at the near wall cells and the corresponding quantities on the wall. A simplified pressure correction equation is derived and proper under-relaxation factors are used so that computational cost is reduced without adversely affecting the convergence rate. The numerical procedure is validated by comparing the computed pressure distribution on the surface of NACA 0012 and NACA 4412 hydrofoils for different angles of attack with experimental data. The grid dependency of the solution is studied by varying the number of cells of the C-type structured mesh. The computed lift coefficients of NACA 4412 hydrofoil at different angles of attack are also compared with experimental results to further substantiate the validity of the proposed methodology.


Author(s):  
A. F. Najafi ◽  
M. H. Saidi ◽  
M. S. Sadeghipour ◽  
M. Souhar

In this study the developing turbulent swirling pipe flow is investigated both numerically and analytically. Governing equations are derived accompanying the boundary layer assumptions. Uniform and solid body rotation (SBR) distributions are taken into account for the axial and tangential velocities at the inlet of the pipe, respectively. Beyond the boundary layers, the flow pattern is considered to be the potential flow. Making use of the fourth-order Runge-Kutta scheme, the numerical solution of the differential equations is obtained. Further more, by simplifying the governing equations for large Rossby number, the analytical solution is performed. The results of numerical and analytical swirl intensity have been compared showing reasonable agreement. As an alternative solution, a CFD analysis has been done as well, having applied FLUENT software to support the ability of our methodology.


2001 ◽  
Vol 436 ◽  
pp. 283-320 ◽  
Author(s):  
T. HAWA ◽  
Z. RUSAK

Bifurcation analysis, linear stability study, and direct numerical simulations of the dynamics of a two-dimensional, incompressible, and laminar flow in a symmetric long channel with a sudden expansion with right angles and with an expansion ratio D/d (d is the width of the channel inlet section and D is the width of the outlet section) are presented. The bifurcation analysis of the steady flow equations concentrates on the flow states around a critical Reynolds number Rec(D/d) where asymmetric states appear in addition to the basic symmetric states when Re [ges ] Rec(D/d). The bifurcation of asymmetric states at Rec has a pitchfork nature and the asymmetric perturbation grows like √Re − Rec(D/d). The stability analysis is based on the linearized equations of motion for the evolution of infinitesimal two-dimensional disturbances imposed on the steady symmetric as well as asymmetric states. A neutrally stable asymmetric mode of disturbance exists at Rec(D/d) for both the symmetric and the asymmetric equilibrium states. Using asymptotic methods, it is demonstrated that when Re < Rec(D/d) the symmetric states have an asymptotically stable mode of disturbance. However, when Re > Rec(D/d), the symmetric states are unstable to this mode of asymmetric disturbance. It is also shown that when Re > Rec(D/d) the asymmetric states have an asymptotically stable mode of disturbance. The direct numerical simulations are guided by the theoretical approach. In order to improve the numerical simulations, a matching with the asymptotic solution of Moffatt (1964) in the regions around the expansion corners is also included. The dynamics of both small- and large-amplitude disturbances in the flow is described and the transition from symmetric to asymmetric states is demonstrated. The simulations clarify the relationship between the linear stability results and the time-asymptotic behaviour of the flow. The current analyses provide a theoretical foundation for previous experimental and numerical results and shed more light on the transition from symmetric to asymmetric states of a viscous flow in an expanding channel. It is an evolution from a symmetric state, which loses its stability when the Reynolds number of the incoming flow is above Rec(D/d), to a stable asymmetric equilibrium state. The loss of stability is a result of the interaction between the effects of viscous dissipation, the downstream convection of perturbations by the base symmetric flow, and the upstream convection induced by two-dimensional asymmetric disturbances.


Dispersion in laminar flow was studied by a method proposed by Aris based on the use of moments. The first four moments for different types of pulses are calculated in the present paper. The solution was applied to the determination of the diffusivities of hydrogen and nitrobenzoic acid in water, and the results obtained are compared with those found by Taylor’s method of analysis. Some criteria for reliable determination of diffusion coefficients are given. The experiments were performed in coiled tubes, and the influence of the curvature of the tube is discussed


Sign in / Sign up

Export Citation Format

Share Document