Nonmodal amplification of disturbances as a route to elasticity-induced turbulence
|
|
- Tabitha West
- 5 years ago
- Views:
Transcription
1 Under consideration for publication in J. Fluid Mech. 1 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence By M I H A I L O R. J O V A N O V I Ć 1 A N D S A T I S H K U M A R 2 1 Department of Electrical and Computer Engineering, University of Minnesota, Minneapolis, MN 55455, USA 2 Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN 55455, USA (Received 3 June 2009) Nonmodal amplification of stochastic disturbances in elasticity-dominated channel flows of Oldroyd-B fluids is analyzed in this work. For streamwise-constant flows with high elasticity numbers µ and finite Reynolds numbers Re, we show that the linearized dynamics can be decomposed into slow and fast subsystems, and establish analytically that the steady-state variance amplification scales as O(µRe 3 ). This demonstrates that large velocity variances can be sustained in weakly inertial/strongly elastic stochastically driven channel flows. The underlying physical mechanism involves polymer stretching that introduces an effective lift-up of flow fluctuations similar to vortex-tilting in inertiadominated flows. The mechanism examined here provides a possible route for a bypass transition to elasticity-induced turbulence and might be exploited to enhance mixing in microfluidic devices. 1. Introduction The classical approach to transition to turbulence examines the linearized governing equations for exponentially growing normal modes. The existence of these unstable modes implies exponential growth of infinitesimal perturbations to the laminar flow, and the corresponding eigenfunctions identify flow patterns that are expected to dominate early stages of transition. This approach agrees with experiments in many relevant flows (e.g., those driven by thermal and centrifugal forces (Trefethen et al. 1993)) but it comes up short in matching experimental observations in wall-bounded shear flows (flows in channels, pipes, and boundary layers). Failure of classical hydrodynamic stability analysis is attributed to the nonnormal nature of the linearized equations which may manifest itself by transient growth of perturbations (Butler & Farrell 1992), protrusion of pseudospectra to the unstable regions (Trefethen et al. 1993; Trefethen & Embree 2005), and large receptivity to ambient disturbances (Farrell & Ioannou 1993; Bamieh & Dahleh 2001; Jovanović & Bamieh 2005). Even in stable regimes owing to nonnormality perturbations that grow transiently before decaying due to viscosity can be configured, irregularities in laboratory design can lead to instability, and background disturbances (such as freestream turbulence or surface imperfections) can be amplified by orders of magnitude. These conclusions can be reached by performing transient growth, pseudospectra, or variance amplification analyses (Grossmann 2000; Schmid & Henningson 2001; Schmid 2007). All of these methods demonstrate the importance of streamwise-elongated flow patterns of high and low streamwise velocity in transitional wall-bounded shear flows of
2 2 M. R. Jovanović and S. Kumar Newtonian fluids; this is at odds with classical modal stability results, but in agreement with experiments conducted in noisy environments (Matsubara & Alfredsson 2001). Transition to turbulence in viscoelastic fluids is important from both fundamental and technological perspectives (Larson 1992). The observation that transition can occur even when the effects of fluid elasticity dominate those of fluid inertia which is a primary cause of transition in Newtonian fluids is particularly intriguing (Larson 2000; Groisman & Steinberg 2000, 2004; Berti et al. 2008). Improved understanding of transition mechanisms in viscoelastic fluids has broad applications, ranging from deeper insight into order-disorder transitions in spatially extended nonlinear dynamical systems to enhanced mixing in microfluidic devices through the addition of polymers (Groisman & Steinberg 2001, 2004). Amplification of stochastic disturbances in channel flows of viscoelastic fluids has recently been investigated using linear systems theory (Hoda, Jovanović & Kumar 2008). For the Oldroyd-B constitutive model one of the simplest describing viscoelastic effects computations reported in Hoda et al. (2008) demonstrated that streamwise-constant disturbances can be considerably amplified in the weakly inertial/strongly elastic regime. As in Newtonian fluids, this amplification is fundamentally nonmodal in nature: it cannot be described using the standard normal mode decomposition of classical hydrodynamic stability analysis (Grossmann 2000; Schmid & Henningson 2001; Schmid 2007). Rather, it arises due to an energy exchange involving the fluctuations in the streamwise/wallnormal polymer stress and the wall-normal gradient of the streamwise velocity (Hoda, Jovanović & Kumar 2009). The key parameters in viscoelastic fluids are: the viscosity ratio, β = η s /(η s + η p ), where η s and η p are the solvent and polymer viscosities; the Reynolds number, Re = ρu o L/(η s + η p ), which represents the ratio of inertial to viscous forces; and the elasticity number, µ = We/Re, which quantifies the strength of elastic forces relative to inertial forces. Here, We = λu o /L is the Weissenberg number, λ is the fluid relaxation time, U o is the largest base flow velocity, L is the channel half-height, and ρ is the fluid density. By modeling ambient disturbances as stochastic forcing to streamwise-constant channel flows of Oldroyd-B fluids (with spanwise wavenumber k z ), an explicit scaling of the variance (energy) amplification with the Reynolds number Re was developed in Hoda et al. (2009), E(k z ; Re, β, µ) = f(k z ; β, µ)re + g(k z ; β, µ)re 3, where f and g denote Re-independent functions. It is worth noting that E quantifies the ensemble-average energy density (associated with the velocity field) of the statistical steady-state (Farrell & Ioannou 1993, 1994). In this work, considering the flows with µ 1, we apply singular perturbation techniques to establish that the steady-state variance scales as E(k z ; Re, β, µ) ˆf(k z ; β)re + ĝ(k z ; β)µre 3. This scaling was hypothesized by Hoda et al. (2009) on the basis of numerical data; here, we furnish an analytical proof of its validity. The last expression should be compared to the expression for the variance amplification in Newtonian fluids (Bamieh & Dahleh 2001), E N (k z ; Re) = f N (k z )Re + g N (k z )Re 3. At low Re the k z -dependence of E N is governed by f N (k z ), E N (k z ; Re) f N (k z )Re, and at high Re it is governed by g N (k z ), E N (k z ) Re 3 g N (k z ) (Bamieh & Dahleh 2001). We note that the k z -dependence of both ˆf and f N is characterized by viscous dissipation, which clearly indicates that the behavior of Newtonian fluids with low Re is dominated
3 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 3 by diffusion. On the other hand, the g-functions exhibit peaks at k z O(1); the values of k z where these peaks take place identify the spanwise length scales of the most energetic response to stochastic forcing in Newtonian fluids with high Re, and in Oldroyd-B fluids with high µre 3. Our results thus show that in viscoelastic fluids even at arbitrarily low (but non-vanishing) Reynolds numbers the steady-state level of velocity variance may be dominated by the Re 3 -term owing to the elastic amplification of disturbances. This reveals the subtle interplay between inertial and elastic forces in Oldroyd-B fluids with low Reynolds/high elasticity numbers. In addition to furnishing an analytical proof for the scaling of the variance amplification at high µ, the other major contribution of this paper is the discovery that the linearized dynamics can be decomposed into slow and fast subsystems. This observation is used to cast the equations into a standard singularly perturbed form for which existing methodology (Kokotović, Khalil & O Reilly 1999) can be applied. The decomposition of the linearized dynamics at high µ is not obvious a priori, and it takes advantage of the intrinsic time scale (λ) in the Oldroyd-B constitutive equation. In addition, it facilitates derivation of an explicit analytical expression for the steady-state velocity variance. Our success with uncovering the hitherto unknown dependence of the energy amplification on elasticity number points to the scaling and modeling steps as prerequisites for applying standard singular perturbation techniques. Our subsequent developments are laid out next. In 2, we describe the streamwiseconstant linearized model and identify a slow-fast decomposition for different components of the frequency response operator. In 3, we employ singular perturbation techniques to determine the scaling of the steady-state velocity variance with elasticity number in strongly elastic polymer solutions. In 4, we provide an explicit analytical expression for the variance amplification and discuss physical mechanisms leading to amplification from different forcing to different velocity components. We also determine the spanwise length scales of flow structures that contribute most to the steady-state variance and show that these structures assume the form of high and low speed streaks. The major contributions are summarized in 5 and the mathematical developments are relegated to the appendices. 2. Streamwise-constant linearized model We consider incompressible channel flows of Oldroyd-B fluids with ɛ = 1/µ 1. The equations governing the dynamics (up to first order) of velocity (v = [ v 1 v 2 v 3 ] T ), pressure (p), and polymer stress tensor (τ ) fluctuations around base flow condition (v, τ ) are brought to a non-dimensional form by scaling length with L, velocity with U o, polymer stresses with η p U o /L, pressure with (η s + η p )U o /L, and time with λ (see Appendix A) ɛ v = Re ( v v + v v) p + (1 β) τ + β 2 v + ɛre d, 0 = v, τ = v + ( v) T τ + Re ɛ ( τ v + τ v + (τ v) T + (τ v) T v τ v τ ). (2.1) Here, a dot signifies a partial derivative with respect to time t, is the gradient, v = v, and v 1, v 2, and v 3 are the velocity fluctuations in the streamwise (x), wallnormal (y), and spanwise (z) directions, respectively. If t denotes time normalized by the convective time scale L/U o, then t = t/we. The linearized momentum equation is driven by the body force fluctuation vector d = [ d 1 d 2 d 3 ] T, which is considered to be purely harmonic in the horizontal directions,
4 4 M. R. Jovanović and S. Kumar d 1 ɛre F1 d 2 ɛre v 2 ω2 v1 F2 J2 ReCp J1 G 1 d 3 ɛre F3 v2 G2 G3 v3 Figure 1. Block diagram of the linearized streamwise-constant model. and stochastic in the wall-normal direction and in time. This spatio-temporal forcing will in turn yield velocity and polymer stress fluctuations of the same nature. Our goal is to study the steady-state level of variance maintained in v by assuming that d is a temporally stationary white Gaussian process with zero mean and unit variance. It is noteworthy that if d denotes the forcing in the standard non-dimensional momentum equation (with time scale L/U o ), then d(r, t) = We d(r, t), where r is the position vector, r = [ x y z ] T. This body-force scaling is introduced to guarantee the same autocorrelation operators of d(r, t) and d(r, t) (Kokotović et al. 1999) Frequency response representation We study the linearized model for streamwise-constant three-dimensional fluctuations, which means that the dynamics evolve in the (y, z)-plane, but flow fluctuations in three spatial directions are considered. This model is analyzed since the largest variance in stochastically forced channel flows of viscoelastic fluids is maintained by streamwiseconstant fluctuations (Hoda et al. 2008). The linearized equations can be brought to an evolution form by removing the pressure from the equations and by expressing v in terms of the wall-normal velocity/vorticity (v 2, ω 2 ) fluctuations (see Appendix B). Furthermore, application of the Fourier transform in t and z allows for elimination of the polymer stresses from the evolution model, which results in an equivalent frequency response representation of the linearized system (see Appendix C). The block diagram in figure 1 provides a systems-level view of the streamwise-constant model. The boxes represent different parts of the system and the circles denote summation of signals. Inputs into each box/circle are represented by lines with arrows directed toward the box/circle, and outputs of each box/circle are represented by lines with arrows leading away from the box/circle. The inputs specify the signals affecting subsystems, and the outputs designate the signals of interest or signals affecting other parts of the system (Åström & Murray 2008). All signals in figure 1 are functions of the wall-normal coordinate y, the spanwise wave-number k z, and the temporal frequency ω, e.g. v 2 = v 2 (y, k z, ω), with the following boundary conditions on v 2 and ω 2, {v 2 (±1, k z, ω) = y v(±1, k z, ω) = ω 2 (±1, k z, ω) = 0}. The capital letters in figure 1 denote the Reynolds-number-independent operators. These operators act in the wall-normal direction and some of them are parameterized by k z (F j and G j, j = 1, 2, 3), while the others depend on k z, ω, β, and ɛ (J 1, J 2, and C p ). The operators F j and G j are given by F 1 = ik z, F 2 = k 2 z 1, F 3 = ik z 1 y, G 1 = (i/k z ), G 2 = I, G 3 = (i/k z ) y, and they, respectively, describe the way the forcing enters into the evolution model, and
5 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 5 the way the velocity fluctuations depend on the wall-normal velocity and vorticity. The operator C p captures the coupling from v 2 to ω 2, and it is defined as C p = C p1 + 1 ɛ(1 + iω) 2 C p2, where C p1 = ik z U (y) denotes the vortex tilting term (Butler & Farrell 1992), and C p2 = ik z (1 β) C p2, Cp2 = U (y) + 2U (y) y, denotes the term arising due to the work done by the polymer stresses on the flow (Hoda et al. 2009). Finally, J 1 and J 2 govern the internal dynamics of the wall-normal vorticity and velocity fluctuations, respectively. These two operators are given by J j = (1 + iω)k j, K j = ( ɛ(iω) 2 I (βt j ɛi)iω T j ) 1. Here, T 1 = and T 2 = 1 2, respectively, stand for the Squire and Orr-Sommerfeld operators in the streamwise-constant model of Newtonian fluids with Re = 1 (Schmid & Henningson 2001), I is the identity operator, = yy k 2 z is a Laplacian with homogeneous Dirichlet boundary conditions, 1 is the inverse of the Laplacian, 2 = yyyy 2k 2 z yy + k 4 z with homogeneous Cauchy (both Dirichlet and Neumann) boundary conditions, i = 1, U(y) = y in Couette flow, U(y) = 1 y 2 in Poiseuille flow, and U (y) = du(y)/dy. In the frequency domain, the velocity and forcing components are related by v i = H ij d j, where H ij denotes the ijth component of the frequency response operator H, v(y, k z, ω; Re, β, ɛ) = [H(k z, ω; Re, β, ɛ) d(, k z, ω; Re, β, ɛ)] (y). The variance maintained in v is given by (Farrell & Ioannou 1994) E(k z ; Re, β, ɛ) = 1 2π trace (H(k z, ω; Re, β, ɛ) H (k z, ω; Re, β, ɛ)) dω, where H is the adjoint of operator H. Using the properties of the trace operator we have E = 3 i,j = 1 E ij, where E ij is the variance maintained in v i by stochastically forcing the linearized model with d j. From the definitions of E and operators H ij, the following Re-scaling of variance amplification was recently obtained in Hoda et al. (2009) E(k z ; Re, β, ɛ) = f(k z ; β, ɛ)re + g(k z ; β, ɛ)re 3. (E) This expression for the steady-state variance of v is valid for all k z, Re, β, and ɛ, and it follows from the fact that the H ij are determined by (see Appendix C) H 1j (k z, ω; Re, β, ɛ) = Re 3 H 1j (k z, ω; β, ɛ), j = 2, 3, H ij (k z, ω; Re, β, ɛ) = Re H ij (k z, ω; β, ɛ), i, j = 2, 3; i = j = 1, H i1 (k z, ω; Re, β, ɛ) = 0, i = 2, 3, where the H ij represent the Re-independent operators, Furthermore, H 11 = ɛ (1 + iω) G 1 K 1 F 1, Hij = ɛ (1 + iω) G i K 2 F j, i, j = 2, 3, H 1j = (1/ ɛ)g 1 K 1 (ɛ(1 + iω) 2 C p1 + C p2 )K 2 F j, j = 2, 3. f = f 11 + f 22 + f 23 + f 32 + f 33, g = g 12 + g 13,
6 6 M. R. Jovanović and S. Kumar with f ij (k z ; β, ɛ) = 1 trace ( Hij (k z, ω; β, ɛ) 2π H ij(k z, ω; β, ɛ) ) dω, and similarly for g 1j (k z ; β, ɛ). Our objective is to establish how the Re-independent functions f and g depend on ɛ in flows with ɛ = 1/µ 1 in order to gain insight into the conditions under which strong elasticity amplifies stochastic disturbances. We next determine evolution equations for each H ij. We show that these equations admit a standard singularly perturbed form which is convenient for uncovering dependence of the variance amplification on elasticity number Evolution equations for H ij The Reynolds-number-independent contributions to the steady-state variance can be determined by recasting each H ij in the evolution form ẋ ij (y, k z, t) = A ij (k z )x ij (y, k z, t) + B j (k z )d j (y, k z, t), v i (y, k z, t) = C i (k z )x ij (y, k z, t), where x ij is a vector of state variables, and (d j, v i ) is the input-output pair for frequency response operator H ij. Note that x ij and operators A ij, B j, and C i will, in general, be different for each H ij. It is a standard fact (Farrell & Ioannou 1993) that the variance of v i sustained by d j is determined by trace (P ij C i C i), where P ij denotes the steady-state auto-correlation operator of x ij, which is found by solving the Lyapunov equation, A ij P ij + P ij A ij = B j B j. Evolution equations for H ij with i = j = 1 and i, j = 2, 3 The operators { H 11 = ɛ (1 + iω) G 1 K 1 F 1 ; H ij = ɛ (1 + iω) G i K 2 F j, i, j = 2, 3} can be represented by the following set of equations (with the left-hand side denoting relations in the frequency domain, and the right-hand side denoting relations in the time domain): σ kj = K k F j d j ɛ σ kj = T k σ kj + (βt k ɛi) σ kj + F j d j, v i = ɛ (1 + iω) G i σ kj v i = ɛ G i (σ kj + σ kj ), with {k = 1 for j = 1; k = 2 for j = 2, 3}. By selecting x kj = σ kj and z kj = σ kj we obtain the following evolution equation [ ] [ ] [ ] [ ] ẋij 0 I xij 0 = + d ɛż ij T k βt k ɛi z j, ij Fj v i = ɛ [ ] [ ] (2.2) x G i G ij i, z ij with {k = 1 for i = 1; k = 2 for i = 2, 3}, homogeneous Dirichlet boundary conditions on x 11 and z 11, and homogeneous Cauchy boundary conditions on x ij and z ij for i, j = 2, 3. Evolution equations for H 1j with j = 2, 3 From Appendix C it follows that H 1j = (1/ ɛ)g 1 K 1 (ɛ(1 + iω) 2 C p1 + C p2 )K 2 F j, j = 2, 3, which yields an equivalent block diagram representation of H12 and H 13 in figure 2. Now, using figure 2 and the definitions of operators K 1 and K 2 we recast H 1j, with
7 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 7 d j ψ F j K2 ɛ(1 + iω) 2 C p1 ϕ K 1 φ G 1/ ɛ v1 Cp2 Figure 2. Block diagram of H 1j(k z, ω; β, ɛ), j = 2, 3. j = 2, 3, into the following set of equations (with the left-hand side denoting relations in the frequency domain, and the right-hand side denoting relations in the time domain): ψ = K 2 F j d j ɛ ψ = T 2 ψ + (βt 2 ɛi) ψ + F j d j, ϕ = ( ɛ(1 + iω) 2 C p1 + C p2 ) ψ ϕ = Cp1 (ɛ ψ + 2ɛ ψ + ɛψ) + C p2 ψ, φ = K 1 ϕ ɛ φ = T 1 φ + (βt 1 ɛi) φ + ϕ, v 1 = (1/ ɛ) G 1 φ v 1 = (1/ ɛ) G 1 φ, with homogeneous Dirichlet boundary conditions on φ, and homogeneous Cauchy boundary conditions on ψ. By selecting x = [ ψ φ ] T, z = [ ψ φ ] T, we obtain a singularly perturbed realization of H 1j, j = 2, 3, [ ẋ ɛż ] = [ 0 I A 21 (ɛ) A 22 (ɛ) v 1 = 1/ɛ [ C 1 0 ] [ x z ] [ ] [ x 0 + z B ] 2, ] d j, where all operators are partitioned conformably with the elements of x and z, [ ] [ ] T A 21 (ɛ) = 2 0 F, B C p2 + C p1 (T 2 + ɛi) T 2 = j, 1 C p1 F j [ ] βt A 22 (ɛ) = 2 ɛi 0, C C p1 (βt 2 + ɛi) βt 1 ɛi 1 = [ ] 0 G 1. (2.3) Equations (2.2) and (2.3) are in the standard singularly perturbed form (Kokotović et al. 1999) as the time-derivative of the second part of the state is multiplied by a small positive parameter ɛ and the lower-right-hand-corner-blocks of dynamical generators in both (2.2) and (2.3) at ɛ = 0 are invertible. Furthermore, this representation gives evolution equations for different components of the frequency response operator with a lower number of states compared to the original evolution model (B 1). In particular, there are two state variables in the evolution equations for operators H ij with {i = j = 1; i, j = 2, 3} (cf. (2.2)), and four state variables in the evolution equations for operators H 12 and H 13 (cf. (2.3)). (In comparison, there are eight states in the evolution model (B 1).) We exploit the structure of equations (2.2) and (2.3) in 3 to uncover a slow-fast decomposition of each H ij and provide an explicit analytical expression for the steady-state velocity variance in flows of strongly elastic polymer solutions. These analytical developments are then utilized in 4 to clearly identify important physical mechanisms leading to amplification from different forcing to different velocity components. 3. Singular perturbation analysis of variance amplification Considering the case of high µ, ɛ = 1/µ 1, we now apply singular perturbation methods to show how functions f and g in the expression for the steady-state velocity
8 8 M. R. Jovanović and S. Kumar variance (E) scale with the elasticity number. From 2.2 it follows that the evolution equations of operators H ij assume the form [ ] [ ] [ ] [ ] ẋ 0 I x 0 = + d ɛż A 21 (ɛ) A 22 (ɛ) z j, B2 v i = r(ɛ) [ ] [ ] (3.1) x C 1 C 2, z with appropriate boundary conditions on x and z, r(ɛ) = ɛ for {i = j = 1; i, j = 2, 3}, and r(ɛ) = 1/ ɛ for {i = 1; j = 2, 3}. To simplify notation we have omitted i and j indices in (3.1); it is to be noted, however, that x, z, r, and the A-operators depend on both i and j, the B-operators depend on j, and the C-operators depend on i. The following coordinate transformation (Kokotović et al. 1999) [ ] [ ] [ ] ξ I ɛq(ɛ)l(ɛ) ɛq(ɛ) x =, (3.2) ζ L(ɛ) I z can be utilized to fully separate the slow and fast dynamics of system (3.1). Namely, if operators L(ɛ) and Q(ɛ) satisfy A 21 (ɛ) A 22 (ɛ)l(ɛ) ɛl(ɛ)l(ɛ) = 0, I Q(ɛ) (A 22 (ɛ) + ɛl(ɛ)) ɛl(ɛ)q(ɛ) = 0, then the change of coordinates (3.2) brings system (3.1) to the following equivalent representation [ ] [ ] [ ] [ ] ξ As (ɛ) 0 ξ Bs (ɛ) = 1 ζ 0 ɛ A + 1 f (ɛ) ζ ɛ B d j, f v i = r(ɛ) [ C s (ɛ) C f (ɛ) ] [ ] (3.4) ξ, ζ with A s (ɛ) = L(ɛ), A f (ɛ) = A 22 (ɛ) + ɛl(ɛ), B s (ɛ) = Q(ɛ)B 2, B f = B 2, C s (ɛ) = C 1 C 2 L(ɛ), C f (ɛ) = C 2 + ɛ(c 1 C 2 L(ɛ))Q(ɛ). (L) (Q) (3.5) The ξ- and ζ-subsystems in (3.4) identify the slow and fast subsystems, respectively. It is noteworthy that ζ penalizes the deviation of z from the slow invariant manifold z = L(ɛ)x, i.e. ζ = z z = z + L(ɛ)x. (Note that z = L(ɛ)x is an invariant manifold of unforced system (3.1) if z(t 0, ɛ) = L(ɛ)x(t 0, ɛ) implies z(t, ɛ) = L(ɛ)x(t, ɛ) for every t t 0.) On the other hand, the variable ξ is introduced to fully decouple the influence of ζ on the slow subsystem. We refer the reader to Kokotović et al. (1999) for additional details on slow-fast decomposition of singularly perturbed systems. It is now easy to show that the steady-state auto-correlation operator of [ ξ T ζ T ] T takes the form [ ] X(ɛ) Y P(ɛ) = (ɛ), Y(ɛ) (1/ɛ)Z(ɛ) where components of P are to be determined from the following system of equations A s (ɛ)x(ɛ) + X(ɛ)A s(ɛ) = B s (ɛ)b s(ɛ), A f (ɛ)y(ɛ) + ɛy(ɛ)a s(ɛ) = B f B s(ɛ), A f (ɛ)z(ɛ) + Z(ɛ)A f (ɛ) = B f B f. (3.6)
9 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 9 This implies that the steady-state variance of operator H ij (cf. 2.1) is given by E ij = r 2 (ɛ) trace ( X(ɛ)C s(ɛ)c s (ɛ) + (1/ɛ)Z(ɛ)C f (ɛ)c f (ɛ) ) + r 2 (ɛ) trace ( Y(ɛ)C s(ɛ)c f (ɛ) + C f (ɛ)c s(ɛ)y (ɛ) ). (3.7) Now, using the fact that A 21 (ɛ) and A 22 (ɛ) in equations (2.2) and (2.3) are given by we represent L and Q as A 21 (ɛ) = A 21,0 + ɛa 21,1, A 22 (ɛ) = A 22,0 + ɛa 22,1, L(ɛ) = ɛ n L n, Q(ɛ) = n = 0 ɛ n Q n, n = 0 and employ (regular) perturbation analysis to render (L) and (Q) into the following set of conveniently coupled equations: ɛ 0 : ɛ 1 : { L0 = A 1 22,0 A 21,0, = A 1 Q 0 22,0, { L1 = A 1 22,0 (A 21,1 A 22,1 L 0 L 0 L 0 ), Q 1 = (L 0 Q 0 + Q 0 (A 22,1 + L 0 )) A 1 22,0,. Based on (3.5) we conclude that Π(ɛ) = n = 0 ɛn Π n where Π stands for A s, A f, B s, C s, or C f ; this implies that a similar procedure can be used to simplify (3.6) and determine coefficients in the power series expansions of operators X, Y, and Z A s,0 X 0 + X 0 A s,0 = B s,0 B s,0, ɛ 0 : A f,0 Y 0 = B f B s,0, A f,0 Z 0 + Z 0 A f,0 = B f B f, A s,0 X 1 + X 1 A s,0 = ( ( A s,1 X 0 + X 0 As,1) Bs,0 B s,1 + B s,1 Bs,0), ɛ 1 : A f,0 Y 1 + B f B s,1 = ( A f,1 Y 0 + Y 0 As,0), A f,0 Z 1 + Z 1 A f,0 = ( A f,1 Z 0 + Z 0 Af,1),. The above results are used in what follows to reveal the scaling of functions f and g in (E) with elasticity number Scaling of functions f ij with µ We next examine how functions f ij (k z ; β, ɛ) depend on ɛ in flows with ɛ = 1/µ 1. The Re-independent operators H 11 = ɛ (1 + iω) G 1 K 1 F 1, Hij = ɛ (1 + iω) G i K 2 F j, i, j = 2, 3, admit the evolution equations given by (2.2). As shown in 2.1, f(k z ; β, ɛ) = i,j f ij(k z ; β, ɛ), {i = j = 1; i, j = 2, 3}, where f ij denotes the steady-state variance of system (2.2). A direct comparison of equations (2.2) and (3.1) yields r(ɛ) = ɛ, A 21 = T k, A 22 = βt k ɛi, B 2 = F j, C 1 = C 2 = G i,
10 10 M. R. Jovanović and S. Kumar which in combination with (3.7) and the above perturbation analysis can be used to obtain (see Appendix D.1) f(k z ; β, ɛ) = ˆf 0 (k z ; β) + ɛ n ˆfn (k z ; β). Here, ˆf n are functions independent of ɛ. Thus, in streamwise-constant channel flows with ɛ = 1/µ 1, the functions contributing to the Re-scaling of the steady-state velocity variance in (E) approximately become µ-independent, i.e. n = 1 f(k z ; β, ɛ) = ˆf 0 (k z ; β) + O(ɛ), where ˆf 0 (k z ; β) is determined by the terms of the form trace (Z 0 C f,0 C f,0), with A f,0 Z 0 + Z 0 A f,0 = B f B f ; see Appendix D.1 for details Scaling of functions g 1j with µ By examining H 12 and H 13 we can determine the µ-dependence of terms responsible for the Re 3 -scaling of the steady-state variance. As shown in 2.1, g(k z ; β, ɛ) = 3 j = 2 g 1j(k z ; β, ɛ), where g 1j denotes the steady-state variance of system (2.3). By comparing equations (2.3) and (3.1) we see that r(ɛ) = 1/ ɛ and C 2 0. The latter observation implies that C s = C 1 = [ 0 G 1 ] is ɛ-independent and that Cf (ɛ) = ɛc 1 Q(ɛ), which together with (3.7) can be used to obtain g 1j = (1/ɛ) trace (X j,0 C 1C 1 ) + O(1). In fact, a closer examination of (3.7) in conjunction with the above perturbation analysis shows that (see Appendix D.2) g(k z ; β, ɛ) = 1 ɛ ɛ n ĝ n (k z ; β), ɛ 1, n = 0 where ĝ 0 (k z ; β) = j=2,3 trace (X j,0c 1C 1 ), with A s,0 X j,0 + X j,0 A s,0 = B js,0 B js,0. Hence, in streamwise-constant channel flows with ɛ = 1/µ 1, the functions contributing to the Re 3 -scaling of the steady-state velocity variance approximately scale linearly with µ. 4. Variance amplification in elasticity-dominated flows Based on the above developments, it follows that in streamwise-constant Poiseuille and Couette flows of Oldroyd-B fluids with sufficiently large µ, the variance maintained in v is given by E(k z ; Re, β, µ) Re ˆf 0 (k z ; β) + µre 3 ĝ 0 (k z ; β), where ˆf 0 and ĝ 0 are functions independent of Re and µ. Therefore, in elasticity-dominated flows, the terms responsible for the Re- and Re 3 -scaling of the steady-state variance are, respectively, µ-independent and linearly dependent on µ. As noted in 1, this µ-scaling was suggested by numerical data from the frequency response analysis of Hoda et al. (2009). Here, we have shown how it follows directly from the governing equations by decomposing the linearized dynamics into slow and fast subsystems. It is readily shown that ˆf 0 (k z ; β) = f 0 (k z )/β,
11 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 11 (a) (b) (c) Figure 3. Plots of: (a) f 0(k z); (b) g 0(k z) in both Couette (solid curve) and Poiseuille (circles) flows; (c) f 0(k z) + M 2 (1 β) 2 g 0(k z) in Poiseuille flow with β = 0.1 and M = 0 (dashed curve), 10, 15, 20, 25, 30, 35. where the base-flow-independent function f 0 (k z ) is given by f 0 (k z ) = f N (k z ) = 0.5(trace (T 1 1 ) + trace (T 1 2 )), with f N (k z ) being the function that arises in the expression for variance amplification in Newtonian fluids as noted in 1. Furthermore, with Thus, for µ 1, ĝ 0 (k z ; β) = g 0 (k z )(1 β) 2 /β, g 0 (k z ) = (k 2 z/4) trace (T 1 1 E(k z ; Re, β, M) Re β C p2 T 2 2 C p2t 1 1 ). ( f0 (k z ) + M 2 (1 β) 2 g 0 (k z )), where M denotes the viscoelastic Mach number (Joseph 2001), M 2 = µre 2 = WeRe = (U o /c) 2, and c is the shear-wave speed, c 2 = (η s + η p )/(ρλ). This shows that the variance depends affinely on M 2 and that it increases monotonically with a decrease in the ratio of the solvent viscosity to the total viscosity. The analytical expressions for trace (T 1 j ), j = 1, 2, were derived in Bamieh & Dahleh (2001); these can be used to explicitly calculate f 0 (k z ) = f N (k z ), which is illustrated in figure 3(a). In Couette flow, the expression for g 0 (k z ) simplifies to (kz/4) 2 trace (T 2 2 T 1 1 ), and an explicit k z -dependence of g 0 can be derived after some manipulation. The resulting expression for g 0 (k z ) is used to generate the plot in figure 3(b); from this plot we observe the non-monotonic character of g 0 (k z ), with peak values at k z 2.07 (in Couette flow) and k z 2.24 (in Poiseuille flow). (In Poiseuille flow, determination of the expression for g 0 (k z ) is considerably more involved than in Couette flow; however, the method developed by Jovanović & Bamieh (2006) can be used to compute this quantity efficiently without resorting to spatial discretization of the underlying operators.). Plots illustrating influence of the Mach number on the k z -dependence of f 0 (k z ) + M 2 (1 β) 2 g 0 (k z ), in Poiseuille flow with β = 0.1, are shown in figure 3(c). Clearly, the behavior at low M is governed by viscous dissipation, f 0 (k z ); as M increases, additional physical effects become important which leads to appearance of the peaks at k z O(1). As mentioned earlier, the values of k z where these peaks emerge determine the spanwise length scales of the most energetic response to stochastic forcing in Oldroyd-B fluids with high Mach numbers. We next discuss the physical mechanisms leading to large amplification. As evident
12 12 M. R. Jovanovic and S. Kumar (a) (b) Figure 4. Streamwise velocity fluctuations v1 (z, y) containing the most variance in strongly elastic: (a) Couette flow; and (b) Poiseuille flow. from both figure 2 and the expression for g 0 (kz ), the coupling term Cp2 plays an essential role in variance amplification; if this term was zero, the dynamics of weakly inertial/strongly elastic flows would be dominated by viscous dissipation (cf. figure 3(c) with M = 0). We note that Cp2 depends on the background shear U 0 (y) and the spanwise/wallnormal variations in the flow fluctuations. A careful analysis of the governing equations (see Appendix C) demonstrates that this term emerges due to: (a) the wall-normal and T spanwise velocity (v2, v3 ) gradients in the equation for φ2 = [ τ22 τ23 τ33 ] ; (b) stretcht ing of φ2 by the background shear in the equation for φ4 = [ τ12 τ13 ] ; and (c) the φ4 -gradients in the equation for ω2. All of these give rise to polymer stretching, leading to a transfer of energy from the base flow to fluctuations which results in large steadystate variance in weakly inertial/strongly elastic flows. The importance of this mechanism was recently recognized by Doering, Eckhardt & Schumacher (2006), where several examples of solutions to the two-dimensional equations for Couette flow of Oldroyd-B fluids that display transient growth were provided. Such energy transfer has also been observed experimentally in elastic turbulence of a swirling flow between two parallel disks (Groisman & Steinberg 2000, 2004). However, the present work demonstrates that this transfer can be initiated even in rectilinear flows. Streamwise velocity flow structures that contain the most variance in strongly elastic flows with kz = 2.07 (Couette) and kz = 2.24 (Poiseuille) are shown in figure 4. These structures are purely harmonic in z and their wall-normal shapes are determined by 2 1 the principal eigenfunctions of operators (kz2 /4)T 1 1 C p2 T2 C p2 T1 (Farrell & Ioannou 1993). The most amplified sets of fluctuations are given by high (hot colors) and low (cold colors) speed streaks. In Couette flow the streaks occupy the entire channel width, and in Poiseuille flow they are antisymmetric with respect to the channel s centerline. These flow structures have striking resemblance to the initial conditions responsible for the largest transient growth in Newtonian fluids (Butler & Farrell 1992). Despite similarities, the fluctuations shown in figure 4 and in Butler & Farrell (1992) arise due to fundamentally different physical mechanisms: in high Re-flows of Newtonian fluids, vortex tilting is the main driving force for amplification; in high µ/low Re-flows of viscoelastic fluids, it is the polymer stretching mechanism mentioned earlier. This suggests that polymer stretching in elasticity-dominated flows effectively introduces the lift-up of flow fluctuations in a similar manner as vortex tilting does in inertia-dominated flows (Landahl 1975). 5. Conclusions In this work, we have analyzed nonmodal amplification of stochastic disturbances in streamwise-constant channel flows of Oldroyd-B fluids. For large elasticity numbers, the linearized governing equations can be decomposed into slow and fast subsystems, allowing
13 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 13 application of singular perturbation methods to obtain an analytical expression for the variance amplification associated with the velocity field. This expression agrees with that hypothesized by Hoda et al. (2009) on the basis of numerical data. For sufficiently large elasticity number, the variance amplification shows a peak at O(1) spanwise wavenumber, and the corresponding streamwise velocity fluctuations have a structure similar to that seen in high-reynolds-number flows of Newtonian fluids. The mechanism of the energy amplification involves polymer stretching, which gives rise to an energy transfer from the base flow to fluctuations. This transfer can be interpreted as an effective lift-up of velocity fluctuations, similar to the role vortex tilting plays in inertia-dominated flows. The present results further confirm our earlier observations (Hoda et al. 2008, 2009) that stochastic disturbances can be considerably amplified by elasticity even when inertial effects are weak. This amplification can serve as an initial stage of the development of high and low speed streaks, which upon reaching a finite amplitude may undergo secondary instability and thereby provide a bypass transition to elasticity-induced turbulence (Larson 2000; Groisman & Steinberg 2000, 2004; Berti et al. 2008). It is important to point out that although we have considered a particular class of disturbances in this work, our results raise the possibility that other types of disturbances might also be significantly amplified in elasticity-dominated flows. We note that in the present problem, energy amplification does not require the presence of curved streamlines in the base flow, which can give rise to linear instabilities in other geometries when the effects of elasticity dominate those of inertia (Larson 1992). (In the present problem, the base flow is linearly stable (Larson 1992).) However, the finite-amplitude flow structures created by the energy amplification explored here may well contain curved streamlines and be subject to further instabilities that lead to a disordered flow. Elasticity-induced turbulence may find use in promoting mixing in microfluidic devices, where inertial effects are weak due to the small geometries (Groisman & Steinberg 2001, 2004). In polymer processing applications, however, elasticity-induced turbulence is generally undesired (Larson 1992, 2000), and our work may aid the development of control strategies to maintain ordered flows. Finally, we point out that the slow-fast decomposition of the linearized dynamics we have uncovered here does not follow a priori from the governing equations in their original form. Identification of this decomposition was a necessary step in the application of the singular perturbation methods that were used to develop an analytical expression for the variance amplification. The approach taken in this work may be helpful in examining the asymptotic structure of other flows at high elasticity number, especially if such flows are subject to disturbances and have linearized governing equations which are nonnormal in nature. This work was partially supported by the National Science Foundation under CAREER Award CMMI (to M. R. J.) and by the Department of Energy under Award Number DE-FG02-07ER46415 (to S. K.). Appendix A. Governing equations The non-dimensional momentum conservation, mass conservation, and constitutive equations for an Oldroyd-B fluid are given by (Bird et al. 1987; Larson 1999) V t = 1 Re (β 2 V + (1 β) T P ) V V + F, 0 = V, T t = 1 We ( V + ( V)T T) V T + T V + (T V) T, (A 1)
14 14 M. R. Jovanović and S. Kumar where V is the velocity vector, P is the pressure, T is the polymeric contribution to the stress tensor, and F is the spatio-temporal body force per unit mass. These equations have been brought to a dimensionless form by scaling length with the channel half height L, velocity with the largest base velocity U o, time with L/U o, polymer stresses with η p U o /L, pressure with (η s + η p )U o /L, and body force with U 2 o /L. In flows with high elasticity numbers, ɛ = 1/µ 1, it is more convenient to rescale time as t = t/we, which leads to the following set of equations: ɛv t = β 2 V + (1 β) T P Re V V + ɛre F, 0 = V, T t = V + ( V) T T + Re ɛ (T V + (T V)T V T). (A 2) Since we are interested in considering problems with stochastic spatio-temporal excitation, the body-forces in equations (A 2) and (A 1) are related to each other by F(r, t) = We F(r, twe), where r denotes the vector of spatial coordinates; this scaling is introduced to guarantee the same auto-correlation operators of F(r, t) and F(r, t) (Kokotović et al. 1999). Linearized dynamics are given by (2.1) and they are obtained by decomposing each field in (A 2) into the sum of the base flow and fluctuations (i.e., V = v + v, T = τ + τ, P = p + p, F = 0 + d), and keeping terms only to first order in fluctuations. Appendix B. The streamwise-constant evolution model The evolution model of the linearized system is obtained by a standard conversion to the wall-normal velocity/vorticity (v 2, ω 2 ) formulation. The procedure described in Hoda et al. (2008) in combination with the Fourier transform in z converts system (2.1) with streamwise-constant fluctuations ( x ( ) 0) to: ɛ φ 1 = β S 11 φ 1 + (1 β) S 12 φ 2 + ɛre (F 2 d 2 + F 3 d 3 ), (B 1a) φ 2 = φ 2 + S 21 φ 1, (B 1b) ɛ φ 3 = β S 33 φ 3 + Re S 31 φ 1 + (1 β) S 34 φ 4 + ɛre F 1 d 1, (B 1c) φ 4 = φ 4 + (Re/ɛ) (S 41 φ 1 + S 42 φ 2 ) + S 43 φ 3, (B 1d) φ 5 = φ 5 (Re/ɛ) 2 S 51 φ 1 + (Re/ɛ) (S 53 φ 3 + S 54 φ 4 ), (B 1e) v 1 0 G 1 [ ] v 2 = G 2 0 φ1, (B 1f ) φ v 3 G where φ 1 = v 2, φ 2 = [ τ 22 τ 23 τ 33 ] T, φ 3 = ω 2, φ 4 = [ τ 12 τ 13 ] T, φ 5 = τ 11. The S- operators are given by: S 11 = 1 2, S 33 =, S 31 = ik z U (y), S 12 = [ ( ) ] 1 kz 2 y ik z yy + kz 2 kz 2 y, S34 = [ ] ik z y kz 2, S 21 = [ 2 y (i/k z ) ( ) ] yy + kz 2 T 2 y, S43 = (1/kz)S 2 T [ ] [ ] 34, U S 41 = (y) y U (y) U (i/k z )U, S (y) 42 = (y) 0 0 yy 0 U, (y) 0 S 51 = 4U (y)u (y), S 53 = (2i/k z )U (y) y, S 54 = [ 2U (y) 0 ].
15 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 15 The expressions for operators, 2, F j, and G j are given in 2.1. Note that (B 1) represents a system of PDEs in the wall-normal direction and in time parameterized by k z, Re, β, and ɛ. Appendix C. Frequency response operators The frequency response operators H ij, relating the forcing and velocity components, v i = H ij d j, can be obtained by applying the temporal Fourier transform to (B 1). Equation (B 1b) can be used to express φ 2 = [ τ 22 τ 23 τ 33 ] T in terms of the wall-normal velocity φ 1 = v 2 φ2 = S iω φ 1, (C 1) where ω denotes the temporal frequency. Substitution of (C 1) into the temporal Fourier transform of (B 1a) yields φ 1 = ( ɛre ɛiωi βs 11 1 β ) iω S 12S 21 (F 2 d 2 + F 3 d 3 ). (C 2) Based on this and equation (B 1f), it follows that, for streamwise constant fluctuations, streamwise forcing does not influence the wall-normal and spanwise velocities, i.e. H i1 (k z, ω; Re, β, ɛ) 0, i = 2, 3. Moreover, using (B 1f) and the fact that S 12 S 21 = S 11 = 1 2 =: T 2, the operators H ij (k z, ω; Re, β, ɛ), {i, j = 2, 3}, can be written as H ij (k z, ω; Re, β, ɛ) = Re H ij (k z, ω; β, ɛ), i, j = 2, 3, where the Re-independent operators H ij are given by H ij (k z, ω; β, ɛ) = ɛ (1 + iω) G i K 2 F j, i, j = 2, 3, K 2 = (ɛ(iω) 2 I (βt 2 ɛi)iω T 2 ) 1. The following relation between the wall-normal velocity/vorticity (φ 1, φ 3 ) and polymer stresses φ 4 = [ τ 12 τ 13 ] T can be established by substituting (C 1) into the temporal Fourier transform of (B 1d) φ 4 = (Re/ɛ) 1 + iω ( S 41 + S ) 42S 21 φ iω iω S 43φ 3. Substitution of this equation in the temporal Fourier transform of (B 1c) yields ( ) 1 ( (1 β) φ 3 = Re ɛiωi βs iω S 34S 43 S 31 + (1 β)s ( 34 S 41 + S )) 42S 21 φ 1 ɛ(1 + iω) 1 + iω + ( ) 1 (1 β) ɛre ɛiωi βs iω S 34S 43 F 1 d 1. (C 3) Now, since v 1 = G 1 φ 3, by substituting (C 2) into (C 3) and using the fact that in Couette and Poiseuille flows {S 34 S 43 = S 33 = =: T 1, S 34 S 41 = 0, S 34 S 42 S 21 = ik z (U (y) + 2U (y) y )}, it follows that operators H 1j (k z, ω; Re, β, ɛ), {j = 1, 2, 3}, are given by H 11 (k z, ω; Re, β, ɛ) = Re H 11 (k z, ω; β, ɛ), H 1j (k z, ω; Re, β, ɛ) = Re 3 H 1j (k z, ω; β, ɛ), j = 2, 3.
16 16 M. R. Jovanović and S. Kumar Here, the Re-independent operators H 1j are determined by H 11 (k z, ω; β, ɛ) = ɛ (1 + iω) G 1 K 1 F 1, H 1j (k z, ω; β, ɛ) = 1 ( ) G 1 K 1 ɛ(1 + iω) 2 C p1 + C p2 K2 F j, ɛ j = 2, 3, with K 1 = (ɛ(iω) 2 I (βt 1 ɛi)iω T 1 ) 1, C p1 = S 31 = ik z U (y), and C p2 = ik z (1 β) C p2, Cp2 = U (y) + 2U (y) y. Clearly, ik z Cp2 = S 34 S 42 S 21, which indicates that operator C p2 arises from: (a) the wall-normal and spanwise velocity (v 2, v 3 ) gradients, S 21, in the equation for φ 2 ; (b) stretching of φ 2 by the background shear, S 42, in the equation for φ 4 ; and (c) the φ 4 -gradients, S 34, in the equation for φ 3. Appendix D. Scaling of functions f ij and g 1j with µ In this appendix, we briefly describe how equation (3.7) can be combined with the perturbation analysis of 3 to uncover the scaling of variance amplification with µ in elasticity-dominated flows. The key observation here is that each operator in (3.7) exhibits a power series expansion, e.g. X(ɛ) = n = 0 ɛn X n. Thus, equation (3.7) can be rewritten as where Ψ(ɛ) = Φ(ɛ) = n = 0 E ij = r 2 (ɛ) trace ( Ψ(ɛ) ) + (r 2 (ɛ)/ɛ) trace ( Φ(ɛ) ), (D 1) ɛ n Ψ n = X(ɛ)C s(ɛ)c s (ɛ) + Y(ɛ)C s(ɛ)c f (ɛ) + C f (ɛ)c s (ɛ)y (ɛ), ɛ n Φ n = Z(ɛ)C f (ɛ)c f (ɛ). n = 0 In what follows, f ij = E ij for {i = j = 1; i, j = 2, 3} and g ij = E ij for {i = 1; j = 2, 3}. D.1. Scaling of functions f ij with µ We next examine the variance amplification of operators H ij with {i = j = 1; i, j = 2, 3}. For each of these operators r(ɛ) = ɛ and, thus, equation (D 1) simplifies to = ɛ ɛ n trace ( ) Ψ n + ɛ n trace ( ) Φ n. f ij n = 0 The expression for function f = i,j f ij given in 3.1 is a direct consequence of this equation. In particular, we have n = 0 f(k z ; β, ɛ) = i,j f ij (k z ; β, ɛ) = i,j ˆf ij,0 (k z ; β) + O(ɛ), where ɛ-independent functions ˆf ij,0 (k z ; β) are determined by ˆf ij,0 (k z ; β) = trace (Φ 0 ) = trace (Z 0 C f,0 C f,0), A f,0 Z 0 + Z 0 A f,0 = B f B f, Note that the Lyapunov equation for Z 0 follows from the perturbation analysis of 3 and that C f,0 is indexed by i, whereas A f,0 and B f (and consequently Z 0 ) are indexed by j A f,0 = βt k, {k = 1 for j = 1; k = 2 for j = 2, 3}, B f = F j, C f,0 = G i.
17 Nonmodal amplification of disturbances as a route to elasticity-induced turbulence 17 A direct comparison of the above expressions with their Newtonian fluid equivalents (see Bamieh & Dahleh 2001; Jovanović & Bamieh 2005) yields ˆf ij,0 (k z ; β) = f ij,0 (k z )/β. Here, fij,0 (k z ) denotes variance amplification from d j to v i, {i = j = 1; i, j = 2, 3}, in Newtonian fluids with Re = 1. Thus, in flows with high elasticity number, the term responsible for the Re-scaling in (E) is inversely proportional to β and approximately µ- independent. Finally, as shown in Bamieh & Dahleh (2001); Jovanović & Bamieh (2005), f 0 (k z ) = i,j f ij,0 (k z ) = 0.5(trace (T 1 1 ) + trace (T 1 2 )). The k z -dependence of this base-flow-independent function is shown in figure 3(a). D.2. Scaling of functions g 1j with µ Analysis of H 12 and H 13 can be used to determine the µ-dependence of terms responsible for the Re 3 -scaling of the steady-state velocity variance. As observed in 3.2, for both these operators C s = C 1 is ɛ-independent, C f (ɛ) = ɛc 1 Q(ɛ), and r(ɛ) = 1/ ɛ. Thus, operators Ψ(ɛ) and Φ(ɛ) in (D 1) can be written as Ψ(ɛ) = X(ɛ)C 1C 1 + ɛ ( Y(ɛ)C 1C 1 Q(ɛ) + Q (ɛ)c 1C 1 Y (ɛ) ), Φ(ɛ) = ɛ 2 Z(ɛ)Q (ɛ)c 1C 1 Q(ɛ). These two expressions in conjunction with r(ɛ) = 1/ ɛ and equation (D 1) can be used to obtain g 1j = (1/ɛ) trace ( X(ɛ)C 1C 1 ) + trace ( Z(ɛ)Q (ɛ)c 1C 1 Q(ɛ) ) + trace ( Y(ɛ)C 1C 1 Q(ɛ) + Q (ɛ)c 1C 1 Y (ɛ) ), where X(ɛ), Y(ɛ), and Z(ɛ) are indexed by j. As noted before, each of the operators in the equation for g 1j can be decomposed into a power series expansion, which in turn yields, 3 g(k z ; β, ɛ) = g 1j (k z ; β, ɛ) = 1 ɛ n ĝ n (k z ; β), ɛ 1, ɛ j = 2 n = 0 where each ĝ n (k z ; β) is ɛ-independent. In particular, at high µ we have g(k z ; β, ɛ) = (1/ɛ)ĝ 0 (k z ; β) + O(1) = (1/ɛ) 3 trace (X j,0 C 1C 1 ) + O(1). j = 2 A careful analysis of the Lyapunov equation for X j,0, A s,0 X j,0 + X j,0 A s,0 = B js,0 B js,0, j = 2, 3, in combination with the structure of operator C 1 can be used to finally obtain ĝ 0 (k z ; β) = g 0 (k z )(1 β) 2 /β, g 0 (k z ) = (k 2 z/4) trace (T 1 1 C p2 T 2 2 C p2t 1 1 ). An in-depth study of function g 0 (k z ) and its importance in early stages of transition to elasticity-induced turbulence is provided in 4.
18 18 M. R. Jovanović and S. Kumar REFERENCES Åström, K. J. & Murray, R. M Feedback systems: an introduction for scientists and engineers. Princeton University Press. Bamieh, B. & Dahleh, M Energy amplification in channel flows with stochastic excitations. Phys. Fluids 13, Berti, S., Bistagnino, A., Boffetta, G., Celani, A. & Musacchio, S Twodimensional elastic turbulence. Phys. Rev. E 77, (R). Bird, R. B., Curtiss, C. F., Armstrong, R. C. & Hassager, O Dynamics of Polymeric Liquids. Wiley. Butler, K. M. & Farrell, B. F Three-dimensional optimal perturbations in viscous shear flow. Phys. Fluids A 4, Doering, C., Eckhardt, B. & Schumacher, J Failure of energy stability in Oldroyd- B fluids at arbitrarily low Reynolds numbers. J. Non-Newtonian Fluid Mech. 135 (2-3), Farrell, B. F. & Ioannou, P. J Stochastic forcing of the linearized Navier-Stokes equations. Phys. Fluids A 5, Farrell, B. F. & Ioannou, P. J Variance maintained by stochastic forcing of nonnormal dynamical systems associated with linearly stable flows. Phys. Rev. Lett. 72 (8), Groisman, A. & Steinberg, V Elastic turbulence in a polymer solution flow. Nature 405, Groisman, A. & Steinberg, V Efficient mixing at low Reynolds numbers using polymer additives. Nature 410, Groisman, A. & Steinberg, V Elastic turbulence in curvilinear flows of polymer solutions. New J. Phys. 6, 29. Grossmann, S The onset of shear flow turbulence. Rev. Mod. Phys. 72 (2), Hoda, N., Jovanović, M. R. & Kumar, S Energy amplification in channel flows of viscoelastic fluids. J. Fluid Mech. 601, Hoda, N., Jovanović, M. R. & Kumar, S Frequency responses of streamwise-constant perturbations in channel flows of Oldroyd-B fluids. J. Fluid Mech. 625, Joseph, D. D Fluid dynamics of viscoelastic liquids. Springer-Verlag. Jovanović, M. R. & Bamieh, B Componentwise energy amplification in channel flows. J. Fluid Mech. 534, Jovanović, M. R. & Bamieh, B A formula for frequency responses of distributed systems with one spatial variable. Syst. Control Lett. 55, Kokotović, P., Khalil, H. K. & O Reilly, J Singular perturbation methods in control: analysis and design. SIAM. Landahl, M. T Wave breakdown and turbulence. SIAM J. Appl. Math. 28, Larson, R. G Instabilities in viscoelastic flows. Rheol. Acta 31, Larson, R. G The Structure and Rheology of Complex Fluids. Oxford University Press. Larson, R. G Turbulence without inertia. Nature 405, Matsubara, M. & Alfredsson, P. H Disturbance growth in boundary layers subjected to free-stream turbulence. J. Fluid Mech. 430, Schmid, P. J Nonmodal stability theory. Annu. Rev. Fluid Mech. 39, Schmid, P. J. & Henningson, D. S Stability and transition in shear flows. Springer- Verlag. Trefethen, L. N. & Embree, M Spectra and pseudospectra: the behavior of nonnormal matrices and operators. Princeton University Press. Trefethen, L. N., Trefethen, A. E., Reddy, S. C. & Driscol, T. A Hydrodynamic stability without eigenvalues. Science 261,
Energy amplification in channel flows of viscoelastic fluids
J. Fluid Mech. (2008), vol. 601, pp. 407 424. c 2008 Cambridge University Press doi:10.1017/s0022112008000633 Printed in the United Kingdom 407 Energy amplification in channel flows of viscoelastic fluids
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 informationNonmodal amplification of disturbances in inertialess flows of viscoelastic fluids
Nonmodal amplification of disturbances in inertialess flows of viscoelastic fluids Mihailo Jovanović www.umn.edu/ mihailo Shaqfeh Research Group Meeting; July 5, 2010 M. JOVANOVIĆ, U OF M 1 Turbulence
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 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 informationDynamics of inertialess flows of viscoelastic fluids: the role of uncertainty
Dynamics of inertialess flows of viscoelastic fluids: the role of uncertainty Mihailo Jovanović www.umn.edu/ mihailo joint work with: Satish Kumar, Binh Lieu Workshop on Complex Fluids and Flows in Industry
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 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 informationEnergy Transfer Analysis of Turbulent Plane Couette Flow
Energy Transfer Analysis of Turbulent Plane Couette Flow Satish C. Reddy! and Petros J. Ioannou2 1 Department of Mathematics, Oregon State University, Corvallis, OR 97331 USA, reddy@math.orst.edu 2 Department
More informationOn using the streamwise traveling waves for variance suppression in channel flows
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 WeC19.4 On using the streamwise traveling waves for variance suppression
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 informationSystems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well
Systems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well Dedicated to J. Boyd Pearson Bassam Bamieh Mechanical & Environmental Engineering University of California
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 informationDynamics and control of wall-bounded shear flows
Dynamics and control of wall-bounded shear flows Mihailo Jovanović www.umn.edu/ mihailo Center for Turbulence Research, Stanford University; July 13, 2012 M. JOVANOVIĆ 1 Flow control technology: shear-stress
More informationTransient growth on boundary layer streaks. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Transient growth on boundary layer streaks Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Primary instability: TS vs TG Secondary instability of the streaks Movie of streaks
More informationUnderstanding and Controlling Turbulent Shear Flows
Understanding and Controlling Turbulent Shear Flows Bassam Bamieh Department of Mechanical Engineering University of California, Santa Barbara http://www.engineering.ucsb.edu/ bamieh 23rd Benelux Meeting
More informationÉCOLE DE PHYSIQUE DES HOUCHES. Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at
ÉCOLE DE PHYSIQUE DES HOUCHES STABILITY OF FLUID FLOWS Carlo COSSU Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at Département de Mécanique, École Polytechinique http://www.imft.fr/carlo-cossu/
More informationStochastic excitation of streaky boundary layers. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Stochastic excitation of streaky boundary layers Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Boundary layer excited by free-stream turbulence Fully turbulent inflow
More informationUsing Statistical State Dynamics to Study the Mechanism of Wall-Turbulence
Using Statistical State Dynamics to Study the Mechanism of Wall-Turbulence Brian Farrell Harvard University Cambridge, MA IPAM Turbulence Workshop UCLA November 19, 2014 thanks to: Petros Ioannou, Dennice
More informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow Fulvio Martinelli, Maurizio Quadrio, John McKernan and James F. Whidborne Abstract Subcritical flows may experience large
More informationELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY. David O. Olagunju. University of Delaware. Newark, DE 19716
ELASTIC INSTABILITIES IN CONE{AND{PLATE FLOW: SMALL GAP THEORY David O. Olagunju Department of Mathematical Sciences University of Delaware Newark, DE 19716 June 15, 1995 Abstract Consider the axisymmetric,
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 informationFrequency domain analysis of the linearized Navier-Stokes equations
Frequency domain analysis of the linearized Navier-Stokes equations Mihailo R. JovanoviC jrnihailoqengineering.ucsb.edu Bassam Bamieh barnieh@engineering.ucsb.edu Department of Mechanical and Environmental
More informationNonmodal Growth and the Unstratified MRI Dynamo
MPPC general meeting, Berlin, June 24 Nonmodal Growth and the Unstratified MRI Dynamo Jonathan Squire!! and!! Amitava Bhattacharjee MRI turbulence Observed accretion rate in disks in astrophysical disks
More informationReduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD
3rd AIAA Flow Control Conference, 5 8 June 26, San Francisco Reduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD M. Ilak and C. W. Rowley Dept. of Mechanical and Aerospace Engineering,
More informationModeling flow statistics using convex optimization
Modeling flow statistics using convex optimization Abstract A method is proposed to estimate the covariance of disturbances to a stable linear system when its state covariance is known and a dynamic model
More informationFeedback Control of Transitional Channel Flow using Balanced Proper Orthogonal Decomposition
5th AIAA Theoretical Fluid Mechanics Conference 3-6 June 8, Seattle, Washington AIAA 8-3 Feedback Control of Transitional Channel Flow using Balanced Proper Orthogonal Decomposition Miloš Ilak Clarence
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 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 informationModels of Turbulent Pipe Flow
Models of Turbulent Pipe Flow Thesis by Jean-Loup Bourguignon In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California 213
More informationTransient growth for streak-streamwise vortex interactions
Physics Letters A 358 006) 431 437 www.elsevier.com/locate/pla Transient growth for streak-streamwise vortex interactions Lina Kim, Jeff Moehlis Department of Mechanical Engineering, University of California,
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 informationPerturbation Energy Production in Pipe Flow over a Range of Reynolds Numbers using Resolvent Analysis
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition - 8 January 29, Orlando, Florida AIAA 29-3 Perturbation Energy Production in Pipe Flow over a Range of Reynolds
More informationSG2221 Wave Motion and Hydrodinamic Stability. MATLAB Project on 2D Poiseuille Flow Alessandro Ceci
SG2221 Wave Motion and Hydrodinamic Stability MATLAB Project on 2D Poiseuille Flow Alessandro Ceci Base Flow 2D steady Incompressible Flow Flow driven by a constant pressure gradient Fully developed flow
More information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationSUBCRITICAL INSTABILITIES IN PLANE COUETTE FLOW OF VISCO-ELASTIC FLUIDS
SUBCRITICAL INSTABILITIES IN PLANE COUETTE FLOW OF VISCO-ELASTIC FLUIDS Alexander N. Morozov and Wim van Saarloos Instituut-Lorentz for Theoretical Physics, LION, Leiden University, Postbus 956, 23 RA
More informationLes Houches School of Foam: Rheology of Complex Fluids
Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter
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 informationLinear stability of channel entrance flow
Linear stability of channel entrance flow Damien Biau DIAM, University of Genova, Via Montallegro 1, 16145 Genova, Italy Abstract The spatial stability analysis of two dimensional, steady channel flow
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 information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationAmplification of coherent streaks in the turbulent Couette flow: an input output analysis at low Reynolds number
J. Fluid Mech. (2), vol. 643, pp. 333 348. c Cambridge University Press 2 doi:.7/s22299925 333 Amplification of coherent streaks in the turbulent Couette flow: an input output analysis at low Reynolds
More informationNumerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section
Korea-Australia Rheology Journal Vol. 19, No. 2, August 2007 pp. 67-73 Numerical study of flow of Oldroyd-3-Constant fluids in a straight duct with square cross-section Mingkan Zhang, Xinrong Shen, Jianfeng
More informationFlow Instability Matlab Tutorial Bangalore, January 2010 Matthew Juniper, Dan Henningson, Peter Schmid 1
Flow Instability Matlab Tutorial Bangalore, January Matthew Juniper, Dan Henningson, Peter Schmid Simple differentiation matrices Imagine a vector that represents the values of a function y at N equispaced
More informationJournal of Non-Newtonian Fluid Mechanics
J. Non-Newtonian Fluid Mech. 165 (21) 925 931 Contents lists available at ScienceDirect Journal of Non-Newtonian Fluid Mechanics journal homepage: www.elsevier.com/locate/jnnfm Velocity fluctuations in
More informationInertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes
J. Non-Newtonian Fluid Mech., 78 (1998) 27 45 Inertial effect on stability of cone-and-plate flow Part 2: Non-axisymmetric modes Yuriko Renardy *, David O. Olagunju 1 Department of Mathematics and ICAM,
More informationNUMERICAL ANALYSIS OF THE VISCOELASTIC FLUID IN PLANE POISEUILLE FLOW
NUMERICAL ANALYSIS OF THE VISCOELASTIC FLUID IN PLANE POISEUILLE FLOW N. Khorasani and B. Mirzalou 1 Department of Mechanical and Aerospace engineering, Science and Research Branch, Islamic Azad University,
More informationNon-modal stability analysis of stratified two-phase channel flows
Non-modal stability analysis of stratified two-phase channel flows I. Barmak a), A. Yu. Gelfgat, A. Ullmann, and N. Brauner School of Mechanical Engineering, Tel Aviv University, Tel Aviv 69978, Israel
More informationEstimating large-scale structures in wall turbulence using linear models
J. Fluid Mech. (218), vol. 842, pp. 146 162. c Cambridge University Press 218 doi:1.117/jfm.218.129 146 Estimating large-scale structures in wall turbulence using linear models Simon J. Illingworth 1,,
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 informationThe Reynolds experiment
Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).
More informationModal and non-modal stability of particle-laden channel flow. Abstract
Modal and non-modal stability of particle-laden channel flow Joy Klinkenberg,, 2 H. C. de Lange, and Luca Brandt 2 TU/e, Mechanical Engineering, 56 MB, Eindhoven, The Netherlands 2 Linné Flow Centre, KTH
More informationPhysics of irregularity flow and tumultuous conversion in shear flows
African Journal of Physics Vol. 1 (3), 15 pages, November, 2013. Available online at www.internationalscholarsjournals.org International Scholars Journals Full Length Research Paper Physics of irregularity
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationToward low order models of wall turbulence using resolvent analysis
Center for Turbulence Research Proceedings of the Summer Program 2016 305 Toward low order models of wall turbulence using resolvent analysis By K. Rosenberg, T. Saxton-Fox, A. Lozano-Durán, A. Towne AND
More informationLINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP
MAGNETOHYDRODYNAMICS Vol. 48 (2012), No. 1, pp. 147 155 LINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP N. Vetcha, S. Smolentsev, M. Abdou Fusion Science and Technology Center, UCLA,
More informationPerturbation of System Dynamics and the Covariance Completion Problem
2016 IEEE 55th Conference on Decision and Control (CDC) ARIA Resort & Casino December 12-14, 2016, Las Vegas, USA Perturbation of System Dynamics and the Covariance Completion Problem Armin Zare, Mihailo
More informationChapter 13. Eddy Diffusivity
Chapter 13 Eddy Diffusivity Glenn introduced the mean field approximation of turbulence in two-layer quasigesotrophic turbulence. In that approximation one must solve the zonally averaged equations for
More informationThe effective slip length and vortex formation in laminar flow over a rough surface
The effective slip length and vortex formation in laminar flow over a rough surface Anoosheh Niavarani and Nikolai V. Priezjev Movies and preprints @ http://www.egr.msu.edu/~niavaran A. Niavarani and N.V.
More informationNonmodal stability analysis
Nonmodal stability analysis Peter Schmid Laboratoire d Hydrodynamique (LadHyX) CNRS Ecole Polytechnique, Palaiseau Motivation: bring together expertise in combustion physics, fluid dynamics, stability
More informationINTERFACIAL WAVE BEHAVIOR IN OIL-WATER CHANNEL FLOWS: PROSPECTS FOR A GENERAL UNDERSTANDING
1 INTERFACIAL WAVE BEHAVIOR IN OIL-WATER CHANNEL FLOWS: PROSPECTS FOR A GENERAL UNDERSTANDING M. J. McCready, D. D. Uphold, K. A. Gifford Department of Chemical Engineering University of Notre Dame Notre
More informationTransient and asymptotic growth of two-dimensional perturbations in viscous compressible shear flow
PHYSICS OF FLUIDS VOLUME 12, NUMBER 11 NOVEMBER 2000 Transient and asymptotic growth of two-dimensional perturbations in viscous compressible shear flow B. F. Farrell Division of Engineering and Applied
More informationProject Topic. Simulation of turbulent flow laden with finite-size particles using LBM. Leila Jahanshaloo
Project Topic Simulation of turbulent flow laden with finite-size particles using LBM Leila Jahanshaloo Project Details Turbulent flow modeling Lattice Boltzmann Method All I know about my project Solid-liquid
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 informationLinear Instability of asymmetric Poiseuille flows
Report no. 7/ Linear Instability of asymmetric Poiseuille flows Dick Kachuma & Ian J. Sobey We compute solutions for the Orr-Sommerfeld equations for the case of an asymmetric Poiseuille-like parallel
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
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 informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More information(1) Transition from one to another laminar flow. (a) Thermal instability: Bernard Problem
Professor Fred Stern Fall 2014 1 Chapter 6: Viscous Flow in Ducts 6.2 Stability and Transition Stability: can a physical state withstand a disturbance and still return to its original state. In fluid mechanics,
More informationTransient growth and minimal defects: Two possible initial paths of transition to turbulence in plane shear flows
PHYSICS OF FLUIDS VOLUME 16, NUMBER 10 OCTOBER 2004 Transient growth and minimal defects: Two possible initial paths of transition to turbulence in plane shear flows Damien Biau a) and Alessandro Bottaro
More informationEXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES
ICMAR 2014 EXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES Introduction A.V. Ivanov, Y.S. Kachanov, D.A. Mischenko Khristianovich Institute of
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 informationDirect Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids
Direct Simulation of the Motion of Solid Particles in Couette and Poiseuille Flows of Viscoelastic Fluids by P. Y. Huang 1, J. Feng 2, H. H. Hu 3 and D. D. Joseph 1 1 Department of Aerospace Engineering
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationReduced-Order Models for Feedback Control of Transient Energy Growth
9th AIAA Flow Control Conference, 5 9 June, Atlanta, Georgia. -9 Reduced-Order Models for Feedback Control of Transient Energy Growth Aniketh Kalur and Maziar S. Hemati Aerospace Engineering and Mechanics,
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 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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationCorrection of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves
PHYSICS OF FLUIDS 19, 082105 2007 Correction of Lamb s dissipation calculation for the effects of viscosity on capillary-gravity waves J. C. Padrino and D. D. Joseph Aerospace Engineering and Mechanics
More informationTowards low-order models of turbulence
Towards low-order models of turbulence Turbulence in Engineering Applications Long Programme in Mathematics of Turbulence IPAM, UCLA Ati Sharma Beverley McKeon + Rashad Moarref + Paco Gomez Ashley Willis
More informationROUTES TO TRANSITION IN SHEAR FLOWS
ROUTES TO TRANSITION IN SHEAR FLOWS Alessandro Bottaro with contributions from: D. Biau, B. Galletti, I. Gavarini and F.T.M. Nieuwstadt ROUTES TO TRANSITION IN SHEAR FLOWS Osborne Reynolds, 1842-1912 An
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationA Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
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 informationMagnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation
Center for Turbulence Research Annual Research Briefs 006 363 Magnetic waves in a two-component model of galactic dynamo: metastability and stochastic generation By S. Fedotov AND S. Abarzhi 1. Motivation
More informationJ. Non-Newtonian Fluid Mech. 91 (2000) Received in revised form 1September1999
J. Non-Newtonian Fluid Mech. 91 (2000 85 104 An experimental/theoretical investigation of interfacial instabilities in superposed pressure-driven channel flow of Newtonian and well-characterized viscoelastic
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 informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationChapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers)
Chapter 5 Phenomena of laminar-turbulent boundary layer transition (including free shear layers) T-S Leu May. 3, 2018 Chapter 5: Phenomena of laminar-turbulent boundary layer transition (including free
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationSTABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET. N. Vetcha, S. Smolentsev and M. Abdou
STABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET N. Vetcha S. Smolentsev and M. Abdou Fusion Science and Technology Center at University of California Los Angeles CA
More informationStability of two-layer viscoelastic plane Couette flow past a deformable solid layer
J. Non-Newtonian Fluid Mech. 117 (2004) 163 182 Stability of two-layer viscoelastic plane Couette flow past a deformable solid layer V. Shankar Department of Chemical Engineering, Indian Institute of Technology,
More informationSingular Value Decomposition of the frequency response operator
Measures of input-output amplification (across frequency) Largest singular value Hilbert-Schmidt norm (power spectral density) Systems with one spatial variable Two point boundary value problems EE 8235:
More informationCompletion of partially known turbulent flow statistics
2014 American Control Conference (ACC) June 4-6, 2014. Portland, Oregon, USA Completion of partially known turbulent flow statistics Armin Zare, Mihailo R. Jovanović, and Tryphon T. Georgiou Abstract Second-order
More informationCompletion of partially known turbulent flow statistics via convex optimization
Center for Turbulence Research Proceedings of the Summer Program 2014 345 Completion of partiall known turbulent flow statistics via convex optimization B A. Zare, M. R. Jovanović AND T. T. Georgiou Second-order
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationInterfacial waves in steady and oscillatory, two-layer Couette flows
Interfacial waves in steady and oscillatory, two-layer Couette flows M. J. McCready Department of Chemical Engineering University of Notre Dame Notre Dame, IN 46556 Page 1 Acknowledgments Students: M.
More informationCharacteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-1103 Characteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
More information