Spectral Element Direct Numerical Simulation of Sodium Flow over a Backward Facing Step
|
|
- Martina Barnett
- 5 years ago
- Views:
Transcription
1 Spectral Element Direct Numerical Simulation of Sodium Flow over a Backward Facing Step ABSTRACT Jure Oder, Jernej Urankar, Iztok Tiselj Jožef Stefan Institute Jamova cesta 39 SI-1000, Ljubljana, Slovenia jure.oder@ijs.si, iztok.tiselj@ijs.si In this paper we present the preliminary results of direct numerical simulations of a turbulent flow of a liquid metal past the backward-facing step (BFS) with finite dimensions. The BFS geometry can be visualised as a channel where one of the walls has a shape of a step. The flow is flowing from the narrower part to the wider part. The simulations are performed in three dimensions. For the inflow boundary condition over the BFS, a fully developed turbulent velocity field is used. To obtain this fully developed turbulent inflow, a recycling boundary condition is used with which a plane of velocity values from the middle of domain is copied and used to set the inflow boundary condition. Simulations are performed with the NEK5000 code. The most notable feature of this code is the use of spectral elements to solve for velocity, temperature and any other passive scalar. It is an open source code developed by the Argonne National Laboratory. Spectral element method is a hybrid method between finite element method and a collocation spectral method. The method divides the computational domain into finite elements, within which a spectral method is used to solve for variables. This method allows for the use of spectral method in irregularly shaped geometries and to perform direct numerical simulations in such geometries. We performed some preliminary analysis of the calculated friction Reynolds number and the velocity profile before the step. Some analysis of the reattachment of flow was also performed. The main purpose of this work is to test the numerical set-up to later perform calculations with temperature field as a passive scalar. Dimensional walls with internal heating will be added to simulate the heat production in the walls. This work is part of our contribution to the SESAME project of Horizon2020 research programme. 1 INTRODUCTION In the last decades Direct Numerical Simulation (DNS) became an important research tool of the turbulent heat transfer [1]. Results of the DNS studies attained a similar role as experimental results and are typically used for verification and validation of the CFD codes and turbulent models in these codes. Particular attention was paid to the DNS of the fully 511.1
2 511.2 developed turbulent channel flow, as it reveals the basic mechanisms of the convective heat transfer between the fluid and the solid wall. In the recent decade, this topic is becoming more important also in the field of nuclear engineering, where it is used to study thermal fatigue related problems [2], [3]. The first DNSs of heat transfer in the channel flow were made by Kim and Moin [4] in 1989 and Kasagi et al. [5] in They were followed by a significant number of simulations performed at various Reynolds and Prandtl numbers [6]. Various papers of simulating the BFS geometry can be found in the literature. The more interesting for our case are the simulations performed with the spectral element code performed by Patera in 1984 [7]. The simulations performed therein were two dimensional and in a laminar flow regime. More recently similar simulations were performed in 2009 by Barri et al. [8]. Biswas et al. [9] performed simulations of a three dimensional BFS geometry, but at lower Reynolds numbers and for various expansion ratios. This paper presents the middle stage in the development of simulations of turbulent heat transfer in the BFS geometry. The simulations with turbulent channel flow with heat transfer were performed and presented by Oder et al. [10]. The goal of that paper was to test the code and to learn how to use it to obtain a fully developed turbulent flow in an infinite channel. In this regard, the BFS geometry in this paper differs significantly from the previous case. Only the velocity part was simulated in this paper and the periodic boundary condition cannot be used in this geometry. To develop the turbulent inflow, we used the so called recycling boundary condition. The values of velocity from previous time-step in the middle of domain were saved and then used as the boundary condition for the inflow. Because only one domain was used for the whole BFS geometry, the forcing, which was used in the previous paper, cannot be used in a similar way. Instead, we prescribed the volumetric inflow rate through the narrow part of the BFS geometry. Simulations presented in this work were performed with the numerical code NEK5000 [11], which is based on spectral element method. The spectral element code was already validated with results from our previous simulations [10][12] obtained with the spectral method code. Implementation of the spectral element code is part of our efforts within the Horizon2020 project SESAME, where we intend to use spectral elements to analyse the conjugate heat transfer in the backward facing step geometry, where spectral methods limited to very simple geometries cannot be applied. 2 MATHEMATICAL MODEL A sketch of the BFS geometry is presented in Figure 1. The geometry is chosen so as to conform to the geometry of the BFS experiment that is scheduled to be performed within the SESAME project. The experimental setup is being developed by the Karlsruhe Institute of Technology (KIT). The expansion ratio of the BFS geometry is fixed and set to 2. This means that the height of the step is equal to the height of the channel before the step. The outflow plane has a shape of a square. Computational domain of the test section is shown in Figure 1. Dimensionless lengths are obtained with normalisation with the half-width of the channel in front of the step. The plane from which the velocities are copied to develop the fully turbulent inflow is also marked in Figure 1. The velocities in the previous time-step are copied (and distributed between parallel processes), scaled to enforce the constant volumetric flow and then set as a
3 511.3 boundary condition at the inflow. The outflow is set to have a pressure equal to zero, but we also check the velocity so as to not allow for any inflow at that boundary. The other boundary conditions are set to no-slip walls. Copy velocity field Origin Outflow 2 z y x 6 6 Flow direction Figure 1: A sketch of the BFS geometry with coordinate system and dimensionless lengths. The inflow boundary condition is recycled from the plane that is marked The dimensionless governing equations of the fluid, normalized with channel half width h and friction velocity u τ are similar to the equations found in the paper of Kim et al. [13]: u + = 0, (1) u + t = (u+ u + ) + 1 Re τ 2 u + p, (2) where u + = (u +, v +, w + ) is the dimensionless velocity vector, p the dimensionless pressure and Re τ is the friction Reynolds number. No gravity is present in the equation. At this moment, the gravity is not planned to be implemented and simulated in the BFS geometry. 2.1 SPECTRAL ELEMENTS To simulate the BFS geometry, the open source implementation of spectral element method NEK5000 is used. This code is developed by the Argonne National Laboratory. The spectral element method is a hybrid method between finite element method and a spectral method. Similar to the finite element method, the domain is divided into finite elements, but on the other hand, similar to the spectral method, the solution within each element is obtained with a spectral method. The method was introduced by Patera in the year 1984 [7]. The spectral element method combines the benefits of both parent methods. The spectral method contributes the quick convergence rate and because of the finite elements, irregularly shaped geometries can be simulated. NEK5000 uses the collocation spectral method to obtain the solution within the element. These specially chosen points in which the solution is calculated within the element, are in the case of NEK5000 the Gauss-Lobatto-Legendre (GLL) points. These points are the N + 1 roots of the Legendre polynomial of degree N. The points are not uniformly distributed, but are rather clustered near the border of the element. Similar to finite element method, the spectral element method solved differential equations by constructing the stiffness and mass matrix.
4 511.4 The BFS domain in our case was divided first into regions and then into around 100,000 elements. The spectral method used 7 collocation points within each element and in each direction and the whole domain used around collocation points. Figure 2 is showing the distribution of spectral elements in the z = 0 plane. The size of the spectral elements after the step was varied in such a way, that the lower part of the domain had bigger elements in the y direction. The details of the sizes and number of elements in each region of BFS are given in Table 1. The region marked as Box 1 covers all elements from inlet till the length 6. Region Box 2 entails all elements from length 6 till the step. In region Box 3 are all the elements from the step onwards that are in the upper part of the domain and region Box 4 all the elements below region Box 3. Figure 2: The distribution and sizes of elements in the plane where z = 0. The elements after the step in the lower part of domain are larger in y direction. Table 1: Details of the distribution of the elements and the size of regions. Region Dimensions No. of elements Box Box Box Box
5 RESULTS The initial conditions for the simulation of the turbulent flow were set so that the average dimensionless velocity in the channel was around 20 with sinusoidal perturbations. As mentioned above, the velocity at the inflow was copied from other part of domain and scaled with the average velocity in Box 1 region and then multiplied by This velocity is roughly the velocity obtained from previous simulations of channel flow. The initial perturbations evolved into a stationary turbulent flow in whole domain in about 10 5 time steps. The average velocity component in each direction was calculated in each box (see Figure 2) to determine if steady state was reached. After the initial time, the simulation was running for another time steps before averaging was started. The averaging was performed during the time steps or roughly, from simulation time 29 to 49. A slice of average velocity component in x direction through the middle of the domain is shown in Figure 3. The simulation was running in parallel on 256 processors for about a week. Figure 3: The slice at z = 2 of averaged velocity in x direction. Our attention is first focused on the channel before the step. This is necessary to establish the simulated Re τ value. Figure 4 shows the average x velocity component and RMS of velocity components. These variables are calculated by averaging over time and over x spatial coordinate from inlet up to length 6 of the channel, with z fixed to 2 (the middle of the channel). Figure 4(a) shows the average component in the direction of the flow. It can be seen that the distribution of velocity is not completely symmetric yet, thus our averaging interval must be extended. The dashed line shows the average velocity component mirrored across the middle of the channel. From the average velocity Re τ can be calculated. To obtain a value for it, a parabola (which was required to be zero at the edges of the channel) was fitted to the first and last few values of averaged velocity component in Figure 4(a). On the lower wall (left side in Figure 4(a)) of the channel (Box 1) Re τ was equal to 219. On the upper plate Re τ was equal to 218. This is a slight asymmetry which can also be observed in Figure 4(a). In Figure 4(b), (c) and (d) the RMS calculated in similar way as the average velocity is shown. Again, the dashed line shows the values mirrored across the middle of the channel. Here also a slight asymmetry can be observed due to the short averaging time.
6 511.6 (a) (b) (c) Figure 4: Figures show the average velocity component in the x direction (a), and three components of RMS values. The dashed line shows the values mirrored through the middle of the channel. There are interesting structures in the flow behind the step. There are two vortices below the channel. This can be observed in Figure 5 which is showing the streamlines of the averaged flow. The streamlines are of course three-dimensional, but the figure shows the projection onto a plane. Two vortices can be observed. Right after the step the vortex rotating clockwise can be observed. In the lower edge of the step a smaller, weaker vortex rotating in the counter-clockwise direction can be seen. Figure 6 shows time averaged profile in the middle of the channel behind the step. We can observe that the maximum velocity at the step is not exactly in the middle of the channel. This might also be the reason, why profiles in Figure 4 are not symmetric. (d) Figure 5: Streamlines of the flow behind the step in the BFS geometry. The figure shows the two-dimensional projection of three-dimensional streamlines. The color marks the velocity of the flow.
7 511.7 Figure 6: x velocity component profile in the middle of the channel behind the step (plane z = 0). In the large vortex in the lower part of the domain the component of velocity in the x direction is directed opposite to the main flow. After long enough distance behind the BFS the flow must again attain the same general direction as before the step. Between these two regions there lies a point where the flow reattaches to the lower wall behind the step. Figure 7 shows the time averaged velocity component in the x direction at the lower wall behind the step (the first computational point above the lowest wall). From it we can conclude that the reattachment happens at about 15 units behind the step. Figure 7: The velocity component in the x direction in the plane y = (first point above y = 2 which is the lowest point in the domain). The figure shows, where the velocity in the general flow direction changes the sign.
8 CONCLUSIONS We presented some preliminary results of DNS of flow in BFS geometry. The fully turbulent flow was developed using the recycling boundary condition technique in which data from the middle of the domain is taken, copied, scaled and imposed as a boundary condition at the inflow. The analysis shows that the friction Reynolds number Re τ calculated in the simulation is higher than aimed for. From the presented results we are making a conclusion, that the buffer domain Box 2 might not be large enough for this case. Our spatial resolution might not be entirely sufficient for accurate DNS since our local friction Reynolds numbers are higher than the friction Reynolds number of the channel flow, where this resolution was verified. Grid sensitivity analysis remains one of our future tasks. In the future we are planning to add the solid walls and the temperature equations. The solid wall, which will be located after the step at the lowest wall, will have a fixed thickness and will have a distributed constant volumetric heat source. ACKNOWLEDGMENTS This project has received funding from the Euratom research and training programme under grant agreement No and from the Ministry of Higher Education, Science and Technology, Republic of Slovenia, Research programme P REFERENCES [1] Moin, P., Mahesh, K., Direct Numerical Simulation: A Tool in Turbulence Research. Annual Review of Fluid Mechanics, Vol. 30, pp [2] Brillant, G., Husson, S., Bataille, F., Subgrid-scale diffusivity: wall behaviour and dynamic methods, J. Appl. Mech., Trans ASME, 73(3): , [3] Aulery, F., Toutant, A., Brillant, G., Monod, R., Bataille, F., Numerical simulations of sodium mixing in a T-junction, Appl. Therm. Eng., 37:38-43, [4] Kim, J., Moin, P., Transport of Passive Scalars in a Turbulent Channel Flow. Turbulent Shear Flows VI, Springer-Verlag, Berlin, pp. 85. [5] Kasagi, N., Tomita, Y., Kuroda, A., Direct Numerical Simulation of Passive Scalar Field in a Turbulent Channel Flow. Jornal of Heat Transfer -Transactions of ASME, Vol. 114, [6] Kawamura, H., Ohsaka, K, Abe, H., Yamamoto, K., DNS of Turbulent Heat Transfer in Channel Flow with low to medium-high Prandtl number fluid. International Journal of Heat and Fluid Flow, Vol. 19, [7] Patera, A.T A spectral element method for fluid dynamics: laminar flow in a channel expansion. Journal of computational Physics 54, 468. [8] Barri, M.; El Khoury, G. K.; Andersson, H. I. & Pettersen, B. DNS of backward-facing step flow with fully turbulent inflow. Int. J. Numer. Meth. Fluids, 2009, 64,
9 511.9 [9] Biswas, G.; Breuer, M.; Durst, F. Backward-Facing Step Flows for Various Expansion Ratios at Low and Moderate Reynolds Numbers. Journal of Fluids Engineering, 2004, 126, [10] Oder, J.; Urankar, J.; Iztok T. Jenčič, I. (Ed.) Spectral Element Direct Numerical Simulation of Heat Transfer in Turbulent Channel Sodium Flow 24th International Conference Nuclear Energy for New Europe - NENE 2015, Portorož, Slovenia, September 14-17, 2015 [11] Fischer, P.F., Lottes, J.W., Kerkemeier, S.G., Nek5000 web page. (accessed 27. August 2015) [12] Tiselj, I., Bergant, R., Mavko, B., Bajsič, I., Hetsroni, G., DNS of Turbulent Heat Transfer in Channel Flow With Heat Conduction in the Solid Wall. Journal of Heat Transfer - ASME, 123, [13] Kim, J.; Moin, P. & Moser, R. Turbulence statistics in fully developed channel flow at low Reynolds number. Journal of Fluid Mechanics 1987, 177,
INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT FLOW NEAR THE HEATED WALL
International Conference Nuclear Energy for New Europe 22 Kranjska Gora, Slovenia, September 9-2, 22 www.drustvo-js.si/gora22 INFLUENCE OF THE BOUNDARY CONDITIONS ON A TEMPERATURE FIELD IN THE TURBULENT
More informationNumerical simulations of heat transfer in plane channel flow
Numerical simulations of heat transfer in plane channel flow Najla EL GHARBI 1, 3, a, Rafik ABSI 2, b and Ahmed BENZAOUI 3, c 1 Renewable Energy Development Center, BP 62 Bouzareah 163 Algiers, Algeria
More 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 informationCHARACTERIZATION OF THE SECONDARY FLOW IN TURBULENT RECTANGULAR DUCTS WITH VARYING ASPECT RATIO
June - July, 21 Melbourne, Australia 9 8A-1 CHARACTERIZATION OF THE SECONDARY FLOW IN TURBULENT RECTANGULAR DUCTS WITH VARYING ASPECT RATIO Ricardo Vinuesa Linné FLOW Centre, KTH Mechanics SE-1 Stockholm,
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationDIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING
DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio
More informationInfluence of high temperature gradient on turbulence spectra
NIA CFD Conference, Hampton, August 6-8, 2012 Future Directions in CFD Research, A Modeling and Simulation Conference Influence of high temperature gradient on turbulence spectra Sylvain Serra Françoise
More informationnek5000 massively parallel spectral element simulations
nek5000 massively parallel spectral element simulations PRACE Scientific Seminar HPC Boosts Science, 22th February 2011 P. Schlatter & D. S. Henningson Linné Flow Centre, KTH Mechanics Fluid flows Tornado,
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 informationFALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
Proceedings of the International Conference on Mechanical Engineering 2 (ICME2) 8-2 December 2, Dhaka, Bangladesh ICME-TH-6 FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
More informationVery large-scale structures observed in DNS of turbulent channel flow with passive scalar transport
5 th Australasian Fluid Mechanics Conference 3-7 December 4 Very large-scale structures observed in DNS of turbulent channel flow with passive scalar transport H. Kawamura, H. Abe and Y. Matsuo Department
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationDIRECT NUMERICAL SIMULATIONS OF VARIABLE-ASPECT-RATIO TURBULENT DUCT FLOWS AT LOW TO MODERATE REYNOLDS NUMBERS
DIRECT NUMERICAL SIMULATIONS OF VARIABLE-ASPECT-RATIO TURBULENT DUCT FLOWS AT LOW TO MODERATE REYNOLDS NUMBERS Ricardo Vinuesa MMAE Department Illinois Institute of Technology Chicago, IL 60616, USA rvinuesa@hawk.iit.edu
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationNUMERICAL SIMULATIONS OF FLOW, HEAT TRANSFER AND CONJUGATE HEAT TRANSFER IN THE BACKWARD-FACING STEP GEOMETRY
Mecánica Computacional Vol XXXII, págs. 3265-3278 (artículo completo) Carlos G. García Garino, Aníbal E. Mirasso, Mario A. Storti, Miguel E. Tornello (Eds.) Mendoza, Argentina, 19-22 Noviembre 2013 NUMERICAL
More informationNumerical analysis of a fully developed non-isothermal particle-laden turbulent channel flow
Arch. Mech., 3, 1, pp. 77 91, Warszawa 11 Numerical analysis of a fully developed non-isothermal particle-laden turbulent channel flow M. JASZCZUR Department of Fundamental Research in Energy Engineering
More 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 informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationFLOW CHARACTERISTICS IN A VOLUTE-TYPE CENTRIFUGAL PUMP USING LARGE EDDY SIMULATION
FLOW CHARACTERISTICS IN A VOLUTE-TYPE CENTRIFUGAL PUMP USING LARGE EDDY SIMULATION Beomjun Kye Keuntae Park Department of Mechanical & Aerospace Engineering Department of Mechanical & Aerospace Engineering
More informationHeat Transfer and Fluid Flow Characteristics of Separated Flows Encountered in a Backward-Facing Step Under the Effect of Suction and Blowing
E. Abu-Nada e-mail: eiyad@hu.edu.jo A. Al-Sarkhi B. Akash I. Al-Hinti Department of Mechanical Engineering, Hashemite University, Zarqa, 13115, Jordan Heat Transfer and Fluid Flow Characteristics of Separated
More informationROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 ROLE OF LARGE SCALE MOTIONS IN TURBULENT POISEUILLE AND COUETTE FLOWS Myoungkyu Lee Robert D. Moser
More informationSTABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET. N. Vetcha, S. Smolentsev and M. Abdou
STABILITY ANALYSIS FOR BUOYANCY-OPPOSED FLOWS IN POLOIDAL DUCTS OF THE DCLL BLANKET N. Vetcha S. Smolentsev and M. Abdou Fusion Science and Technology Center at University of California Los Angeles CA
More informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationColloquium FLUID DYNAMICS 2013 Institute of Thermomechanics AS CR, v.v.i., Prague, October 23-25, 2013 p.1
Colloquium FLUID DYNAMICS 2013 Institute of Thermomechanics AS CR, v.v.i., Prague, October 23-25, 2013 p.1 ON THE REYNOLDS NUMBER ROLE IN STRUCTURE OF RECIRCULATION ZONE BEHIND BACKWARD FACING STEP IN
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationDNS of a buoyancy-dominated turbulent flow using an immersed boundary method
Center for Turbulence Research Annual Research Briefs 2008 231 DNS of a buoyancy-dominated turbulent flow using an immersed boundary method By S. Kang, G. Iaccarino AND F. Ham 1. Motivation and objectives
More informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More 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 informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationEffects of Downstream Wall Heating on the Flow Characteristics of Inclined Backward Facing Step
International Journal of Pure and Applied Mathematics Volume 119 No. 12 2018, 37-41 ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Effects of Downstream Wall Heating on the Flow Characteristics
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationEffect on heat transfer for laminar flow over Backward Facing Step with square cylinder placed inside using Higher Order Compact Scheme
ISSN (e): 2250 3005 Volume, 06 Issue, 01 January 2016 International Journal of Computational Engineering Research (IJCER) Effect on heat transfer for laminar flow over Backward Facing Step with square
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationSG Turbulence models for CFD
SG2218 2012 Turbulence models for CFD Stefan Wallin Linné FLOW Centre Dept of Mechanics, KTH Dept. of Aeronautics and Systems Integration, FOI There are no simple turbulent flows Turbulent boundary layer:
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 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 informationS. Varela, G. Usera, A. Vernet & J. A. Ferré. Abstract
Advances in Fluid Mechanics X 313 Numerical simulation of heat transfer and temperature distribution in a printed circuit board enclosure model for different geometries and Reynolds numbers 1 2 1 1 S.
More informationCHALMERS UNIVERSITY OF TECHNOLOGY
Report No. /3 CHALMERS UNIVERSITY OF TECHNOLOGY Department of Thermo and Fluid Dynamics Large-eddy simulations of turbulent flows in stationary and rotating channels and in a stationary square duct by
More 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 informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationHelical Coil Flow: a Case Study
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Helical Coil Flow: a Case Study Marco Cozzini Renewable Energies and Environmental Technologies (REET) Research Unit, Fondazione Bruno Kessler
More informationA new numerical approach for Soret effect on mixed convective boundary layer flow of a nanofluid over vertical frustum of a cone
Inter national Journal of Pure and Applied Mathematics Volume 113 No. 8 2017, 73 81 ISSN: 1311-8080 printed version); ISSN: 1314-3395 on-line version) url: http://www.ijpam.eu ijpam.eu A new numerical
More informationInvestigation of reattachment length for a turbulent flow over a backward facing step for different step angle
MultiCraft International Journal of Engineering, Science and Technology Vol. 3, No. 2, 2011, pp. 84-88 INTERNATIONAL JOURNAL OF ENGINEERING, SCIENCE AND TECHNOLOGY www.ijest-ng.com 2011 MultiCraft Limited.
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationHISTORY EFFECTS FOR CAMBERED AND SYMMETRIC
HISTORY EFFECTS FOR CAMBERED AND SYMMETRIC WING PROFILES A. Tanarro, R. Vinuesa and P. Schlatter Linné FLOW Centre, KTH Mechanics and Swedish e-science Research Centre (SeRC), SE-1 44 Stockholm, Sweden
More informationGeneration of initial fields for channel flow investigation
Generation of initial fields for channel flow investigation Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (Mai 2) In the framework of the DFG-funded research
More informationOptimization and control of a separated boundary-layer flow
Optimization and control of a separated boundary-layer flow Journal: 2011 Hawaii Summer Conferences Manuscript ID: Draft lumeetingid: 2225 Date Submitted by the Author: n/a Contact Author: PASSAGGIA, Pierre-Yves
More informationTurbulent momentum and heat transfer in ducts of rhombic cross section
Turbulent momentum and heat transfer in ducts of rhombic cross section Naoya Fukushima and Nobuhide Kasagi Department of Mechanical Engineering, The University of Tokyo Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656,
More 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 informationTurbulent Heat Transfer in A Channel Flow with Arbitrary. Directional System Rotation
Turbulent Heat Transfer in A Channel Flow with Arbitrary Directional System tation Haibin Wu * and Nohuhide Kasagi Department of mechanical engineering, The University of Tokyo, Hongo 7-3-1, Bunkyo ku,
More informationSimulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code
1 st. Annual (National) Conference on Industrial Ventilation-IVC2010 Feb 24-25, 2010, Sharif University of Technology, Tehran, Iran IVC2010 Simulation and improvement of the ventilation of a welding workshop
More informationPROOF COPY JHR. Three-Dimensional Simulations of Enhanced Heat Transfer in a Flat Passage Downstream From a Grooved Channel
M. Greiner Professor of Mechanical Engineering, University of Nevada, Reno, Nevada 89557 e-mail: greiner@unr.edu P. F. Fischer Mathematics and Computer Science Division, Argonne National Laboratory, Argonne,
More information3D Numerical Study on Laminar Forced Convection in V-Baffled Square Channel
American Journal of Applied Sciences 10 (10): 1287-1297, 2013 ISSN: 1546-9239 2013 Boonloi and Jedsadaratanachai, This open access article is distributed under a Creative Commons Attribution (CC-BY) 3.0
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationApplication of wall forcing methods in a turbulent channel flow using Incompact3d
Application of wall forcing methods in a turbulent channel flow using Incompact3d S. Khosh Aghdam Department of Mechanical Engineering - University of Sheffield 1 Flow control 2 Drag reduction 3 Maths
More informationINSTITUTTET FOR BYGNINGSTEKNIK DEPT. OF BUILDING TECHNOLOGY AND STRUCTURAL ENGINEERING AALBORG UNIVERSITET AAU AALBORG DANMARK
INSTITUTTET FOR BYGNINGSTEKNIK DEPT. OF BUILDING TECHNOLOGY AND STRUCTURAL ENGINEERING AALBORG UNIVERSITET AAU AALBORG DANMARK Lars Davidson and Peter V. Nielsen A Study of Laminar Backward-Facing Step
More informationSECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES
SECONDARY MOTION IN TURBULENT FLOWS OVER SUPERHYDROPHOBIC SURFACES Yosuke Hasegawa Institute of Industrial Science The University of Tokyo Komaba 4-6-1, Meguro-ku, Tokyo 153-8505, Japan ysk@iis.u-tokyo.ac.jp
More 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 informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
More informationEnergy fluxes in turbulent separated flows
Journal of Physics: Conference Series PAPER OPEN ACCESS Energy fluxes in turbulent separated flows To cite this article: J.-P. Mollicone et al 2016 J. Phys.: Conf. Ser. 759 012003 View the article online
More informationCranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom. Cranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom
MultiScience - XXX. microcad International Multidisciplinary Scientific Conference University of Miskolc, Hungary, 21-22 April 2016, ISBN 978-963-358-113-1 NUMERICAL INVESTIGATION OF AN INCOMPRESSIBLE
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationThermal Nanofluids flow in Corrugated Facing Step - A Review
Thermal Nanofluids flow in Corrugated Facing Step - A Review Kafel Azeez Mohammed Department of Mechanical and Manufacturing Engineering, UPM, MALAYSIA. Abstract In this paper, the studies of thermal fluid
More informationEffects of Free-Stream Vorticity on the Blasius Boundary Layer
17 th Australasian Fluid Mechanics Conference Auckland, New Zealand 5-9 December 2010 Effects of Free-Stream Vorticity on the Boundary Layer D.A. Pook, J.H. Watmuff School of Aerospace, Mechanical & Manufacturing
More informationEnhancement of Heat Transfer by an Electric Field for a Drop Translating at Intermediate Reynolds Number
Rajkumar Subramanian M. A. Jog 1 e-mail: milind.jog@uc.edu Department of Mechanical, Industrial, and Nuclear Engineering, University of Cincinnati, Cincinnati, OH 45221-0072 Enhancement of Heat Transfer
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationNumerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Numerical Solutions for the Lévêque Problem of Boundary Layer Mass or Heat Flux Ekkehard Holzbecher Weierstrass Institute for Applied
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationCOMPUTACIONAL SIMULATION OF TWO-DIMENSIONAL TRANSIENT NATURAL CONVECTION IN VOLUMETRICALLY HEATED SQUARE ENCLOSURE
COMPUTACIONAL SIMULATION OF TWO-DIMENSIONAL TRANSIENT NATURAL CONVECTION IN VOLUMETRICALLY HEATED SQUARE ENCLOSURE Camila Braga Vieira, camila@lasme.coppe.ufrj.br Jian Su, sujian@lasme.coppe.ufrj.br Universidade
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationCurriculum Vitae of Sergio Pirozzoli
Curriculum Vitae of Sergio Pirozzoli Address University of Rome La Sapienza Department of Mechanical and Aerospace Engineering Via Eudossiana 18 00184, Roma Contact tel.: +39 06 44585202 fax : +39 06 4881759
More informationFlow field in the compressor blade cascade NACA
Flow field in the compressor blade cascade NACA 65-100 Tomáš Turek Thesis Supervisor: Ing. Tomáš Hyhlík, Ph.D. Abstract An investigation is made into the effects of a flow field in the compressor blade
More informationSYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL
ISSN 1726-4529 Int j simul model 7 (2008) 4, 186-197 Original scientific paper SYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL Ternik, P. University
More informationPeriodic planes v i+1 Top wall u i. Inlet. U m y. Jet hole. Figure 2. Schematic of computational domain.
Flow Characterization of Inclined Jet in Cross Flow for Thin Film Cooling via Large Eddy Simulation Naqavi, I.Z. 1, Savory, E. 2 and Martinuzzi, R. J. 3 1,2 The Univ. of Western Ontario, Dept. of Mech.
More informationComputational Fluid Dynamics Modelling of Natural Convection in Copper Electrorefining
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Abstract Computational Fluid Dynamics Modelling of Natural Convection in Copper Electrorefining A computational
More informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
More informationRayleigh-Taylor Unstable Flames
Rayleigh-Taylor Unstable Flames Elizabeth P. Hicks 1,2 and Robert Rosner 2 CIERA, Northwestern University 1 University of Chicago 2 CIERA Conference: September 2, 2011 1 Type Ia Supernovae Image: NASA
More informationEffect of near-wall treatments on airflow simulations
Proceedings of 29 International Conference on Computational Methods for Energy Engineering and Environment: ICCM3E. Sousse, Tunisia, 2-22 November, 29, pp. Effect of near-wall treatments on airflow simulations
More informationCONVECTIVE HEAT TRANSFER OVER A WALL MOUNTED CUBE AT DIFFERENT ANGLE OF ATTACK USING LARGE EDDY SIMULATION
S301 CONVECTIVE HEAT TRANSFER OVER A WALL MOUNTED CUBE AT DIFFERENT ANGLE OF ATTACK USING LARGE EDDY SIMULATION by Habibollah HEIDARZADEH, Mousa FARHADI *, and Kurosh SEDIGHI Faculty of Mechanical Engineering,
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 informationOptimum algorithm for channel flow analysis in direct numerical simulation method
International Journal of Civil Engineering Optimum algorithm for channel flow analysis in direct numerical simulation method M. R. Kavianpour,*, E. Rajabi Received: May 0, Revised: February 0,Accepted:
More informationVortex wake and energy transitions of an oscillating cylinder at low Reynolds number
ANZIAM J. 46 (E) ppc181 C195, 2005 C181 Vortex wake and energy transitions of an oscillating cylinder at low Reynolds number B. Stewart J. Leontini K. Hourigan M. C. Thompson (Received 25 October 2004,
More informationContribution of Reynolds stress distribution to the skin friction in wall-bounded flows
Published in Phys. Fluids 14, L73-L76 (22). Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Koji Fukagata, Kaoru Iwamoto, and Nobuhide Kasagi Department of Mechanical
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
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 informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
More informationApplied Thermal and Fluid Engineering. Energy Engineering (Thermal Engineering Laboratory)
Applied Thermal and Fluid Engineering Energy Engineering (Thermal Engineering Laboratory) Professor Assoc. Professor Hajime Nakamura Shunsuke Yamada Outline of Research In our laboratory, we have been
More informationNumerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell Method
Send Orders for Reprints to reprints@benthamscience.net 30 The Open Transport Phenomena Journal, 2013, 5, 30-41 Open Access Numerical Simulation of Rotating Turbulent Channel Flow by the Vortex in Cell
More informationPatterns of Turbulence. Dwight Barkley and Laurette Tuckerman
Patterns of Turbulence Dwight Barkley and Laurette Tuckerman Plane Couette Flow Re = U gap/2 ν Experiments by Prigent and Dauchot Re400 Z (Spanwise) 40 24 o Gap 2 Length 770 X (Streamwise) Examples: Patterns
More informationCFD STUDIES IN THE PREDICTION OF THERMAL STRIPING IN AN LMFBR
CFD STUDIES IN THE PREDICTION OF THERMAL STRIPING IN AN LMFBR K. Velusamy, K. Natesan, P. Selvaraj, P. Chellapandi, S. C. Chetal, T. Sundararajan* and S. Suyambazhahan* Nuclear Engineering Group Indira
More informationFlow characteristics of curved ducts
Applied and Computational Mechanics 1 (007) 55-64 Flow characteristics of curved ducts P. Rudolf a *, M. Desová a a Faculty of Mechanical Engineering, Brno University of Technology,Technická,616 69 Brno,
More information