DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
|
|
- Felicity Byrd
- 5 years ago
- Views:
Transcription
1 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 of Mechanical Engineering University of Manitoba Winnipeg, MB, R3T 5V6, Canada BingChen.Wang@umanitoba.ca Xingjun Fang Department of Mechanical Engineering University of Manitoba Winnipeg, MB, R3T 5V6, Canada fangx@myumanitoba.ca ABSTRACT In this paper, direct numerical simulation (DNS) is performed to study the effects of the Coriolis force on turbulent heat transfer within a square duct subjected to spanwise system rotation. In order to investigate the rotation effects on the heat transfer process, a wide range of rotation numbers have been tested. In response to the system rotation, secondary flows appear as large streamwise counterrotating vortices, which interact intensely with the four boundary layers of the square duct and have a significant impact on the statistics of the coupled velocity and temperature fields. It is observed that at sufficiently high rotation numbers, a Taylor-Proudman region appears and the flow approaches a complete laminarization state near the top and side walls. The influence of large organized secondary flows on the turbulent kinetic energy of the velocity field and turbulent scalar energy of the temperature field have been investigated through a spectral analysis. The effects of system rotation on the production of the temperature variance have been examined through a budget analysis of its transport equation. INTRODUCTION Turbulent heat transfer in a closed duct subjected to spanwise system rotation has many important industrial applications, such as turbine blade cooling and rotating heat exchangers. Pallares & Davidson (2002) conducted large-eddy simulation (LES) to study turbulent heat transfer in stationary and spanwise rotating square ducts under different wall heating conditions and Grashof numbers. They observed significant variation of Nusselt numbers with different boundary conditions (constant wall temperature or heat flux) and showed that turbulent intensity was strongly influenced by centrifugal buoyancy effect. Later, Pallares et al. (2005) further extended their LES studies to higher rotation number cases. They investigated the influences of secondary flow on the overall heat transfer performance and observed increased local Nusselt numbers on the unstable wall compared to the non-rotation case. Thus far, transient simulation of turbulent heat transfer in a spanwise rotating square duct has been limited to only several LES studies as reviewed above. Application of DNS to this type of research has so far focused only on the turbulent flow field in the literature. Dai et al. (2015) conducted a DNS study of a spanwise rotating square duct flow. Through a careful examination of the budget balance of the mean momentum equations, they demonstrated the influence of secondary flows on turbulence statistics obtained at different rotation numbers. Very recently, Fang et al. (2017) performed DNS of spanwise rotating square duct flows at very high rotation numbers. They observed that the Coriolis force dominates the transport of Reynolds stresses and turbulent kinetic energy (TKE), and forces the spectra of streamwise and vertical velocities to synchronize within a wide range of turbulence scales. For turbulent heat transfer within a spanwise rotating channel, Ω y T w x Flow 20πδ -2Ωu dt w /dx = constant 2Ωv x. Heat flux q w Figure 1. Side and front views of the computational domain and coordinate system. The two Coriolis force components in the rotating reference frame are 2Ωv and 2Ωu. The flow in the duct is cooled along the streamwise direction by imposing a constant wall heat flux on all four walls. the transport processes of momentum and thermal energy are coupled, and the physical mechanism is complicated not only by the two Coriolis force components induced by the system rotation but also by the presence of four solid side walls of the square duct. Currently, detailed DNS study of this special subject is still lacking in literature. In view of this, we aim at conducting a comprehensive comparative DNS study of the turbulent heat and fluid flow confined within a square duct subjected to spanwise system rotations, with a wide range of rotation numbers varying from = 2Ωδ/u τ = 0 (corresponding to the non-rotating state) to = 18 (which is the highest rotation number of the current literature). Here, u τ is the wall friction velocity, δ represents half the channel height, and Ω denotes the angular speed of the spanwise system rotation. Also as our research objectives, the first- and second-order statistical moments of the velocity and temperature fields will be examined, the spectral characteristics of the TKE and temperature fluctuations will be studied, and the physical mechanisms underlying the interactions of the large secondary flow structures with the four boundary layers developed along the four walls of the square duct will be investigated. TEST CASES AND NUMERICAL ALGORITHM Figure 1 shows the computational domain and coordinate system of the test case. The flow is confined within a closed square duct with the cross-sectional area L y L z = 2δ 2δ. Following Fang et al. (2017), in order to ensure the TKE associated with the most energetic eddies to be well captured in the streamwise direction, a very long square duct is used, with its streamwise computational domain size set to L x = 20πδ. Both the flow and temperature fields are assumed to be fully developed and the periodic boundary conditions are applied to the 2δ y Ω 2δ z
2 inlet and outlet of the duct. No-slip boundary condition is imposed on all four walls for the velocity field. The temperature field is assumed to be a passive scalar. A constant wall heat flux is imposed on all four walls of the square duct. The Prandtl number of the fluid is assumed to be 1.0 and the Reynolds number is fixed to Re τ = u τ δ/ν = 150. In order to examine the rotation effect, in total, eight rotation numbers are tested for = 0,,,, 4.5, 9.0, 12.0 and In the spanwise rotating reference frame, the continuity, momentum and thermal energy equations for an incompressible flow can be expressed as u i t u i x i = 0, (1) + (u i u j )= δ 1i Π 1 p + ν 2 u i + 2ε i j3 Ωu j, (2) x j ρ x i x j x j T t + dt w (u j T)= u 1 x j dx + α 2 T, (3) x j x j respectively. Here, coordinates x 1, x 2 and x 3 correspond to x, y and z, respectively, as in Fig. 1, and velocity components u 1, u 2 and u 3 are also denoted as u, v and w, respectively. In Eq. (3), T w = T w (x) represents the local mean peripheral temperature over the wall at a given streamwise location. In the above equation, ρ, p, ν and α represent the density, pressure, kinematic viscosity and thermal diffusivity of the fluid, respectively, and ε i jk and δ i j are the Levi- Civita and Kronecker symbols, respectively. In Eqs. (2) and (3), Π and dt w /dx are constant, representing the streamwise pressure and wall-temperature gradients, respectively. Given the fact that the duct flow is both hydraulically and thermally fully developed and is cooled down steadily under the constant surface heat flux condition, the discharge of thermal energy from the flow to the ambient environment occurs steadily along the streamwise direction. In consequence, both the bulk mean temperature and the local peripheral temperature of the wall decrease linearly in the streamwise direction, i.e. dt b /dx = dt w /dx = 2 q w /(ρc p U b δ). This further allows us to apply a periodical boundary condition to the temperature field in the streamwise direction, if the linear streamwise distribution of the bulk mean temperature is explicitly subtracted from the temperature field. In effect, in Eq. (3), T denotes the temperature deviation from the local bulk mean temperature, and T+(dT w /dx)x is the actual local temperature. This treatment method for the temperature field is similar to that used in the LES study of Pallares et al. (2005). DNS is conducted using an open-source spectral-element method code (so-called Semtex ) originally developed by Blackburn & Sherwin (2004) for calculation of turbulent velocity fields. In order to perform heat transfer calculation, an additional thermal energy equation is implemented in Semtex by the authors with the same spectral accuracy. The computer code is developed using C++ and FORTRAN programming languages and parallelized using message passing interface libraries. Fourier series are used in the homogeneous streamwise direction, and the quadrilateral spectral-element method based on high-order Gauss-Lobatto- Legendre Lagrange interpolants is employed in the cross-stream directions. The time integration is achieved by utilizing the high-order splitting method proposed by Karniadakis et al. (1991). In total, the current DNS uses grid points in the streamwise, vertical and spanwise directions, respectively. The cross-section of the duct is divided into rectangular elements. The values of y + and z + of the first grid off the walls are below 0.4 in all tested cases. All simulations are conducted using a 240-core local computer cluster. In this paper, the temperature field is non-dimensionalized using the wall friction temperature defined as T τ = q w /ρc p u τ, k x φ θθ = =4.5 =12.0 = λ x /δ Figure 2. Pre-multiplied temperature spectra k x φ θθ of different rotation numbers at position (y/δ, z/δ)=(0, -0.9). where C p represents the specific heat of the fluid. The nondimensionalized temperature is expressed as θ = (T T w )/T τ. In presenting the results, a pair of angular brackets are used to denote averaging both in time and over the homogeneous streamwise direction. As such, the mean temperature is denoted as θ, and the instantaneous temperature can be decomposed as θ = θ +θ. RESULTS AND DISCUSSIONS In order to validate that the streamwise computational domain size (L x = 20πδ) is long enough to resolve turbulent coherent structures of the temperature field, Fig. 2 compares the pre-multiplied spectra k x φ θθ at five rotation numbers. Here, k x and λ x denote the wavenumber and wavelength in the streamwise direction, respectively. As seen in the figure, in all the plotted cases, the magnitude of k x φ θθ at the largest resolved streamwise length scale is much smaller than its peak value. This indicates that the streamwise computational domain in current DNS is sufficiently long such that the turbulent scalar energy associated with different scales of flow motions are very well captured. Figure 3 compares the streamlines of the mean velocity field and contours of θ in the y-z plane at four different rotation numbers. From the figure, secondary flow appears as four pairs of counter-rotating vortices in each corner of the duct in the nonrotating case. In contrast, the pattern of secondary flow appears as a pair of counter-rotating vortices occupying the entire cross-section at = 18. The peak value of θ in the rotating case occurs at the vertical location below the duct center around y/δ =. This is due to the downwash of the mean flow caused by the secondary flow associated with the system rotation. For detailed discussion of the velocity field, the reader is referred to Fang et al. (2017). In order to compare the heat transfer performance, the Nusselt number is calculated as Nu = 1 θ θ b Γ n dγ, where n and Γ denote the wall-normal direction and four peripheries of the crosssection of the square duct, respectively, and θ b represents the nondimensional bulk mean temperature. Figure 4 shows the variations of the friction velocity and Nusselt number at different walls with the rotation number. In the figure, superscripts ( ) t, ( ) s and ( ) b denote the value of a quantity on the top, side and bottom walls, respectively. It is well-known that for a zero-pressure-gradient heated turbulent boundary layer developed over a solid smooth still flat plate, the following analogy between velocity and temperature wall laws hold in the viscous sublayer, i.e., u + = y + and θ = Pr y + (cor-
3 1.6 u τ t / u τ 1.4 u s τ / u τ u b τ / u τ 1.2 u τ (a) The = 0 and cases (a) Friction velocity 1 Nu t Nu s Nu b 1 Nu 1 1 (b) The = 4.5 and 18.0 cases Figure 3. Time averaged cross-stream streamlines superimposed with contours of θ for different rotation numbers. (b) Nusselt number Figure 4. The variation of the friction velocity and Nusselt number at different duct walls with the rotation number. respondingly, θ / n=pr u / n and Nu Pr u 2 τ ). This indicates that a large (small) Nusselt number relates to a large (small) wall friction velocity. The appearance the Coriolis forces as a result of system rotation has a significant impact on the transport processes of the momentum and thermal energy. As a result, such analogy between molecular diffusion processes of the momentum and thermal energy within the viscous sublayer holds only when the Coriolis force effects are negligible, and breaks as the rotation number increases. Indeed, based on their LES study of spanwise rotating square duct flow, Pallares & Davidson (2002) showed that this analogy holds approximately at low rotation numbers for 0< < The rotation number range tested here and in our previous study (Fang et al., 2017) is the widest in literature. Our DNS results show that as increases larger than 0.75, the trend of the Nusselt number becomes inconsistent with that of the wall friction velocity. For instance, as seen in Fig. 4, the friction velocity at the bottom wall (u b τ) peaks between = and 1.5, however, the value of Nu b increases monotonically within the rotation number range of 0 9. It is interesting to see from Fig. 4(b) that the value of Nu t shows a sharply decreasing trend for 0, whereas Nu t remains almost unchanged for This can be reasoned from the fact that the flow near the top wall becomes laminarized for (Fang et al., 2017) and, consequently, turbulent heat transfer is diminished. Figure 5 compares the vertical profiles of u + and θ at dif- ferent rotation numbers in the central plane (z/δ = 0), respectively. As seen in Fig. 5(a), the value of u + decreases monotonically as the rotation number increases. In contrast, from Fig. 5(b), it is interesting to observe that the variation of θ with the rotation number exhibits a non-symmetrical and non-monotonic pattern. More specifically, for 0, the magnitude of θ within the upper half duct (for 0 < y/δ < 1) decreases as the rotation number increases. In contrast, the level of θ increases as further increases beyond 4.5. It is known that a Taylor-Proudman (TP) region (in which the mean streamwise velocity does not vary along the axis of rotation, i.e., u / z 0) exists in both turbulent and laminar duct flows at a sufficiently high spanwise rotation number (Speziale, 1982; Fang et al., 2017). The presence of the TP region is evident at higher rotation numbers. For instance, as shown vividly in Fig. 3(b), at = 4.5, the TP region (as indicated by the temperature contours) dominates the center of the duct for z/δ < 0.7. It should be indicated here that although a passive scalar is useful for indicating the fluid flow structures, its transport process is governed separately by the convection-diffusion equation (i.e., Eq. (3)) rather than the momentum equation itself. Therefore, the dispersion pattern of the passive thermal energy may differ from the dominant flow pattern resulted purely from the momentum system. As shown in Fig. 3(b), as continues to increase to 18, the TP region can no loner be effectively demonstrated by using the contours of θ. In order to fur-
4 u + 1 u + 1 θ (a) Mean streamwise velocity u + = (b) Mean temperature θ Figure 5. Vertical profiles of the non-dimensional mean streamwise velocity u and the mean temperature θ of different rotation numbers in the central plane (z/δ = 0). Arrow shows the monotonic variation of = 0,,,, 4.5, 9.0, 12.0 and ther investigate the effect of system rotation on the transport of momentum and temperature, Fig. 6 shows the spanwise distributions of the mean velocity and temperature at different rotation numbers. As seen in Fig. 6(a), a TP region occurs at = and its spanwise extension increases monotonically as the rotation number further increases. In contrast, it is interesting to see in Fig. 6(b) that as the rotation number increases, the variation of θ is non-monotonic. Specifically, the TP effect on the distribution of θ is the strongest at = 4.5 (which extends from z/δ = 0.8 to 0.8), and weakens as further increases from 4.5 to 18. Figures 7 (a) and (b) show the vertical profiles of TKE (k = u i u i /2) and temperature variance θ θ for different rotation numbers in the central plane (z/δ = 0), respectively. It is evident from Fig. 7(a) that the turbulent intensity generally reduces as increases. Particularly, for cases of, the value of TKE in the upper half duct (0 y/δ 1) is almost zero, which is an indication of laminarization over the top wall. Consequently, as seen in Fig. 7(b), the value of θ θ approaches zero in the laminarized region. It is also interesting to observe that the peak magnitude of θ θ increases as the rotation number increases beyond. Figure 8 compares the vertical profiles of turbulent heat fluxes u θ and v θ at different rotation numbers in the central plane (z/δ = 0). From Fig. 8(a), the value of u θ monotonically decreases as the rotation number increases in general. In contrast, it θ z / δ (a) Mean streamwise velocity u z / δ (b) Mean temperature θ Figure 6. Spanwise profiles of the non-dimensional mean streamwise velocity u + and the mean temperature θ of different rotation numbers at half duct height (y/δ = 0). is interesting to see in Fig. 8(b) that the peak magnitude of v θ increases as the rotation number increases with the peak location fixed around y/δ = 0.9. Although both velocity and temperature fluctuations are suppressed individually by system rotation close to the bottom wall (see Fig. 7), Fig. 8(b) suggests that the correlation between v and θ near the bottom wall is actually enhanced as the rotation number increases. To further understand this interesting observation, the transport equation of turbulent heat flux u k θ can be studied, which reads u j u k θ x j } {{ } H k + u 1 u k dt w T τ dx 1 }{{} Pk s + u ju k θ + u x jθ u k j x j }{{} Pk c + u j u k θ x j } {{ } D t k 2ε k j3 Ω u jθ } {{ } C k α 2 θ u k x j x ν 2 u k θ + p θ = 0. j x j x j x k }{{}}{{} V k P k Here, H k, D t k and V k are the convection, turbulent diffusion and molecular diffusion-dissipation terms, respectively. P c k and Ps k represent the turbulent production and turbulent production term due to the streamwise mean temperature gradient, respectively. P k de- (4)
5 = k + θ θ (a) TKE k (b) Temperature variance θ θ Figure 7. Vertical profiles of TKE (k) and temperature variance θ θ of different rotation numbers in the central plane (z/δ = 0). u θ + v θ + 1 notes the pressure diffusion-redistribution term. The Coriolis production term C k is introduced by the system rotation. It is worthy to note that it relates u θ to v θ because C 1 = 2Ω v θ and C 2 = 2Ω u θ. From Fig. 8(b), for cases of, the value of v θ is negative in most of the turbulent regions near the bottom wall. Therefore, C 1 is positive and, consequently, destructs u θ. In other words, system rotation tends to suppress the level of u θ. In contrast, positively-valued u θ near the bottom wall leads to positively-valued C 2 and, consequently, enhances v θ. Figure 9 compares the vertical profiles of the budget balance of v θ for the non-rotating ( = 0) and rotating ( = ) cases at z/δ = 0. It is interesting to see that the Coriolis production term C 2 becomes dominant in the budget balance of v θ at = and its peak location is consistent with that of v θ shown in Fig. 8(b). In other words, system rotation tends to increase the magnitude of v θ. This explains the observation in Fig. 8(b) that the magnitude of v θ increases near the bottom wall as the rotation number increases. It is well known that the Taylor-Görtler (TG) vortices appear as streamwise-elongated counter-rotating roll cells in spanwiserotating turbulent plane-channel flows (Johnston et al., 1972; Kristoffersen & Andersson, 1993). As demonstrated previously in Fang et al. (2017) by using streamwise velocity spectra and here in Fig. 2 by using the streamwise temperature spectra, we confirm that such streamwise elongated structures also exist in a spanwise (a) Streamwise turbulent heat flux u θ (b) Vertical turbulent heat flux v θ + Figure 8. Vertical profiles of turbulent heat fluxes u θ + and v θ + of different rotation numbers in the central plane (z/δ = 0). rotating square duct flow. To further investigate, Fig. 10 displays a typical instantaneous fluctuating temperature fields around the roll cells in a cross-section plane at = To clearly demonstrate the flow structures, only a portion of the domain ( 1 z/δ 0 and 1 z/δ 0.4) is plotted. From Fig. (10), two pairs of counter-rotating roll cells can be observed and, ejection events (featuring positively-valued v ) occur between vortices. As a consequence, near-wall fluid with low temperature is convected upwards and negatively-valued θ is generated. As such, a large magnitude of v θ is maintained at high rotation numbers, a pattern that can be clearly observed in Fig. 8(b). CONCLUSIONS Direct numeric simulation (DNS) is conducted to investigate turbulent heat transfer in a closed square duct subjected to spanwise rotation. The rotation numbers tested here are among the highest in literature. The effects of system rotation on the turbulent flow and temperature fields are examined in terms of the mean temperature, Nusselt number, turbulent heat fluxes, streamwise temperature spectra, and the budget balance of the transport equation of turbulent heat fluxes. Induced by the Coriolis forces, secondary flows appear as streamwise-elongated counter-rotating vortex pairs in the crosssection direction. A Taylor-Proudman region appears in the center of the duct, which has a profound influence on the dispersion of
6 Budget balance of v θ 0.6 H 2 c s D t 2 C 2 V 2 effects are negligible, and breaks as the rotation number increases. Both the velocity and temperature fluctuations (as indicated by TKE and temperature variance θ θ ) become suppressed as the rotation number increases in general. However, it is interesting to observe that the magnitude of the vertical turbulent heat flux ( v θ ) increases as the rotation number increases. By investigating the budget balance of the transport equation of v θ, it is concluded that system rotation tends to enhance v θ through the Coriolis production term which relates to the streamwise turbulent heat flux u θ. Furthermore, by analyzing the turbulent structures in the cross-stream plane, it is observed that the ejection events pump near-wall cold fluid upwards to sustain the level of negatively-valued vertical turbulent heat flux v θ. Budget balance of v θ x / δ (c) Budget balance of v θ for = H 2 c s D t 2 C 2 V x / δ (d) Budget balance of v θ for = Figure 9. Budget balance of v θ in the central plane (z/δ = 0) at = 0 and. All the budget terms are non-dimensionalized using u 2 τ /δ. REFERENCES Blackburn, H. M. & Sherwin, S. J Formulation of a Galerkin spectral element-fourier method for three-dimensional incompressible flows in cylindrical geometries. J. Comp. Phys. 197(2), Dai, Y. J., Huang, W. X., Xu, C. X. & Cui, G. X Direct numerical simulation of turbulent flow in a rotating square duct. Phys. Fluids 27, Fang, X., Yang, Z., Wang, B.-C. & Bergstrom, D. J Direct numerical simulation of turbulent flow in a spanwise rotating square duct at high rotation numbers. Int. J. Heat Fluid Flow, in press. Johnston, J. P., Halleen, R. M. & Lezius, D. K Effects of spanwise rotation on the structures of two-dimensional fully developed turbulent channel flow. J. Fluid Mech. 56(3), Karniadakis, G. E., Israeli, M. & Orszag, S. A High-order splitting methods for the incompressible Navier-Stokes equations. J. Comp. Phys. 97(2), Kristoffersen, R. & Andersson, H. I Direct simulations of low-reynolds-number turbulent flow in a rotating channel. J. Fluid Mech. 256, Pallares, J. & Davidson, L Large-eddy simulations of turbulent heat transfer in stationary and rotating square ducts. Phys. Fluids 14(8), Pallares, J., Grau, F. X. & Davidson, L Pressure drop and heat transfer rates in forced convection rotating square duct flows at high rotation rates. Phys. Fluids 17, Speziale, C. G Numerical study of viscous flow in rotating rectangular ducts. J. Fluid Mech. 122, Figure 10. Contours of instantaneous temperature fluctuation θ superimposed with vectors of fluctuating velocities (w, v ) in a portion (for 1 z/δ 0 and 1 y/δ 0.4) of cross-section (z-y) plane at = the thermal energy. The investigation of Nusselt number at different walls shows that the analogy between the molecular diffusion processes of the momentum and thermal energy within the viscous sublayer (typical of a zero-pressure-gradient boundary layer developed over a solid smooth wall) holds only when the Coriolis force
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 information6.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 informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More 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 informationCHALMERS UNIVERSITY OF TECHNOLOGY
Report No. /3 CHALMERS UNIVERSITY OF TECHNOLOGY Department of Thermo and Fluid Dynamics Large-eddy simulations of turbulent flows in stationary and rotating channels and in a stationary square duct by
More informationTurbulent Heat Transfer in A Channel Flow with Arbitrary. Directional System Rotation
Turbulent Heat Transfer in A Channel Flow with Arbitrary Directional System tation Haibin Wu * and Nohuhide Kasagi Department of mechanical engineering, The University of Tokyo, Hongo 7-3-1, Bunkyo ku,
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 informationSECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES
SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES Yosuke Hasegawa Institute of Industrial Science The University of Tokyo Komaba 4-6-1, Meguro-ku, Tokyo 153-8505, Japan ysk@iis.u-tokyo.ac.jp
More informationWALL 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 informationDIRECT 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 informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
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 Rotating Turbulent Channel Flow by the Vortex in Cell Method
Send Orders for Reprints to reprints@benthamscience.net 30 The Open Transport Phenomena Journal, 2013, 5, 30-41 Open Access Numerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 ROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS Myoungkyu Lee Robert D. Moser
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationThe 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 informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationTurbulence Laboratory
Objective: CE 319F Elementary Mechanics of Fluids Department of Civil, Architectural and Environmental Engineering The University of Texas at Austin Turbulence Laboratory The objective of this laboratory
More informationProblem 4.3. Problem 4.4
Problem 4.3 Problem 4.4 Problem 4.5 Problem 4.6 Problem 4.7 This is forced convection flow over a streamlined body. Viscous (velocity) boundary layer approximations can be made if the Reynolds number Re
More informationINFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL
International Conference Nuclear Energy for New Europe 22 Kranjska Gora, Slovenia, September 9-2, 22 www.drustvo-js.si/gora22 INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationPeriodic 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 informationLecture 30 Review of Fluid Flow and Heat Transfer
Objectives In this lecture you will learn the following We shall summarise the principles used in fluid mechanics and heat transfer. It is assumed that the student has already been exposed to courses in
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationMYcsvtu Notes HEAT TRANSFER BY CONVECTION
www.mycsvtunotes.in HEAT TRANSFER BY CONVECTION CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More 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 informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP013672 TITLE: Direct Numerical Simulation of Turbulent Flow and Heat Transfer in a Square Duct at Low Reynolds Number DISTRIBUTION:
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationHEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1
HEAT TRANSFER BY CONVECTION Dr. Şaziye Balku 1 CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in the
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
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 informationContribution of Reynolds stress distribution to the skin friction in wall-bounded flows
Published in Phys. Fluids 14, L73-L76 (22). Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Koji Fukagata, Kaoru Iwamoto, and Nobuhide Kasagi Department of Mechanical
More 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 informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationRANS simulations of rotating flows
Center for Turbulence Research Annual Research Briefs 1999 257 RANS simulations of rotating flows By G. Iaccarino, A. Ooi, B. A. Pettersson Reif AND P. Durbin 1. Motivation and objectives Numerous experimental
More informationSecond-moment turbulence closure in prediction of channel flow with system rotation
Proceedings of the International Gas Turbine Congress 003 Tokyo November -7, 003 IGTC003Tokyo TS-061 Second-moment turbulence closure in prediction of channel flow with system rotation Masaaki ICHIKAWA,
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 informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationRoughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman
Roughness Sub Layers John Finnigan, Roger Shaw, Ned Patton, Ian Harman 1. Characteristics of the Roughness Sub layer With well understood caveats, the time averaged statistics of flow in the atmospheric
More informationfluctuations based on the resolved mean flow
Temperature Fluctuation Scaling in Reacting Boundary Layers M. Pino Martín CTR/NASA Ames, Moffett Field, CA 94035 Graham V. Candler Aerospace Engineering and Mechanics University of Minnesota, Minneapolis,
More informationOPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW
OPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW Shingo Motoki, Genta Kawahara and Masaki Shimizu Graduate School of Engineering Science Osaka University 1-3 Machikaneyama, Toyonaka, Osaka
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
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 informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationCHARACTERIZATION 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 informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More informationChuichi 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 informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationDirect and Large Eddy Simulation of stably stratified turbulent Ekman layers
Direct and Large Eddy Simulation of stably stratified turbulent Ekman layers Stimit Shah, Elie Bou-Zeid Princeton University 64 th APS DFD Baltimore, Maryland Nov 21, 211 Effect of Stability on Atmospheric
More information3D Numerical Study on Laminar Forced Convection in V-Baffled Square Channel
American Journal of Applied Sciences 10 (10): 1287-1297, 2013 ISSN: 1546-9239 2013 Boonloi and Jedsadaratanachai, This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
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 informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More informationHeat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to A Heated Flat Plate
1 nd International Conference on Environment and Industrial Innovation IPCBEE vol.35 (1) (1) IACSIT Press, Singapore Heat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
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 informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationvector H. If O is the point about which moments are desired, the angular moment about O is given:
The angular momentum A control volume analysis can be applied to the angular momentum, by letting B equal to angularmomentum vector H. If O is the point about which moments are desired, the angular moment
More informationPerformance of Elliptical Pin Fin Heat Exchanger with Three Elliptical Perforations
www.cfdl.issres.net Vol. 3 (2) June 2011 Performance of Elliptical Pin Fin Heat Exchanger with Three Elliptical Perforations Monoj Baruah 1, Anupam Dewan 2c and P. Mahanta 1 1 Department of Mechanical
More informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
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 informationVery large-scale structures observed in DNS of turbulent channel flow with passive scalar transport
5 th Australasian Fluid Mechanics Conference 3-7 December 4 Very large-scale structures observed in DNS of turbulent channel flow with passive scalar transport H. Kawamura, H. Abe and Y. Matsuo Department
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More 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 informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More 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 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 informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationChapter 3 NATURAL CONVECTION
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 Chapter 3 NATURAL CONVECTION Mehmet Kanoglu Copyright The McGraw-Hill Companies,
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER
THERMAL SCIENCE: Year 2018, Vol. 22, No. 2, pp. 963-972 963 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER by Jitesh RANA, Anshuman SILORI, Rajesh MAITHANI *, and
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 informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
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 information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationTurbulent momentum and heat transfer in ducts of rhombic cross section
Turbulent momentum and heat transfer in ducts of rhombic cross section Naoya Fukushima and Nobuhide Kasagi Department of Mechanical Engineering, The University of Tokyo Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656,
More informationy * x * Shumpei HARA
th International Symposium on Turbulence and Shear Flow Phenomena (TSFP), Chicago, USA, July, 7 ANALYSIS OF THE TURBULENT KINETIC ENERGY BUDGET FOR MEANDERING MOTION APPEARING IN A BACKWARD-FACING STEP
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