Oscillating cylinder in viscous fluid: calculation of flow patterns and forces
|
|
- Bertina Wilson
- 5 years ago
- Views:
Transcription
1 J Eng Math (2011) 70: DOI /s Oscillating cylinder in viscous fluid: calculation of flow patterns and forces Farah Rashid Magnus Vartdal John Grue Received: 9 April 2010 / Accepted: 12 July 2010 / Published online: 5 August 2010 The Author(s) This article is published with open access at Springerlink.com Abstract Laminar and large-eddy-simulation (LES) calculations with the dynamic Smagorinsky model evaluate the flow and force on an oscillating cylinder of diameter D = 2R in otherwise calm fluid, for β = D 2 /νt in the range and Keulegan Carpenter number K = U m T/D in the range (ν kinematic viscosity, T oscillation period, U m maximal velocity). Calculations resolving the streakline patterns of the Honji instability exemplify the local flow structures in the cylinder boundary layer (β , K 2) but show that the drag and inertia force are not affected by the instability. The present force calculations conform with the classical Stokes Wang solution for all cases below flow separation corresponding to K<2 (with β<61400). The LES calculations of flow separation and vortical flow resolve the flow physics containing a large range of motion scales; it is shown that the energy in the temporal turbulent fluctuations (in fixed points) are resolved. Accurate calculation of the flow separation occurring for K>2has strong implication for the force on the cylinder. Present calculations of the force coefficients for K up to 4 and β = are in agreement with experiments by Otter (Appl Ocean Res 12: , 1990). Drag coeffients when flow separation occurs are smaller than found in U-tube experiments. Inertia coefficients show strong decline for large K (up to 8) and moderate β = 1035 but is close to unity for K = 4 and β = The finest grid has cells, finest radial r/r = , number of points along the cylinder circumference of 180, z/r = and a time step of T. Keywords Flow physics Flow separation Honji instability LES calculation Smagorinsky 1 Introduction Many papers have considered the forces on a long cylinder exposed to oscillatory flow, including, particularly, a group of classical papers by Stokes [1], Wang [2], Honji [3] and Hall [4], experimental papers by e.g. Sarpkaya [5], Otter [6], Chaplin [7], Sarpkaya [8], Johanning et al. [9] and Chaplin and Mouazé [10], and (quite) recent computational papers by Lu et al. [11], Lu and Ling [12] and Wu et al. [13]. A combined experimental and computational study was given by Dütsch et al. [14]. In the very recent paper by An et al. [15] detailed computations and flow visualizations of Honji vortices below flow separation are provided. The present paper falls into the category of computations. We compute the flow and forces on a long cylinder that oscillates in a fluid otherwise at rest. The F. Rashid M. Vartdal J. Grue (B) Mechanics Division, Department of Mathematics, University of Oslo, Oslo, Norway johng@math.uio.no
2 282 F. Rashid et al. calculations are performed for ranges of the β-number, Keulegan Carpenter number (K) and Reynolds number (Re). β is up to 61400, K up to 8 and Re up to The Navier Stokes equations are fully resolved in the laminar flow range, while large-eddy simulation (LES) is employed when the flow is turbulent. LES obtains filtered variants of the velocity and pressure fields, assuming that the filter width captures the motion scales of the inertia range of the turbulent motion [16]. Pioneering works on LES include the original Smagorinsky model of the sub-grid scale (SGS) motion [17,18]. Notable are the development and success of the dynamic Smagorinsky model, as proposed by Germano et al. [19] and further extended by Lilly [20] and Meneveau et al. [21]. Thanks to the steadily increasing computer power, the computational challenge with LES is merely how it best can be performed, rather than how it can be afforded. For recent references on SGS closure and LES, see the book by Pope [22] and the review by Pope [23]. In this paper we describe LES and laminar-flow computations solving the plain or filtered Navier Stokes equations. The computations are carried out using an extension of the CDP code (v2.5.1), which is an unstructured finite-volume based CFD code developed at Stanford University s Centre for Integrated Turbulence Simulations (CITS) over a period of 10 years. It has been applied to a variety of complex fluid-mechanical problems, including: staggered jet in cross-flow, flow in jet-engine fuel-nozzle rig and full wheel combustor simulation [24,25]. The code emphasizes good kinetic energy-conservation properties, which is a necessity for accurate turbulence simulations [26]. Implementation details can be found in [27]. The present effort represents a first step towards a long-term goal, which is to extend the LES code to free-surface motion and the interaction with marine bodies. An important example includes nonlinear wave diffraction from slender vertical cylinders (offshore wind turbines) where a typical diameter may be 5 m, wave amplitude up to 18.5 m (corresponding to the crest height of the Draupner wave), implying a Keulegan Carpenter number up to about 20. The estimate indicates a need to resolve both flow separation and wave breaking at the geometry. 1.1 Circular cylinder in sinusoidal flow Sinusoidally oscillating flow past a fixed circular cylinder has acquired a great deal of attention from fluid dynamicists and engineers because of its theoretical significance and practical importance. The forces acting on a general body are assumed to be the sum of the inertia force, viscous force and drag. The total force can be expressed in terms of Morison s equation which for a long circular cylinder of diameter D takes the form F = 4 1 πρd2 C m U ρdc d U U. (1) Here, ρ denotes the density of the fluid and U = U m sin(2πt/t) the velocity of the flow, with T the period. By (1) the force is represented by the inertia coefficient C m and drag coefficient C d. The oscillatory flow can be characterised by two parameters, viz. the Keulegan Carpenter number defined by K = U m T/D, indicating the ratio between the convective and local acceleration, and the Reynolds number defined by Re = U m D/ν where ν denotes the kinematic viscosity. Alternatively, Re is replaced by a frequency parameter β defined by β = D 2 /νt = Re/K. Reference to the Stokes number St in place of the β-number is sometimes used (St = β); see below. Stokes [1] was the first to show that the force on a fixed cylinder exposed to a sinusoidal flow of a viscous fluid is dependent on both K and Re. Wang [2], using the method of inner and outer expansions, extended Stokes analysis. The Stokes Wang solution valid for K 1, and β 1, can be reduced to the following form [ (πβ) (πβ) (πβ) 2 C d = 3π 3 2K ], (2) C m = 2 + 4(πβ) (πβ) 2 3. (3) If the cylinder is oscillating in a fluid otherwise at rest, like in the present calculations, the inertia coefficient is expressed by Ĉ m = 1 + 4(πβ) (πβ) 2 3. In an experimental glass-tank facility with fluid otherwise at rest, Honji [3] oscillated back and forth a vertical circular cylinder of height H = 30 cm and diameter D, the latter ranging between 2 and 3.8 cm, such that H/D 1.
3 Oscillating cylinder in viscous fluid 283 He described a three-dimensional instability of the flow appearing in the form of an oscillatory cell pattern in the boundary layer of the oscillating cylinder. The cylinder motion was given by 1 2 d(t) = 1 2 d 0 sin(2πf t), with 1 2 d 0 the amplitude and f = 1/T the frequency. In his original work Honji expressed the appearance of the instability in terms of the non-dimensional amplitude d 0 /D and Stokes number St = fd 2 /ν. (The St-number equals the β-number.) Investigating the range 0.4 <d 0 /D < 1.2 and 70 < St < 700 he identified the domain where the cell pattern appeared: for relatively large St, exceeding about 400, the cell pattern appeared for amplitudes in the range, about 0.4 <d 0 /D < 0.7. For smaller Stokes number, upper and lower bounds of d 0 /D both increased, however, with the same difference (d 0 /D) max (d 0 /D) min, approximately. The behavior of the cells exhibited locally a mushroom-like streakline pattern. Honji found no regular patters for d 0 /D above the upper bound where the flow exhibited flow separation and turbulence. While the present work was in progress, An et al. [15] have, by means of a finite-element method, provided detailed simulations and flow visualizations of the Honji vortices, for K = 2 and β They were able to identify steady-streaming structures induced by the Honji rolls. Performing U-tube measurements with a cylinder at rest, but oscillating the water, for β-numbers in the range between 1035 and 11240, Sarpkaya [5] has maintained that presence of the Honji instability along the cylinder boundary causes a drag force that is higher than predicted by the Stokes Wang solution (2). This has, however, not been found in experiments with a smooth cylinder oscillating in still water. With β = 61400, Otter [6] showed that the damping forces are close to the Stokes Wang solution for K<1.88 and that the drag coefficient then increases for larger K because of flow separation at the cylinder. In Otter s experiment, theoretical instability value of transition to Taylor Görtler instability [4] occurred at a Keulegan Carpenter number of The present numerical calculations of C d and C m agree with Otter s experimental values. Moreover, we find in the regime of the Honji rolls, that their appearance do not change the force coefficient C d. This is the same both prior to and after their appearance. Chaplin and Mouazé [11] have presented a collection of available experimental measurements of the drag coefficient for 100 <β< and K (about) below the theoretical stability bound of Taylor Görtler vortices given by Hall [4]. Experiments by Otter [6] and one of the data sets by Johanning [28] exhibit drag coefficients corresponding to the Stokes Wang solution (JR Chaplin (2010), personal communication). This is also true for parts of the data set of Sarpkaya [5]. An unpublished experimental data set for β = showed C d /W (W meaning the Stokes Wang C d ), and computational attempts showed similar forces (JR Chaplin (2010), personal communication). Other experimental, sets included in Chaplin and Mouazé s summary [10] showed drag coefficients twice the Stokes Wang solution [7,8,28 30]. A problem with the U-tube experiments, where the cylinder is fixed and the water oscillating, is that the input flow may become polluted by boundary-layer effects and turbulence. The scatter in C d at high values of β and small K remains an experimental and computational challenge. The paper is organized as follows. Following the introduction, Sect. 2 describes the numerical code, its implementation, discretization aspects and organization of the computations. In Sect. 3 detailed computations of the Honji instability and the particular mushroom-shaped flow patterns in the cylinder boundary layer are presented. Calculations of the forces, obtaining the inertia and drag coefficients are discussed in Sect. 4 as well as comparison to other force calculations and measurements. A short Sect. 5 discusses the highly local and three-dimensional flow pattern occurring during flow separation, exemplifying also the frequency components of the velocity field. Finally, a conclusion is given in Sect. 6. A description of the particular implementation of the dynamic Smagorinsky model is found in the Appendix. 2 Numerical approach In this investigation the motion is modeled in a frame of reference O x,y,z fixed to the oscillating cylinder. The z-axis coincides with the cylinder axis, x-axis is along the oscillatory cylinder motion, y-axis in the lateral direction, and the plane z = 0 midways on the cylinder vertical. Coordinates (x,y,z) or alternatively (x 1,x 2,x 3 ) are used.
4 284 F. Rashid et al. For the present calculations a modified version of the CDP code is employed [27]. It is an unstructured finite-volume node-based code designed for large-eddy simulation (LES) for potentially very large grids using parallel computers. For the present calculations, the incompressible Navier Stokes version of CDP is modified with a source term in order to incorporate the moving frame of reference of the cylinder. Additionally, for LES, we are using the dynamic Smagorinsky model (see Appendix). This results in the following set of equations U j = 0, (4) x j DU j Dt = ν 2 U j tr ij 1 p + Sδ 1j, (5) x i x i x i ρ x j where U j denotes the filtered velocity, p the filtered pressure and tij r the residual-stress tensor. The latter is expressed in terms of the dynamic Smagorinsky model; see the Appendix. Further, S is the source term applying along the x-direction given by S = ωu m cos(2πt/t), (6) where U m is the maximum velocity of the cylinder as defined in the Sect. 1.1 above (ω = 2π/T ). All the discretizations of the terms above are found in [27] except for the source term which is discretized according to SdV P = SV P, (7) V P where V P is the dual control volume of the node P. 2.1 Discretization aspects and organization of computations We have made fourteen runs, a k, some with different resolutions (runs b1,2 and j1,2,3), see Tables 1 and 2. The cylinder has radius R and height H ;thisis6r, however, an additional run c of the Honji instability is made with H = 8R (Table 1), as well as a convergence run k with H = 4R for high values of β and K (Table 2). The computational domain is represented by a box of dimensions (60R,20R,6R) in (x,y,z), with periodic extensions in all directions. In runs a,b,c of the Honji instability, for β = 197 and 300 (Table 1), the discretization along the vertical is z/r = 0.1. The total number of cells in these runs is In runs d-j z/r = and run k z/r = (Table 2). In the angular direction 128 points are used in runs a-j2,k and 180 in run j3. In the radial direction the grid cell nearest to the cylinder boundary is r/r = (runs a c), r/r = (d g,i,j1) and r/r = (j2,j3,k) (Table 2). The discretization along the plane z = 0 is shown in Fig. 1 (runs j2,k). Plot a shows the section of the entire computational domain, with totally sub-volumes; color plot b shows a magnification of the discretization at the cylinder boundary, including computational results of the velocities and velocity shear in the boundary layer, at (x,y,z) = (0,R,0), and plot c shows a stronger magnification. The computation shown in color Fig. 1b has parameters K = 4, β = Convergence tests are performed for K = 4, β = using a coarser grid with cells (run j1) and finer one with cells (run j3). In the runs with K = 4 the time step is T, in all other runs it is 0.001T. Table 1 Wavelength of Honji mushroom λ D Run H (β,k) Present Honji Lu and Ling An et al. a 6R (300, 2) b1,2 6R (197, 2.07) c 8R (197, 2.07)
5 Oscillating cylinder in viscous fluid 285 Table 2 Force coefficients for different values of β and K W means Stokes Wang solution given in (2 3) Run β K Ĉ m (Ĉ m + 1)/W C d C d /W a b1, d e f g h i j1,2, k The simulations were carried out at the University of Oslo on a computer with 16 cores operating at 2.13 GHz. The simulation time for each case was in the order of 1 week. 3 Simulations of the Honji rolls In the laminar range, secondary streaming is generated in the cylinder s boundary layer. For sufficiently large oscillation amplitude (and Keulegan Carpenter number), roll pattern becomes generated, superposed on this streaming, as experimentally found by Honji [3]; see Sect. 1.1, and also in recomputations by Lu and Ling [12] and An et al. [15]. Below we shall compare the present computations to those found in [11,12,15]. The roll pattern extends vertically along the cylinder boundary and has a maximum in the plane x = 0. The roll pattern is periodic with period 2 1 T, with maxima at 2πt/T = 0, π, 2π etc. The cycle of the z-velocity w is visualized in color Fig. 2 for K = 2.07, β = 197 (finest run b in Table 1), showing the rolls for 2πt/T = 0, π/4, π/2. At the latter time the rolls almost disappear before gaining again, keeping the polarity. This computation is well within the domain where the cell pattern appeared in Honji s experiments, with range 0.4 <d 0 /D < 1.2 (corresponding to 1.26 <K<3.77) and 70 <β= St < 700. Putting the perturbation velocities in the plane x = 0ontheform(v, w) = ( ψ/ z, ψ/ y) the stream function ψ may be integrated, where ψ = constant gives streak lines of the mushroom-shaped instability as visualized in color Fig. 3. The calculations compare well to the experimental measurements by Honji [3] (Fig. 4). The height of the mushroom is h = 0.65R; the outer width is W o = 0.62R and inner width W i = 0.4R. With a number of five mushrooms appearing along the cylinder length of 6R, the wavelength of the instability becomes λ/d = 0.6 inthis case. We have carried out four runs simulating the Honji instability as indicated in Table1 by runs a, b1, b2 and c, for (β, K) values of (300, 2) and (197, 2.07). The difference between runs b1 and b2 is that two different resolutions are used. In run c the cylinder of height 8R is used; in this case the wavelength of the instability becomes λ/d = Both computations agree with Honji s observation of the wavelength, obtained for two cylinders of diameters 2 and 3.8 cm, both of height 30 cm; see his Fig. 14. The calculations by Lu and Ling [12] show a wavelength of 0.59 when their cylinder becomes tall, while An et al. [15] have found λ/d , see their Fig Force calculations The pressure and shear force are integrated over the cylinder surface during the simulation time which is up to 23T, as illustrated in Fig. 5a forrunb1with(β, K) = (197, 2.07). After about t = 2T the force is almost repeatable before this time small transients are present. After a time t = 14T the small instabilities appear at the surface
6 286 F. Rashid et al. Fig. 1 Example of discretization. Plain cut of mid-section. Total number of cells Arrows indicate cylinder motion. a Complete lateral extension. b Cells at the cylinder in moderate close-up. Also shown are velocities (u, v, w)/u along the radial direction, at the widest part of the cylinder. c Cells at the cylinder in strong close-up. Innermost cell has r/r = of the body; in the present calculations they first appear at the intersection between the cylinder surface and the periodic boundaries at z =±3R. The instabilities grow along the cylinder surface towards z = 0, from both sides, developing into the rolls and mushroom-shaped instabilities discussed above. An important result of the simulations is that the presence of the Honji instability does not change the drag coefficient. Apart from some initial transients, the value of C d is the same both before and after the instability appears. Figure 5b illustrates the force history for runs j1 and j2 for time up to 13T, both with (β, K) = (11240, 4).Run j1 marked with the dashed line has cells, while j2 marked with the solid line has cells. Apart from a small difference during the initial phase, the force histories are identical. Due to the separation taking place at the cylinder, the flow is not entirely symmetric, causing a lateral force acting on the cylinder. This is shown in
7 Oscillating cylinder in viscous fluid 287 Fig. 2 The cycle of the z-velocity w in x = 0plane, at the total z-direction of the cylinder. K = 2.07, β = 197. t = 0(a), 2πt/T = π/4(b) and 2πt/T = π/2 (c) (a) (b) (c) Fig. 5c for runs j1 and j2. It is noted that the transverse force is much smaller than the in-line force. More precisely, F y /W has rms of and relative to the max F x /W in runs j1 and j2, respectively, for (β, K) = (11240, 4). Altogether, we have performed fourteen different computational runs evaluating the flow and force characteristics and thereby extracting the force coefficients according to the Morison equation (1); see Table 2. We have already noted one surprise, namely for the runs in the Honji instability range, where the drag force remained unchanged and close to the Stokes Wang solution. This has also been true for the other runs, for the larger values of β, but with K small enough to avoid flow separation (K <2). In this range present computed drag coefficients are close to the Stokes Wang solution, as visualized in Fig. 5d. For the runs when flow separation is dominant, three runs e,f,g, all with β = 1035, and K = 4.86, 6.28, 8show drag coefficients C d = 0.919, , and inertia coefficients Ĉ m = 0.856, 0.780, 0.566, respectively (Re up to 8280), i.e., a strong increase in C d and a strong decline in Ĉ m (Table 2; Fig.5e). In runs j1,2,3, with the much higher β = and K = 4 (and Re = 45000), the present calculation of the drag coefficient becomes C d = This C d -value is in excellent agreement with the measurements of C d performed by Otter [6] of cylinders oscillating in water otherwise at rest, finding C d = 0.31 for (β, K) = (61400, 3.63) and C d = 0.44 for (β, K) = (61400, 4.22). In runs j1,2,3 with (β, K) = (11240, 4) there is not a big drop in the inertia coefficient, which becomes Ĉ m = Comparison to other force calculations and measurements Results of the drag coefficient, particularly in two publications, are relevant to compare with and discuss. These are the LES simulations by Lu et al. [11] and the U-tube experiments by Sarpkaya [5]. In Lu et al. they simulated
8 288 F. Rashid et al. Fig. 3 Mushroom vortices corresponding to the Honji instability in x = 0 plane, along the total z-direction of the cylinder. K = 2.07, β = 197. t = 0 (a), 2πt/T = π/4(b) and 2πt/T = π/2 (c) (a) (b) (c) Fig. 4 a Expansion of a single roll in Fig. 3a. Number of grid cells b Close-up view of streaked flow on cylinder. K = 2.07, β = 197. Adapted from [3, Fig. 10] (a) (b)
9 Oscillating cylinder in viscous fluid 289 Fig. 5 Force history (along x-direction) normalized by Stokes Wang solution (inertia coefficient) a runb1with(β, K) = (197, 2.07); b run j1 (dashed); run j2 (solid) with(β, K) = (11240, 4); c lateral force along y-direction in run j1 (dashed); run j2 (solid); d drag coefficient divided by Stokes Wang solution denoted by W versus K. Present with β = (large circle); β = (triangle); β = 1035 (square); β = 300 and 197 (dot); Otter [6] β = (crosses); e Ĉ m versus K. Same symbols as in d
10 290 F. Rashid et al. particularly the data set from Sarpkaya with β = 1035 finding good agreement between the simulations and the measurements. Apart from results fitting to the Stokes Wang solution, we do not find here the same agreement with the results in Sarpkaya and thus neither with Lu et al. An important disagreement is that the present computations conform to the Stokes Wang solution until separation takes place. In the flow-separation range we compare the present computed C d = 0.36 to Sarkaya s measured C d = 0.5 for (β, K) = (11240, 4), his Fig. 3. We also compare present C d = 0.919, , for K = 4.86, 6.28, 8, respectively, with β = 1035, to Sarpkaya s measured C d s in his Fig. 1. It is evident that his C d -values exceed ours by 40% (i.e., present values multiplied by 1.4 give Sarpkaya s). We note that the inertia coefficients calculated here fit rather well to Sarpkaya s measurements, however. We note a fundamental difference between U-tube experiments such as those performed by Sarpkaya [5] oscillating the water while keeping the cylinder fixed, and the experiments by Otter [6] keeping the water at rest, oscillating the cylinder: it may be speculated if the input oscillatory channel flow gets polluted by turbulence in the former experiments. Present simulations fit well to the experiments by Otter. We may compare the present computations of the lateral force for β = and K = 4 to the computations by Lu et al. for β = 1035 and K = 3, 5, their Figs. 3 and 7, respectively. While they find a lateral force up to 16% of F x for K = 3, and up to and even exceeding 50% of F x for K = 5, we find here a rms of F y of only 2% of F x. This is another poor comparison between the present LES calculations and those by Lu et al. We note that in the latter computations the original Smagorinsky model was used with the Smagorinsky constant put equal to 0.1, as also used in a lot of similar calculations. However, it is evident from several publications that the dynamic Smagorinsky model quite successfully improves the representation of the turbulent motion which is the case when flow separation occurs; see e.g. [23]. Another difference between our and Lu et al. is that here a finer resolution is employed; presently grid cells (runs d k) compared to Lu et al. s grid cells. 5 Local flow structures Obtaining the Honji rolls exemplifies the calculation of the local flow structures. The LES calculations at moderate and high K obtain the velocity field during flow separation. There then is strong vortical motion in a large domain surrounding the cylinder. For K = 4 this domain extends 9 by 4 diameters along the x- and y-directions, respectively. For K = 2 the corresponding extension is 5 by 1 diameters (Fig. 6). The vorticity close to the cylinder is about five times U m /R (Fig. 7), while in the boundary-layer, the vorticity the velocity shear is up to 200 times U m /R. The boundary-layer thickness in this calculation for K = 4,β = is 0.01R while the grid cell nearest to the wall extends R radially (Fig. 1b). During flow separation, the velocity and vorticity fields contain a large range of motion scales which are important to represent accurately. The flow separation involving processes on the small scales has to be computed accurately in order to obtain the force well, where the latter is an example of a global quantity. We exemplify the frequency content of the velocity field by computing the time history of w-velocity in a fixed points (x,y,z) = (x 0,y 0,z 0 ). The auto-correlation function of w(t; x 0,y 0,z 0 )/U m and its Fourier transform are evaluated as exemplified in Fig. 8. Results are obtained in two points, for y 0 = 1.02R and y 0 = 3R with (x 0,z 0 ) = (0, 0) is both cases. This clearly show that the energy in the temporal turbulent fluctuations are resolved in the computations. 6 Conclusions Laminar and LES calculations of the flow and force on an oscillating cylinder in otherwise calm fluid have been performed for β-number ranging from 197 to The LES code is a modification of CDP, the latter used also for other purposes; it employs the dynamic Smagorinsky model. The computational domain was 30 diameters along the oscillation direction and 10 diameters laterally. The cylinder height H which equals the height of the computational
11 Oscillating cylinder in viscous fluid 291 Fig. 6 Vorticity v /U /R. Cut along mid-section z = 0. β = KC = 2(a) and4(b) (a) (b) Fig. 7 Same as Fig. 6 but x = 0 plane, and 3R <z<3r, R<y<4R (a) (b) domain was 6R; in two separate runs H is 8R and 4R. The computational domain extended periodically in all directions. We have found as main result that the force coefficients follow (or are very close to) the analytical Stokes Wang solution for Keulegan Carpenter number K up to 2 where flow separation occurs. The drag coefficient then departs from the Stokes Wang solution. Our calculations of the force coefficients are in agreement with the experiments by Otter [6]. He obtained the results for β = while computations here are performed for a range in β,upto 61400, as indicated above. Present calculations of flow separation, performed for K up to 4 and β = 11240, with corresponding Re = 45000, exhibit a drag coefficient of C d = 0.36, comparing well to the experiments by Otter [6], finding C d = 0.31 for K = 3.63 and β = 61400, and C d = 0.44 for K = 4.22 and β = Convergence calculations were performed for K = 4 and β = as well as for the laminar range. Supplementary computations with smaller β = 1035 and K = 4.86, 6.28, 8 (i.e., Re up to 8280) exhibit C d = 0.919, ,
12 292 F. Rashid et al. Fig. 8 Temporal spectrum of w-velocity in (x 0,z 0 ) = (0, 0) and a y 0 = 1.02R, b y 0 = 3R All present drag coefficient calculations are about 40% smaller than C d s measured by Sarpkaya [5] when flow separation occurs. Inertia coefficients calculated here fit rather well to Sarpkaya s measurements, however. A surprising result of the laminar calculations is that the force coefficients are the same both before and after the appearance of the Honji instability. They agree with the Stokes Wang solution (they are 10% higher) but disagree with U-tube experiments obtaining a doubling of the drag coefficient as compared to the Stokes Wang solution. Particularly Sarpkaya [5,8] has pointed to the Honji instability being the cause of this increase in the drag coefficient. In the present calculations with (β, K) = (197, 2.07) and (300, 2) it takes some time until the Honji instability becomes present; it grows from small instabilities at the intersection between the cylinder geometry and the (periodic) walls at the top and bottom (of the cylinder). The instability then grows along the cylinder side, and becomes fully developed after a number of 14 periods, in the present calculations. Wavelength, height and width of the mushroom-shaped streak patterns induced by the Honji rolls fit well to the original experiments [3]. The periodic appearance and disappearance of the rolls are visualized. Fine and moderate grid resolutions show the same streak pattern. Several experimental sets at high β and K below Hall s stability limit exhibit drag coefficients corresponding to the Stokes Wang solution. These include [6] and one of the data sets by Johanning [28]. The same has been found in an unpublished experimental data set for β = (JR Chaplin (2010), personal communication), and also in other computational attempts than those presented here (JR Chaplin (2010), personal communication). However, Otter s experiments [6] conform also to the Stokes Wang solution well beyond the stability limit, up to where flow separation occurs, at K = 2. These experiments physical and numerical conform with the present calculations. Other experimental sets at high β and below the stability limit [7,8,28 30] show drag coefficients twice the Stokes Wang solution; see the summary by Chaplin and Mouazé s [10]. We note that a problem with the (physical) U-tube experiments, where the cylinder is fixed and the water oscillating, is that the input flow may become polluted by boundary-layer effects and turbulence. The scatter in C d at high values of β and small K remains an experimental and computational challenge. Calculations of the Honji rolls exemplify the local flow structures. The present calculations show that the pattern has no consequence for the force on the cylinder, however. The LES calculations of flow separation and vortical flow resolve the flow physics at small scale. The calculations exhibit velocity and vorticity fields containing a large range of motion scales. Evaluation of the frequency contents of the velocity field shows that the energy of the temporal turbulent fluctuations (in fixed points) is resolved in the computations. It is evidenced that the flow separation occurring for K>2 has strong implication for the force on the cylinder. An accurate numerical representation of the flow physics during the separation is needed in order to accurately calculate the force. This implies a relevant modeling of the scales of the flow separation process, that also is three-dimensional.
13 Oscillating cylinder in viscous fluid 293 The finest grid used had a number of cells; it had a finest radial r/r = , number of points along the cylinder circumference of 180, z/r = and time step of T, for the case K = 4, β = The simulations presented here were carried out on a computer with 16 cores operating at 2.13 GHz. The simulation time for each case was of the order of 1 week. Acknowledgements This paper is presented in memory of Professor E. O. Tuck who, with his brilliant analytical approach to a large variety of problems, particularly has inspired the international environments of Applied Mathematics and Hydrodynamics. This research was funded by the Research Council of Norway through NFR /V30. Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited. Appendix: The dynamic Smagorinsky model In this section we describe how the dynamic Smagorinsky model (DSM) is implemented in the present numerical code. The theoretical descriptions of the DSM is given in e.g. [19,20] and [22]. As our starting point for the description of the dynamic Smagorinsky model we consider the plain Smagorinsky model in which the deviatoric part of the residual stress tensor is approximated by t r ij = 2c S 2 S S ij, (8) where is the filterwidth S = (2S ij S ij ) 1/2 is the characteristic rate of strain, S ij is the filtered rate of strain and c S is the square of the Smagorinsky constant. In the Smagorinsky model c S is, as the name implies, a constant, e.g. in theworkofluetal.[11] the Smagorinsky constant is put equal to 0.1, as also used in a lot of similar calculations. In general, however, we must expect that this is not the case; c S must be viewed as a function of both space and time. The main idea of the DSM is to apply a second (test) filter, of width ˆ, to extract the information about the spatial variability of c S. Usually ˆ 2. Applying the test filter to the filtered N S equations results in a similar equation, but with a different residual stress expression, T ij = ÛiU j ˆ U i ˆ U j, (9) We approximate the deviatoric part of T ij in the Smagorinsky manner so that we get, Tij r = 2c ˆ 2 S Ŝ Ŝij, (10) If we now consider the difference between (10) and the test filtered (8) we obtain the following if we, quite inconsistently, assume c S to be uniform [22], ( L S ij = T ij r t ij r = c S 2 2 Ŝ S ij 2 ˆ 2 ) Ŝ Ŝij = c S M ij. (11) Here L S ij is an approximation to the deviatoric part of the Leonard stresses, where the Leonard stresses are the difference between the two residual stress tensors, L ij = T ij t ij = Ûi U j ˆ U i ˆ U j. As we can see the terms in (12) are known in terms of the filtered velocities and can thus be computed. This means that we can calculate a value of c S that minimizes the difference between L S ij and the deviatoric part of L ij.this is done by minimizing the mean-square error [20] according to, (12) c S = M ij L ij M kl M kl. (13)
14 294 F. Rashid et al. In the present numerical code the eddy viscosity ν T = c S 2 S is computed using the filtered velocity field from the previous timestep. The coefficient c S is computed according to c S 2 = M ij L ij, M kl M kl where L ij is as given above. In this particular application M ij is modified into M ij = 2 (Ŝ Sij ˆ 2 ) / 2 Ŝ Ŝij. (14) (15) The local averaging obtained by the filtering in (14) is necessary for stability [19]. As an additional stabilising step all negative values of eddy viscosity are set to zero. Finally, to complete the description of the model the filters and their widths need to be specified. For the width filter the filtering is implicit and the filterwidth is taken to be proportional to the mesh spacing, i.e., = (V ) 1/3, where V is the control volume. The ˆ filter is a weighted average of the value at the node and its neighbours. The filter width ratio is taken to be ˆ / = 6. References 1. Stokes GG (1851) On the effect of the internal friction of fluids on the motion of pendulums. Trans Camb Philos Soc 9: Wang C-Y (1968) On high-frequency oscillatory viscous flows. J Fluid Mech 32: Honji H (1981) Streaked flow around an oscillating circular cylinder. J Fluid Mech 107: Hall P (1984) On the stability of the unsteady boundary layer on a cylinder oscillating transversely in a viscous fluid. J Fluid Mech 146: Sarpkaya T (1986) Force on a circular cylinder in viscous oscillatory flow at low Keulegan-Carpenter numbers. J Fluid Mech 165: Otter A (1990) Damping forces on a cylinder oscillating in a viscous fluid. Appl Ocean Res 12: Chaplin JR (2000) Hydrodynamic damping of a cylinder at β J Fluids Struct 14: Sarpkaya T (2001) Hydrodynamic damping and quasi-coherent structures at large Stokes numbers. J Fluids Struct 15: Johanning L, Bearman PW, Graham JMR (2001) Hydrodynamic damping of a large scale surface piercing circular cylinder in planar oscillatory motion. J Fluids Struct 15: Chaplin JR, Mouazé D (2004) Disturbed laminar flow over an oscillating cylinder. In: XXI international congress on theoretical and applied mechanics, Aug 2004, Warsaw, Poland. pdf 11. Lu X, Dalton C, Zhang J (1997) Application of large eddy simulation to an oscillating flow past a circular cylinder. J Fluids Eng 119: Lu X-Y, Ling GC (2003) Three-dimensional instability of an oscillating viscous flow past a circular cylinder. Appl Math Mech (English edition) 24: Wu J-Z, Lu X-Y, Zhuang L-X (2007) Integral forces acting on a body due to local flow structures. J Fluid Mech 576: Dütsch H, Durst F, Becker S, Lienhart H (1998) Low-Reynolds-number flow around an oscillating circular cylinder at low Keulegan-Carpenter numbers. J Fluid Mech 360: An H, Cheng L, Zhao M (2010) Direct numerical simulation of oscillatory flow around a circular cylinder at low KC number. J Fluid Mech (in press) 16. Leonard A (1975) Energy cascade in large-eddy simulations of turbulent flows. Adv Geophys 18(1): Smagorinsky J (1963) General circulation experiments with the primitive equations. Mon Weather Rev 91(3): Lilly DK (1967) The representation of small-scale turbulence in numerical simulation experiments. In: Goldstine HH (ed) Proceedings of IBM scientific computing symposium on environmental sciences. IBM, Yorktown Heights, NY, pp Germano M, Piomelli U, Moin P, William HC (1991) A dynamic subgrid-scale eddy viscosity model. Phys Fluids A Fluid Dyn 3: Lilly DK (1992) A proposed modification of the Germano subgrid-scale closure method. Phys Fluids A Fluid Dyn 4: Meneveau C, Lund TS, Cabot WH (1996) A Lagrangian dynamic subgrid-scale model of turbulence. J Fluid Mech 319: Pope SB (2000) Turbulent flows. Cambridge University Press, Cambridge 23. Pope SB (2004) Ten questions concerning the large-eddy simulation of turbulent flows. New J Phys 6:35. doi: / / 6/1/ CITS (2004) Annual technical report. Center for Integrated Turbulence Simulations, Stanford. pdf/reports/cits2004.pdf
15 Oscillating cylinder in viscous fluid CITS (2005) Annual technical report. Center for Integrated Turbulence Simulations, Stanford. pdf/reports/cits2005.pdf 26. Ham F, Iaccarino G (2004) Energy conservation in collocated discretization schemes on unstructured meshes. Centre for Turbulence Research, Stanford Ham F, Mattasson K, Iaccarino G (2006) Accurate and stable finite volume operators for unstructured flow solvers. Centre for Turbulence Research, Stanford. Annual Research Briefs Johanning L (2003) PhD thesis. Imperial College, London 29. Bearman PW, Russel MP (1996) Measurements of hydrodynamic damping of bluff bodies with application to the prediction of viscous damping of TLP hulls. In: 21st symposium on naval hydrodynamics, Trondheim, Norway, pp Chaplin JR, Subbiah K (1998) Hydrodynamic damping of a cylinder in still water and in a transverse current. Appl Ocean Res 20:
Submitted for publ. in Journal of Engineering Mathematics, Tuck Special Issue, 9 April 2010
Submitted for publ. in Journal of Engineering Mathematics, Tuck Special Issue, 9 April 2010 Oscillating cylinder in viscous fluid: calculation of flow patterns and forces Farah Rashid, Magnus Vartdal and
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationTransactions on Modelling and Simulation vol 16, 1997 WIT Press, ISSN X
Numerical and experimental investigation of oscillating flow around a circular cylinder P. Anagnostopoulos*, G. Iliadis* & S. Kuhtz^ * University of Thessaloniki, Department of Civil Engineering, Thessaloniki
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationLES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE
LES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE By Chin-Tsau Hsu, 1 Xiyun Lu, and Man-Kim Kwan 3 ABSTRACT: Oscillatory flows over a flat plate are studied by using Large Eddy Simulation (LES)
More informationDirect Numerical Simulations on the Uniform In-plane Flow around an Oscillating Circular Disk
Proceedings of the Twenty-third (2013) International Offshore and Polar Engineering Anchorage, Alaska, USA, June 30 July 5, 2013 Copyright 2013 by the International Society of Offshore and Polar Engineers
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationModeling 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 informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationPredicting natural transition using large eddy simulation
Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationVortex wake and energy transitions of an oscillating cylinder at low Reynolds number
ANZIAM J. 46 (E) ppc181 C195, 2005 C181 Vortex wake and energy transitions of an oscillating cylinder at low Reynolds number B. Stewart J. Leontini K. Hourigan M. C. Thompson (Received 25 October 2004,
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationPrediction of unsteady heat transfer from a cylinder in crossflow
Center for Turbulence Research Proceedings of the Summer Program 202 07 Prediction of unsteady heat transfer from a cylinder in crossflow By S. T. Bose, B. C. Wang AND M. Saeedi The accuracy of a tensorial
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationModule 3: Velocity Measurement Lecture 15: Processing velocity vectors. The Lecture Contains: Data Analysis from Velocity Vectors
The Lecture Contains: Data Analysis from Velocity Vectors Velocity Differentials Vorticity and Circulation RMS Velocity Drag Coefficient Streamlines Turbulent Kinetic Energy Budget file:///g /optical_measurement/lecture15/15_1.htm[5/7/2012
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 informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationFlow Structure Investigations in a "Tornado" Combustor
Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas
More informationStatistical 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 informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationSimulation 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 informationApplied Mathematics and Mechanics (English Edition)
Appl. Math. Mech. -Engl. Ed., 39(9), 1267 1276 (2018) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-018-2364-7 Direct numerical simulation of turbulent flows through
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
More informationMultiscale 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 informationLarge-eddy simulation of an industrial furnace with a cross-flow-jet combustion system
Center for Turbulence Research Annual Research Briefs 2007 231 Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system By L. Wang AND H. Pitsch 1. Motivation and objectives
More informationS. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives
Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationLES and unsteady RANS of boundary-layer transition induced by periodically passing wakes
Center for Turbulence Research Proceedings of the Summer Program 2 249 LES and unsteady RANS of boundary-layer transition induced by periodically passing wakes By F. E. Ham, F. S. Lien, X. Wu, M. Wang,
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 informationABSTRACT INTRODUCTION
Numerical simulation of pulsating flow around a cube C. Dargent, D. Dartus, J. George Institut de Mecanique des Fluides de Toulouse, Avenue du Professeur Camille Soula, 31400 Toulouse, France ABSTRACT
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationIMPACT OF LES TURBULENCE SUBGRID MODELS IN THE JET RELEASE SIMULATION
IMPACT OF LES TURBULENCE SUBGRID MODELS IN THE JET RELEASE SIMULATION E. S. FERREIRA JÚNIOR 1, S. S. V. VIANNA 1 1 State University of Campinas, Faculty of Chemical Engineering E-mail: elmo@feq.unicamp.br,
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationJet pumps for thermoacoustic applications: design guidelines based on a numerical parameter study
arxiv:158.5119v1 [physics.flu-dyn] 2 Aug 215 Jet pumps for thermoacoustic applications: design guidelines based on a numerical parameter study Jet pump design guidelines Joris P. Oosterhuis a), Simon Bühler,
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationA NEW MODEL FOR ESTIMATING NEUTRAL PLANE IN FIRE SITUATION
A NEW MODEL FOR ESTIMATING NEUTRAL PLANE IN FIRE SITUATION JY Zhang¹,*, Jane WZ Lu² and R Huo¹ 1 PhD student, State Key Laboratory of Fire Science, University of Science and Technology of China, Hefei,
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More 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 informationNumerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation
, pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
More informationLES ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT
2th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics ANALYSIS ON CYLINDER CASCADE FLOW BASED ON ENERGY RATIO COEFFICIENT Wang T.*, Gao S.F., Liu Y.W., Lu Z.H. and Hu H.P. *Author
More informationNumerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder
Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More information6. Basic basic equations I ( )
6. Basic basic equations I (4.2-4.4) Steady and uniform flows, streamline, streamtube One-, two-, and three-dimensional flow Laminar and turbulent flow Reynolds number System and control volume Continuity
More informationarxiv: v1 [physics.flu-dyn] 11 Oct 2012
Low-Order Modelling of Blade-Induced Turbulence for RANS Actuator Disk Computations of Wind and Tidal Turbines Takafumi Nishino and Richard H. J. Willden ariv:20.373v [physics.flu-dyn] Oct 202 Abstract
More informationMODELLING THE INTERACTION BETWEEN WATER WAVES AND THE OSCILLATING WATER COLUMN WAVE ENERGY DEVICE. Utku Şentürk, Aydoğan Özdamar
Mathematical and Computational Applications, Vol. 16, No. 3, pp. 630-640, 2011. Association for Scientific Research MODELLING THE INTERACTION BETWEEN WATER WAVES AND THE OSCILLATING WATER COLUMN WAVE ENERGY
More informationLES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2
MAGNETOHYDRODYNAMICS Vol. 00 (1964), No. 00, pp. 1 5 LES modeling of heat and mass transfer in turbulent recirculated flows E. Baake 1, B. Nacke 1, A. Umbrashko 2, A. Jakovics 2 1 Institute for Electrothermal
More informationStudy on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time
World Journal of Mechanics, 2018, 8, 301-310 http://www.scirp.org/journal/wjm ISSN Online: 2160-0503 ISSN Print: 2160-049X Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration
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 informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationCHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM
CHARACTERISTIC OF VORTEX IN A MIXING LAYER FORMED AT NOZZLE PITZDAILY USING OPENFOAM Suheni and Syamsuri Department of Mechanical Engineering, Adhi Tama Institute of Technology Surabaya, Indonesia E-Mail:
More informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
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):
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 informationComputers and Mathematics with Applications. Investigation of the LES WALE turbulence model within the lattice Boltzmann framework
Computers and Mathematics with Applications 59 (2010) 2200 2214 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Investigation
More informationResistance in Unsteady Flow Search for a Model
Resistance in Unsteady Flow Search for a Model T. Sarpkaya Department of Mechanical Engineering 700 Dyer Road, Rm: 339 (ME-SL) Naval Postgraduate School Monterey, CA 93943-5100 phone: (408) 656-3425 fax:
More informationVortex Induced Vibrations
Vortex Induced Vibrations By: Abhiroop Jayanthi Indian Institute of Technology, Delhi Some Questions! What is VIV? What are the details of a steady approach flow past a stationary cylinder? How and why
More informationNatural frequency analysis of fluid-conveying pipes in the ADINA system
Journal of Physics: Conference Series OPEN ACCESS Natural frequency analysis of fluid-conveying pipes in the ADINA system To cite this article: L Wang et al 2013 J. Phys.: Conf. Ser. 448 012014 View the
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 informationNUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven
More informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationDesign and Modeling of Fluid Power Systems ME 597/ABE Lecture 7
Systems ME 597/ABE 591 - Lecture 7 Dr. Monika Ivantysynova MAHA Professor Fluid Power Systems MAHA Fluid Power Research Center Purdue University Content of 6th lecture The lubricating gap as a basic design
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 informationNumerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow
Fluid Structure Interaction VII 147 Numerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow H. G. Sung 1, H. Baek 2, S. Hong 1 & J.-S. Choi 1 1 Maritime and Ocean
More informationO. A Survey of Critical Experiments
O. A Survey of Critical Experiments 1 (A) Visualizations of Turbulent Flow Figure 1: Van Dyke, Album of Fluid Motion #152. Generation of turbulence by a grid. Smoke wires show a uniform laminar stream
More informationDynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh
40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationThe effect of top tension on VIV model analysis of a vertical flexible riser
The Second Conference of Global Chinese Scholars on Hydrodynamics The effect of top tension on VIV model analysis of a vertical flexible riser Muyu Duan 1,2, Bowen Fu 1, Decheng Wan 1* 1 State Key Laboratory
More informationSimulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems
Center for Turbulence Research Annual Research Briefs 27 241 Simulation of a lean direct injection combustor for the next high speed civil transport (HSCT) vehicle combustion systems By H. El-Asrag, F.
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationAutomatic Eddy Viscosity Assignment for 2-D Hydrodynamic Model of Szczecin Bay
PUBLS. INST. GEOPHYS. POL. ACAD. SC., E-6 (390), 006 Automatic Eddy Viscosity Assignment for -D Hydrodynamic Model of Szczecin Bay Ryszard EWERTOWSKI 1, 1 Technical University of Szczecin, Building and
More informationLarge-Eddy Simulation of the Lid-Driven Cubic Cavity Flow by the Spectral Element Method
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2006) DOI: 10.1007/s10915-005-9039-7 Large-Eddy Simulation of the Lid-Driven Cubic Cavity Flow by the Spectral Element Method Roland Bouffanais,
More informationA Comparative Analysis of Turbulent Pipe Flow Using k And k Models
A Comparative Analysis of Turbulent Pipe Flow Using k And k Models 1 B. K. Menge, 2 M. Kinyanjui, 3 J. K. Sigey 1 Department of Mathematics and Physics. Technical University of Mombasa, P.O BOX 90420-80100,Mombasa,
More informationNatural convection adjacent to a sidewall with three fins in a differentially heated cavity
ANZIAM J. 48 (CTAC2006) pp.c806 C819, 2007 C806 Natural convection adjacent to a sidewall with three fins in a differentially heated cavity F. Xu 1 J. C. Patterson 2 C. Lei 3 (Received 31 August 2006;
More informationDatabase analysis of errors in large-eddy simulation
PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)
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 informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationNUMERICAL 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 informationINTRODUCTION 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 informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS
ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS Karsten Lindegård Jensen 1, B. Mutlu Sumer 1, Giovanna Vittori 2 and Paolo Blondeaux 2 The pressure field in an oscillatory boundary layer
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationContents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1
Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors
More information