INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL
|
|
- Colleen Andrews
- 5 years ago
- Views:
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 ABSTRACT Jure Oder, Jernej Urankar, Iztok Tiselj Jožef Stefan Institute Jamova cesta 39 SI-1000, Ljubljana, Slovenia
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 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 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 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 informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationNumerical 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 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 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 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 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 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 informationThe 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 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 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 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 informationDNS 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 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 informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationGeneration 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 informationCOMPRESSIBLE 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 informationBoundary 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 informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More information6.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 informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationInterpretation 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 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 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 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 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 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 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 informationThe 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 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 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 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 informationJournal 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 informationWall 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 informationMass 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 informationNEW 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 informationAn 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 informationMEASUREMENTS 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 informationThe 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 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 informationTurbulence 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 information7.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 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 informationVelocity 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 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 informationApplied 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 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 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 informationENGR 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 informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationLaminar 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 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 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 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 informationExercise 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 information4.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 informationSOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG
SOME DYNAMICAL FEATURES OF THE TURBULENT FLOW OF A VISCOELASTIC FLUID FOR REDUCED DRAG L. Thais Université de Lille Nord de France, USTL F9 Lille, France Laboratoire de Mécanique de Lille CNRS, UMR 817
More informationEddy 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 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 informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationChapter 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 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: 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 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 informationNumerical 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 informationDAY 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 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 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 informationDirect 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 informationTheoretical 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 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 informationChapter 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 informationTurbulent 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 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 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 informationTransport 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 informationOutlines. 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 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 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 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 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 informationENGR 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 informationCompressible 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 informationDifferential 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 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 informationWall 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 informationRAYLEIGH-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 informationEffect 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 information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationMODELLING 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 informationFLOW-NORDITA Spring School on Turbulent Boundary Layers1
Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on
More informationLectures 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 informationA 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 informationMIXED 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 informationConvection 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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information