On the role of pressure in elasto-inertial turbulence
|
|
- Buddy Franklin
- 6 years ago
- Views:
Transcription
1 Journal of Turbulence ISSN: (Print) (Online) Journal homepage: On the role of pressure in elasto-inertial turbulence Vincent E. Terrapon, Yves Dubief & Julio Soria To cite this article: Vincent E. Terrapon, Yves Dubief & Julio Soria (215) On the role of pressure in elasto-inertial turbulence, Journal of Turbulence, 16:1, 26-43, DOI: 1.18/ To link to this article: The Authors. Published by Taylor & Francis. Published online: 2 Sep 214. Submit your article to this journal Article views: 674 View related articles View Crossmark data Citing articles: 7 View citing articles Full Terms & Conditions of access and use can be found at Download by: [ ] Date: 1 December 217, At: 11:23
2 Journal of Turbulence, 215 Vol. 16, No. 1, 26 43, On the role of pressure in elasto-inertial turbulence Vincent E. Terrapon a, Yves Dubief b and Julio Soria c,d a Department of Aerospace and Mechanical Engineering, University of Liege, Liège, Belgium; b School of Engineering, University of Vermont, Burlington, VT, USA; c Department of Mechanical and Aerospace Engineering, Monash University, Melbourne, Australia; d Department of Aeronautical Engineering, King Abdulaziz University, Jeddah, Kingdom of Saudi Arabia (Received 7 January 214; accepted 3 August 214) Downloaded by [ ] at 11:23 1 December 217 The dynamics of elasto-inertial turbulence is investigated numerically from the perspective of the coupling between polymer dynamics and flow structures. In particular, direct numerical simulations of channel flow with Reynolds numbers ranging from 1 to 6 are used to study the formation and dynamics of elastic instabilities and their effects on the flow. Based on the splitting of the pressure into inertial and polymeric contributions, it is shown that the polymeric pressure is a non-negligible component of the total pressure fluctuations, although the rapid inertial part dominates. Unlike Newtonian flows, the slow inertial part is almost negligible in elasto-inertial turbulence. Statistics on the different terms of the Reynolds stress transport equation also illustrate the energy transfers between polymers and turbulence and the redistributive role of pressure. Finally, the trains of cylindrical structures around sheets of high polymer extension that are characteristics of elasto-inertial turbulence are shown to be correlated with the polymeric pressure fluctuations. Keywords: elasto-inertial turbulence; FENE-P; maximum drag reduction; DNS; pressure split 1. Introduction Polymer additives are known for producing upward of 8% drag reduction in turbulent wall-bounded flows through strong alteration and reduction of turbulent activity [1]. The changes in flow dynamics induced by polymers do not lead to flow relaminarisation but, atmost, to a universal asymptotic state called maximum drag reduction (MDR) [2]. At the same time, polymer additives have also been shown to promote transition to turbulence [3], or even lead to a chaotic flow at very low Reynolds number as in elastic turbulence [4,5]. These seemingly contradictory effects of polymer additives can be explained by the interaction between elastic instabilities and the flow s inertia characterising elasto-inertial turbulence, hereafter referred to as EIT [6,7]. EIT is a state of small-scale turbulence driven by elastic instabilities that exists by either creating its own extensional flow patterns or by exploiting extensional flow topologies. EIT could provide answers to phenomena that current understanding of MDR cannot, such as the absence of a log-law in finite-reynolds numbers MDR flows [8,9], and the phenomenon of early turbulence. Moreover, it supports De Gennes picture [1] that drag reduction derives from two-way energy transfers between Corresponding author. vincent.terrapon@ulg.ac.be C 214 The Authors. Published by Taylor & Francis. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The moral rights of the named authors have been asserted.
3 Journal of Turbulence C f Downloaded by [ ] at 11:23 1 December 217 Figure 1. Skin friction coefficient C f as a function of the Reynolds number Re H for two Weissenberg numbers, Wi H = 4( )andwi H = 3 ( ). Lines indicate correlations for laminar (, C f = 12/Re H ) and turbulent (, C f =.73Re 1/4 H [11]) Newtonian channel flow, and for MDR (, C f =.42Re.55 H [12]); Newtonian solutions are also included ( ). Unlike for Newtonian flows, the skin friction of viscoelastic flows at sufficiently high Weissenberg number departs from the laminar correlation at subcritical Reynolds number and then smoothly transitions to the MDR asymptote. turbulent kinetic energy of the flow and elastic energy of polymers at small scales, resulting in an overall modification of the turbulence energy cascade at high Reynolds numbers. As shown by viscoelastic pipe experiments and direct numerical simulations [6,7], an elastic instability can occur at a Reynolds number smaller than the transition in Newtonian pipe flow if the polymer concentration and Weissenberg number are sufficiently large. Moreover, it is observed that the skin friction coefficient C f = τ w /(1/2ρUb 2 ), where ρ is the density, U b the bulk velocity and τ w the wall shear stress, then follows the characteristic MDR friction law, as shown in Figure 1. Visualisations of numerical simulations indicate that thin sheets of locally high polymer stretch, tilted away from the wall and elongated in the flow direction, create trains of spanwise cylindrical structures of alternating sign (see Figure 2). In low-elasticity flows, EIT is drown out by the canonical Newtonian near-wall vortices. At higher elasticity, the polymer drag reduction mechanism inhibits Newtonian vortices and EIT dominates. It is hypothesised that EIT may be an asymptotic state of parallel wall-bounded flows over a large range of Reynolds numbers. It is suggested that the formation of sheets of polymer stretch results from the unstable nature of the nonlinear advection of low-diffusivity polymers [7]. These sheets, hosting a significant increase in extensional viscosity, create a strong local anisotropy, with a formation of local low-speed jet-like flow. The response of the flow is through pressure, whose role is to redistribute energy across components of momentum, resulting in the formation of waves, or trains of alternating rotational and straining motions. Once triggered, EIT is self-sustained since the elastic instability creates the very velocity fluctuations it feeds upon. The underlying mechanism driving EIT is here further investigated through an analysis of the pressure, and its interaction with topological structures of the flow and polymer stress. The approach relies on the splitting of the pressure into inertial and polymeric contributions [13 15]. Additionally, energy transfers between polymers and turbulence are Re H
4 28 V.E. Terrapon et al. Downloaded by [ ] at 11:23 1 December 217 Figure 2. Contour of the normalised polymer extension C ii /L on selected planes and instantaneous isosurface of the second invariant Q A of the velocity gradient tensor in the lower half of the channel for the subcritical case Re H = 1 and Wi H = 4; Q A =.1 (red) and Q A =.1 (cyan). Trains of mostly spanwise cylindrical Q A structures of alternating sign form around thin sheets of large polymer extension. further analysed through statistics on the different terms of the Reynolds stress transport equation. The paper is organised as follows. Section 2 describes the numerical simulations and the corresponding parameters, and summarises the main theoretical aspects regarding the splitting of the pressure. Results are then presented in Section 3. Finally, the key findings are summarised and discussed in Section Method 2.1. Direct numerical simulations Channel flow simulations are performed in a Cartesian domain, where x, y and z are the streamwise, wall-normal and spanwise directions, respectively. For a polymer solution, the non-dimensional equations governing the flow are the continuity equation, u =, where u = (u, v, w) is the velocity vector, and the momentum equation, t u (u )u = p β 2 u 1 β T. (1) Re H Re H
5 Journal of Turbulence 29 Table 1. Parameters used for the three viscoelastic and the Newtonian cases considered here. For all simulations, the FENE-P model is used with the maximum polymer extension parameter L = 2 and β =.9. Re H a Re h b h c Wi H d Wi h e Wi f DM [%] g Colour h Downloaded by [ ] at 11:23 1 December 217 a Reynolds number based on the bulk velocity U b and the total channel height H. b Reynolds number based on the bulk velocity U b and the channel half-height h (for comparison with previous publications). c Channel half-height h normalised by the friction velocity u 2 τ = ν γ and the viscosity ν. d Weissenberg number based on the solution relaxation time λ and the integral time scale H/U b. e Weissenberg number based on the solution relaxation time λ and the integral time scale h/u b (for comparison with previous publications). f Weissenberg number based on the solution relaxation time λ and the wall-shear rate γ of the corresponding Newtonian case. g Drag modification (reduction or increase) measured as the relative change of the skin friction coefficient with the comparable Newtonian flow friction coefficient (laminar or turbulent). h Colour scheme used throughout the paper for each specific case. The Reynolds number is based on the bulk velocity U b, the full channel height H = 2h and the kinematic viscosity ν of the solution, Re H = U b H/ν. The polymer stress tensor T is computed using the FENE-P model [16]: T = 1 ( ) C Wi H 1 tr(c)/l I 2, (2) where I is the unit tensor and C is the polymer conformation tensor, whose transport equation is t C (u )C = ( u) C C ( u) T T. (3) The properties of the polymer solution are the ratio β of solvent viscosity to the zero-shear viscosity of the solution, the maximum polymer extension L (so that C ii < L 2 ) and the Weissenberg number Wi H = λu b /H based on the solution relaxation time λ and the integral time scale H/U b. The Weissenberg number can also be defined as Wi = λ γ, where γ is the wall shear-rate of the corresponding Newtonian flow at each Reynolds number. 1 Equations (1) (3) are solved using finite differences on a staggered grid and a semi-implicit time advancement scheme [17,18]. After a thorough resolution study, a domain size of 1h 2h 5h with computational nodes was chosen. All results have been verified on domains with a factor 2 in horizontal dimensions and a factor 2 in resolution in each direction. The CFL number was set to.15 to guarantee the boundedness of C. Three different viscoelastic and one Newtonian cases are considered here, as summarised in Table 1. The lower Reynolds number corresponds to a subcritical 2 flow (Re H < Re H,c, where Re H,c = 1719 defines the intersection between the laminar and MDR friction drag lines as shown in Figure 1), while the larger Reynolds number corresponds to a value for which the Newtonian flow is turbulent. The Weissenberg numbers at Re H = 6 are chosen to achieve high drag reduction (HDR) and MDR, respectively. For
6 3 V.E. Terrapon et al. all three viscoelastic cases, L = 2 and β =.9 were used. The corresponding statistics can be found in [7] Inertial and polymeric contributions to pressure In order to investigate the role of pressure in the mechanism underlying EIT, a similar approach to Mansour et al. [13], Kim [14] and Ptasinski et al. [15] is followed. Taking the divergence of the momentum equation (1) leads to a Poisson equation for the pressure: 2 p = 2Q A 1 β ( T), (4) Re Downloaded by [ ] at 11:23 1 December 217 where Q A = ( 1/2) i u j j u i = ( 1/2) i j u i u j is the second invariant of the velocity gradient tensor. In contrast to a Newtonian flow, a second term appears on the right-hand side, which represents the contribution from the polymeric stress. For a periodic channel flow, the pressure satisfies Equation (4) with the boundary conditions p y = y=±h β 2 v Re y 2 1 β y=±h Re 2 T yy y 2 (5) y=±h at the walls and periodicity in x and z. By splitting the right-hand side of Equation (4) into different terms and separating the effect of the wall boundary condition, their respective contributions to the total pressure can be isolated. In particular, we consider here following splitting for the pressure fluctuations p (x) = p(x) P (x): p (x) = p r (x) p s (x) p p (x) p St (x), (6) where fluctuations are denoted by and mean quantities by. The first three contributions on the right-hand side of Equation (6) are solutions of [19]: 2 p r = 2dU v dy x, (7) 2 p s = u i x j u j x i d2 v 2 dy 2, (8) 2 p p = 1 β Re 2 T ij x i x j, (9) with the homogeneous wall boundary conditions p r y = p s y=±h y = p p y=±h y =. (1) y=±h The rapid inertial part, p r, is linear in the velocity fluctuations and represents the immediate response to a change imposed on the mean field, while the slow inertial part, p s, feels this change through nonlinear interactions [14]. In addition to these two inertial terms, the pressure in a viscoelastic flow also has an elastic contribution, p p, originating in the
7 Journal of Turbulence 31 polymeric stress. Finally, the effect of the wall boundary condition is represented by the Stokes pressure, p St, which satisfies with the inhomogeneous boundary conditions p St y = y=±h 2 p St =, (11) β 2 v Re y 2 1 β 2 T yy y=±h Re y 2. (12) y=±h Downloaded by [ ] at 11:23 1 December 217 As the Stokes pressure is typically much smaller than the other contributions, it is not considered here. This pressure split is applied to the simulation results for the four cases mentioned above and statistics are computed in order to identify the relative contributions from inertia and elasticity, as shown below. 3. Results Equations (7) (9) have been solved with homogeneous wall boundary conditions to obtain the three contributions p r, p s and p p. About 25 fields have been collected for each of the four cases considered, on which statistics have been performed as a function of the wall-normal coordinate.the exponent indicates wall units, i.e., a normalisation by the friction velocity u τ = (ν γ ) 1/2 and viscosity ν Statistics Pressure fluctuations The root-mean-square (rms) values of the different pressure contributions and of their sum across the channel height are shown in Figure 3. 3 The most striking observation is that the level of pressure fluctuations at Re H = 6 remains almost constant, or even slightly increases at the wall and at the centre of the channel, when the Weissenberg number is increased, i.e., when the turbulence is damped (Figure 3(a)). This differs markedly from the velocity fluctuations, which strongly decrease with drag reduction (see Figures 1 and 2 in [7]). Pressure fluctuations are, however, much lower in the subcritical case (Re H = 1). The qualitative behaviour is different between Newtonian and viscoelastic flows. In particular, the rms profile is much flatter in the latter case. The splitting of the pressure shows that the marked peak in the pressure fluctuations seen in the buffer layer for the Newtonian case originates in the slow pressure (see Figure 3(b)). Figure 3(b) and 3(c) demonstrates that polymers strongly damp the slow pressure contribution while the rapid part slightly increases, which explains the flatter profiles in the viscoelastic cases. 4 The dominating contribution comes thus from the rapid pressure in contrast to a Newtonian turbulent channel flow, where the slow part dominates almost everywhere [19]. This result differs from the analysis of Ptasinski et al. [15], who found that slow and rapid contributions are more or less equal at MDR. However, their simulation was performed at a much lower Weissenberg number and with a lower maximum extensibility parameter than the present calculations. The elastic pressure, absent in the Newtonian case, has a non-negligible and increasing contribution when the Weissenberg number is increased. If the Weissenberg number is high
8 32 V.E. Terrapon et al rms 1 rms Downloaded by [ ] at 11:23 1 December 217 rms (a) All cases (c) Re H = 6, Wi H =4(HDR) rms (b) Newtonian, Re H = 6 (d) Re H = 6, Wi H =3(MDR) Figure 3. Pressure rms as a function of. (a) Total pressure for the Newtonian case at Re H = 6 ( ), and the viscoelastic cases Re H = 1, Wi H = 4( ); Re H = 6, Wi H = 4 (HDR, ); and Re H = 6, Wi H = 3 (MDR, ). (b) (d) Pressure contributions for selected cases: p r ( );p s ( ); p p ( ); p rs = p r p s ( ); p rsp = p r p s p p ( ). enough, this polymeric contribution is larger than the slow part. This is not the case here at HDR (Figure 3(c)) where both rapid and slow parts are larger than for the MDR case, owing to the stronger turbulence. Figure 4(a) shows the ratio of the rms of the elastic pressure p p to the inertial pressure, p rs = p r p s, indicating that the polymer contribution is about 15% of p rs at Wi H = 4 and increases to about 35% 4% at the larger Weissenberg number. The ratio of the polymeric to inertial pressure contributions also remains quite constant in the near-wall region. The ratio of the rms of the slow pressure p s to the inertial pressure, p rs is shown in Figure 4(b). This ratio also remains rather flat with only a slight increase towards the channel centre. The slow part is only about 6% and 2% of the inertial pressure at Re H = 1 and Re H = 6, respectively. It is interesting to notice that the ratio of the polymeric part to the inertial part seems to depend mostly on the Weissenberg number, while the ratio of the slow part to the inertial
9 Journal of Turbulence p, rms / rs, rms.3.2 s, rms / rs, rms Downloaded by [ ] at 11:23 1 December 217 (a) Figure 4. (a) Ratio of the polymeric p p to the inertial p rs pressure rms as a function of. (b) Ratio of the slow p s to the inertial p rs pressure rms as a function of. Same colour labels as in Figure 3. part appears to depend on the Reynolds number. This is, however, only speculative and more cases need to be analysed. Moreover, the splitting into an inertial and elastic contribution is a simplified view as the polymer stress in Equation (1) can create velocity fluctuations and, thus, indirectly contributes to the inertial pressure and conversely, as discussed below Reynolds stress transport Dubief et al. [7] suggested that the role of pressure is to redistribute turbulent kinetic energy across components of momentum, resulting in the formation of waves, or trains of alternating rotational and straining motions. To illustrate this, various terms in the transport equation for the Reynolds stress u i u j, and, in particular, for the components of the turbulent kinetic energy u α u α (no summation is implied on the subscript α) are further analysed: (b) Du i u j Dt = u i u j t ( ) u i U u j k = u U j i x u k u U i j k x u k k x k }{{} ( ) u i x u j u k β u i u j β k Re x k Re 2 u u i j x k x }{{}}{{ k } ɛ T ij ij ( ) ( ) u p j u p i 1 β u T T ik j u jk i, (13) x i x j Re x k x k }{{}}{{} ij P ij u j f p,i u i f p,j
10 34 V.E. Terrapon et al. where P ij is the production term, T ij the transport term, ɛ ij the dissipation term and f p = T the polymer body force. The last term in Equation (13) represents the creation of Reynolds stress due to polymer stress fluctuations. It can be separated into two contributions: 1 β Re ( u j T ik u i x k ) T jk = 1 β x k Re ( ) u j x T ik u i T jk k ( ) 1 β u T ik j T u i jk. (14) Re x k x k }{{} P p ij Downloaded by [ ] at 11:23 1 December 217 The first term on the right-hand side is in the divergence form and thus corresponds to the transport of Reynolds stress by fluctuating polymer stress. The second term, P p ij,isthe polymer stress work. In particular, P p ii represents the transfer of energy between the mean turbulent kinetic energy and the mean elastic energy of the polymers [7,2]. This term can be either positive or negative, and a positive value corresponds to an energy transfer from the polymers to the turbulence. Similarly, the velocity pressure gradient correlation ij = u p i u p j (15) x j x i can be separated into two contributions, ij = ij d (p) ij, representing the pressure strain and the pressure diffusion, respectively [19]. They are defined as ( u ) ij = p i u j = 2p x j x S ij, (16) i d (p) ij = p u i x j p u j x i. (17) Taking half the trace of Equation (13) leads to the transport equation for the turbulent kinetic energy k. In the case of a fully developed channel flow, the only contribution to P ii is the streamwise component P xx. The pressure diffusion term, d (p) ii, is in the divergence form and thus represents a transport term. On the other hand, ii is deviatoric and thus does not contribute to changes in k, but merely redistributes the turbulent kinetic energy across the components. In particular, in a Newtonian channel flow, turbulent kinetic energy is produced from the mean shear through P xx ; this energy is then redistributed from u x u x to u y u y and u z u z through the terms αα. This well-known picture changes, however, in the viscoelastic case. In particular, the terms Pαα p offer a new path for the production of k, i.e., the transfer of energy from the polymers to the turbulence. Moreover, this term does not solely contribute to u x u x,butto all three components. Figure 5 shows the production terms for the diagonal components of the Reynolds stress tensor. Compared to a Newtonian flow, the production term P αα is much lower in a drag-reduced flow (see Figure 5(a) and 5(b)). Nevertheless, the most interesting aspect is
11 Journal of Turbulence P xx.3.2 p P xx, P xx.2.1 Downloaded by [ ] at 11:23 1 December 217 (a) Streamwise, Newtonian p, P zz P p yy (c) Wall-normal and spanwise, viscoelastic -.2 (b) Streamwise, viscoelastic Figure 5. (a) and (b) Streamwise production terms P xx ( )andpp xx ( ) as a function of. The polymer contribution is positive, indicating a transfer of energy from the polymers to the turbulence. (c) Wall-normal and spanwise production terms Pyy p ( )andpp zz ( ) as a function of. In this case, energy is transferred from the turbulence to the polymers, but the levels are much lower than in the streamwise direction. Same colour labels as in Figure 3. that the streamwise transfer of energy from the polymers to the turbulence, P p xx, dominates. In other words, the turbulence is mainly sustained by the polymers, as already recognised in previous studies [7,18,2]. On the other hand, the polymer production for the two other components of the diagonal is negative (see Figures 5(c)), indicating a transfer of energy from the flow to the polymers. The transfer rates are, however, 1 times lower than for the streamwise component. It is also interesting to see that P p xx has almost the same level for both HDR and MDR, while all other production terms are much lower at MDR. The analysis above indicates that energy is transferred from the polymers to the flow through the streamwise component, and that a small portion of this energy goes back to the polymers through the wall-normal and spanwise components. The redistribution of energy from u x u x to u y u y and u z u z is achieved through the pressure strain term ij. The different components of the diagonal, αα,areshowninfigures 6 8. Their qualitative behaviour is
12 36 V.E. Terrapon et al xx xx Downloaded by [ ] at 11:23 1 December 217 xx (a) Newtonian, Re H = (c) Re H = 6, Wi H =4(HDR) xx (b) Re H =1,Wi H = (d) Re H =6,Wi H =3(MDR) Figure 6. Streamwise pressure strain component xx as a function of obtained from the different pressure contributions: p r ( );p s ( ); p p ( ); p rs = p r p s ( ); p rsp = p r p s p p ( ). Same colour labels as in Figure 3. xx is negative in most of the domain indicating an energy transfer away from the streamwise component. very similar to a Newtonian flow at a similar Reynolds number [19], but with lower levels owing to the turbulence reduction. One can observe that xx is negative while yy and zz are positive for 1, which demonstrates the redistribution of energy from the streamwise component to both the wall-normal and spanwise components. Note, however, that xx is negative and yy positive very close to the wall. Similar conclusions can be drawn from ij and u j f p,i u i f p,j (not shown here). Unlike for Newtonian flows, the largest contribution to the pressure-strain term comes from the rapid part, and a non-negligible contribution from the polymer part is observed. The polymeric contribution is comparatively larger for yy and lower for zz, which tends to indicate that the elastic contribution slightly favours a two-dimensional flow, the threedimensionality being mostly driven by the inertial part. On the other hand, the slow part is very weak, except for yy at HDR in the centre of the channel. Interestingly, the behaviour of the rapid contribution to yy differs between the Newtonian and the viscoelastic cases (see
13 Journal of Turbulence yy yy (a) Newtonian, Re H = (b) Re H = 1, Wi H =4 Downloaded by [ ] at 11:23 1 December 217 yy (c) Re H = 6, Wi H =4(HDR) yy (d) Re H = 6, Wi H =3(MDR) Figure 7. Wall-normal pressure strain component yy as a function of obtained from different pressure contributions. Same line style as in Figure 6 and same colour labels as in Figure 3. yy is positive in most of the domain indicating an energy transfer to the wall-normal component. Figure 7). In the Newtonian case, it is positive in the near-wall region and then becomes negative for 1; in the viscoelastic case, it is exactly the opposite. Actually, the rapid contribution in the viscoelastic case behaves exactly like the slow contribution in the Newtonian case. The Reynolds shear stress production P xy is shown in Figure 9. The same qualitative behaviour is observed for Newtonian and viscoelastic flows, but with a lower magnitude in drag-reduced flows due to the reduction of turbulence intensity by polymers. P xy is negative, indicating that Reynolds shear stress is produced. However, this production of Reynolds shear stress is compensated by a positive polymer contribution through Pxy.The p lower magnitude of the production term P xy and this cancelling effect between P xy and Pxy p explains the very low levels of Reynolds shear stress seen in MDR flows (see Figure 2 of [7]). A difference between Newtonian and viscoelastic flows is also observed for the x y component of the velocity pressure gradient correlation. In particular, xy in the viscoelastic cases is negative (i.e., production of Reynolds shear stress) around 1, which is not observed in the Newtonian case (see Figure 1). This behaviour is directly caused by the polymer pressure contribution p p.
14 38 V.E. Terrapon et al zz.2 zz (a) Newtonian, Re H = (b) Re H = 1, Wi H =4 Downloaded by [ ] at 11:23 1 December 217 zz (c) Re H = 6, Wi H =3(HDR) zz (d) Re H = 6, Wi H =3(MDR) Figure 8. Spanwise pressure strain component zz as a function of obtained from the different pressure contributions. Same line style as in Figure 6 and same colour labels as in Figure 3. zz is positive everywhere indicating an energy transfer to the spanwise component P xy -.4 P xy, P p xy (a) Newtonian, Re H =6 -.1 (b) Viscoelastic Figure 9. Production terms P x ( )andpp xy ( ) of Reynolds shear stress u x u y as a function of. Same colour labels as in Figure 3. Both terms cancel each other in the viscoelastic cases, explaining the very low levels of Reynolds shear stress observed.
15 Journal of Turbulence xy.3.2 xy Downloaded by [ ] at 11:23 1 December 217 xy (a) Newtonian, Re H = (c) Re H = 6, Wi H =4(HDR) -.5 xy (b) Re H = 1, Wi H =4 -.5 (d) Re H = 6, Wi H =3(MDR) Figure 1. Pressure strain component xy as a function of obtained from the different pressure contributions: p r ( - );p s ( ); p p ( ); p rs = p r p s ( ); p rsp = p r p s p p ( ). Same colour labels as in Figure Instantaneous fields As shown in Figure 2, the flow is characterised by trains of cylindrical Q A structures of alternating sign around thin sheets of large polymer extension. This is further illustrated in Figure 11, which shows the contour of the polymer stress component T xx and isolines of Q A in an x y plane for the lower Reynolds number case. The long thin sheets of large polymer extension surrounded by cylindrical structures are clearly visible. A small region of the x y plane around one of these sheets (indicated by the dashed box in Figure 11)isshowninFigures 12 and 13. The first figure shows the instantaneous contour of the polymeric pressure contribution p p and isolines of Q A, and illustrates the strong correlation between these two quantities. The largest fluctuations of p p are of alternating sign and located on each side of the sheet, mostly in between cylindrical Q A structures. The wavelength of these structures is the same as the wavelength of the Q A structures. On the other hand, the inertial pressure p rs does not show such a clear correlation. Large fluctuations of p rs are rather correlated with the centre of strong Q A structures, as shown in
16 4 V.E. Terrapon et al. Figure 11. Instantaneous contour of the polymer stress component T xx and isolines of Q A (dash lines represent negative values) in an x y plane for Re H = 1 and Wi H = 4. The flow is from left to right and the dashed box represents the region plotted in the next two figures. Downloaded by [ ] at 11:23 1 December 217 Figure 12. Instantaneous contour of the polymeric pressure contribution p p and isolines of Q A (dash lines represent negative values) in an x y plane for Re H = 1 and Wi H = 4. The region plotted corresponds to the dashed box in Figure 11. A clear correlation between p p and Q A is visible. Figure 13. Instantaneous contour of the inertial pressure contribution p rs and isolines of Q A (dash lines represent negative values) in an x y plane for Re H = 1 and Wi H = 4. The region plotted corresponds to the dashed box in Figure 11. Unlike the previous figure, no clear correlation between p rs and Q A is visible.
17 Journal of Turbulence 41 Downloaded by [ ] at 11:23 1 December 217 Figure 14. Schematics of the typical structures observed around thin sheets of large polymer extension (black line) in the near-wall region: second invariant of the velocity gradient tensor Q A (blue), pressure p (red) and polymer body force f p (green). Dotted lines indicate a negative value. Figure 13. In particular, as in inertial turbulence, positive values of Q A correspond to low pressure regions, and thus, to vortical regions [21]. The analysis of the polymer body force f p = T (not shown here) indicates that the polymer body force is mostly parallel to the sheet with opposite sign on each side. Additionally, f p also alternates direction along the sheet with the same wavelength as the Q A structures and opposite sign, indicating that the polymers are most likely the driving force that creates these structures. The overall picture gained from this analysis is schematically summarised in Figure 14. It is believed that the combined effect of advection at low diffusivity and existing flow perturbations leads to the formation of sheets of high polymer extension and, in turn, to a large increase in the extensional viscosity, and thus in the polymer stress. Small perturbations of the sheets cause the polymer body force to alternate direction, and thereby, create these cylindrical Q A structures. The pressure adapts to ensure zero divergence of the velocity and redistribute part of the turbulent kinetic energy from the streamwise to the other components. Once triggered, this process appears to be self-sustained, at least over the hundreds of flow through-time simulated here. Therefore, it can be conjectured that these characteristic trains of cylindrical structures are mostly driven by the polymers. Nonetheless, an indirect effect from inertia probably also contributes to the dynamics. Finally, note that, since the present analysis is based on correlations, these causal connections are only speculative. Further analysis is still required to ascertain this. 4. Conclusion and future work A new state of small-scale turbulence, elasto-inertial turbulence, has been recently discovered in polymeric solutions. EIT has not only provided answers to previously unexplained phenomena such as early turbulence [6,22], it is also a unique window into the creation of turbulence by means other than those known in Newtonian turbulence. The current understanding of EIT is that it is driven by both inertial and elastic instabilities, leading to strong alterations of the flow dynamics compared to the Newtonian case. In particular, the flow departs from the laminar solution at a Reynolds number lower than the transitional Reynolds number observed for a Newtonian fluid. The topology of coherent structures identified by the second invariant Q A of the velocity gradient tensor is dramatically different
18 42 V.E. Terrapon et al. Downloaded by [ ] at 11:23 1 December 217 than the coherent structures observed in Newtonian near-wall turbulence. Interestingly, our simulations show that EIT coherent structures extend to larger Reynolds numbers where the skin friction reaches the MDR asymptotic law. The spatial organisation of regions of highly stretch polymers is also found to be similar between Re = 1 and 6. In this context, the present analysis aimed at better characterising the dynamics of EIT and, in particular, to study the role of pressure in the relationship between the Q A structures and the organisation of polymer stretch. The splitting of the pressure into different components, in particular, into an elastic and an inertial part, has provided a tool to assess the relative contributions of both components to the overall dynamics. It is shown that the elastic, or polymer, pressure is a non-negligible component of the total pressure fluctuations, although the rapid part dominates. Unlike Newtonian flows, the slow part is much lower in elasto-inertial turbulence. Statistics also demonstrate the redistributive role of pressure, transferring energy from the streamwise to the other components. Finally, a schematic description of the typical structures encountered in EIT is proposed. It is postulated that those small-scale structures, associated with thin sheets of large polymer extension, are directly driven by the polymers. Nonetheless, an indirect inertial contribution is still possible, and most likely required for self-sustaining dynamics. In a broader context, elasto-inertial turbulence offers a new perspective on the interaction between polymers and flow. First, its existence at Re = 1, for which the Newtonian flow is laminar, clearly establishes a pathway of energy from polymers to flow, lending support to De Gennes theory of energy transfers between polymers and flow [1]. Then, the existence of EIT at larger Reynolds number raises some questions as to its role in MDR and the nature of MDR, which will be addressed in future publications. Many questions remain unanswered, such as the mechanism by which small perturbations of the sheet lead to a body force of alternating sign, the role of those elastic instabilities, the validity of the above conclusions at larger Reynolds numbers, the two-dimensional nature of these instabilities. Finally, the possible existence of EIT in other flows, specifically non-wall-bounded flows, remains a speculation based on our current definition: EIT occurs in parallel flows with a mean shear. Consequently, one may anticipate that EIT may live in shear flows (mixing layers, jets) or in a Kolmogorov flow. Further analysis is thus required to obtain a definite description of the physical mechanisms underlying EIT, which will be part of future work. Funding V.E. Terrapon acknowledges the financial support of a Marie Curie FP7 Career Integration Grant [grant number PCIG1-GA ) and of a grant from the Fonds Spéciaux pour la Recherche of the University of Liege (Crédits classiques 213) [grant number C-13/19]. The Vermont Advanced Computing Center is gratefully acknowledged for providing the computing resources. Y. Dubief acknowledges the support from the National Institutes of Health [grant number P1HL4673 (project 1]. J. Soria acknowledges the support of the Australian Research Council [grant number DP19562]. This work was initiated during the Center for Turbulence Research 212 Summer Program. Notes 1. In laminar flows, γ = 6U b /H,sothatWi = 6Wi H. 2. The Newtonian flow is laminar. 3. Note that p rsp 2 = (p r p s p p )2 p r 2 p s 2 p p The subcritical case (not shown here) is qualitatively similar to the MDR case but with lower values.
19 Journal of Turbulence 43 Downloaded by [ ] at 11:23 1 December 217 References [1] C.M. White and M.G. Mungal, Mechanics and prediction of turbulent drag reduction with polymer additives, Annu. Rev. Fluid Mech. 4 (28), pp [2] P.S. Virk, H.S. Mickley, and K.A. Smith, The ultimate asymptote and mean flow structure in Toms phenomenon, J. Appl. Mech. 37 (197), pp [3] J.W. Hoyt, Laminar-turbulent transition in polymer solutions, Nature 27 (1977), pp [4] A. Groisman and V. Steinberg, Elastic turbulence in a polymer solution flow, Nature 45 (2), pp [5] V. Steinberg, Elastic stresses in random flow of a dilute polymer solution and the turbulent drag reduction problem, C. R. Phys. 1 (29), pp [6] D. Samanta, Y. Dubief, M. Holzner, C. Schšfer, A.N. Morozov, C. Wagner, and B. Hof, Elasto-inertial turbulence, Proc. Natl. Acad. Sci. 11 (213), pp [7] Y. Dubief, V.E. Terrapon, and J. Soria, On the mechanism of elasto-inertial turbulence, Phys. Fluids 25 (213), [8] C.M. White, Y. Dubief, and J. Klewicki, Re-examining the logarithmic dependence of the mean velocity distribution in polymer drag reduced wall-bounded flow, Phys. Fluids 24 (212), [9] B.R. Elbing, M. Perlin, D.R. Dowling, and S.L. Ceccio, Modification of the mean near-wall velocity profile of a high-reynolds number turbulent boundary layer with the injection of drag-reducing polymer solutions, Phys. Fluids 25 (213), [1] P.G. De Gennes, Introduction to Polymer Dynamics, Cambridge University Press, Cambridge, 199. [11] R.D. Dean, Reynolds number dependence of skin friction and other bulk flow variables in two-dimensional rectangular duct flow, J. Fluids Eng. 1 (1978), pp [12] P.S. Virk, E.W. Merrill, H.S. Mickley, K.A. Smith, and E.L. Mollo-Christensen, The Toms phenomenon: Turbulent pipe flow of dilute polymer solutions, J. Fluid Mech. 3 (1967), pp [13] N.N. Mansour, J. Kim, and P. Moin, Reynolds-stress and dissipation-rate budgets in a turbulent channel flow, J. Fluid Mech. 194 (1988), pp [14] J. Kim, On the structure of pressure fluctuations in simulated turbulent channel flow, J. Fluid Mech. 25 (1989), pp [15] P.K. Ptasinski, B.J. Boersma, F.T.M. Nieuwstadt, M.A. Hulsen, B.H.A.A. van den Brule, and J.C.R. Hunt, Turbulent channel flow near maximum drag reduction: Simulations, experiments and mechanisms, J. Fluid Mech. 49 (23), pp [16] R.B. Bird, R.C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids. Vol. 2: Kinetic Theory, Wiley-Interscience, [17] Y. Dubief, C.M. White, V.E. Terrapon, E.S.G. Shaqfeh, P. Moin, and S.K. Lele, On the coherent drag-reducing and turbulence-enhancing behaviour of polymers in wall flows, J. Fluid Mech. 514 (24), pp [18] Y. Dubief, V.E. Terrapon, C.M. White, E.S.G. Shaqfeh, P. Moin, and S.K. Lele, New answers on the interaction between polymers and vortices in turbulent flows, Flow Turbul. Combust. 74 (25), pp [19] G.A. Gerolymos, D. Sénéchal, and I. Vallet, Wall effects on pressure fluctuations in turbulent channel flow, J. Fluid Mech. 72 (213), pp [2] V. Dallas, J.C. Vassilicos, and G.F. Hewitt, Strong polymer-turbulence interactions in viscoelastic turbulent channel flow, Phys. Rev. E 82 (21), [21] Y. Dubief and F. Delcayre, On coherent-vortex identification in turbulence, J. Turbul. 1 (2) pp [22] Y. Dubief, C.M. White, E.S.G. Shaqfeh, and V.E. Terrapon, Polymer Maximum Drag Reduction: A Unique Transitional State, Center for Turbulence Research Annual Research Briefs, Stanford, CA, 21.
Events of High Polymer Activity in Drag Reducing Flows
Flow Turbulence Combust (27) 79:23 32 DOI.7/s494-7-975-5 Events of High Polymer Activity in Drag Reducing Flows C. Tesauro B. J. Boersma M. A. Hulsen P. K. Ptasinski F. T. M. Nieuwstadt Received: May 26
More informationPolymer maximum drag reduction: A unique transitional state
Center for Turbulence Research Annual Research Briefs 47 Polymer maximum drag reduction: A unique transitional state By Y. Dubief, C. M. White, E. S. G. Shaqfeh AND V. E. Terrapon. Motivation and objectives
More informationSOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG
SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG L. Thais Université de Lille Nord de France, USTL F9 Lille, France Laboratoire de Mécanique de Lille CNRS, UMR 817
More informationModel-based analysis of polymer drag reduction in a turbulent channel flow
2013 American Control Conference ACC Washington, DC, USA, June 17-19, 2013 Model-based analysis of polymer drag reduction in a turbulent channel flow Binh K. Lieu and Mihailo R. Jovanović Abstract We develop
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationINTENSE FOCAL AND REYNOLDS STRESS STRUCTURES OF A SELF-SIMILAR ADVERSE PRESSURE GRADIENT TURBULENT BOUNDARY LAYER
INTENSE FOCAL AND REYNOLDS STRESS STRUCTURES OF A SELF-SIMILAR ADVERSE PRESSURE GRADIENT TURBULENT BOUNDARY LAYER Atsushi Sekimoto, Vassili Kitsios, Callum Atkinson Laboratory for Turbulence Research in
More informationMelbourne, Victoria, 3010, Australia. To link to this article:
This article was downloaded by: [The University Of Melbourne Libraries] On: 18 October 2012, At: 21:54 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationCentro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal
DEVELOPING CLOSURES FOR TURBULENT FLOW OF VISCOELASTIC FENE-P FLUIDS F. T. Pinho Centro de Estudos de Fenómenos de Transporte, FEUP & Universidade do Minho, Portugal C. F. Li Dep. Energy, Environmental
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationarxiv: v2 [physics.flu-dyn] 20 Jun 2017
Exceeding the asymptotic limit of polymer drag reduction George H. Choueiri, Jose M. Lopez, and Björn Hof Institute of Science and Technology Austria, 3400 Klosterneuburg, Austria (Dated: June 21, 2017)
More informationEstimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows
Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 140 ISSN: 2369-8136 DOI: 10.11159/htff17.140 Estimation of Turbulent
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationHomogeneous isotropic turbulence in dilute polymers
J. Fluid Mech. (5), vol. 531, pp. 1 1. c 5 Cambridge University Press doi:1.117/s1153666 Printed in the United Kingdom 1 Homogeneous isotropic turbulence in dilute polymers By E. DE ANGELIS, 1 C. M. C
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
More informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationRole of polymers in the mixing of Rayleigh-Taylor turbulence
Physics Department University of Genova Italy Role of polymers in the mixing of Rayleigh-Taylor turbulence Andrea Mazzino andrea.mazzino@unige.it Guido Boffetta: University of Torino (Italy) Stefano Musacchio:
More informationContribution of Reynolds stress distribution to the skin friction in wall-bounded flows
Published in Phys. Fluids 14, L73-L76 (22). Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Koji Fukagata, Kaoru Iwamoto, and Nobuhide Kasagi Department of Mechanical
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationSkin-Friction Drag Reduction by Dilute Polymer Solutions in Turbulent Channel Flow
Skin-Friction Drag Reduction by Dilute Polymer Solutions in Turbulent Channel Flow by Dong-Hyun Lee A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationMaximum drag reduction asymptotes and the cross-over to the Newtonian plug
J. Fluid Mech. (26), vol. 551, pp. 185 195. c 26 Cambridge University Press doi:1.117/s221125795 Printed in the United Kingdom 185 Maximum drag reduction asymptotes and the cross-over to the Newtonian
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationWhat is Turbulence? Fabian Waleffe. Depts of Mathematics and Engineering Physics University of Wisconsin, Madison
What is Turbulence? Fabian Waleffe Depts of Mathematics and Engineering Physics University of Wisconsin, Madison it s all around,... and inside us! Leonardo da Vinci (c. 1500) River flow, pipe flow, flow
More informationInterpretation of the mechanism associated with turbulent drag reduction in terms of anisotropy invariants
J. Fluid Mech. (27), vol. 577, pp. 457 466. c 27 Cambridge University Press doi:1.117/s221127583 Printed in the United Kingdom 457 Interpretation of the mechanism associated with turbulent drag reduction
More informationSECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES
SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES Yosuke Hasegawa Institute of Industrial Science The University of Tokyo Komaba 4-6-1, Meguro-ku, Tokyo 153-8505, Japan ysk@iis.u-tokyo.ac.jp
More informationarxiv: v1 [physics.flu-dyn] 16 Nov 2018
Turbulence collapses at a threshold particle loading in a dilute particle-gas suspension. V. Kumaran, 1 P. Muramalla, 2 A. Tyagi, 1 and P. S. Goswami 2 arxiv:1811.06694v1 [physics.flu-dyn] 16 Nov 2018
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationApplication of wall forcing methods in a turbulent channel flow using Incompact3d
Application of wall forcing methods in a turbulent channel flow using Incompact3d S. Khosh Aghdam Department of Mechanical Engineering - University of Sheffield 1 Flow control 2 Drag reduction 3 Maths
More informationRegularity diagnostics applied to a turbulent boundary layer
Center for Turbulence Research Proceedings of the Summer Program 208 247 Regularity diagnostics applied to a turbulent boundary layer By H. J. Bae, J. D. Gibbon, R. M. Kerr AND A. Lozano-Durán Regularity
More informationGeneration of initial fields for channel flow investigation
Generation of initial fields for channel flow investigation Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (Mai 2) In the framework of the DFG-funded research
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More informationModifications of the V2 Model for Computing the Flow in a 3D Wall Jet Davidson, L.; Nielsen, Peter Vilhelm; Sveningsson, A.
Aalborg Universitet Modifications of the V2 Model for Computing the Flow in a D Wall Jet Davidson, L.; Nielsen, Peter Vilhelm; Sveningsson, A. Published in: Proceedings of the International Symposium on
More informationCOMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES
COMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES Somnath Ghosh, Rainer Friedrich Lehrstuhl für Aerodynamik, Technische Universität München Boltzmannstrasse 15, D-85748 Garching,
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
More informationCoherent structures in stably stratified plane Couette flow
Coherent structures in stably stratified plane Couette flow D. Olvera * & R. R. Kerswell School of Mathematics, University of Bristol, Bristol, UK. * do2542@bristol.ac.uk Abstract A large body of recent
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationReynolds stress budgets in Couette and boundary layer flows
Reynolds stress budgets in Couette and boundary layer flows By Jukka Komminaho and Martin Skote Dept. of Mechanics, KTH, SE-1 44, Stockholm, Sweden Reynolds stress budgets for both Couette and boundary
More informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationDirect numerical simulation (DNS) investigation of turbulent open channel flow of a Herschel Bulkley fluid
Paste 211 R.J. Jewell and A.B. Fourie (eds) 211 Australian Centre for Geomechanics, Perth, ISBN 978--986154-3-2 https://papers.acg.uwa.edu.au/p/114_38_guang/ Direct numerical simulation (DNS) investigation
More informationROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 ROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS Myoungkyu Lee Robert D. Moser
More informationHEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS
HEAT TRANSFER OF SIMPLIFIED PHAN-THIEN TANNER FLUIDS IN PIPES AND CHANNELS Paulo J. Oliveira Departamento de Engenharia Electromecânica, Universidade da Beira Interior Rua Marquês D'Ávila e Bolama, 600
More informationIntermittency in the Transition to Turbulence for a Shear-Thinning Fluid in Hagen-Poiseuille Flow
Journal of Applied Fluid Mechanics, Vol. 7, No. 1, pp. 1-6, 214. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Intermittency in the Transition to Turbulence for a Shear-Thinning
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationShear Turbulence. Fabian Waleffe. Depts. of Mathematics and Engineering Physics. Wisconsin
Shear Turbulence Fabian Waleffe Depts. of Mathematics and Engineering Physics, Madison University of Wisconsin Mini-Symposium on Subcritical Flow Instability for training of early stage researchers Mini-Symposium
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationThe interaction of vorticity and rate-of-strain in homogeneous sheared turbulence
PHYSICS OF FLUIDS VOLUME 12, NUMBER 4 APRIL 2000 The interaction of vorticity and rate-of-strain in homogeneous sheared turbulence K. K. Nomura a) and P. J. Diamessis Department of Mechanical and Aerospace
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationFLOW-NORDITA Spring School on Turbulent Boundary Layers1
Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on
More informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationOPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW
OPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW Shingo Motoki, Genta Kawahara and Masaki Shimizu Graduate School of Engineering Science Osaka University 1-3 Machikaneyama, Toyonaka, Osaka
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationO. A Survey of Critical Experiments
O. A Survey of Critical Experiments 1 (A) Visualizations of Turbulent Flow Figure 1: Van Dyke, Album of Fluid Motion #152. Generation of turbulence by a grid. Smoke wires show a uniform laminar stream
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
More informationA new numerical framework to simulate viscoelastic free-surface flows with the finitevolume
Journal of Physics: Conference Series PAPER OPEN ACCESS A new numerical framework to simulate viscoelastic free-surface flows with the finitevolume method Related content - Gravitational collapse and topology
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationOutline: Phenomenology of Wall Bounded Turbulence: Minimal Models First Story: Large Re Neutrally-Stratified (NS) Channel Flow.
Victor S. L vov: Israeli-Italian March 2007 Meeting Phenomenology of Wall Bounded Turbulence: Minimal Models In collaboration with : Anna Pomyalov, Itamar Procaccia, Oleksii Rudenko (WIS), Said Elgobashi
More informationCONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
Second International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 6-8 December 1999 CONTRIBUTION TO EXTRUDATE SWELL FROM THE VELOCITY FACTOR IN NON- ISOTHERMAL EXTRUSION
More informationvery elongated in streamwise direction
Linear analyses: Input-output vs. Stability Has much of the phenomenology of boundary la ut analysis of Linearized Navier-Stokes (Cont.) AMPLIFICATION: 1 y. z x 1 STABILITY: Transitional v = (bypass) T
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationApplied Computational Fluid Dynamics
Lecture 9 - Kolmogorov s Theory Applied Computational Fluid Dynamics Instructor: André Bakker André Bakker (2002-2005) Fluent Inc. (2002) 1 Eddy size Kolmogorov s theory describes how energy is transferred
More informationApplied Mathematics and Mechanics (English Edition)
Appl. Math. Mech. -Engl. Ed., 39(9), 1267 1276 (2018) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-018-2364-7 Direct numerical simulation of turbulent flows through
More informationNonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto
Nonequilibrium Dynamics in Astrophysics and Material Science 2011-11-02 @ YITP, Kyoto Multi-scale coherent structures and their role in the Richardson cascade of turbulence Susumu Goto (Okayama Univ.)
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationEXACT SOLUTIONS TO THE NAVIER-STOKES EQUATION FOR AN INCOMPRESSIBLE FLOW FROM THE INTERPRETATION OF THE SCHRÖDINGER WAVE FUNCTION
EXACT SOLUTIONS TO THE NAVIER-STOKES EQUATION FOR AN INCOMPRESSIBLE FLOW FROM THE INTERPRETATION OF THE SCHRÖDINGER WAVE FUNCTION Vladimir V. KULISH & José L. LAGE School of Mechanical & Aerospace Engineering,
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationTHE SUBORDINATION OF THE THREE- DIMENSIONAL FLOW INSTALLATION IN THE CONVERGING CHANNEL ON RHEOLOGICAL CHARACTERISTICS OF POLYMER STREAM
International Journal of Civil Engineering and Technology (IJCIET) Volume 9, Issue 13, December 2018, pp. 949-956, Article ID: IJCIET_09_13_095 Available online at http://www.iaeme.com/ijciet/issues.asp?jtype=ijciet&vtype=9&itype=13
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationAdvances in the Analysis and Prediction of Turbulent Viscoelastic Flows
Old Dominion University ODU Digital Commons CCPO Publications Center for Coastal Physical Oceanography 2014 Advances in the Analysis and Prediction of Turbulent Viscoelastic Flows T. B. Gatski Old Dominion
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationFlow Transition in Plane Couette Flow
Flow Transition in Plane Couette Flow Hua-Shu Dou 1,, Boo Cheong Khoo, and Khoon Seng Yeo 1 Temasek Laboratories, National University of Singapore, Singapore 11960 Fluid Mechanics Division, Department
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationVictor S. L vov, Baku, January 21, Turbulent Drag Reduction by Dilute Solution of Polymers
Victor S. L vov, Baku, January 21, 2009 Turbulent Drag Reduction by Dilute Solution of Polymers In collaboration with T.S. Lo (TL), A. Pomyalov (AP), I. Procaccia (IP), V. Tiberkevich (VT). Weizmann Inst.,
More informationBuoyancy Fluxes in a Stratified Fluid
27 Buoyancy Fluxes in a Stratified Fluid G. N. Ivey, J. Imberger and J. R. Koseff Abstract Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More information