Calculations on a heated cylinder case
|
|
- Sydney Patrick
- 5 years ago
- Views:
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+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationExplicit 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 informationValidation 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 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 informationTutorial 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 informationPHYSICAL 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 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 informationOn 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 informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationHeat 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 informationExternal 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 informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationPrediction of unsteady heat transfer from a cylinder in crossflow
Center for Turbulence Research Proceedings of the Summer Program 202 07 Prediction of unsteady heat transfer from a cylinder in crossflow By S. T. Bose, B. C. Wang AND M. Saeedi The accuracy of a tensorial
More informationIntroduction to Turbulence AEEM Why study turbulent flows?
Introduction to Turbulence AEEM 7063-003 Dr. Peter J. Disimile UC-FEST Department of Aerospace Engineering Peter.disimile@uc.edu Intro to Turbulence: C1A Why 1 Most flows encountered in engineering and
More informationNumerical 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 informationForced 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 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 informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More 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 informationLARGE-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 informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationTutorial 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 informationActive 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 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 informationA 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 informationIntroduction 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 informationTURBULENT 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 informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More 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 information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More 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 informationρ 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 informationCFD 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 informationNumerical 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 informationA 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 informationME 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 informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More 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 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 informationDay 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 informationConjugate 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 informationOPTIMAL 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 informationExplicit 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 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 informationA 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 informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
More 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 informationNumerical 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 informationChapter 3 NATURAL CONVECTION
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 Chapter 3 NATURAL CONVECTION Mehmet Kanoglu Copyright The McGraw-Hill Companies,
More informationProceedings 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 informationForced 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 informationTABLE 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 informationWall 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 informationA 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 informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationCFD 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 informationFluid 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 informationDepartment 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 informationStudy 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 informationApplication 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 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 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 informationPART 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 informationElliptic 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 informationChapter 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 informationComparison 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 informationME 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 informationThe 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 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 informationStudies 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 informationMestrado 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 informationRANS simulations of rotating flows
Center for Turbulence Research Annual Research Briefs 1999 257 RANS simulations of rotating flows By G. Iaccarino, A. Ooi, B. A. Pettersson Reif AND P. Durbin 1. Motivation and objectives Numerous experimental
More informationComparison 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 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 informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationRSM 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 informationInitial 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 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 information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationNumerical 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 informationWALL 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 informationLaminar 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 informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationTurbulence 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 informationForced 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 information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationNumerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation
, pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,
More informationDetailed 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 informationExperimental 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 informationConvection 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 informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationInfluence 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 informationNumerical 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 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 informationAnalysis 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