INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL

Size: px
Start display at page:

Download "INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL"

Transcription

1 International Conference Nuclear Energy for New Europe 22 Kranjska Gora, Slovenia, September 9-2, 22 INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL ABSTRACT Robert Bergant, Iztok Tiselj Jožef Stefan Institute Reactor Engineering Division Jamova 39, SI- Ljubljana, Slovenia robert.bergant@ijs.si, iztok.tiselj@ijs.si Direct Numerical Simulation (DNS of the fully developed velocity and temperature fields in the two-dimensional turbulent channel flow was performed for friction Reynolds number Re τ = 5 and Prandtl number Pr =.7. Two thermal boundary conditions (BCs, isothermal and isoflux, were carried out. The main difference between two ideal types of boundary conditions is in temperature fluctuations, which retain a nonzero value on the wall for isoflux BC, and zero for isothermal BC. Very interesting effect is seen in streamwise temperature auto-correlation functions. While the auto-correlation function for isothermal BC decreases close to zero in the observed computational domain, the decrease of the auto-correlation function for the isoflux BC is slower and remains well above zero. Therefore, another DNS at two times longer computational domain was performed, but results did not show any differences larger than the statistical uncertainty. INTRODUCTION Fully developed channel or flume flow is of great importance from the scientific and engineering point of view. Several experimental and numerical studies have been carried out to increase the understanding of the mechanics of the near-wall turbulent flows. Among them, the Direct Numerical Simulation (DNS became an important research tool of the turbulent heat transfer in the last 5 years. DNS means precise solving of Navier-Stoke s equations without any extra turbulent models. It is capable to describe all relevant scales in the velocity and thermal fields (Kolmogorov micro scales []. Some reviews of DNS results are presented by Kim&Moin [2] and Kasagi et al. [3] at moderate Reynolds number and Prandtl numbers of about one or less. Later Kawamura, et al. [4], Na and Hanratty [5] performed the DNS of the turbulent channel at Prandtl numbers up to ten. All DNSs mentioned above were performed at isothermal wall boundary condition for dimensionless temperature. Fig. shows fluid section bounded by bottom and top walls. Walls are heated by a constant heat source and cooled by a flow, which is assumed to be incompressible. In this paper two different thermal boundary conditions on the walls, isothermal and isoflux BCs, are considered. If a flow of air (low density, low heat capacity and low thermal conductivity is heated by a thick metal wall, the influence of the turbulence on the wall is practically negligible due to the high density, high heat capacity and high conductivity of the wall [6]. In such case the first ideal boundary condition (isothermal BC is approached. If water flow is 26.

2 26.2 heated by a tiny metal foil, the influence of the turbulence on the wall is strong because of the small wall thickness and relatively high density, high heat capacity and high thermal conductivity of water in comparison with density, thermal conductivity and heat capacity of the foil. The second boundary condition (isoflux BC is almost reached in such case. 2 EQUATIONS AND NUMERICAL PROCEDURE The flow in the channel is assumed to be fully developed. Both walls of the channel are heated by constant heat source, while the fluid flows between them, as shown in Fig.. L 3 - HEATED WALL 2 FLOW L2=2h=2 L Z HEATED WALL X Y, Figure. Computational domain. The dimensionless equations normalized with flume width h, friction velocity u τ kinematic viscosity ν, and friction temperature T τ = qw ( u τ ρ f c pf can be found in the papers of Kasagi [3] or Kawamura [4]: u v = v u v = u t v v 2 v ( u u p lx θ t = Re τ v 2 ux ( u θ θ Re τ Pr u B ( (2 (3 Reτ = uτ h / ν is the friction Reynolds number and Pr is Prandtl number. Terms l v x (unit vector in streamwise direction and u x u B appear in the equations (4 and (5 due to the numerical scheme that requires boundary conditions in streamwise and spanwise directions. u B is bulk mean velocity averaged over the time and entire volume. Therefore, the Reynolds number Re can be calculated as Re = ub Reτ. Dimensionless wall temperature difference is defined as θ ( x y, z, t ( x, y, z t Tw T,, =. (4 Tτ As can be seen from Eqs. (-3 temperature is assumed to be a passive scalar. This assumption introduces two approximations: neglected buoyancy, 2 neglected temperature dependence of the material properties - especially viscosity and heat conductivity. Results of Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

3 26.3 the present study are thus very accurate only for the systems, where the temperature differences are not too large, while some caution is required for the systems, where the temperature differences are not negligible. The outer sides of the walls are adiabatic, which mean that heat released in the walls is entirely transferred to the fluid. As mentioned before, two ideal boundary conditions are considered. The isothermal boundary condition at the wall is θ =, (5 while the isoflux boundary condition is imposed with the mean dimensionless temperature at heated wall fixed to zero and the ( = ± = θ y (6 dθ dy ( y = ± =. (7 It is evident that velocity at the interface wall-fluid is zero. The equations are solved with pseudo-spectral scheme using Fourier series in x and z directions and Chebyshev polynomials in the wall-normal y direction. Numerical procedure and the code of Gavrilakis [7] modified by Lam and Banerjee [8], and Lam [9] is used to solve the continuity and momentum equations. Equations (-3 are periodic in streamwise (x and spanwise (z directions. All our DNS simulations were performed at Re = 458 (Re τ = 5 and Pr =.7 in two different size of computational domains, shorter and longer domains (see table. Fig. 6 shows that the length of the shorter domain might not be long enough. Therefore another DNS with two times longer domain but with the same grid spacing in all three directions was performed to check the differences. It should be stressed that resolution that was chosen, is capable to describe all velocity and thermal structures according to the Kolmogorov theory []. Distances in wall units instead of IS units, are useful because turbulent flows with different Reynolds numbers can be easily compared in the near wall region. In wall units the height of the channel is equal to the two times friction Reynolds number. In our case the height is 3 wall units, which is according to the Fig. equal to 2 units of the computational box (units that are used in the computer code. The averaging was performed over time interval t =36 respectively t =9 after the fully developed flow was achieved (6 respectively time steps. Table : Summary of DNS calculations. time simulation grid comp. size resolution averaging Pr=.7 Re τ =5 (x, z, y x y z step time shorter domain 28x97x64 5π x 2 x π longer domain 256x97x64 π x 2 x π Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

4 26.4 RESULTS AND DISCUSSION The mean temperature profiles are given in Fig. 2. As can be seen, the type of thermal boundary condition does not affect the mean temperature profiles. Comparison of shorter and longer domain shows that two times longer length of the computational size does not give us different results. θ short_q long_q short_t long_t 2, y Figure 2: Profiles of mean temperature (Q - isoflux BC, T - isothermal BC. The differences between isothermal and isoflux boundary conditions can be clearly seen in Fig. 3. The thermal boundary condition affects only the temperature fluctuations, so the velocity fluctuations are not shown in this paper. In the case of isoflux thermal boundary condition the θ RMS remains constant across the viscous sublayer. In both cases there are no differences due to the different computational domain. Differences, which are not larger than 2% can be treated as statistical uncertainties. 3 RMS_θ 2,5 2,5 short_q short_t long_q long_t,5, y Figure 3: Profiles of RMS temperature fluctuations. Fig. 4 shows the turbulent axial heat flux versus the dimensionless distance from the wall. The location of the peak turbulent axial heat flux values is approximately at 5 wall units. The minor differences between thermal BCs appear near the wall (see Tiselj et. al. [] for details. The location of a peak is at the same distance from wall, the only visible difference is the height of heat flux peak. The peak in the case of isoflux BC is about 5% Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

5 26.5 higher than in the case of isothermal BC. If we compare heat flux for shorter and longer computational domain, one can see that there are no significant differences between the DNSs of different streamwise lengths θu short_q short_t long_q long_t y Figure 4: Profiles of the turbulent axial heat flux. The profiles of the normal heat fluxes are plotted in Fig. 5. The distinct feature is that the thermal wall boundary condition does not affect turbulent normal heat flux. The differences between different computational lengths are also negligible. However, the longer computational domain does not give any improvements in the meaning of heat fluxes.,8,7,6 θw,5,4,3,2 short_q short_t long_q long_t, y Figure 5: Profiles of the turbulent wall normal heat flux. Fig. 6 shows two-point auto-correlation functions of temperature (for two boundary conditions and streamwise velocity in streamwise direction calculated at distance y =5. from the heated wall. Auto-correlation function in the numerical simulation is used to find out if the dimensions of the channel are suitable. Because of the periodical boundary conditions the fluid must be well mixed in the first half of the channel, which means that the situation in the middle of the channel has nothing to do with the situation at the entrance (or end of the channel. The points at the end of the channel are neighbors to the points at the entrance of the channel due to the periodic boundary conditions. This is the reason why the auto-correlation functions are calculated only in the first half of the channel. The influence of the periodic Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

6 26.6 boundary conditions on the DNS are avoided, when the periodicity length ensures that autocorrelation functions of all fields fall to zero in the directions in which the periodic boundary conditions are imposed. An important difference is seen between two different boundary conditions (Fig. 6. The auto-correlation function for isothermal BC decreases close to zero in the observed computational domain, while the auto-correlation function for the isoflux BC decreases slower and remains well above zero. Therefore, the goal of the present study is to verify if the length of the computational domain for the isoflux BC should be longer to satisfy the generally accepted criteria for DNS about the sufficient length of the computational domain. For this reason auto-correlation function of isoflux BC at two times longer computational domain is added. In this concrete example (longrtt_q in Fig. 6 function decreases to zero and satisfies the sufficient length of the computational domain. The physical background of the differences between the auto-correlation functions (Fig. 6 at different thermal boundary conditions is not completely clear yet. We attempted to get a clearer physical picture with a new DNS study, which has been performed at Pr=.7 and with two times longer computational domain (2359 vs. 472 wall units but with the same grid resolution. Besides thermal BC comparison, also DNSs at "standard" and two times longer domain was performed (see Figs. -5. These results lead to next conclusions: The differences between the DNS in the "standard" length computational domain and the DNS in the extended computational domain did not show any differences larger than the statistical uncertainty for first-order statistics and also for the power spectra. 2 The periodicity length, which is long enough for the velocity field as an origin of the turbulence (Fig. 6, is long enough also for the passive scalar fields, despite the behavior of the streamwise two-point correlation at isoflux BC shown in Fig. 6. RR,2,8,6,4,2 shortruu shortrtt_q shortrtt_t longrtt_q -, x Figure 6: Two-point auto-correlation in streamwise direction at y =5.: temperature (isothermal, isoflux and streamwise velocity. Fig. 7 shows the spanwise two-point auto-correlations for temperatures at two different BCs and for streamwise velocity. The position of the minimum in the streamwise velocity correlation at approximately 5 wall units shows the presence of the low speed streaks. This minimum shows that certain flow structures - coherent structures exist close to the wall. These structures are the streamwise vortices that can be measured with the position of the first minimum. The first minimum shows the average distance between the high speed and low speed regions. Temperature minimums in fig. 7 show that they are slightly shifted to the right comparing to the velocity minimum. The temperature, which is a passive scalar do not have influence on the turbulence and the coherent structures close to the walls. Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

7 26.7,2,8,6 Ruu Rtt_Q Rtt_T RR,4,2 -,2 -, z Figure 7: Two-point auto-correlation in spanwise direction at y =5.: temperature (isothermal, isoflux and streamwise velocity. 3 CONCLUSIONS A direct numerical simulation of the fully developed channel flow was performed at friction Reynolds number Re τ = 5, and Prandtl number Pr =.7 (air. The fluid temperature was a passive scalar and therefore did not have the influence on the turbulence. Two different thermal boundary conditions at the heated wall were carried out. The auto-correlation functions are used to verify if the dimensions of the computational domain are suitable. Because of the periodical boundary conditions the fluid must be well mixed in the first half of the channel, where the auto-correlation functions must fall to zero. The temperature streamwise auto-correlation function at isoflux BC remained well above zero, therefore another DNS at two times longer domain was performed. The results did not show any differences larger than statistical error. On the other side, spanwise auto-correlation functions are used also to show the presence of the coherent structures of the near wall flow. These structures are streamwise vortices with the characteristic diameter of approximately wall units. Thermal boundary conditions did not affect the mean temperature profiles. The main difference in the first order statistics was seen in temperature fluctuations near the wall. Temperature fluctuations retained a nonzero value on the wall for isoflux boundary condition, and zero for isothermal wall boundary condition. The minor differences were seen in the turbulent axial heat flux profiles, where the peak in the case of isoflux BC is about 5% higher than in the case of isothermal BC. NOMENCLATURE h k L,L 3 p Pr q w Re τ Rtt channel half height wave number streamwise and spanwise length of turbulent box pressure Prandtl number wall-to-fluid heat flux friction Reynolds number auto-correlation function for temperature Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

8 26.8 Ruu auto-correlation function for streamwise velocity t time T τ = qw /( u τ ρ f c pf friction temperature u,v,w velocity components in x, y and z directions u τ = τ w / ρ friction velocity u B bulk mean velocity x streamwise distance y distance from the wall z spanwise distance Greek α = λ / ρc p thermal diffusivity ϑ = ( T w T / Tτ dimensionless temperature difference λ thermal conductivity ν kinematic viscosity ρ density r x unit vector in x direction (,, Subscripts and superscripts ( f fluid ( normalized by u τ, T τ, ν REFERENCES [] Tennekes, H., Lumley, J.L., 972, "A First Course in Turbulence", MIT Press, Cambridge, MA. [2] Kim, J., Moin, P., 989, "Transport of Passive Scalars in a Turbulent Channel Flow", Turbulent Shear Flows VI, Springer-Verlag, Berlin, pp. 85. [3] Kasagi, N., Tomita, Y., Kuroda, A., 992, "Direct Numerical Simulation of Passive Scalar Field in a Turbulent Channel Flow", Journal of Heat Transfer -Transactions of ASME, Vol. 4, [4] Kawamura, H., Ohsaka, K, Abe, H., Yamamoto, K., 998, "DNS of Turbulent Heat Transfer in Channel Flow with low to medium-high Prandtl number fluid", International Journal of Heat and Fluid Flow, Vol. 9, [5] Na, Y., Hanratty, T.J., 2, "Limiting Behavior of Turbulent Scalar Transport Close to a Wall", International Journal of Heat and Mass Transfer, Vol. 43, pp [6] Tiselj, I., Bergant, R., Mavko, B., Bajsič, I., Hetsroni, G., October 2, ýdns of turbulent heat transfer in Channel Flow with Heat Conduction in the Solid Wallý, Journal of Heat Transfer, Vol. 23, pp.: [7] Gavrilakis S., Tsai, H.M., Voke, P.R., Leslie, D.C., 986, "Direct and Large Eddy Simulation of Turbulence", Notes on Numerical Fluid Mechanics Vol. 5, ed. U. Schumann, R. Friedrich, Vieweg, Braunschweig, D.B.R., pp. 5. [8] Lam, K.L., Banerjee, S., 988, "Investigation of Turbulent Flow Bounded by a Wall and a Free Surface", Fundamentals of Gas-Liquid Flows, Edited by Michaelides, Sharma, Vol. 72, 29-38, ASME, Washington DC. [9] Lam, K.L., 989, "Numerical Investigation of Turbulent Flow Bounded by a Wall and a Free-Slip Surface", Ph.D. Thesis, Univ. Calif. Santa Barbara. [] Tiselj, I., Pogrebnyak, E., C. Li, Mosyak, A., Hetsroni, G., 2, "Effect of wall boundary condition on scalar transfer in a fully developed turbulent flume", Physics of Fluids, 3 (4, pp Proceedings of the International Conference Nuclear Energy for New Europe, Kranjska Gora, Slovenia, Sept. 9-2, 22

Spectral Element Direct Numerical Simulation of Sodium Flow over a Backward Facing Step

Spectral Element Direct Numerical Simulation of Sodium Flow over a Backward Facing Step Spectral Element Direct Numerical Simulation of Sodium Flow over a Backward Facing Step ABSTRACT Jure Oder, Jernej Urankar, Iztok Tiselj Jožef Stefan Institute Jamova cesta 39 SI-1000, Ljubljana, Slovenia

More information

Numerical simulations of heat transfer in plane channel flow

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

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

LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,

More information

Very large-scale structures observed in DNS of turbulent channel flow with passive scalar transport

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

The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface

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

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

DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT 10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department

More information

Numerical analysis of a fully developed non-isothermal particle-laden turbulent channel flow

Numerical analysis of a fully developed non-isothermal particle-laden turbulent channel flow Arch. Mech., 3, 1, pp. 77 91, Warszawa 11 Numerical analysis of a fully developed non-isothermal particle-laden turbulent channel flow M. JASZCZUR Department of Fundamental Research in Energy Engineering

More information

The JHU Turbulence Databases (JHTDB)

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

SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES

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

2.3 The Turbulent Flat Plate Boundary Layer

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

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

The JHU Turbulence Databases (JHTDB)

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

More information

The Effect of Schmidt Number on Air-Water Interface Mass Transfer

The Effect of Schmidt Number on Air-Water Interface Mass Transfer Proceeding of The th International Conference on Multiphase Flow 1-New Orleans May 7 to June 1, 1 New Orleans, ouisiana, U. S. A. The Effect of Schmidt Number on Air-Water Interface Mass Transfer Y. Hasegawa

More information

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

Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Published in Phys. Fluids 14, L73-L76 (22). Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Koji Fukagata, Kaoru Iwamoto, and Nobuhide Kasagi Department of Mechanical

More information

Principles of Convection

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

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

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

More information

DNS of turbulent channel flow with conjugate heat transfer: effect of thermal boundary conditions on the second moments and budgets

DNS of turbulent channel flow with conjugate heat transfer: effect of thermal boundary conditions on the second moments and budgets DNS of turbulent channel flow with conjugate heat transfer: effect of thermal boundary conditions on the second moments and budgets Cédric Flageul, Sofiane Benhamadouche, Eric Lamballais, Dominique Laurence

More information

BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW

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

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and

More information

Generation of initial fields for channel flow investigation

Generation of initial fields for channel flow investigation Generation of initial fields for channel flow investigation Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (Mai 2) In the framework of the DFG-funded research

More information

COMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES

COMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES COMPRESSIBLE TURBULENT CHANNEL AND PIPE FLOW: SIMILARITIES AND DIFFERENCES Somnath Ghosh, Rainer Friedrich Lehrstuhl für Aerodynamik, Technische Universität München Boltzmannstrasse 15, D-85748 Garching,

More information

Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity

Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Boundary layer flows The logarithmic law of the wall Mixing length model for turbulent viscosity Tobias Knopp D 23. November 28 Reynolds averaged Navier-Stokes equations Consider the RANS equations with

More information

Available online at ScienceDirect. Procedia Engineering 90 (2014 )

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

6.2 Governing Equations for Natural Convection

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

More information

Computation of turbulent Prandtl number for mixed convection around a heated cylinder

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

Interpretation of the mechanism associated with turbulent drag reduction in terms of anisotropy invariants

Interpretation of the mechanism associated with turbulent drag reduction in terms of anisotropy invariants J. Fluid Mech. (27), vol. 577, pp. 457 466. c 27 Cambridge University Press doi:1.117/s221127583 Printed in the United Kingdom 457 Interpretation of the mechanism associated with turbulent drag reduction

More information

Turbulence - Theory and Modelling GROUP-STUDIES:

Turbulence - Theory and Modelling GROUP-STUDIES: Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence

More information

CHALMERS UNIVERSITY OF TECHNOLOGY

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

Applied Mathematics and Mechanics (English Edition)

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

Lecture 30 Review of Fluid Flow and Heat Transfer

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

ROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS

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

6. Laminar and turbulent boundary layers

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

The Effect Of MHD On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel

The Effect Of MHD On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel The Effect Of MH On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel Rasul alizadeh,alireza darvish behanbar epartment of Mechanic, Faculty of Engineering Science &

More information

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS

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

Turbulent Heat Transfer in A Channel Flow with Arbitrary. Directional System Rotation

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

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

An evaluation of a conservative fourth order DNS code in turbulent channel flow Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical

More information

Journal of Fluid Science and Technology

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

More information

Wall turbulence with arbitrary mean velocity profiles

Wall turbulence with arbitrary mean velocity profiles Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the

More information

Mass Transfer in Turbulent Flow

Mass Transfer in Turbulent Flow Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.

More information

NEW SCALING LAWS FOR TURBULENT POISEUILLE FLOW WITH WALL TRANSPIRATION

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

More information

An investigation of turbulent oscillatory heat transfer in channel flows by large eddy simulation

An investigation of turbulent oscillatory heat transfer in channel flows by large eddy simulation International Journal of Heat and Mass Transfer 47 (4) 6 7 www.elsevier.com/locate/ijhmt An investigation of turbulent oscillatory heat transfer in channel flows by large eddy simulation Lei Wang, Xi-Yun

More information

MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM

MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM Hajime Nakamura Department of Mechanical Engineering, National Defense Academy 1-10-0 Hashirimizu,

More information

The mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):

The mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient): 8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation

More information

Turbulent momentum and heat transfer in ducts of rhombic cross section

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

Turbulence Modulation by Micro-Particles in Boundary Layers

Turbulence Modulation by Micro-Particles in Boundary Layers 1 Turbulence Modulation by Micro-Particles in Boundary Layers Maurizio Picciotto, Andrea Giusti, Cristian Marchioli and Alfredo Soldati Centro Interdipartimentale di Fluidodinamica e Idraulica and Dipartimento

More information

7.6 Example von Kármán s Laminar Boundary Layer Problem

7.6 Example von Kármán s Laminar Boundary Layer Problem CEE 3310 External Flows (Boundary Layers & Drag, Nov. 11, 2016 157 7.5 Review Non-Circular Pipes Laminar: f = 64/Re DH ± 40% Turbulent: f(re DH, ɛ/d H ) Moody chart for f ± 15% Bernoulli-Based Flow Metering

More information

CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW

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

Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step

Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step Copyright c 2004 Tech Science Press CMC, vol.1, no.3, pp.275-288, 2004 Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step B. Wang 1, H.Q. Zhang 1, C.K. Chan 2 and X.L.

More information

Defense Technical Information Center Compilation Part Notice

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

Applied Computational Fluid Dynamics

Applied Computational Fluid Dynamics Lecture 9 - Kolmogorov s Theory Applied Computational Fluid Dynamics Instructor: André Bakker André Bakker (2002-2005) Fluent Inc. (2002) 1 Eddy size Kolmogorov s theory describes how energy is transferred

More information

Turbulent boundary layer

Turbulent boundary layer Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic

More information

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

Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data 35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander

More information

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value

More information

Turbulent drag reduction by streamwise traveling waves

Turbulent drag reduction by streamwise traveling waves 51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For

More information

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition

Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,

More information

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

Publication 97/2. An Introduction to Turbulence Models. Lars Davidson,   lada ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature

More information

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

Before we consider two canonical turbulent flows we need a general description of turbulence. Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent

More information

OPTIMIZATION OF HEAT TRANSFER ENHANCEMENT IN PLANE COUETTE FLOW

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

Exercise 5: Exact Solutions to the Navier-Stokes Equations I

Exercise 5: Exact Solutions to the Navier-Stokes Equations I Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel

More information

4.2 Concepts of the Boundary Layer Theory

4.2 Concepts of the Boundary Layer Theory Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very

More information

SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG

SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG L. Thais Université de Lille Nord de France, USTL F9 Lille, France Laboratoire de Mécanique de Lille CNRS, UMR 817

More information

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):

Eddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr): AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress

More information

MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow

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

Computational Fluid Dynamics 2

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

Chapter 5. The Differential Forms of the Fundamental Laws

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

More information

CONVECTIVE HEAT TRANSFER

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

Turbulence: Basic Physics and Engineering Modeling

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

FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE

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

Numerical Simulation of Rocket Engine Internal Flows

Numerical Simulation of Rocket Engine Internal Flows Numerical Simulation of Rocket Engine Internal Flows Project Representative Masao Furukawa Authors Taro Shimizu Nobuhiro Yamanishi Chisachi Kato Nobuhide Kasagi Institute of Space Technology and Aeronautics,

More information

DAY 19: Boundary Layer

DAY 19: Boundary Layer DAY 19: Boundary Layer flat plate : let us neglect the shape of the leading edge for now flat plate boundary layer: in blue we highlight the region of the flow where velocity is influenced by the presence

More information

CONVECTIVE HEAT TRANSFER

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

HEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1

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

Direct numerical simulation of self-similar turbulent boundary layers in adverse pressure gradients

Direct numerical simulation of self-similar turbulent boundary layers in adverse pressure gradients Direct numerical simulation of self-similar turbulent boundary layers in adverse pressure gradients By Martin Skote, Dan S. Henningson and Ruud A. W. M. Henkes Direct numerical simulations of the Navier-Stokes

More information

Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder

Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder 326 Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder Qiusheng LIU, Katsuya FUKUDA and Zheng ZHANG Forced convection transient

More information

Numerical Heat and Mass Transfer

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

Chapter 5 Principles of Convection heat transfer (Text: J. P. Holman, Heat Transfer, 8 th ed., McGraw Hill, NY)

Chapter 5 Principles of Convection heat transfer (Text: J. P. Holman, Heat Transfer, 8 th ed., McGraw Hill, NY) hapter 5 Principles of onvection heat transfer (Tet: J. P. Holman, Heat Transfer, 8 th ed., McGra Hill, NY) onsider a fluid flo over a flat plate ith different temperatures (Fig 5-1) q A ha( T T ) since

More information

Turbulent Flows. g u

Turbulent Flows. g u .4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at

More information

Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions

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

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

Transport of a passive scalar at a shear-free boundary in fully developed turbulent open channel flow

Transport of a passive scalar at a shear-free boundary in fully developed turbulent open channel flow PHYSICS OF FLUIDS VOLUME 11, NUMBER 9 SEPTEMBER 1999 Transport of a passive scalar at a shear-free boundary in fully developed turbulent open channel flow R. A. Handler, a) J. R. Saylor, R. I. Leighton,

More information

Outlines. simple relations of fluid dynamics Boundary layer analysis. Important for basic understanding of convection heat transfer

Outlines. simple relations of fluid dynamics Boundary layer analysis. Important for basic understanding of convection heat transfer Forced Convection Outlines To examine the methods of calculating convection heat transfer (particularly, the ways of predicting the value of convection heat transfer coefficient, h) Convection heat transfer

More information

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

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

More information

ROLE OF THE VERTICAL PRESSURE GRADIENT IN WAVE BOUNDARY LAYERS

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

Problem 4.3. Problem 4.4

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

A combined application of the integral wall model and the rough wall rescaling-recycling method

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

ENGR Heat Transfer II

ENGR Heat Transfer II ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks

More information

Compressible turbulent channel flow with impedance boundary conditions

Compressible turbulent channel flow with impedance boundary conditions Center for Turbulence Research Proceedings of the Summer Program 2014 305 Compressible turbulent channel flow with impedance boundary conditions By C. Scalo, J. Bodart, S. K. Lele AND L. Joly We have performed

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

Analysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell

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

Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes

Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Excerpt from the Proceedings of the COMSOL Conference 9 Boston Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Daoyun Song *1, Rakesh K. Gupta 1 and Rajendra P. Chhabra

More information

RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8

RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical

More information

Effect of near-wall treatments on airflow simulations

Effect of near-wall treatments on airflow simulations Proceedings of 29 International Conference on Computational Methods for Energy Engineering and Environment: ICCM3E. Sousse, Tunisia, 2-22 November, 29, pp. Effect of near-wall treatments on airflow simulations

More information

2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;

2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ   v  = 0; Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received

More information

MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION

MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION MODELLING TURBULENT HEAT FLUXES USING THE ELLIPTIC BLENDING APPROACH FOR NATURAL CONVECTION F. Dehoux Fluid Mechanics, Power generation and Environment Department MFEE Dept.) EDF R&D Chatou, France frederic.dehoux@edf.fr

More information

FLOW-NORDITA Spring School on Turbulent Boundary Layers1

FLOW-NORDITA Spring School on Turbulent Boundary Layers1 Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on

More information

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6 Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture

More information

A SEAMLESS HYBRID RANS/LES MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS

A SEAMLESS HYBRID RANS/LES MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS A SEAMS HYBRID RANS/ MODEL WITH DYNAMIC REYNOLDS-STRESS CORRECTION FOR HIGH REYNOLDS NUMBER FLOWS ON COARSE GRIDS P. Nguyen 1, J. Uribe 2, I. Afgan 1 and D. Laurence 1 1 School of Mechanical, Aerospace

More information

MIXED CONVECTION OF NEWTONIAN FLUID BETWEEN VERTICAL PARALLEL PLATES CHANNEL WITH MHD EFFECT AND VARIATION IN BRINKMAN NUMBER

MIXED CONVECTION OF NEWTONIAN FLUID BETWEEN VERTICAL PARALLEL PLATES CHANNEL WITH MHD EFFECT AND VARIATION IN BRINKMAN NUMBER Bulletin of Engineering Tome VII [14] ISSN: 67 389 1. Rasul ALIZADEH,. Alireza DARVISH BAHAMBARI, 3. Komeil RAHMDEL MIXED CONVECTION OF NEWTONIAN FLUID BETWEEN VERTICAL PARALLEL PLATES CHANNEL WITH MHD

More information

Convection Heat Transfer. Introduction

Convection Heat Transfer. Introduction Convection Heat Transfer Reading Problems 12-1 12-8 12-40, 12-49, 12-68, 12-70, 12-87, 12-98 13-1 13-6 13-39, 13-47, 13-59 14-1 14-4 14-18, 14-24, 14-45, 14-82 Introduction Newton s Law of Cooling Controlling

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

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