Calculations on a heated cylinder case

Size: px
Start display at page:

Download "Calculations on a heated cylinder case"

Transcription

1 Calculations on a heated cylinder case J. C. Uribe and D. Laurence 1 Introduction In order to evaluate the wall functions in version 1.3 of Code Saturne, a heated cylinder case has been chosen. The case has been treated experimentally by Scholten and Murray [6] and numerically by Szczepanik et al. [7]. This is essentially an unsteady flow. In the case of the heated cylinder, plots of the Nusselt number show an impingement zone at the front of the cylinder where the heat transfer is high and again at the rear of the cylinder where the two vortices that are shed from the cylinder meet. The first part of the cylinder has a laminar region which extends up to around α = 8 9 (α = corresponds to the impingement point). The effects of increasing the Reynolds number are to bring forward the separation point and increase the heat transfer [4]. In this laminar region, the Nusselt number depends on the square root of the Reynolds number. Sanitjai and Goldstein [5] showed that the heat transfer is strongly dependent on the Prandtl number, by examining different gases and liquids with.7 < P r < 176. Another important parameter that influences the flow is the turbulence intensity at the inlet. The Reynolds number chosen for the present study is Re = 36 based on the diameter of the cylinder and bulk velocity at the inlet. The turbulence intensity is set.34% to match the experiments in Scholten and Murray [6]. The fluid is air at a temperature of 293K with P r =.7. 2 Computational Parameters Three meshes have been generated for this study. They all have the same domain, the axes of which coincide with the centre of the cylinder. The inlet plane is located at x = 8D where D is the cylinder diameter. The outlet plane is located at x = 2D. The top and bottom boundaries are set as inlet planes at y = ±8D. The first mesh (Mesh A) has 268 cells in total, with 164 cells around the cylinder. The first cell was found to have a minimum y + 1, and therefore another mesh was created to ensure a minimum value of y Mesh B has cells, with 192 cells around the cylinder periphery and with a minimum y + = 17.7 (see figure 1). The cells are clustered in the rear part of the cylinder as can be seen in figure 2. The cells are equally spaced in the direction normal to the cylinder surface, and extend to a radius of 5D and then are stretched towards the outer 1

2 boundaries. Mesh C has been created to obtain a low Reynolds model solution having y Figure 1: Distribution of y + for Mesh B Figure 2: Close-up of Mesh B y + < 1 in the first cell along the cylinder wall. This has 2916 cells with 192 cells around the cylinder wall and clustered in the wall normal direction. All solutions using Code Saturne have been obtained using a centred scheme for velocity, turbulent variables and temperature. The standard k ε model is used for the high Reynolds meshes and the SST model is used for the low Reynolds mesh. The version of Code Saturne used was v1.3g. The effects of different numerical parameters such as convection scheme, inlet conditions and solution algorithm have been studied. Air was used in all calculations, with the following properties: ρ = 1.25 Kg/m 3, µ =.1511x1 4 Kg/ms, Cp = 15 J/KgK, P r =.7. A turbulent Prandltl number of.9 was used. There main calculations were performed on this case. The k ε model with standard wall function in both StarCd and Code Saturne, the modified turbulence wall function in Code Saturne (see section 3.2) and the SST model in both codes. Additionally some other studies were conducted in Code Saturne such as the effect of numerical parameters including inlet conditions, numerical schemes and the introduction of scalable wall functions and linear production modifications. 3 Boundary Conditions 3.1 Inlet The inlet velocity was calculated from the Reynolds number and the air viscosity yielding U ref = m/s. The freestream Temperature was set to T = K. The inlet turbulence intensity (T u) was set to.34% according to the experiment therefore obtaining a value for the kinetic energy from k = 1.5T u 2 Uref 2. In order to obtain a complete set of inlet conditions, a value for ε needs to be prescribed. This can be done by prescribing a 2

3 length scale l and computing the dissipation as: ε = C3/4 µ k 3/2 l Is then up to the user to impose the turbulent length scale l. Another way is to impose a turbulence viscosity ratio (ν t /ν) and obtain the dissipation as: (1) ε = ν ν t C µ k 2 (2) The two methods present a problem since there is no reference to either the length scale or turbulence viscosity ratio from the experiment. The only length scale available is the cylinder diameter, D. It is known that in most wind tunnels the turbulent viscosity ration can vary from 5 to 5 depending on the level of turbulence. As it will be shown, the effect of the choice of value of dissipation greatly influences the solution. Whereas in StarCd equation (1) is used with the value of l prescribed by the user, Code Saturne computes the dissipation as: ε = 1 C3/4 µ k 3/2 (3) κl Equating equation (1) and (3) after inserting lengths l 1 and l 3 respectively leads to: l 1 = κl 3 /1. We see that in equation (3) the length scale corresponds to the edge of the Log layer when l 3 corresponds to the diameter of a pipe (the mixing-length increases as κy until about 1% of the diameter after which it stays constant). Similarly for a bluff body, the integral scale in the wake is about 5% of the body diameter. For the present case, setting l 1 = D (i.e. using equation (1)) yields a value of ε about twenty five times larger than setting l 3 = D (using equation (3)) which would then say that eddies are very small compared to the cylinder. There is not enough information on the experimental set up to determine which choice is the best, but by using equation (3) the turbulent viscosity ratio is ν t /ν = 2.84 which is a more realistic value for the low turbulence case treated here. The effect of the choice of dissipation value at the inlet can be seen in section Wall Standard Wall Function. At the wall three different approaches were tested. First the standard wall function, which is implemented differently in Code Saturne and in StarCD. In Code Saturne the value at the face has to be imposed, for the kinetic energy a value of k = u 2 k /C1/2 µ is prescribed, whereas for ε a constant flux condition is imposed so that the gradient calculated by the code at the first cell is equal to the theoretical value resulting in: u 3 k ε F = ε I + d (4) κ(d/2) 2 3

4 where F is the face centre, I is projection of the cell centre on the normal vector of the wall face and d is the distance between these two points. In StarCD the treatment of the wall boundary condition is different. According to the users manual the velocity is calculated from: u + = u u F (τ w /ρ) 1/2 (5) and the turbulence variables are treated as follows: k is solved at near-wall cells, after appropriate modifications. The latter consist of: setting the diffusive flux of k at the wall to zero (i.e. assuming k n = (6) evaluating the turbulence generation on the assumptions that the dominant strain rate is u/ y, evaluated from (5) and the dominant shear stress is the tangential component, which is uniform in the cell and equal to τ w (these assumptions are, of course, part of the underlying set on which the wall functions are based); evaluating the total turbulence dissipation by volume integration of the wall function distribution of ε, given by ε + = C 3/4 µ /κ By contrast, the transport equation for ε or ω is not employed at boundary cells; rather, the nodal values of these quantities are obtained directly from the algebraic wall-function relation given by equations ε + = εy k 3/2 (7) ω + = ωy u τ (8) This procedure whereby the value of epsilon is modified at the cell centre instead of applying a boundary condition at the cell face is in fact very common Modified Wall Function in Code Saturne. In order to compare the wall function implementation, an attempt was made to implement the boundary conditions for k and ε in a similar way as in StarCD. Although the manual does not describe the exact way the kinetic energy equation terms are represented at the first cell, a simple approach was taken here (see [2] for several possible ways to do this). The production term is calculated as: P k = τ w C1/4 µ k 1/2 κy (9) 4

5 and the boundary conditions are set as: k n = (1) ε = C3/4 µ k 3/2 κy note that the ε equation is still solved inside the wall boundary cell but the boundary conditions on ε is of Dirichlet type with a face value of ε = u 3 k /κy.the value of y has been chosen as the distance from the cell centre to the boundary face. This is arbitrary and could also be chosen as 1/2 or 1/4 of that value as in some cell vertex discretisation codes (N3S, ULYSSE). However it must be noted that a real Neumann condition on k as in (1) can now be applied instead of k = u 2 k /C1/2 µ which is what is currently coded. Finally the scalable wall function approach was also tested. 4 Results 4.1 Velocity Field The pressure coefficient around the cylinder can be seen in figure 3 and the friction coefficient in figure 4. From the figures it can be seen how the standard wall function approach (11) 2 1 Cp Saturne Cp StarCD Cp SST Saturne Cp Saturne Modified WF C p Figure 3: Pressure Coefficient Figure 4: Friction coefficient produces late separation compared to the low Reynolds solution (SST). By imposing a modified boundary condition (see section 3.2, equations (9), (1) and (11)), the separation point is still the same as with the standard wall function, but the pressure coefficient is closer to the SST solution. It is important to note that due to the limitations on the StarCD user subroutines, the friction coefficient was calculated by approximating τ w ν(u I /d) and although when using wall functions is it not a reliable approximation, it gives an idea of where the separation point lays. The averaged streamlines are shown in figures 5, 6, 7 and 8. It can be seen how the vortices predicted by the SST model on the fine mesh are larger since the separation point 5

6 is near = 9. It is also interesting to note that although the difference on the separation point between Code Saturne and StarCD is not very large (about 1 ) the recirculation zone is different not only on length (Code Saturne: 1.1D, StarCD:.7D) but also in the angle of the separated stream lines. The effect of applying a modified boundary condition is not very noticeable in terms of recirculation length which is only.1d larger than with the standard boundary condition in Code Saturne. In contrast to the velocity field, by Figure 5: k ε, standard wall functions, Figure 6: Code Saturne. StarCD. k ε, standard wall functions, Figure 7: SST fine mesh, Code Saturne. Figure 8: k ε, Modified WF, Code Saturne. changing the turbulent boundary condition, the behaviour for k and ε at the first cell around the cylinder is similar to the StarCD solution, as it can be seen from figures 9 and 1. The k profile around the cylinder is higher and the ε lower (when compared to StarCD) which in turn produces a large value for the turbulent viscosity. This difference could be due to the different implementation of the ε equation as explain in sections 3.2. It is also important to note than in this study, only the wall boundary condition for the turbulent variables has been modified, no attempt has been made to try to change the velocity wall function. 6

7 k StarCD k Saturne k Saturne Modified WF.25 ε StarCD ε Saturne ε Saturne Modified WF Figure 9: k at the first cell around the cylinder. 4.2 Temperature Field Figure 1: ε at the first cell around the cylinder In figure 11 the Nusselt number around the cylinder is shown. Both codes produce very similar SST-low-Re-mesh results, close to the experimental value except after separation where the plateau is over predicted. The use of wall functions yields different predictions. StarCD predicts a maximum which is not located at the impingement points at the front or the back but close to the 9 angle. This is a feature difficult to explain from the physical point of view. On the other hand, Code Saturne predicts correct levels at the impingement but the minimum (separation point) is wrongly estimated. The use of the modified wall function reproduces the shape of predicted by StarCD but at higher level. This could be because of the implementation of the turbulent wall functions. In StarCD the ε equation is not solved at the first cell, whereas in Code Saturne the face value is imposed as u 3 k /κy. Another explanation could be the formulation of the production and dissipation terms in the k equation which is not know in StarCD. It is clear that the treatment of the turbulent variables at the first cell is very important for this case. 4.3 Effects of the convection scheme To test the sensitivity of the solution to the convetion scheme, another calculation was computed with an upwind scheme for the turbulent variables. Even though the velocities had a centered scheme, the resulting Nusselt number is greatly affected. This is not the case for the SST calculation on the fine mesh where the use of different schemes for the turbulent variables practically does not affect the results. 4.4 Effect of the inlet length scale choice This particular case is sensitive to the length scale used to compute the dissipation at the inlet as it can be seen from figure 13. Here the results of using l 1 = D, that is equation (1) or l 3 = D in (3) can be seen. By using equation (1) the value of ε at the inlet is smaller 7

8 4 3 Exp k-ε Saturne k-ε Star CD SST Saturne SST StarCD k-ε Saturned Modified WF Nu Figure 11: Nusselt number distribution 3 25 Exp k-ε Saturne k-ε Saturne upwind k-ε StarCD k-ε StarCD upwind 2 Nu Figure 12: Effect of the convection scheme on Nusselt number distribution. 8

9 therefore is producing a higher viscosity ratio. At the impingement point, the use of a higher length scale (equation (1)) produces a value for the Nusselt number 25% higher for StarCd and 62% for Code Saturne. This illustrates the high sensitivity of a key parameter to a basically unknown inlet condition, something users should be made aware of in e.g. tutorials Exp k-ε Saturne k-ε Star CD StarCD hi length scale Saturne hi length scale 2 Nu Figure 13: Effect of the inlet length scale on Nusselt number distribution. Using equation (3): StarCD, Code Saturne. Using equation (1): Star CD Code Saturne 4.5 Other computations Another two possibilities available in Code Saturne have been studied, namely the use of the scalable wall function in the fine mesh with the k ε and the linear production model on the coarse mesh. The scalable wall function was design to allow the high Reynolds models to use a fine with y + 11 The linear production model [3] reduces the production of k at stagnation points by making it proportional to the strain instead of the quadratic proportionality of the traditional k ε models. The effect of these two approaches can be seen in figure 14. The corresponding streamlines can be seen in figures 15 and 16. The scalable wall function produces a much larger over prediction of the Nusselt number, similar to what is expected from a low Reynolds k ε model without corrections. The recirculation length is smaller than the one produced by the standard wall function. The linear production model reduces the Nusselt number at the impingement point as it would be expected, therefore laminarising the flow around the front of the cylinder. This laminar region is present until the separation point which is earlier than the standard wall function approach. 9

10 8 6 Exp k-ε Saturne standard k-ε Saturne Scalable WF k-ε Saturne Linear production Nu Figure 14: Nusselt number with Scalable wall function and linear production. Figure 15: k ε, scalable wall functions. Figure 16: k ε, linear production. 1

11 4.6 Higher Reynolds number The same case was treated at a higher Reynolds number (Re = 89) to compare with the experiments of Sanitjai and Goldstein [5]. The same mesh was used and with the nondimensional distance varying as 48 < y + < 18. In the experiment the measure turbulence intensity was T u =.3% which is comparable wit the previous case at lower Reynolds number (T u =.34%). Equation (3) was used to compute the dissipation at the inlet with a length scale of 1D. The resulting Nusselt number predictions can be seen in figure 17. When using the presently modified or more universal boundary conditions for k and ε as explained in section 3.2 both Code Saturne and StarCD have similar patterns as in the previous case but this time there is a large overprediction of the Nusselt number. The original Code Saturne boundary condition now also gives a severe overestimation of the heat transfer, in particular at the stagnation point. This also shows that the occasional agreement at this stagnation point in previous graphs was coincidental and due to a fortunate chose of value of epsilon at the inlet. Here again the low Reynolds number model captures the impingement and separation quite well but overpredicts the heat transfer at the rear of the cylinder. 6 5 Exp k-ε StarCD SST Saturne k-ε Saturne k-ε Saturne modified WF Nu Figure 17: Nusselt number with Scalable wall function and linear production. 5 Conclusions The use of the wall function in this case yields unsatisfactory results. Mesh requirements have been followed to ensure that the value of y + of the first cell is inside the log law region, and the circumferential resolution has been demonstrated to be sufficient. By comparing the results with StarCD, the question of the behaviour of the wall function cannot be answered, since both codes yield different results. By resesting the bound- 11

12 ary conditions on k and ε in a similar way to what is the commonly agreed norm for wall functions, better agreement between the codes can be achieved. The results with the scalable wall function differ from the standard wall function, implying that it might be a low Reynolds number effect, therefore needing more points after the large first cell. This case has been proven to be very sensitive to different parameters. The choice of lengthscale at the inlet boundary is crucial for the prediction of flow and thermal fields but since there is no information about the experimental set up it remains an open question. The original boundary conditions for Code Saturne seem to be more sensitive to the inlet values or the Re number than show StarCD and Code Saturne with modified boundary conditions. The numerical parameters such as the convection scheme and resolution method also play an important role and need to be carefully taken into account. Although the results with standard wall function in Code Saturne differ from the StarCD, it can be concluded that the main difference arise from the implementation of the k and ε wall functions. Unfortunately this is not a case that can yield to a decision of which way of implementing the these boundary conditions is better since both results are not acceptable and only a low Reynolds model seems to be able to capture the main important features such as separation point and heat transfer distribution. Another inconvenience of this case is the lack of data, such as reliable inlet conditions and data about the flow field such as friction coefficient and nature of the flow in the wake of the cylinder. [1] presents an interesting set of heat transfer measurements for a small tube bundle in cross-flow (4 rows of tubes) which might be interesting to consider. This could also be accompanied by a refined LES with heat transfer to provide combined velocity, turbulence and temperature data, thus resolving a number of unknowns. References [1] E. Buyuk. Heat transfer and flow structures around circular cylinders in cross flow. Journal of Engineering and environmental science, 23: , [2] S. Gant. Development and Application of a New Wall Function for Complex Turbulent Flows. PhD thesis, UMIST, 22. [3] V. Guimet and D. Laurence. A linearised turbulent production in the k ε model for engineering applications. In Engineering Turbulence and Measurements 5, page 11, 22. [4] H. Nakamura and T. Igarashi. Unsteady heat transfer from a circular cylinder for Reynolds numbers from 3 to 15. International Journal of Heat and Fluid Flow, 25: , 24. [5] S. Sanitjai and R. J. Goldstein. Forced convection heat transfer from a circular cylinder in crossflow to air and liquids. International Journal of Heat and Mass Transfer, 47: ,

13 [6] J. W. Scholten and D. B. Murray. Unsteady heat transfer and velocity of a cylinder in cross flow I. low freestream turbulence. International Journal of Heat and Mass Transfer, 41(1): , [7] K. Szczepanik, A. Ooi, L. Aye, and G. Rosengarten. A numerical study of heat transfer from a cylinder in cross flow. In 15 th Australasian Fluids Mechanics Conference. The university of Sydney,

Tutorial for the supercritical pressure pipe with STAR-CCM+

Tutorial for the supercritical pressure pipe with STAR-CCM+ Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with

More information

Turbulent Boundary Layers & Turbulence Models. Lecture 09

Turbulent Boundary Layers & Turbulence Models. Lecture 09 Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects

More information

Explicit algebraic Reynolds stress models for internal flows

Explicit algebraic Reynolds stress models for internal flows 5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence

More information

Validation 3. Laminar Flow Around a Circular Cylinder

Validation 3. Laminar Flow Around a Circular Cylinder Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous

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

Tutorial for the heated pipe with constant fluid properties in STAR-CCM+

Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ Tutorial for the heated pipe with constant fluid properties in STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities

More information

PHYSICAL MECHANISM OF CONVECTION

PHYSICAL MECHANISM OF CONVECTION Tue 8:54:24 AM Slide Nr. 0 of 33 Slides PHYSICAL MECHANISM OF CONVECTION Heat transfer through a fluid is by convection in the presence of bulk fluid motion and by conduction in the absence of it. Chapter

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

On the transient modelling of impinging jets heat transfer. A practical approach

On the transient modelling of impinging jets heat transfer. A practical approach Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,

More information

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling

Numerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline

More information

Heat Transfer F12-ENG Lab #4 Forced convection School of Engineering, UC Merced.

Heat Transfer F12-ENG Lab #4 Forced convection School of Engineering, UC Merced. 1 Heat Transfer F12-ENG-135 - Lab #4 Forced convection School of Engineering, UC Merced. October 23, 2012 1 General purpose of the Laboratory To gain a physical understanding of the behavior of the average

More information

External Forced Convection :

External Forced Convection : External Forced Convection : Flow over Bluff Objects (Cylinders, Spheres, Packed Beds) and Impinging Jets Chapter 7 Sections 7.4 through 7.8 7.4 The Cylinder in Cross Flow Conditions depend on special

More information

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds. Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,

More information

Prediction of unsteady heat transfer from a cylinder in crossflow

Prediction of unsteady heat transfer from a cylinder in crossflow Center for Turbulence Research Proceedings of the Summer Program 202 07 Prediction of unsteady heat transfer from a cylinder in crossflow By S. T. Bose, B. C. Wang AND M. Saeedi The accuracy of a tensorial

More information

Introduction to Turbulence AEEM Why study turbulent flows?

Introduction to Turbulence AEEM Why study turbulent flows? Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and

More information

Numerical study of forced convection around heated horizontal triangular ducts

Numerical study of forced convection around heated horizontal triangular ducts Advanced Computational Methods and Experiments in Heat Transfer XI 0 Numerical study of forced convection around heated horizontal triangular ducts O. Zeitoun, M. E. Ali & A. Nuhait King Saud University,

More information

Forced Convection: Inside Pipe HANNA ILYANI ZULHAIMI

Forced Convection: Inside Pipe HANNA ILYANI ZULHAIMI + Forced Convection: Inside Pipe HANNA ILYANI ZULHAIMI + OUTLINE u Introduction and Dimensionless Numbers u Heat Transfer Coefficient for Laminar Flow inside a Pipe u Heat Transfer Coefficient for Turbulent

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

Numerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers

Numerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An

More 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

LARGE-EDDY SIMULATIONS OF FLUID FLOW AND HEAT TRANSFER AROUND A PARABOLIC TROUGH SOLAR COLLECTOR

LARGE-EDDY SIMULATIONS OF FLUID FLOW AND HEAT TRANSFER AROUND A PARABOLIC TROUGH SOLAR COLLECTOR LARGE-EDDY SIMULATIONS OF FLUID FLOW AND HEAT TRANSFER AROUND A PARABOLIC TROUGH SOLAR COLLECTOR Ahmed Amine Hachicha, Ivette Rodríguez and Assensi Oliva 1 1 Heat and Mass Transfer Technological Center

More information

Boundary-Layer Theory

Boundary-Layer Theory Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22

More information

Tutorial 1. Where Nu=(hl/k); Reynolds number Re=(Vlρ/µ) and Prandtl number Pr=(µCp/k)

Tutorial 1. Where Nu=(hl/k); Reynolds number Re=(Vlρ/µ) and Prandtl number Pr=(µCp/k) Tutorial 1 1. Explain in detail the mechanism of forced convection. Show by dimensional analysis (Rayleigh method) that data for forced convection may be correlated by an equation of the form Nu = φ (Re,

More information

Active Control of Separated Cascade Flow

Active Control of Separated Cascade Flow Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.

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

A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM

A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM 206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk

More information

Introduction to Heat and Mass Transfer. Week 9

Introduction to Heat and Mass Transfer. Week 9 Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional

More information

TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS

TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,

More information

Convective Mass Transfer

Convective Mass Transfer Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface

More 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

1. Introduction, tensors, kinematics

1. Introduction, tensors, kinematics 1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and

More information

Fundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.

Fundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D. Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer

More 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

ρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ

ρ t + (ρu j ) = 0 (2.1) x j +U j = 0 (2.3) ρ +ρ U j ρ Chapter 2 Mathematical Models The following sections present the equations which are used in the numerical simulations documented in this thesis. For clarity, equations have been presented in Cartesian

More information

CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2

CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2 CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration G.S.T.A. Bangga 1*, W.A. Widodo 2 1,2 Department of mechanical engineering Field of study energy conversion Institut

More information

Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder

Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:3 32 Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder Ressan Faris Al-Maliky Department of

More information

A COMPARISON OF HEAT TRANSFER AROUND A SINGLE SERRATED FINNED TUBE AND A PLAIN FINNED TUBE

A COMPARISON OF HEAT TRANSFER AROUND A SINGLE SERRATED FINNED TUBE AND A PLAIN FINNED TUBE IJRRAS 2 (2) February A COMPARISON OF HEAT TRANSFER AROUND A SINGLE SERRATED FINNED TUBE AND A PLAIN FINNED TUBE S. R. MCILWAIN School of Engineering, University of the West of Scotland, Hamilton, Scotland

More information

ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics

ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics For Friday, October 26 th Start reading the handout entitled Notes on finite-volume methods. Review Chapter 7 on Dimensional Analysis

More information

Table of Contents. Foreword... xiii. Preface... xv

Table of Contents. Foreword... xiii. Preface... xv Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...

More 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

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

Day 24: Flow around objects

Day 24: Flow around objects Day 24: Flow around objects case 1) fluid flowing around a fixed object (e.g. bridge pier) case 2) object travelling within a fluid (cars, ships planes) two forces are exerted between the fluid and the

More information

Conjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane

Conjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane Conjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane V. Esfahanian 1, A. Shahbazi 1 and G. Ahmadi 2 1 Department of Mechanical Engineering, University of Tehran, Tehran, Iran 2 Department

More information

OPTIMAL DESIGN OF CLUTCH PLATE BASED ON HEAT AND STRUCTURAL PARAMETERS USING CFD AND FEA

OPTIMAL DESIGN OF CLUTCH PLATE BASED ON HEAT AND STRUCTURAL PARAMETERS USING CFD AND FEA International Journal of Mechanical Engineering and Technology (IJMET) Volume 9, Issue 5, May 2018, pp. 717 724, Article ID: IJMET_09_05_079 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=9&itype=5

More information

Explicit algebraic Reynolds stress models for boundary layer flows

Explicit algebraic Reynolds stress models for boundary layer flows 1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both

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

A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case

A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case A Discussion of Low Reynolds Number Flow for the Two-Dimensional Benchmark Test Case M. Weng, P. V. Nielsen and L. Liu Aalborg University Introduction. The use of CFD in ventilation research has arrived

More information

UNIT II CONVECTION HEAT TRANSFER

UNIT II CONVECTION HEAT TRANSFER UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid

More information

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders

Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute

More 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

Numerical modeling of complex turbulent flows

Numerical modeling of complex turbulent flows ISTP-16, 5, PRAGUE 16 TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA Numerical modeling of complex turbulent flows Karel Kozel Petr Louda Jaromír Příhoda Dept. of Technical Mathematics CTU Prague, Karlovo

More information

Chapter 3 NATURAL CONVECTION

Chapter 3 NATURAL CONVECTION Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 Chapter 3 NATURAL CONVECTION Mehmet Kanoglu Copyright The McGraw-Hill Companies,

More information

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

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

More information

Forced Convection Around Obstacles

Forced Convection Around Obstacles Chapter 4 Forced Convection Around Obstacles 4.1. Description of the flow This chapter is devoted to heat transfer on bodies immersed in a stream. We consider a solid characterized by the length scale

More information

TABLE OF CONTENTS CHAPTER TITLE PAGE

TABLE OF CONTENTS CHAPTER TITLE PAGE v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review

More information

Wall treatments and wall functions

Wall treatments and wall functions Wall treatments and wall functions A wall treatment is the set of near-wall modelling assumptions for each turbulence model. Three types of wall treatment are provided in FLUENT, although all three might

More information

A concept for the integrated 3D flow, heat transfer and structural calculation of compact heat exchangers

A concept for the integrated 3D flow, heat transfer and structural calculation of compact heat exchangers Advanced Computational Methods and Experiments in Heat Transfer XIII 133 A concept for the integrated 3D flow, heat transfer and structural calculation of compact heat exchangers F. Yang, K. Mohrlok, U.

More information

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force

More information

CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer

CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer Accepted for publication in International Journal of Heat and Mass Transfer CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer Thijs Defraeye a, *,

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

Department of Mechanical Engineering

Department of Mechanical Engineering Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible

More information

Study of flow and heat transfer in a 3D pin-fin channel

Study of flow and heat transfer in a 3D pin-fin channel Study of flow and heat transfer in a 3D pin-fin channel Jean-Eloi Lombard Swiss Institute of Technology Lausanne Supervisor : Dr. Mark Sawley Assistants : Laura Augello, Thierry Cornu December 27, 2011

More information

Application of COMSOL Multiphysics in Transport Phenomena Educational Processes

Application of COMSOL Multiphysics in Transport Phenomena Educational Processes Application of COMSOL Multiphysics in Transport Phenomena Educational Processes M. Vasilev, P. Sharma and P. L. Mills * Department of Chemical and Natural Gas Engineering, Texas A&M University-Kingsville,

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

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

PART 1B EXPERIMENTAL ENGINEERING. SUBJECT: FLUID MECHANICS & HEAT TRANSFER LOCATION: HYDRAULICS LAB (Gnd Floor Inglis Bldg) BOUNDARY LAYERS AND DRAG

PART 1B EXPERIMENTAL ENGINEERING. SUBJECT: FLUID MECHANICS & HEAT TRANSFER LOCATION: HYDRAULICS LAB (Gnd Floor Inglis Bldg) BOUNDARY LAYERS AND DRAG 1 PART 1B EXPERIMENTAL ENGINEERING SUBJECT: FLUID MECHANICS & HEAT TRANSFER LOCATION: HYDRAULICS LAB (Gnd Floor Inglis Bldg) EXPERIMENT T3 (LONG) BOUNDARY LAYERS AND DRAG OBJECTIVES a) To measure the velocity

More information

Elliptic relaxation for near wall turbulence models

Elliptic relaxation for near wall turbulence models Elliptic relaxation for near wall turbulence models J.C. Uribe University of Manchester School of Mechanical, Aerospace & Civil Engineering Elliptic relaxation for near wall turbulence models p. 1/22 Outline

More information

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University

Chapter 7: External Forced Convection. Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Chapter 7: External Forced Convection Dr Ali Jawarneh Department of Mechanical Engineering Hashemite University Objectives When you finish studying this chapter, you should be able to: Distinguish between

More information

Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2

Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,

More information

ME 331 Homework Assignment #6

ME 331 Homework Assignment #6 ME 33 Homework Assignment #6 Problem Statement: ater at 30 o C flows through a long.85 cm diameter tube at a mass flow rate of 0.020 kg/s. Find: The mean velocity (u m ), maximum velocity (u MAX ), and

More information

The effect of geometric parameters on the head loss factor in headers

The effect of geometric parameters on the head loss factor in headers Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract

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

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Studies on flow through and around a porous permeable sphere: II. Heat Transfer Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India

More information

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13

Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13 Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :

More information

RANS simulations of rotating flows

RANS simulations of rotating flows Center for Turbulence Research Annual Research Briefs 1999 257 RANS simulations of rotating flows By G. Iaccarino, A. Ooi, B. A. Pettersson Reif AND P. Durbin 1. Motivation and objectives Numerous experimental

More information

Comparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building

Comparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building May 26, Volume 3, Issue 5 Comparison of Numerical Prediction of Pressure Coefficient on Rectangular Small Building Neel M. Patel, 2 Satyen D. Ramani PG Student, Department of Civil Engineering, 2 Asst.

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

Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.

Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative

More information

RSM AND v 2 -f PREDICTIONS OF AN IMPINGING JET IN A CROSS FLOW ON A HEATED SURFACE AND ON A PEDESTAL

RSM AND v 2 -f PREDICTIONS OF AN IMPINGING JET IN A CROSS FLOW ON A HEATED SURFACE AND ON A PEDESTAL 16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia -7 December 007 RSM AND v -f PREDICTIONS OF AN IMPINGING JET IN A CROSS FLOW ON A HEATED SURFACE AND ON A PEDESTAL D. Rundström

More information

Initial and Boundary Conditions

Initial and Boundary Conditions Initial and Boundary Conditions Initial- and boundary conditions are needed For a steady problems correct initial conditions is important to reduce computational time and reach convergence Boundary conditions

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

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

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable

More information

WALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS

WALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue

More information

Laminar Flow. Chapter ZERO PRESSURE GRADIENT

Laminar Flow. Chapter ZERO PRESSURE GRADIENT Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible

More information

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral

More information

Turbulence Solutions

Turbulence Solutions School of Mechanical, Aerospace & Civil Engineering 3rd Year/MSc Fluids Turbulence Solutions Question 1. Decomposing into mean and fluctuating parts, we write M = M + m and Ũ i = U i + u i a. The transport

More information

Forced Convection from Heated Rectangular Cylinder in Turbulent Regime using PANS Method

Forced Convection from Heated Rectangular Cylinder in Turbulent Regime using PANS Method Proceedings of the 6 th International and 43 rd National Conference on Fluid Mechanics and Fluid Power December 15-17, 2016, MNNITA, Allahabad, U.P., India Forced Convection from Heated Rectangular Cylinder

More information

7. Basics of Turbulent Flow Figure 1.

7. Basics of Turbulent Flow Figure 1. 1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds

More information

Turbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.

Turbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics. Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena

More information

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer 1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic

More information

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation

Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation , pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,

More information

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

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

More information

Experimental and Numerical Investigations of Flow Characteristics Over Half Sphere Positioned at Unfavorable Angle to The Flow

Experimental and Numerical Investigations of Flow Characteristics Over Half Sphere Positioned at Unfavorable Angle to The Flow International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:12 No:05 35 Experimental and Numerical Investigations of Flow Characteristics Over Half Sphere Positioned at Unfavorable Angle

More information

Convection in Three-Dimensional Separated and Attached Flow

Convection in Three-Dimensional Separated and Attached Flow Convection in Three-Dimensional Separated and Attached Flow B. F. Armaly Convection Heat Transfer Laboratory Department of Mechanical and Aerospace Engineering, and Engineering Mechanics University of

More information

Simulation of Aeroelastic System with Aerodynamic Nonlinearity

Simulation of Aeroelastic System with Aerodynamic Nonlinearity Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School

More information

Influence of the Heat Transfer on the Pressure Field in Radial Diffusers Flows

Influence of the Heat Transfer on the Pressure Field in Radial Diffusers Flows Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2010 Influence of the Heat Transfer on the Pressure Field in Radial Diffusers Flows Jose

More information

Numerical and Experimental Study on the Effect of Guide Vane Insertion on the Flow Characteristics in a 90º Rectangular Elbow

Numerical and Experimental Study on the Effect of Guide Vane Insertion on the Flow Characteristics in a 90º Rectangular Elbow Numerical and Experimental Study on the Effect of Guide Vane Insertion on the Flow Characteristics in a 90º Rectangular Elbow Sutardi 1, Wawan A. W., Nadia, N. and Puspita, K. 1 Mechanical Engineering

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

Analysis of Heat Transfer in Pipe with Twisted Tape Inserts

Analysis of Heat Transfer in Pipe with Twisted Tape Inserts Proceedings of the 2 nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, April 30 May 1, 2015 Paper No. 143 Analysis of Heat Transfer in Pipe with Twisted Tape Inserts

More information