Journal of Fluid Science and Technology

Size: px
Start display at page:

Download "Journal of Fluid Science and Technology"

Transcription

1 Bulletin of the JSME Vol.11, No.1, 2016 Journal of Fluid Science and Technology Relationships between local topology of vortical flow and pressure minimum feature derived from flow kinematics in isotropic homogeneous turbulence Katsuyuki NAKAYAMA, Yasumasa OHIRA and Hideki HASEGAWA Department of Mechanical Engineering, Aichi Institute of Technology Yakusa-cho, Toyota, Aichi , JAPAN Division of Mechanical Engineering, Graduate School of Engineering, Aichi Institute of Technology Yakusa-cho, Toyota, Aichi , JAPAN Received 27 February 2015 Abstract Relationships between local flow topology (geometry) of vortical flow and pressure minimum feature derived from its flow kinematics are investigated in isotropic homogeneous turbulence, in a new topological point of view with respect to vortical flow symmetry in a swirl plane. Consistency of the pressure minimum plane and the swirl plane, and relationships between the pressure minimum and vortical flow symmetry in their development are also analyzed. The pressure minimum feature is specified by the λ 2 definition and also the ˇλ 2 definition that is the integrated definition of the, Q and λ 2 definitions and specifies the pressure minimum in the swirl plane. The swirlity ϕ that represents the geometrical average of the intensity of swirling flow is applied to the statistical analysis of the flow topology and pressure minimum. It shows that the pressure minimum requires a certain flow symmetry, and that the development or decay of the pressure minimum is associated with that of the vortical flow symmetry especially in the state (process) of a vortex attaining or losing the pressure minimum feature. The vortices with high intensity of ϕ have the feature that the effect of vorticity components parallel to the swirl plane decreases in specifying the pressure minimum plane by the λ 2 definition, and then this pressure minimum plane tends to approach the swirl plane. Key words : Vortical flow, Symmetry, Pressure minimum, Swirlity, Velocity gradient tensor 1. Introduction Vortices are associated with many fluid phenomena and engineering fields in various scales, such as fluid machinery, power plants, wind turbines, and bio-inspired flight objects. Elucidation of vortical flow phenomena in such fields is an important issue to control or suppress the vortices. An important characteristic of a vortex with swirling flow is its stability by its flow geometry (topology). The symmetry of the vortical flow is associated with the stability (Lundgren, 1982), and the vortical flow induces the pressure minimum feature. The pressure minimum supports an inflow motion that prevents the diffusion of the vorticity and induces vortex stretching (Tennekes and Lumley, 1972). Thus the relationships between the detail vortical flow geometry (symmetry) and the feature of the pressure minimum is of great interest and important in vortex dynamics and the flow control. In definition or identification of a vortex, many definitions have been proposed from various physical aspects of vortices (Chong et al., 1990,Hunt et al., 1988,Jeong and Hussain, 1995,Kida and Miura, 1998,Cucitore et al., 1999,Zhou et al., 1999, Chakraborty et al., 2005, Wu et al., 2005, Haller, 2005, Finn and Boghosian, 2006, Zhang and Choudhury, 2006, Kolář, 2007, Nakayama et al., 2007, Nakayama et al., 2014). Although no universal definition has been established, the (Chong et al., 1990), Q (Hunt et al., 1988), and λ 2 (Jeong and Hussain, 1995) definitions in local approach are Paper No

2 popular definitions which are used frequently (Chakraborty et al., 2005). The definition specifies the invariant swirling motion with eigenvalues of the velocity gradient tensor u. The Q definition is associated with the positive second invariant Q of u, and specifies a vortex where the vorticity exceeds the irrotational straining, which is equivalent to the positive pressure Laplacian. The λ 2 definition ascertains the signs of the pressure Hessians discarding the unsteady straining and viscous effects, which identifies the local pressure minimum region induced by vortical (swirling) motion. Because the pressure minimum feature induced by the vortical or swirling flow should be specified in its swirl plane, a vortex definition, say ˇλ 2 definition, has been proposed to specify this feature (Nakayama et al., 2014). The ˇλ 2 definition is a unified definition that integrates these three definitions, which specifies the invariant vortical (swirling) flow and resulting pressure minimum in the swirl plane, including the vorticity exceeding the strain in the plane. As for identification of flow geometry, because the eigenvalues of u can specify the local flow geometry (pattern) which is Galilei invariant, they not only derived the -definition where = (Q/3) 3 + (R/2) 2 (> 0) and R denotes the third invariant of u, but also have contributed in the classification of the local flow and thus the analysis of turbulent flow (Chong et al., 1990,Blackburn et al., 1996). Although the essential physical interpretation of the complex eigenvalues was lacked and the eigenvalues were difficult to identify the detail flow geometry, a mathematical analysis of the local velocity with its decomposition into the radial and azimuthal velocities has derived their interpretation and several properties specifying the detail geometry (Nakayama, 2014). It has been shown that the imaginary part of the complex eigenvalues represents the geometrical average of intensity of the azimuthal flow, and the real part does the arithmetic average of that of the radial flow. Consequently the properties to specify the symmetry of vortical flow have been derived from u, and one symmetry property is associated with the pressure minimum feature in the swirl plane specified by the ˇλ 2 definition (Nakayama et al., 2014). Then this definition and the symmetry property of the vortical flow enable to examine relationships between the flow geometry and the pressure minimum. The present paper investigates the relationships between the vortical flow geometry and the pressure minimum feature, in local flow kinematics. Vortices in isotropic homogeneous turbulence are our subject as a study of primary characteristics of a vortex. The pressure minimum is specified by the criteria in the λ 2 and ˇλ 2 definitions that are associated with the pressure minimum plane and pressure feature in the swirl plane induced by the vortical flow. The relation of these pressure minimums to the symmetry of vortical flow, and relationships between their development/decay are analyzed. The swirlity ϕ (Nakayama, 2014) representing the intensity of swirling flow of a vortex is used as a parameter to examine the features of the flow geometry and pressure minimum. It shows that the pressure minimum is associated with the symmetry of the vortical flow especially in the state or process of a vortex attaining or losing the pressure minimum, in which the pressure minimum feature and symmetry tend to develop or decay each other. The pressure minimum plane specified by the λ 2 definition differs from the swirl plane in general. However, vortices with high ϕ have the feature that these two planes approach. This characteristic of the vortices indicates that the vorticity components parallel to the swirl plane decrease, and the pressure minimum plane specified by the λ 2 definition approaches the swirl plane. 2. Representation of Flow topology and pressure minimum We summarize the symmetry property of the local vortical flow geometry and swirlity (Nakayama, 2014), and the λ 2 and ˇλ 2 definitions with the representation of u in a coordinate system associated with the swirl plane. We consider a local flow geometry around a point in an instantaneous velocity field. In the reference frame with respect to the point, the local velocity v i (i = 1, 2, 3) around the point can be expressed as: v i = dx i /dt = ( v i / x j )x j, where the summation convention is applied and v i / x i = 0 in incompressible fluids. Then the invariant local flow geometry is specified by the eigenvalues and eigenvectors of u (= [ v i / x j ]) (Chong et al., 1990). We assume that u has a pair of complex conjugate eigenvalues ε R ± iψ (ψ > 0) and their eigenvectors ξ pl ± iη pl, where i is the imaginary unit, and a real eigenvalue ε a and its eigenvector ξ axis. Then the flow trajectory can be represented as 2e ε Rt (ξ pl cos ψt η pl sin ψt) + e ε at ξ axis, where the flow swirls in the plane defined by ξ pl and η pl (Chong et al., 1990, Zhou et al., 1999), hereafter referred to as the swirl plane P. Swirlity ϕ represents the geometrical average of the intensity of the azimuthal velocity component in v i, i.e., the eigenvalues of the specific quadratic forms of the azimuthal velocity, and is specified as ϕ = Q + 3ε 2 a/4. The condition that 0 < ϕ is equivalent to the criterion that u has complex conjugate eigenvalues, i.e., > 0 (Chong et al., 1990), and ϕ = ψ in case of the vortical flow, i.e., 0 < ϕ (Nakayama, 2014). We note that both ξ pl and η pl are the degenerate eigenvectors of the matrix (A ε R E) 2 corresponding to the eigenvalue ψ 2, where A is a matrix of u in a coordinate system, i.e., A = [a i j ] = [ v i / x j ] (i, j = 1, 2, 3). The directions of ξ pl and η pl, and the relative length of these vectors, i.e., the ratio c = ξ pl / η pl, are specified by the eigenequations of A in terms 2

3 Fig. 1 flow geometries in swirl plane, with same complex eigenvalues (ε R = 0.5, ψ = 2) and different c. of the complex eigenvalues; A(ξ pl ± iη pl ) = (ε R ± iψ)(ξ pl ± iη pl ). Then we can specify that ξ pl and η pl are orthogonal, i.e., ξ pl η pl. The vortical flow geometry in P specified by the complex eigenvalues and eigenvectors of u, i.e., 2e ε Rt (ξ pl cos ψt η pl sin ψt), changes according to c. The geometry is elliptic in high c, and this symmetry increases as c approaches 1, as shown in Fig. 1. Then c is necessary to specify the geometry uniquely, and considered an invariant property of the symmetry of the swirling motion. c is also associated with the flow symmetry of the radial and azimuthal velocities in the plane, thus it represents the symmetry of whole vortical flow (Nakayama, 2014). We consider an orthonormal coordinate system ˇx i (i = 1, 2, 3) with bases ě i (i = 1, 2, 3) where ě 1 and ě 2 are parallel to ξ pl and η pl, respectively. Importantly, the ˇx 1 - ˇx 2 plane is P and the invariant subspace of u. u in this coordinate system, Ǎ = [ǎ i j ], can be expressed in the form Ǎ = ε R cψ ǎ 13 ψ/c ε R ǎ ε a. It is noted that ǎ 12 ǎ 21 except when ξ pl = η pl (c = 1) and that ǎ 13 = ǎ 23 = 0 when ξ axis P. c specifies the symmetry of the rotational components of u in P. In the ˇx i coordinate system, we consider the formulation of the λ 2 and ˇλ 2 definitions with the above representation of u. The second derivative of the pressure obtained by differentiating the Navier-Stokes equation, neglecting the unsteady and viscous strain terms yields (1) 1 ρ p,i j = s iks k j + ω ik ω k j, (2) where the subscript after the comma denotes the derivative in the directions of indices, and ρ and p denote the density and the pressure discarding the above-mentioned terms, respectively (Jeong and Hussain, 1995). s i j and ω i j denote the rate-of-strain tensor and vorticity tensor, respectively. If two of the eigenvalues λ i (i = 1, 2, 3) of p,i j /ρ are negative, i.e. λ 2 < 0 (here we set λ 1 λ 2 λ 3 ), then p has a local minimum characteristic, which yields the λ 2 definition. The ˇλ 2 definition focuses on the pressure minimum of p in P. The Hesse matrix of p, H = [h i j ] = [ p,i j /ρ] (i, j = 1, 2, 3), can be expressed as H = (AA + t A t A)/2, where the superscript t denotes the transpose of the matrix. Then, from Eq. (1), H in the ˇx i coordinate system, Ȟ, is derived as follows: ε 2 R ψ2 (c 1/c)ε R ψ ι 1 Ȟ = (c 1/c)ε R ψ ε 2 R ψ2 ι 2, (3) ι 1 ι 2 ε 2 a where ι 1 = ( ε R ǎ 13 + cψǎ 23 )/2 and ι 2 = ( ψǎ 13 /c ε R ǎ 23 )/2. The block matrix of Ȟ corresponding to P is expressed as Ȟ pl = [ȟ i j ] (i, j = 1, 2), which specifies the characteristics of p in P. Then the eigenvalues of Ȟ pl, ˇλ i (i = 1, 2; ˇλ 1 ˇλ 2 ), are given as: ˇλ 1, ˇλ 2 = (ε 2 R ψ2 ) (c 1 c )ε R ψ. (4) Because Ȟ pl is symmetric, ˇλ i (i = 1, 2) have the same sign if Ȟ pl > 0. The two conditions that Ȟ pl > 0 and ˇλ 1, ˇλ 2 < 0 can be integrated by the inequality with respect to ψ: c ε R < ψ, (5) 23

4 where c = c (1 c) or 1/c (c < 1). Then Eq. (5) is the condition that ˇλ 2 < 0, which represents the local minimum feature of p in P. Here we rotate the bases ě i (i = 1, 2) in P along ě 3 so that they coincide with the two orthonormal eigenvectors ˇζ i of ˇλ i (i = 1, 2), respectively. H in this orthonormal coordinate system ˇx i (i = 1, 2, 3), Ȟ is expressed as: ˇλ 1 0 ι 1 Ȟ = 0 ˇλ 2 ι 2 ι 1 ι 2 ε 2 a, where ι i (i = 1, 2) denotes the inner product of ˇζ i and ι (ι 1, ι 2 ), i.e., ι i = ˇζ i ι. Then p in P, ˇp, is expressed as ˇp/ρ = {ˇλ 1 ( ˇx 1 )2 + ˇλ 2 ( ˇx 2 )2 }/2. We examine the relationship among the λ 2 and ˇλ 2 definitions. Because p,ii /ρ = a i j a ji = Tr(H) (Jeong and Hussain, 1995), p,ii ( p) has the following relation: 1 ρ p = 3 i=1 λ i = 3 i=1 ˇλ i, where ˇλ 3 ε 2 a. The eigenequation of H, Φ λ2 (λ), can be expressed in terms of Ȟ in Eq. (6) as follows: Φ λ2 (λ) = Ψ(λ) ι 2 2 (ˇλ 1 λ) ι 2 1 (ˇλ 2 λ), (8) where Ψ(λ) (ˇλ 1 λ)(ˇλ 2 λ)(ˇλ 3 λ). Equation (8) is influenced by ι (ι 1, ι 2 ) which does not associate with P, and indicates Φ λ2 (ˇλ 1 ) < 0 and Φ λ2 (ˇλ 2 ), Φ λ2 (ˇλ 3 ) > 0. Then λ i ˇλ i (i = 1, 2) and ˇλ 3 λ 3, which indicates that the ˇλ 2 definition satisfies the λ 2 definition (λ 2 < 0) (Nakayama et al., 2014). ι is composed of ǎ 13 and ǎ 23 that are associated with the vorticity components parallel to P, i.e., ˇω 1 = ǎ 23 and ˇω 2 = ǎ 13, where ˇω i denotes the vorticity vector component in the ˇx i coordinate system. It indicates an important characteristic that λ i is influenced by ˇω 1 and ˇω 2. If ξ axis P, then ι = ι = o, i.e., ˇω 1 = ˇω 2 = 0. In this case, Φ λ2 (λ) is identical to Ψ(λ), which yields λ i = ˇλ i (i = 1, 2, 3). The difference between the λ 2 and ˇλ 2 definitions is the plane where the pressure minimum is specified. The ˇλ 2 definition ascertains the pressure minimum feature in P, on the other hand, λ 2 does the existence of a plane where the pressure has the minimum feature. This difference yields the characteristics of the two definitions that the ˇλ 2 definition is more severe and satisfies the λ 2 definition. 3. Numerical analysis We analyze the vortices in an isotropic homogeneous decaying turbulence by the pseudo-spectral method in the region (2π) 3, composed of nodes (Jiménez et al., 1993,Kida and Miura, 1998,Ishihara et al., 2007). For the wavenumber vector k = (k 1, k 2, k 3 ), k < 121 where k = (k i k i ) 1/2, the phase shifting method is used for dealiasing. The time step is in the fourth order Runge-Kutta method. An energy spectrum E(k) = (k/k p ) 4 exp{ 2(k/k p ) 2 } (k = k, k p = 4) (Kida and Miura, 1998) gives an initial velocity field with random phases of k, where the Taylor Reynolds number Re λ = 311, Taylor microscale λ T = 0.59, Kolmogorov length η = 0.015, and eddy turnover time T eddy = The kinematic viscosity is Hereafter ϕ, λ 2 and ˇλ 2 are nondimensionalized by their root mean square values in the corresponding time, and c is expressed as c (0 < c 1). Figure 2 shows contours of ˇλ 2 and c in vortical regions conditioned on 0 < ϕ, in an instantaneous velocity field where Re λ = 35. The regions occupied only by the contours of c (c = 0.75, 0.85) where ˇλ 2 contours (ˇλ 2 = 3, 2) do not overlap with are vortical regions where negative ˇλ 2 (pressure minimum) is not as high as those of the contours. Figure 3 shows a zoomed vortical region indicated by the arrow in Fig. 2, where the contours of λ 2 = 3, 2 are also shown. In these vortical regions, both symmetry and pressure minimum are higher in the core region. The contours of λ 2, ˇλ 2 = 2 cover those of λ 2, ˇλ 2 = 3 respectively, and of c = 0.75, Then the contours of λ 2, ˇλ 2 = 3 and c = 0.75 overlap in part, and the contour of c = 0.85 is almost inside those of λ 2, ˇλ 2 = 3. Figure 4 shows characteristics of joint probability density functions (JPDFs) of (c, λ 2 ) and (c, ˇλ 2 ) in the vortical region (0 < ϕ) where Re λ = 35, 50 and 104, respectively. It shows that λ 2 and ˇλ 2 have wide ranges in their values for an arbitrary value of c. However, both λ 2 and ˇλ 2 require high c in the high pressure minimum. These characteristics are similar irrespective of Re λ described above, although the distribution of the JPDF becomes diffusive and has a higher correlation coefficient cor in Re λ = 35. As for the development or decay of the pressure minimum, Fig. 5 shows JPDFs of c/ t and λ 2 / t as well as ˇλ 2 / t, where Re λ = 35, 50 and 104, and the time t is nondimensionalized by Kolmogorov time at the moment. cor between c/ t (6) (7) 24

5 Nakayama, Ohira and Hasegawa, Journal of Fluid Science and Technology, Vol.11, No.1 (2016) Fig. 2 Contours of λ 2 and c where c = 0.75, 0.85 and λ 2 = 2, 3 in an instantaneous velocity field where Reλ = 35. Fig. 3 Contours of c = 0.75, 0.85 and λ2, λ 2 = 3, 2 in a zoomed vortical region indicated by the arrow in Fig. 2. (a) λ2 and (b) λ 2. and λ2 / t or λ 2 / t is low in the whole vortical region. We focus on their relationships in the vortical region with low λ2 and λ 2, where the probability density seems to be concentrated as shown in Fig. 4. Figure 6 shows JPDFs of (ϕ, λ2 ) and (ϕ, λ 2 ) in the vortical region (0 < ϕ) where Reλ = 35. The distributions of the JPDFs have a shape of a quadratic curve, and, in the JPDF of (ϕ, λ 2 ), the region where ϕ 1 is the most concentrated region. Thus we extract the vortical region in terms of 0.7 < ϕ < 1.2 where the probability is approximately 0.24 (24%) and the rate of vortices satisfying λ 2 < 0 is approximately 61%. Figure 7 shows JPDFs of ( c/ t, λ2 / t) and ( c/ t, λ 2 / t) in this region. Then a feature is clearly exhibited that both λ2 and λ 2 develop or decay with c. 4. Discussion 4.1. Pressure minimum feature and flow symmetry In Fig. 4, the very weak vortices that do not have the pressure minimum feature, i.e., 0 < λ2, λ 2, have a wide range of the symmetry c. However, as positive λ2 and λ 2 is higher, the range of c becomes narrow and c is lower. Thus vortices with high positive λ2 or λ 2 have a tendency to have low symmetry. In the comparison of the regions where 0 < λ2 and 0 < λ 2, the 0 < λ 2 region (area) is larger, which indicates that vortices where 0 < λ 2 may have a pressure minimum plane except P, i.e., λ2 < 0. Even though a vortex does not have the pressure minimum feature in P, it may have that feature in the other plane, i.e., the eigenplane specified by the λ2 definition, by the effect of the vorticity components parallel to P. Then some vortices that satisfy the λ2 definition do not satisfy the λ 2 definition. As λ2 or λ 2 becomes negative and the vortices have the pressure minimum feature, vortices have a certain symmetry and seem to have a lower limit of c according to λ2 or λ 2. Furthermore, if the vortices have higher pressure minimum (higher negative λ2 or λ 2 ), then the symmetry increases. This reason is associated with the relationship between ϕ and c, as shown in Fig. 8 that is a JPDF of ϕ and c where Reλ = 35. In Fig. 8, ϕ and c have a high correlation, and high ϕ requires 25

6 Nakayama, Ohira and Hasegawa, Journal of Fluid Science and Technology, Vol.11, No.1 (2016) Fig. 4 JPDFs of (c, λ2 ) ((a), (b) and (c)) and (c, λ 2 ) ((d), (e) and (f)) in terms of 0 < ϕ, where Reλ = 35, 50, 104. (note: cor denotes the correlation coefficient.) high c to some extent. Because the high pressure minimum is derived from the high intensity of the swirling motion, i.e., high ϕ, higher pressure minimum needs higher c. In high pressure minimum vortices, because λ2 and λ 2 have a feature that λ2 < λ 2 (< 0), λ2 may seem to indicate higher pressure minimum than λ 2. However, their features of the distributions in JPDFs of (c, λ2 ) and (c, λ 2 ) are similar in these regions. We examine the difference between the pressure minimum plane Λ specified by the λ2 definition and P. Because the normal vector of Λ can be expressed by the eigenvector of λ3, say ζ 3, the angle θ between ζ 3 and the normal vector of P denoted as e 3 ξ pl η pl are analyzed in terms of λ2 and λ 2. Figure 9 shows JPDFs of (θ, λ2 ) and (θ, λ 2 ). It indicates that θ approaches 0, i.e., Λ approaches P, as λ2 and λ 2 become high negative. Then the vorticity components parallel to P become less in the high pressure minimum vortices. This feature is the reason that λ2 and λ 2 have a similar feature in their high negative regions Relationships between flow symmetry and pressure minimum in dynamics Figure 5 shows the feature that, in distributions of the JPDFs of ( c/ t, λ2 / t) and ( c/ t, λ 2 / t), high concentrated regions have downward (negative) slope shapes, although each cor is low. This feature is exhibited clearer in Fig. 7, where λ2 and λ 2 are low. This reason is also associated with the relationship between ϕ and c, that is, high ϕ requires high c. In 26

7 Nakayama, Ohira and Hasegawa, Journal of Fluid Science and Technology, Vol.11, No.1 (2016) Fig. 5 JPDFs of ( c/ t, λ2 / t) ((a), (b) and (c)) and ( c/ t, λ 2 / t ) ((d), (e) and (f)) in terms of 0 < ϕ, where Reλ = 35, 50, 104. Fig. 6 Fig. 7 JPDFs of (a) (ϕ, λ2 ) and (b) (ϕ, λ 2 ) where Reλ = 35. JPDFs of (a) ( c/ t, λ2 / t) and (b) ( c/ t, λ 2 / t ) in terms of 0.7 < ϕ < 1.2, where Reλ = 35 (cor: (a) and (b) 0.648). 27

8 Fig. 8 JPDFs of (c, ϕ) (Re λ = 35, cor = 0.838). Fig. 9 JPDFs of (a) (θ, λ 2 ) and (b) (θ, ˇλ 2 ) (Re λ = 35). addition, Fig. 8 shows that such tendency between ϕ and c is prominent in low ϕ region. Figure 10 shows the correlation coefficients of ( c/ t, λ 2 / t) and ( c/ t, ˇλ 2 / t) in terms of ϕ < 0.4, 0.8, and 1.5 < ϕ, in Re λ = 35. It is clear that the low ϕ region has high correlation with respect to ( c/ t, λ 2 / t) and ( c/ t, ˇλ 2 / t), whereas high ϕ region has lower correlation. In particular, the most concentrated probability density region in JPDF of (ϕ, λ 2 ) or (ϕ, ˇλ 2 ) in Fig. 6 has higher correlation as shown in Fig. 7, where vortices with/without the pressure minimum feature are mixed. Thus, in the state or process of a vortex attaining or losing the pressure minimum feature, c and λ 2 or ˇλ 2 have a tendency to codevelop or codecay. In Fig. 10, each range of ϕ has a higher cor than that in whole vortical region. Thus the reason of the lower cor in the whole vortical region can be derived from the different characteristics of the relationships of ( c/ t, λ 2 / t) or ( c/ t, ˇλ 2 / t) in terms of the range of λ 2 or ˇλ 2. On the other hand, it is clear that ˇλ 2 has the higher correlation with c than λ 2 in Figs. 5, 7, and 10, because the pressure minimum in P, not in other plane searched by the λ 2 definition, is directly associated with c. 5. Conclusion The relationships between the vortical flow symmetry and the pressure minimum feature in terms of the flow kinematics specified by the λ 2 and ˇλ 2 definitions were investigated in the vortices in isotropic homogeneous decaying turbulence. 28

9 Fig. 10 Correlation coefficients of ( c/ t, λ 2 / t) and ( c/ t, ˇλ 2 / t) in terms of ϕ < 0.4, ϕ < 0.8, and 1.5 < ϕ (Re λ = 35). The pressure minimum feature requires a certain flow symmetry, and the symmetry increases according to the degree of the pressure minimum. The development or decay of the pressure minimum is associated with that of the vortical flow symmetry especially in the state or process of attaining or losing the pressure minimum feature. The pressure minimum plane specified by the λ 2 definition approaches the swirl plane as the swirlity increases. It yields a characteristic that the λ 2 definition tends to be similar to the ˇλ 2 definition. References Blackburn, H. M., Mansour, N. N., and Cantwell, B. J., Topology of fine-scale motions in turbulent channel flow, J. Fluid Mech., Vol.310 (1996), pp Chong, M. S., Perry, A. E., and B. J. Cantwell, A general classification of three-dimensional flow fields, Phys. Fluids, Vol.A2 No.5 (1990), pp Cucitore, R., Quadrio, M., and Baron, A., On the effectiveness and limitations of local criteria for the identification of a vortex, E. J. Mech. B/Fluids Vol.18, No.2 (1999), pp Chakraborty, P., Balachandar, S., and Adrian, R. J., On the relationships between local vortex identification schemes, J. Fluid Mech. Vol.535 (2005), No Finn, L. I., and Boghosian, B. M., A global variational approach to vortex core identification, Physica A Vol. 362 (2006), pp Haller, G., An objective definition of a vortex, J. Fluid Mech. Vol.525 (2005), pp Hunt, J. C. R., Wray, A. A., and Moin, P., Eddies, streams, and convergence zones in turbulent flows, Center for Turbulence Research CTR-S88 (1988), pp Ishihara, T., Kaneda, Y., Yokokawa, M., Itakura, K. and Uno. A., Small-scale statistics in high-resolution direct numerical simulation of turbulence: Reynolds number dependence of one-point velocity gradient statistics, J. Fluid Mech., Vol.592 (2007), pp Jeong, J., and Hussain, F., On the identification of a vortex, J. Fluid Mech. Vol.285 (1995), pp Jiménez, J., Moin, P., Moser, R., and Keefe, L., Ejection mechanisms in the sublayer of a turbulent channel, Phys. Fluids Vol.31, No.6 (1988), pp Jiménez, J., Wray, A. A., Saffman, P. G., and Rogallo, R. S., The structure of intense vorticity in isotropic turbulence, J. Fluid Mech. Vol.255 (1993), pp Kida, S., and Miura, H., Identification and analysis of vortical structures, E. J. Mech. B/Fluids Vol.17 No.4 (1998), pp Kolář, V., Vortex identification: new requirements and limitations, Int. J. Heat Fluid Flow, Vol.28 (2007), pp Lundgren, T. S., Strained spiral vortex model for turbulent fine structure, Phys. Fluids, Vol.25 No.12 (1982), pp Nakayama, K., Physical properties corresponding to vortical flow geometry, Fluid Dyn. Res., Vol.46 (2014), Nakayama, K., Sugiyama, K., and Takagi, S., A unified definition of a vortex derived from vortical flow and the resulting pressure minimum, Fluid Dyn. Res. Vol.46 (2014), Nakayama, K., Umeda, K., Ichikawa, T., Nishihara, Y., and Takagi, S., Definition of swirl function and its application to identification method of swirling motion, Transactions of the Japan Society of Mechanical Engineers, Series B, Vol.73, No.725 (2007), pp (in Japanese). 29

10 Strawn, R. C., Kenwright, D. N., and Ahmad, J., Computer visualization of vortex wake systems, AIAA J. Vol.37 No.4 (1999), pp Tennekes, H., and Lumley, J. L., A first course in turbulence, MIT Press (1972), chap.3. Wu, J.-Z., Xiong, A.-K., and Yang, Y.-T., Axial stretching and vortex definition, Phys. Fluids Vol.17 (2005), Zhang, S. and Choudhury, D., Eigen helicity density: a new vortex identification scheme and its application in accelerated inhomogeneous flows, Phys. Fluids Vol.18 (2006), Zhou, J., Adrian, R. J., Balachandar, S., and Kendall, T. M., Mechanisms for generating coherent packets of hairpin vortices in channel flow, J. Fluid Mech. Vol.387 (1999), pp

Brief Notes on Vortex Identification

Brief Notes on Vortex Identification Brief Notes on Vortex Identification VÁCLAV KOLÁŘ Institute of Hydrodynamics Academy of Sciences of the Czech Republic Pod Patankou 30/5, 16612 Prague 6 CZECH REPUBLIC kolar@ih.cas.cz Abstract: An update

More information

Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data

Preliminary 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 information

Chuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa

Chuichi Arakawa Graduate School of Interdisciplinary Information Studies, the University of Tokyo. Chuichi Arakawa Direct Numerical Simulations of Fundamental Turbulent Flows with the Largest Grid Numbers in the World and its Application of Modeling for Engineering Turbulent Flows Project Representative Chuichi Arakawa

More information

Regularity diagnostics applied to a turbulent boundary layer

Regularity 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 information

Flow topology and non-local geometry of structures in a flat-plate turbulent boundary layer

Flow topology and non-local geometry of structures in a flat-plate turbulent boundary layer Center for Turbulence Research Proceedings of the Summer Program 2010 65 Flow topology and non-local geometry of structures in a flat-plate turbulent boundary layer By I. Bermejo-Moreno, C. Atkinson, S.

More information

Melbourne, Victoria, 3010, Australia. To link to this article:

Melbourne, 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 information

Tutorial 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 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 information

A new statistical tool to study the geometry of intense vorticity clusters in turbulence

A new statistical tool to study the geometry of intense vorticity clusters in turbulence Journal of Physics: Conference Series PAPER OPEN ACCESS A new statistical tool to study the geometry of intense vorticity clusters in turbulence To cite this article: Alberto Vela-Martin and Takashi Ishihara

More information

Citation PHYSICS OF FLUIDS (2005), 17(5)

Citation PHYSICS OF FLUIDS (2005), 17(5) Title Energy dissipation in spiral vortex straight vortex tube Author(s) Kawahara, G Citation PHYSICS OF FLUIDS (005), 17(5) Issue Date 005-05 URL http://hdl.handle.net/433/50071 Copyright 005 American

More information

Part II: Lagrangian Visualization

Part II: Lagrangian Visualization Tutorial: Visualization of Time-Varying Vector Fields Part II: Lagrangian Visualization Filip Sadlo VISUS Universität Stuttgart Germany Overview Part I: Vortices Part II: Vector Field Topology Part I:

More information

MULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES

MULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 MULTIDIMENSIONAL TURBULENCE SPECTRA - STATISTICAL ANALYSIS OF TURBULENT VORTICES

More information

Nonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto

Nonequilibrium 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 information

Kinematics of fluid motion

Kinematics of fluid motion Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a three-dimensional, steady fluid flow are dx

More information

Multiscale Computation of Isotropic Homogeneous Turbulent Flow

Multiscale Computation of Isotropic Homogeneous Turbulent Flow Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible

More information

Rortex and comparison with eigenvalue-based vortex identification criteria

Rortex and comparison with eigenvalue-based vortex identification criteria Rortex and comparison with eigenvalue-based vortex identification criteria Yisheng Gao 1 and Chaoqun Liu 1,a) 1 Department of Mathematics, University of Texas at Arlington, Arlington, Texas 76019, USA

More information

CHAPTER 5 KINEMATICS OF FLUID MOTION

CHAPTER 5 KINEMATICS OF FLUID MOTION CHAPTER 5 KINEMATICS OF FLUID MOTION 5. ELEMENTARY FLOW PATTERNS Recall the discussion of flow patterns in Chapter. The equations for particle paths in a three-dimensional, steady fluid flow are dx -----

More information

Fluid Dynamics Exercises and questions for the course

Fluid 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 information

DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT

DNS 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 information

Modeling of turbulence in stirred vessels using large eddy simulation

Modeling of turbulence in stirred vessels using large eddy simulation Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,

More information

The interaction of vorticity and rate-of-strain in homogeneous sheared turbulence

The 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 information

Flow and added small-scale topologies in a turbulent premixed flame

Flow and added small-scale topologies in a turbulent premixed flame Flow and added small-scale topologies in a turbulent premixed flame L. Cifuentes*, A. Kempf* and C. Dopazo** luis.cifuentes@uni-due.de *University of Duisburg-Essen, Chair of Fluid Dynamics, Duisburg -

More information

Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows

Contribution 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 information

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS

LARGE 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 information

Antonia, R. A., Chambers, A. J., Friehe, C. A. & Van Atta, C. W Temperature. ramps in the atmospheric surface layer. J. Atmos. Sci. 36,

Antonia, R. A., Chambers, A. J., Friehe, C. A. & Van Atta, C. W Temperature. ramps in the atmospheric surface layer. J. Atmos. Sci. 36, 129 Bibliography Andreotti, B. 1997 Studying Burgers models to investigate the physical meaning of the alignments statistically observed in turbulence. Phys. Fluids 9, 735 742. Antonia, R. A., Chambers,

More information

Turbulence Modeling I!

Turbulence 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 information

Topology of fine-scale motions in turbulent channel flow

Topology of fine-scale motions in turbulent channel flow J. Fluid Mech. (1996), vol. 310, pp. 269 292 Copyright c 1996 Cambridge University Press 269 Topology of fine-scale motions in turbulent channel flow By H U G H M. B L A C K B U R N 1, NAGI N. MANSOUR

More information

Final abstract for ONERA Taylor-Green DG participation

Final abstract for ONERA Taylor-Green DG participation 1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,

More information

The mean velocity profile in the smooth wall turbulent boundary layer : 1) viscous sublayer

The mean velocity profile in the smooth wall turbulent boundary layer : 1) viscous sublayer The mean velocity profile in the smooth wall turbulent boundary layer : 1) viscous sublayer u = τ 0 y μ τ = μ du dy the velocity varies linearly, as a Couette flow (moving upper wall). Thus, the shear

More information

Max Planck Institut für Plasmaphysik

Max Planck Institut für Plasmaphysik ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately

More information

FLOW-NORDITA Spring School on Turbulent Boundary Layers1

FLOW-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 information

Turbulent velocity fluctuations need not be Gaussian

Turbulent velocity fluctuations need not be Gaussian J. Fluid Mech. (1998), vol. 376, pp. 139 147. Printed in the United Kingdom c 1998 Cambridge University Press 139 Turbulent velocity fluctuations need not be Gaussian By JAVIER JIMÉNEZ School of Aeronautics,

More information

On the Euler rotation angle and axis of a subgrid-scale stress model

On the Euler rotation angle and axis of a subgrid-scale stress model Center for Turbulence Research Proceedings of the Summer Program 212 117 On the Euler rotation angle and axis of a subgrid-scale stress model By M. Saeedi, B.-C. Wang AND S. T. Bose In this research, geometrical

More information

A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS

A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /

More information

INTRODUCTION OBJECTIVES

INTRODUCTION OBJECTIVES INTRODUCTION The transport of particles in laminar and turbulent flows has numerous applications in engineering, biological and environmental systems. The deposition of aerosol particles in channels and

More information

Detailed Outline, M E 521: Foundations of Fluid Mechanics I

Detailed Outline, M E 521: Foundations of Fluid Mechanics I Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic

More information

Direct numerical simulation of homogeneous turbulence with hyperviscosity

Direct numerical simulation of homogeneous turbulence with hyperviscosity PHYSICS OF FLUIDS 17, 015106 (2005) Direct numerical simulation of homogeneous turbulence with hyperviscosity A. G. Lamorgese, D. A. Caughey, and S. B. Pope Sibley School of Mechanical and Aerospace Engineering,

More information

Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders

Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Toshiki Morita 1, Takashi Watanabe 2 and Yorinobu Toya 3 1. Graduate School of Information

More information

DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING

DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio

More information

Vortex Identification in the Wake of a Wind Turbine Array

Vortex Identification in the Wake of a Wind Turbine Array Portland State University PDXScholar Dissertations and Theses Dissertations and Theses Winter 3-25-205 Vortex Identification in the Wake of a Wind Turbine Array Aleksandr Sergeyevich Aseyev Portland State

More information

Turbulent energy density and its transport equation in scale space

Turbulent energy density and its transport equation in scale space PHYSICS OF FLUIDS 27, 8518 (215) Turbulent energy density and its transport equation in scale space Fujihiro Hamba a) Institute of Industrial Science, The University of Toyo, Komaba, Meguro-u, Toyo 153-855,

More information

Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600

Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600 Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation

More information

Lecture 3: The Navier-Stokes Equations: Topological aspects

Lecture 3: The Navier-Stokes Equations: Topological aspects Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without

More information

Axial motion and scalar transport in stretched spiral vortices

Axial motion and scalar transport in stretched spiral vortices PHYSICS OF FLUIDS VOLUME 13, NUMBER 9 SEPTEMBER 1 Axial motion and scalar transport in stretched spiral vortices D. I. Pullin a) Graduate Aeronautical Laboratories 15-5, California Institute of Technology,

More information

The study of vortex ring/wall interaction. for artificial nose improvement

The study of vortex ring/wall interaction. for artificial nose improvement ESAIM : Proceedings, Vol. 7, 1999, 258-269 Third International Workshop on Vortex Flows and Related Numerical Methods http://www.emath.fr/proc/vol.7/ 258 The study of vortex ring/wall interaction for artificial

More information

CHARACTERIZATION OF THE SECONDARY FLOW IN TURBULENT RECTANGULAR DUCTS WITH VARYING ASPECT RATIO

CHARACTERIZATION OF THE SECONDARY FLOW IN TURBULENT RECTANGULAR DUCTS WITH VARYING ASPECT RATIO June - July, 21 Melbourne, Australia 9 8A-1 CHARACTERIZATION OF THE SECONDARY FLOW IN TURBULENT RECTANGULAR DUCTS WITH VARYING ASPECT RATIO Ricardo Vinuesa Linné FLOW Centre, KTH Mechanics SE-1 Stockholm,

More information

Journal of Fluid Science and Technology

Journal of Fluid Science and Technology Science and Technology Effect of Molecular Diffusivities on Countergradient Scalar Transfer in a Strong Stable Stratified Flow (Study on the Linear and Nonlinear Processes by using RDT) Kouji NAGATA, Takashi

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (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 information

Concentration and segregation of particles and bubbles by turbulence : a numerical investigation

Concentration and segregation of particles and bubbles by turbulence : a numerical investigation Concentration and segregation of particles and bubbles by turbulence : a numerical investigation Enrico Calzavarini Physics of Fluids Group University of Twente The Netherlands with Massimo Cencini CNR-ISC

More information

M.A. Aziz, M.A.K. Azad and M.S. Alam Sarker Department of Applied Mathematics, University of Rajshahi-6205, Bangladesh

M.A. Aziz, M.A.K. Azad and M.S. Alam Sarker Department of Applied Mathematics, University of Rajshahi-6205, Bangladesh Research Journal of Mathematics Statistics 2(2): 56-68, 2010 ISSN: 2040-7505 Maxwell Scientific Organization, 2009 Submitted Date: November 27, 2009 Accepted Date: December 21, 2010 Published Date: June

More information

Dynamics and Statistics of Quantum Turbulence in Quantum Fluid

Dynamics and Statistics of Quantum Turbulence in Quantum Fluid Dynamics and Statistics of Quantum Turbulence in Quantum Fluid Faculty of Science, Osaka City University Michikazu Kobayashi May 25, 2006, Kansai Seminar House Contents 1. 2. 3. 4. 5. Introduction - history

More information

Development of high vorticity structures in incompressible 3D Euler equations.

Development of high vorticity structures in incompressible 3D Euler equations. Development of high vorticity structures in incompressible 3D Euler equations. D. Agafontsev (a), A. Mailybaev (b), (c) and E. Kuznetsov (d). (a) - P.P. Shirshov Institute of Oceanology of RAS, Moscow,

More information

Generation of initial fields for channel flow investigation

Generation 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 information

INVESTIGATION OF SWIRLING FLOW IN DIFFUSERS INSTALLED AT THE EXIT OF AN AXIAL-FLOW PUMP

INVESTIGATION OF SWIRLING FLOW IN DIFFUSERS INSTALLED AT THE EXIT OF AN AXIAL-FLOW PUMP TASK QUARTERLY 5 No 4 (2001), 603 610 INVESTIGATION OF SWIRLING FLOW IN DIFFUSERS INSTALLED AT THE EXIT OF AN AXIAL-FLOW PUMP ALEXEY N. KOCHEVSKY Department of Fluid Mechanics, Sumy State University, Rimsky-Korsakov

More information

Kirchhoff s Elliptical Vortex

Kirchhoff s Elliptical Vortex 1 Figure 1. An elliptical vortex oriented at an angle φ with respect to the positive x axis. Kirchhoff s Elliptical Vortex In the atmospheric and oceanic context, two-dimensional (height-independent) vortices

More information

Energy transfer and dissipation in forced isotropic turbulence

Energy transfer and dissipation in forced isotropic turbulence Energy transfer and dissipation in forced isotropic turbulence W. D. McComb, 1 A. Berera, 1 S. R. Yoffe, 2 and M. F. Linkmann 1 1 SUPA, School of Physics and Astronomy, University of Edinburgh, James Clerk

More information

The JHU Turbulence Databases (JHTDB)

The JHU Turbulence Databases (JHTDB) The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay

More information

Before we consider two canonical turbulent flows we need a general description of turbulence.

Before 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 information

NEW SCALING LAWS FOR TURBULENT POISEUILLE FLOW WITH WALL TRANSPIRATION

NEW SCALING LAWS FOR TURBULENT POISEUILLE FLOW WITH WALL TRANSPIRATION NEW SCALING LAWS FOR TURBULENT POISEUILLE FLOW WITH WALL TRANSPIRATION V.S. Avsarkisov, M. Oberlack Chair of Fluid Dynamics Technische Universität Darmstadt Petersenstr. 3, 6487 Darmstadt, Germany avsarkis@fdy.tu-darmstadt.de,

More information

EXACT 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 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 information

O. A Survey of Critical Experiments

O. 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 information

+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):

+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2): International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS

More information

Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size

Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous

More information

On the decay of two-dimensional homogeneous turbulence

On the decay of two-dimensional homogeneous turbulence On the decay of two-dimensional homogeneous turbulence J. R. Chasnov a) The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong Received 19 June 1996; accepted 3 September

More information

Numerical Simulation and Its Possibility of Application to Green Environment

Numerical Simulation and Its Possibility of Application to Green Environment Numerical Simulation and Its Possibility of Application to Green Environment Anna Kuwana 1,a, Tetuya Kawamura 2,b 1 Faculty of Science and Technology, Gunma University, 1-5-1 Tenjin-cho, Kiryu City, Gunma

More information

AER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly

AER1310: 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

Lecture 2. Turbulent Flow

Lecture 2. Turbulent Flow Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small

More information

Publication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada

Publication 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 information

Behavior of heavy particles in isotropic turbulence

Behavior of heavy particles in isotropic turbulence PHYSICAL REVIEW E 77, 637 28 Behavior of heavy particles in isotropic turbulence Jaedal Jung, Kyongmin Yeo,* and Changhoon Lee Department of Mechanical Engineering, Yonsei University, 34 Shinchon-dong,

More information

Intermittent distribution of heavy particles in a turbulent flow. Abstract

Intermittent distribution of heavy particles in a turbulent flow. Abstract Intermittent distribution of heavy particles in a turbulent flow Gregory Falkovich 1 and Alain Pumir 2 1 Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100 Israel 2 I.N.L.N., 1361

More information

IMPLEMENTATION OF POD AND DMD METHODS IN APACHE SPARK FRAMEWORK FOR SIMULATION OF UNSTEADY TURBULENT FLOW IN THE MODEL COMBUSTOR

IMPLEMENTATION OF POD AND DMD METHODS IN APACHE SPARK FRAMEWORK FOR SIMULATION OF UNSTEADY TURBULENT FLOW IN THE MODEL COMBUSTOR ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

MULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS

MULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS MULTISCALE MODEL. SIMUL. Vol. 6, No. 4, pp. 1317 1346 c 2008 Society for Industrial and Applied Mathematics MULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS

More information

Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions

Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi

More information

Proceedings of the 4th Joint US-European Fluids Engineering Division Summer Meeting ASME-FEDSM2014 August 3-7, 2014, Chicago, Illinois, USA

Proceedings of the 4th Joint US-European Fluids Engineering Division Summer Meeting ASME-FEDSM2014 August 3-7, 2014, Chicago, Illinois, USA Proceedings of the 4th Joint US-European Fluids Engineering Division Summer Meeting ASME-FEDSM4 August 3-7, 4, Chicago, Illinois, USA FEDSM4-38 SUPPRESSION OF UNSTEADY VORTEX SHEDDING FROM A CIRCULAR CYLINDER

More information

(U c. t)/b (U t)/b

(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 information

An evaluation of a conservative fourth order DNS code in turbulent channel flow

An 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 information

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran

MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER

More information

Chapter 4: Fluid Kinematics

Chapter 4: Fluid Kinematics Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to

More information

F11AE1 1. C = ρν r r. r u z r

F11AE1 1. C = ρν r r. r u z r F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions

More information

CHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.

CHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1. 1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy

More information

Lagrangian statistics in turbulent channel flow

Lagrangian statistics in turbulent channel flow PHYSICS OF FLUIDS VOLUME 16, NUMBER 3 MARCH 2004 Jung-Il Choi, a) Kyongmin Yeo, and Changhoon Lee b) Department of Mechanical Engineering and Yonsei Center for Clean Technology, Yonsei University, 134

More information

Turbulent Vortex Dynamics

Turbulent Vortex Dynamics IV (A) Turbulent Vortex Dynamics Energy Cascade and Vortex Dynamics See T & L, Section 2.3, 8.2 We shall investigate more thoroughly the dynamical mechanism of forward energy cascade to small-scales, i.e.

More information

Scale interactions and scaling laws in rotating flows at moderate Rossby numbers and large Reynolds numbers

Scale interactions and scaling laws in rotating flows at moderate Rossby numbers and large Reynolds numbers Scale interactions and scaling laws in rotating flows at moderate Rossby numbers and large Reynolds numbers P.D. Mininni NCAR, Boulder, Colorado, USA, and Departamento de Física, Facultad de Cs. Exactas

More information

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.

Periodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain. Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.

More information

6.2 Governing Equations for Natural Convection

6.2 Governing Equations for Natural Convection 6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed

More information

Turbulence - Theory and Modelling GROUP-STUDIES:

Turbulence - 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 information

Self-similar spiral flow structure in low Reynolds number isotropic and decaying turbulence

Self-similar spiral flow structure in low Reynolds number isotropic and decaying turbulence PHYSICAL REVIEW E VOLUME 54, NUMBER 1 JULY 1996 Self-similar spiral flow structure in low Reynolds number isotropic and decaying turbulence J. C. Vassilicos 1 and James G. Brasseur 2 1 Department of Applied

More information

The velocity-scalar cross spectrum of stretched spiral vortices

The velocity-scalar cross spectrum of stretched spiral vortices PHYSICS OF FLUIDS OLUME 15, NUMBER 2 FEBRUARY 23 P. A. O Gorman a) and D. I. Pullin Graduate Aeronautical Laboratories 15-5, California Institute of Technology, Pasadena, California 91125 Received 2 July

More information

WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION

WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science

More information

SOME 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 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 information

TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD

TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD Proceedings of The Second International Conference on Vortex Methods, September 26-28, 21, Istanbul, Turkey TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD Teruhiko Kida* Department of Energy Systems

More information

Turbulent boundary layer

Turbulent 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 information

Entropic Evaluation of Dean Flow Micromixers

Entropic Evaluation of Dean Flow Micromixers COMSOL Conference, Boston, 2013 Brian Vyhnalek, Petru S. Fodor and Miron Kaufman Physics Department Cleveland State University Entropic Evaluation of Dean Flow Micromixers ABSTRACT We study the use of

More information

1. Introduction, tensors, kinematics

1. 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 information

http://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error

More information

LOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT

LOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT LOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT J.-P. Laval (1) CNRS, UMR 8107, F-59650 Villeneuve d Ascq, France (2) Univ Lille Nord de France, F-59000 Lille,

More information

NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD

NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias

More information

Anisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray

Anisotropic 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 information

Quantum vortex reconnections

Quantum vortex reconnections Quantum vortex reconnections A.W. Baggaley 1,2, S. Zuccher 4, Carlo F Barenghi 2, 3, A.J. Youd 2 1 University of Glasgow 2 Joint Quantum Centre Durham-Newcastle 3 Newcastle University 4 University of Verona

More information

Homogeneous Turbulence Dynamics

Homogeneous Turbulence Dynamics Homogeneous Turbulence Dynamics PIERRE SAGAUT Universite Pierre et Marie Curie CLAUDE CAMBON Ecole Centrale de Lyon «Hf CAMBRIDGE Щ0 UNIVERSITY PRESS Abbreviations Used in This Book page xvi 1 Introduction

More information