Dean Instability in Double-Curved Channels

Total Page:16

File Type:pdf, Size:1020Kb

Dean Instability in Double-Curved Channels Dean Instability in Double-Curved Channels J.-D. Debus,∗ M. Mendoza,† and H. J. Herrmann‡ ETH Z¨urich, Computational Physics for Engineering Materials, Institute for Building Materials, Wolfgang-Pauli-Str. 27, HIT, CH-8093 Z¨urich (Switzerland) We study the Dean instability in curved channels using the lattice Boltzmann model for generalized metrics. For this purpose, we first improve and validate the method by measuring the critical Dean number at the transition from laminar to vortex flow for a streamwise curved rectangular channel, obtaining very good agreement with the literature values. Taking advantage of the easy implementation of arbitrary metrics within our model, we study the fluid flow through a double- curved channel, using ellipsoidal coordinates, and study the transition to vortex flow in dependence of the two perpendicular curvature radii of the channel. We observe not only transitions to 2- cell vortex flow, but also to 4-cell and even 6-cell vortex flow, and we find that the critical Dean number at the transition to 2-cell vortex flow exhibits a minimum when the two curvature radii are approximately equal. PACS numbers: 47.11.-j, 47.20.-k I. INTRODUCTION Numerical simulations of curved channel flow have also been popular, see for example Finlay et al. [1, 5, 6], who Curved channel flow has attracted fluid dynamics re- extended Dean’s stability criterion to channels of differ- searchers for over a hundred years because of its practi- ent aspect ratios η = ri/ro, using high-accuracy linear cal applications in engineering and fascinating features at stability analysis. high Reynolds numbers. At low Reynolds numbers, the All the studies mentioned above deal with a very ide- flow through a curved duct is laminar, and its velocity alized channel geometry, namely with channels having profile resembles the parabolic profile of plane channel a rectangular cross section, which are uniformly curved flow. Following Ref. [1], we will refer to this type of along the streamwise direction only. Thinking of practi- flow as “curved channel Poiseuille flow”. For increasing cal engineering applications (e.g. water ducts), however, Reynolds number, the flow pattern is determined by cen- channels do not possess a perfectly rectangular cross sec- trifugal forces until it becomes unstable and bifurcates to tion, but can also be bent perpendicularly to the stream- secondary flow. In 1928, W. R. Dean showed analytically wise direction for instance under the influence of external for the case of a narrow channel that the fully developed forces such as gravity. The resulting cross section be- flow between two concentric cylinders (see Fig. 1) can comes the section of a circular ring, such that the chan- be characterized by a single non-dimensional parameter, nel possesses a second, cross-sectional curvature. In this the Dean number e, which depends on the Reynolds paper, we study the Dean instability in such a double- number e as wellD as on the channel geometry [2]. It is curved channel, which, to the best of our knowledge, has defined asR never been done before. Introducing a second curvature, we expect qualitative deviations from the behaviour of d the Dean flow as compared to the idealized cylindrical e := e , D R · r case, since the translational symmetry along the cylinder R axis is then broken. where d = r r is the channel width, r , r are the o i i o For the double-curved channel, we study the flow at radii of the inner− and outer cylinder respectively and different Dean numbers and observe a bifurcation from denotes the curvature radius of the inner wall (for curved channel Poiseuille flow to vortex flow, which is aR cylinder, = r ). In particular, Dean performed a i similar to the vortex flow in the cylindrical channel, show- stability analysisR to show that flow between two concen- arXiv:1411.3133v1 [physics.flu-dyn] 12 Nov 2014 ing a pair of counter-rotating vortices. We measure the tric cylinders becomes unstable for small perturbations critical Dean number at the bifurcation point in depen- when the Dean number e exceeds a critical value of D dence of the two perpendicular curvatures. ec 36. Above this value, pairs of counter-rotating streamwise-orientedD ≈ vortices develop, which is known as First, we vary the streamwise curvature radius θ while keeping the cross-sectional curvature radius R Dean vortex flow. The theoretical predictions by Dean Rφ have been confirmed in various experiments, see for ex- fixed. Surprisingly, we find that the lowest Dean num- ample H¨ammerlin [3] or Brewster, Grosberg & Nissan [4]. ber at which the Dean instability occurs corresponds to the configuration in which both curvature radii are ap- proximately equal, θ φ. For channels with both weaker or stronger streamwiseR ≈ R curvature, we measure an ∗ [email protected] (almost linear) increase of the critical Dean number with † [email protected] . For strongly-curved channels, we even observe a sec- Rθ ‡ [email protected] ond bifurcation to vortex flow with 4 counter-rotating 2 vortices. II. LATTICE BOLTZMANN METHOD Secondly, we vary the cross-sectional curvature radius while keeping the streamwise curvature fixed. Again, In this section, we will shortly review the method used φ to simulate the motion of the fluid. For further details weR find a minimum of the critical Dean number for the case , consistent with the first study. For in- we refer to [12, 13]. θ φ The method is based on the Boltzmann equation, creasingR ≈ R, the curvature perpendicular to the stream- φ which describes the motion of fluid particles in terms of wise directionR becomes smaller and smaller, such that the a distribution function f. In a curved three-dimensional double-curved channel approaches a cylindrical geome- Riemann space, the Boltzmann equation reads try in the limit . Correspondingly, for increas- Rφ → ∞ ing φ, the critical Dean number for the double-curved ∂f ∂f ∂f R + ξi + F i = [f] (1) channel approaches the value of the critical Dean number ∂t ∂xi ∂ξi C for the cylindrical channel in our simulations. [12], where f = f(xi,pi) denotes the distribution function, which depends on the space coordinates xi = (x1, x2, x3) as well as on the momentum pi = For the simulations, we use the lattice Boltzmann m(ξ1, ξ2, ξ3), where m and ξi denote the mass and ve- (LB) method, which has been developed during the last locity of the fluid particles respectively, and the mass is decades to simulate fluids by means of simple arithmetic set to m = 1. We are using the Einstein sum convention operations instead of discretizing and solving the com- throughout the whole paper, i.e. Latin indices run over plicated macroscopic equations of continuum fluid me- the spatial directions 1 to 3. The spatial metric enters the chanics. The method itself is based on the Boltzmann Boltzmann equation through the force F i := Γi ξj ξk, kinetic equation, which describes the motion of the mi- − jk which depends on the Christoffel symbols Γi and thus croscopic fluid particles instead of the macrosopic con- jk drives the particles along the geodesics of the curved tinuum. Because of its simplicity and straightforward space. Collisions between fluid particles are accounted parallelizability, the LB method has gained more and for by the collision operator [f], for which we use the more popularity among scientists and engineers in the Bhatnagar-Gross-Krook (BGK)C approximation [15] past. A review about the LB method is given in Ref. [7]. While it was originally designed to solve fluid flows, f f eq [f]= − , the LB method has even been applied to electrodynamics C − τ [8] and magnetohydrodynamics [9] as well as relativistic eq [10] and ultra-relativistic flows [11]. Most of the LB ap- where τ denotes the relaxation parameter and f is the plications use standard Cartesian coordinates (e.g. for Maxwell-Boltzmann equilibrium distribution. The latter fluid flow in rectangular cavities), which cannot be ap- is given by plied straightforwardly to more complex curved geome- eq ρ 1 i i j j tries. However, the lattice Boltzmann method has re- f = exp ξ u gij ξ u 3/2 −2θ − − cently been extended to general metrics being defined by (2πθT ) T a metric tensor [12, 13], which offers a variety of inter- i [12], where θT is the normalized temperature, ρ and u esting new applications. With this extension at hand, it denote the macroscopic density and velocity of the fluid becomes possible to simulate fluids in arbitrary geome- (normalized by the speed of sound cs) and gij are the tries, while standard LB methods are restricted to simple components of the metric tensor g. In the following, we geometries. will always assume the isothermal case θT = 1. In this paper, we improve the method of Ref. [12] by The macroscopic density ρ and velocity ~u are given increasing the accuracy of the forcing term in the LB by the zeroth and first order moment of the distribution equation. The improvement is validated for the case of a function, cylindrical channel, for which we measure the dependence ρ = f √g d3ξ, ρui = f ξi√g d3ξ, (2) of the critical Dean number on the channel aspect ratio Z Z η = ri/ro. Comparing our results to the numerical values given in Ref. [6] by Finlay et al., we find very good where √g d3ξ := √det g dξ1dξ2dξ3 denotes the invariant agreement. volume element of the momentum space. Since conser- vation of mass and momentum are intrinsic features of Since our model can handle arbitrary geometries, we the Boltzmann equation, the macroscopic density ρ and can easily introduce a second perpendicular curvature in velocity ui automatically fulfil the hydrodynamic conser- the channel by choosing ellipsoidal coordinates.
Recommended publications
  • Effect of Dean Number on Heating Transfer Coefficients in an Flat Bottom Agitated Vessel
    IOSR Journal of Engineering May. 2012, Vol. 2(5) pp: 945-951 Effect of Dean Number on Heating Transfer Coefficients in an Flat Bottom Agitated Vessel 1*Ashok Reddy K 2 Bhagvanth Rao M 3 Ram Reddy P 1*..Professor Dept. of Mechanical Engineering ,Prof Rama Reddy College of Engineering & Technology Nandigoan(V) Pantacheru(M) Greater Hyderabad A P 2. Director CVSR College of Engineering Ghatkesar(M) R R Dist A P 3 .Director Malla Reddy College of Engineering Hyderbad A P ABSTRACT Experimental work has been reviewed using non- Newtonian and Newtonian fluids like Sodium Carboxymethal Cellulose with varying concentrations 0.05%,0.1%,0.15% and 0.2% with two different lengths (L=2.82m,L=2.362m) and outer diameter do=6.4mm,inner diameter di=4.0mm in an mixing vessel completely submerged helical coil. In our present study, the effect impeller speed, discharge ,flow behavior index(n) and consistency index (K) of test fluids with dimensionless curvature of the helical coil based on characteristic diameter were presented.. Keywords: Dean Number, Curvature, Heating Rates, Agitated Vessel I. Introduction To improve the performance of heat transfer rates, the fundamental aspect of flow patterns is an essential and important criteria in the helical coil design theory. The helical coil silent features consists of a inner diameter di and mean coil diameter Dm which is also known as pitch circle diameter of the coil. The ratio of coil pipe diameter to the mean coil diameter is known as dimensionless curvature .The ratio of pitch of the helical coil to the developed length of one turn is called non dimensional pitch .The angle which makes with plane geometry perpendicular to the axis of the one turn of the coil dimensions ic defined as helix angle.
    [Show full text]
  • Introduction to Biofluid Mechanics Portonovo S
    To protect the rights of the author(s) and publisher we inform you that this PDF is an uncorrected proof for internal business use only by the author(s), editor(s), reviewer(s), Elsevier and typesetter TNQ Books and Journals Pvt Ltd. It is not allowed to publish this proof online or in print. This proof copy is the copyright property of the publisher and is confidential until formal publication. CHAPTER 16 Introduction to Biofluid Mechanics Portonovo S. Ayyaswamy OUTLINE 16.1 Introduction e2 Exercises e71 16.2 The Circulatory System in the Acknowledgments e72 Human Body e2 Literature Cited e72 16.3 Modeling of Flow in Blood Vessels e18 Supplemental Reading e73 16.4 Introduction to the Fluid Mechanics of Plants e65 CHAPTER OBJECTIVES • To properly introduce the subject of biofluid • To review the parametric impact of the prop- mechanics including the necessary language erties of rigid, flexible, branched, and curved fl • To describe the components of the human cir- tubes on blood ow culation system and document their nominal • To provide an overview of fluid transport in characteristics plants • To present analytical results of relevant models of steady and pulsatile blood flow Fluid Mechanics http://dx.doi.org/10.1016/B978-0-12-405935-1.00016-2 e1 Copyright © 2016 Elsevier Inc. All rights reserved. To protect the rights of the author(s) and publisher we inform you that this PDF is an uncorrected proof for internal business use only by the author(s), editor(s), reviewer(s), Elsevier and typesetter TNQ Books and Journals Pvt Ltd. It is not allowed to publish this proof online or in print.
    [Show full text]
  • Extended Stokes Series for Dean Flow in Weakly Curved Pipes
    EXTENDED STOKES SERIES FOR DEAN FLOW IN WEAKLY CURVED PIPES A thesis presented for the degree of Doctor of Philosophy of Imperial College by Florencia Amalia Tettamanti Department of Mathematics Imperial College 180 Queen’s Gate, London SW7 2BZ June 24, 2012 2 I certify that this thesis, and the research to which it refers, are the product of my own work, and that any ideas or quotations from the work of other people, published or otherwise, are fully acknowledged in accordance with the standard referencing practices of the discipline. Signed: 3 Copyright Copyright in text of this thesis rests with the Author. Copies (by any process) either in full, or of extracts, may be made only in accordance with instructions given by the Author and lodged in the doctorate thesis archive of the college central library. Details may be obtained from the Librarian. This page must form part of any such copies made. Further copies (by any process) of copies made in accordance with such instructions may not be made without the permission (in writing) of the Author. The ownership of any intellectual property rights which may be described in this thesis is vested in Imperial College, subject to any prior agreement to the contrary, and may not be made available for use by third parties without the written permission of the University, which will prescribe the terms and conditions of any such agreement. Further information on the conditions under which disclosures and exploitation may take place is available from the Imperial College registry. 4 5 Abstract This thesis considers steady, fully-developed flows through weakly curved pipes using the extended Stokes series method.
    [Show full text]
  • (English Edition) a Review on the Flow Instability of Nanofluids
    Appl. Math. Mech. -Engl. Ed., 40(9), 1227–1238 (2019) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-019-2521-9 A review on the flow instability of nanofluids∗ Jianzhong LIN†, Hailin YANG Department of Mechanics, State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China (Received Mar. 22, 2019 / Revised May 6, 2019) Abstract Nanofluid flow occurs in extensive applications, and hence has received widespread attention. The transition of nanofluids from laminar to turbulent flow is an important issue because of the differences in pressure drop and heat transfer between laminar and turbulent flow. Nanofluids will become unstable when they depart from the thermal equilibrium or dynamic equilibrium state. This paper conducts a brief review of research on the flow instability of nanofluids, including hydrodynamic instability and thermal instability. Some open questions on the subject are also identified. Key words nanofluid, thermal instability, hydrodynamic instability, review Chinese Library Classification O358, O359 2010 Mathematics Subject Classification 76E09, 76T20 1 Introduction Nanofluids have aroused significant interest over the past few decades for their wide applica- tions in energy, machinery, transportation, and healthcare. For example, low concentration of particles causes viscosity changes[1–2], decreases viscosity with increasing shear rate[3], changes friction factors and pressure drops[4–6], and improves heat transfer[7–9]. Moreover, the unique flow properties of nanofluids are determined by the flow pattern. Heat transfer[10] and pressure drop[11] are also much lower and higher, respectively, in laminar flow than in turbulent flow. By adding particles to the fluid, the thermal entropy generation and friction are of the same order of magnitude as in turbulent flow, while the effect of heat transfer entropy generation strongly outweighs that of the friction entropy generation in laminar flow.
    [Show full text]
  • Blood Flow in Arteries
    November 8, 1996 13:58 Annual Reviews KUTEXT.DOC AR23-12 Annu. Rev. Fluid Mech. 1997. 29:399–434 Copyright c 1997 by Annual Reviews Inc. All rights reserved BLOOD FLOW IN ARTERIES David N. Ku G. W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0405 ABSTRACT Blood flow in arteries is dominated by unsteady flow phenomena. The cardiovas- cular system is an internal flow loop with multiple branches in which a complex liquid circulates. A nondimensional frequency parameter, the Womersley num- ber, governs the relationship between the unsteady and viscous forces. Normal arterial flow is laminar with secondary flows generated at curves and branches. The arteries are living organs that can adapt to and change with the varying hemo- dynamic conditions. In certain circumstances, unusual hemodynamic conditions create an abnormal biological response. Velocity profile skewing can create pock- ets in which the direction of the wall shear stress oscillates. Atherosclerotic dis- ease tends to be localized in these sites and results in a narrowing of the artery lumen—a stenosis. The stenosis can cause turbulence and reduce flow by means of viscous head losses and flow choking. Very high shear stresses near the throat of the stenosis can activate platelets and thereby induce thrombosis, which can totally block blood flow to the heart or brain. Detection and quantification of by CAPES on 04/01/09. For personal use only. stenosis serve as the basis for surgical intervention. In the future, the study of arterial blood flow will lead to the prediction of individual hemodynamic flows in any patient, the development of diagnostic tools to quantify disease, and the design of devices that mimic or alter blood flow.
    [Show full text]
  • Secondary Flow Structures Under Stent-Induced Perturbations for Cardiovascular Flow in a Curved Artery Model
    SECONDARY FLOW STRUCTURES UNDER STENT-INDUCED PERTURBATIONS FOR CARDIOVASCULAR FLOW IN A CURVED ARTERY MODEL Autumn L. Glenn, Kartik V. Bulusu, and Michael W. Plesniak Department of Mechanical and Aerospace Engineering The George Washington University 801 22 nd Street, N.W, Washington, D.C. 20052 [email protected] Fangjun Shu Department of Mechanical and Aerospace Engineering New Mexico State University MSC 3450, P.O. Box 30001, Las Cruces, NM 88003-8001 [email protected] ABSTRACT by the complex geometry and inherent forcing functions, as Secondary flows within curved arteries with unsteady well as changes due to disease. Strong evidence linking forcing are well understood to result from amplified cellular biochemical response to mechanical factors such as centrifugal instabilities under steady-flow conditions and are shear stress on the endothelial cells lining the arterial wall has expected to be driven by the rapid accelerations and received considerable interest (Berger and Jou, 2000; Barakat decelerations inherent in such waveforms. They may also and Lieu, 2003; White and Frangos, 2007; Melchior and affect the function of curved arteries through pro-atherogenic Frangos, 2010). Secondary flow structures may affect the wall shear stresses, platelet residence time and other vascular wall shear stress in arteries, which is known to be closely response mechanisms. related to atherogenesis (Mallubhotla et al., 2001; Evegren et Planar PIV measurements were made under multi-harmonic al., 2010). non-zero-mean and physiological carotid artery waveforms at In curved tubes, secondary flow structures characterized various locations in a rigid bent-pipe curved artery model. by counter-rotating vortex pairs (Dean vortices) are well- Results revealed symmetric counter-rotating vortex pairs that understood to result from amplified centrifugal instabilities developed during the acceleration phases of both multi- under steady flow conditions.
    [Show full text]
  • On Dimensionless Numbers
    chemical engineering research and design 8 6 (2008) 835–868 Contents lists available at ScienceDirect Chemical Engineering Research and Design journal homepage: www.elsevier.com/locate/cherd Review On dimensionless numbers M.C. Ruzicka ∗ Department of Multiphase Reactors, Institute of Chemical Process Fundamentals, Czech Academy of Sciences, Rozvojova 135, 16502 Prague, Czech Republic This contribution is dedicated to Kamil Admiral´ Wichterle, a professor of chemical engineering, who admitted to feel a bit lost in the jungle of the dimensionless numbers, in our seminar at “Za Plıhalovic´ ohradou” abstract The goal is to provide a little review on dimensionless numbers, commonly encountered in chemical engineering. Both their sources are considered: dimensional analysis and scaling of governing equations with boundary con- ditions. The numbers produced by scaling of equation are presented for transport of momentum, heat and mass. Momentum transport is considered in both single-phase and multi-phase flows. The numbers obtained are assigned the physical meaning, and their mutual relations are highlighted. Certain drawbacks of building correlations based on dimensionless numbers are pointed out. © 2008 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved. Keywords: Dimensionless numbers; Dimensional analysis; Scaling of equations; Scaling of boundary conditions; Single-phase flow; Multi-phase flow; Correlations Contents 1. Introduction .................................................................................................................
    [Show full text]
  • Numerical Investigation of the Effect of Curvature and Reynold's Number
    mp Co uta & tio d n e a li Mwangi et al., J Appl Computat Math 2017, 6:3 l p M p a Journal of A t h DOI: 10.4172/2168-9679.1000363 f e o m l a a n t r ISSN: 2168-9679i c u s o J Applied & Computational Mathematics Research Article Open Access Numerical Investigation of the Effect of Curvature and Reynold’s Number to Radial Velocity in a Curved Porous Pipe Mwangi DM1* Karanja S1 and Kimathi M2 1Meru University of Science and Technology, Meru, Kenya 2Technical University of Kenya, Nairobi Kenya, Nairobi, Kenya Abstract Different irrigation methods are being used in agriculture. However, due to scarcity of water, irrigation methods that use water efficiently are needed. The Motivation of this study is the increasing use of porous pipes to meet this requirement. The objective of this study is to investigate the effect of curvature and Reynold’s number on radial velocity profile of water across a porous wall of a curved pipe with circular cross-section, constant permeability, k and porosity, φ. The momentum equations of the two dimensional flow are written in toroidal coordinates. The main flow in the pipe is only characterized by δ and Re as the only non-dimensional groups of numbers. We also considered the flow to be fully developed, unsteady, laminar and irrotational. Darcy law is used to analyse the flow across the porous membrane. The main flow was coupled with the flow through the porous wall of the pipe. The equations were solved using finite difference method.
    [Show full text]
  • The Effect of Entrance Flow Development on Vortex Formation
    The effect of entrance flow development on vortex formation and wall shear stress in a curved artery model Christopher Coxa) and Michael W. Plesniakb) Department of Mechanical and Aerospace Engineering, The George Washington University, Washington, DC 20052, USA We numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three- dimensional vortex pairs and wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e. uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross-sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid 180◦ curved pipe with circular cross-section and constant curvature, neglecting effects of taper, torsion and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses—a potentially significant correlation since cardiovascular disease is known to progress along the inner wall of curved arteries under varying degrees of flow development. Keywords: cardiovascular flows, wall shear stress, vortex dynamics, computational fluid dynamics, discontinuous spectral element I.
    [Show full text]
  • Steady Laminar Flow in a 90° Bend DOI: 10.1177/1687814016669472 Aime.Sagepub.Com
    Research Article Advances in Mechanical Engineering 2016, Vol. 8(9) 1–9 Ó The Author(s) 2016 Steady laminar flow in a 90° bend DOI: 10.1177/1687814016669472 aime.sagepub.com Asterios Pantokratoras Abstract This research undertook an assessment of the characteristics of laminar flow in a 90° bend. The research utilised the computational fluid dynamics software tool ANSYS Fluent, adopting a finite volume method. The study focused on pres- sure decreases, velocity profiles, Dean vortices and Dean cells present in the bend. The results indicated that the Dean vortices first appear slightly at the bend, take a clear form at the exit pipe and disappear further downstream. This article contains also an original chart for determining the Darcy–Weisbach friction factor, while the areas of pipe where Dean vortices and Dean cells formed are presented in a bifurcation diagram. Keywords Curved pipe, bend, laminar, friction factor, Dean vortices Date received: 11 May 2016; accepted: 18 August 2016 Academic Editor: Ishak Hashim Introduction analysis of data concerning declining pressure during turbulent flow in smooth-pipe curves. Pipe bends of An investigation of flow through bends has consider- varying angle were attached to linear piping, to form a able relevance and application to a vast array of situa- consistent entrance and exit flow. Laser-Doppler data tions, given their presence in almost all pipelines. were acquired by Enayet et al.,5 by assessing a 90° pipe Interest has arisen, for instance, concerning the charac- curve’s laminar flow characteristics, where the pipe’s teristics of falls in pressure in the emerging and com- cross-sectional profile had a curvature of d = 0.18.
    [Show full text]
  • Bifurcations in Curved Duct Flow Based on a Simplified Dean Model
    Rigo, Biau & Gloerfelt Bifurcations in curved duct flow based on a simplified Dean model Leonardo Rigo, Damien Biau,a) and Xavier Gloerfelt Dynfluid Laboratory, Ecole´ Nationale Superieure d’Arts et M´etiers. 151 Boulevard de l’Hˆopital, 75013 Paris, France (Dated: 10 April 2020) We present a minimal model of an incompressible flow in square duct subject to a slight curvature. Using a Poincar´e-like section we identify stationary, periodic, aperiodic and chaotic regimes, depending on the unique control parameter of the problem: the Dean number (De). Aside from representing a simple, yet rich, dynamical system the present simplified model is also representative of the full problem, reproducing quite accurately the bifurcation points observed in the literature. We analyse the bifurcation diagram from De = 0 (no curvature) to De = 500, observing a periodic segment followed by two separate chaotic regions. The phase diagram of the flow in the periodic regime shows the presence of two symmetric steady states, the system oscillates around these solutions following a heteroclinic cycle. In the appendix some quantitative results are provided for validation purposes, as well as the python code used for the numerical solution of the Navier-Stokes equations. PACS numbers: Valid PACS appear here Keywords: Navier-Stokes equations, deterministic chaos, curved flow, Dean vortices The flow in a curved square duct is governed by for the flow in a pipe or square duct subject to a slight the Navier-Stokes (NS) equations in polar coor- curvature, based on a truncation of the Navier-Stokes dinates. The dynamics of this kind of flow can (NS) equations written in polar coordinates.
    [Show full text]
  • Steady Streaming Motion in Entrance Region of Curved Tubes During Oscillatory Flow
    Avestia Publishing Journal of Fluid Flow, Heat and Mass Transfer Volume 4, Year 2017 ISSN: 2368-6111 DOI: 10.11159/jffhmt.2017.005 Steady Streaming Motion in Entrance Region of Curved Tubes during Oscillatory Flow Masaru Sumida Kindai University, Faculty of Engineering 1 Takaya Umenobe, Higashi-Hiroshima, 739-2116, Japan [email protected] Abstract - An experimental investigation was carried out on Nomenclature the development of oscillatory laminar flow in curved tubes a radius of tube, ad2 (d: tube diameter) with a turn angle of 180 (i.e., U-shaped tubes) and curvature D Dean number, D =Re Rc-1/2 radius ratios Rc of 4 and 10. Experiments were conducted Q Flow rate under Womersley numbers of 18 to 30 and Dean numbers D R curvature radius of tube axis of 150 and 300, giving (2/D) above 0.655, for which the characteristics of the oscillatory flow in the fully developed Rc curvature radius ratio, Rc=R/a region of the curved tubes are roughly divisible into two. The Re Reynolds number, Re =2awm,o / secondary flow induced in the cross section was visualized by a S secondary flow velocity solid tracer method with a stroboscopic light and was recorded T oscillation period, T 2 ( : angular photographically at several stations upstream and frequency) downstream from the entrance of the curved tube. The u,v, w velocities in x, y and z directions, respecttively development of the secondary flow and transitions between the wm axial velocity averaged over cross section straight-tube and curved-tube oscillatory flows were x, y, z coordinate system, see Figure 3 systematically investigated in detail.
    [Show full text]