Simulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code
|
|
- Catherine Heath
- 5 years ago
- Views:
Transcription
1 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 using a Finite volume scheme code Mani Mahdinia M.Sc. Student Mahdinia@mech.sharif.edu Amir Ghasemi M.Sc. Student AGhasemy.v@gmail.com Bijan Farhanieh Professor bfia@sharif.edu Sharif University of Technology, Mechanical Engineering School, Tehran,Iran Abstract A Finite-volume CFD code is developed for the solution of the Navier-Stokes equations in an industrial ventilation domain. The point at the middle of the domain, having a temperature of 500 C, is the place where a welding process is being accomplished. Due to the high temperatures at the vicinity, three fans (one supply and two exhausts) are used. The purpose of this study is to investigate the effect of different parameters on the local ventilation and developing the most effective ventilation scheme for this applied problem. In order to find the velocity components and the temperatures, the equations are solved in a Cartesian coordinate system and a collocated grid. To prevent the checkerboard distribution, velocity components on the faces of the control volumes are obtained by the Rhie and Chow interpolation scheme. SIMPLE algorithm is used for the modeling of the pressure-velocity coupling. Convection fluxes have been discretised, using the upwind scheme and the resulting equations have been transformed into a tri-diagonal matrix with the use of the ADI technique and solved with the Tomas Tridiagonal algorithm. For the sake of maximum accuracy, the grid has been refined at the points of larger velocity gradients. The developed code accuracy has been verified against the cavity flow test case and checked for the grid independency. The effect of different variables like the wall lengths, boundary conditions and flow rates has been studied. Using all the above mentioned techniques and results, an effective ventilation scheme has been proposed for the ventilation of the welding stand. Keywords: Industrial Ventilation, Computational Fluid Dynamics, SIMPLE Algorithm, Quick Scheme, Rhie- Chow Interpolation Scheme 1) Introduction The improvements in the industrial technologies, has lead to a growing need for the more carefully designed and manipulated methods of ventilation. This includes manufacturing industries like car production factories, subway environment ventilation, clean room HVAC and so many others. Amongst these diverse applications, those that involve higher temperatures, like the one in an industrial welding stand may need more attention and care. The temperature in the vicinity of a welding point may become as high as 1000 degrees Celsius, which may result in an enormous temperature rise due to the accumulation of the generated heat, leading to an uncomfortable or even harmful environment for the workers, if not properly cooled. Since different cases of industrial ventilations, involve many kinds of flow configurations and geometries, it is nearly impossible to rely completely on the experimental methods. Regarding this fact, other tools such as Computational Fluid Dynamics techniques, have gained more attention in this field, specifically in the recent decades, thanks to the improvements in the computer technology. Another beneficial factor regarding the CFD methods is that it enables one to analyze the performance of the designed ventilation system for various arrangement and parameters and also select the best. To achieve such an optimization using an experimental approach may demand a huge amount of resources, even if it is possible. The problem at hand pertains to the cooling and ventilation of a welding stand, used in the car manufacturing industries. The geometry of the problem is shown in Fig. (1). A welding process is being accomplished in the middle of the domain and the three surrounding internal walls isolates the welder from its surroundings. The article organization is as follows. After this brief introduction, the governing equations and the boundary conditions are derived for the problem at hand. Then the discretization method and the numerical solution algorithm are investigated in detail. After this part, the code accuracy verification is accomplished by testing its results for the well-known cavity flow. In the next section, the grid independency of the results is investigated. Finally, the simulation results are presented and the effect of various parameters is taken into consideration. The effect of the inlet and outlet relative velocities and Reynolds numbers has been determined. Also the impact of the type of the outer B.C., considering the symmetry, wall and free boundary condition types has been investigated. Other geometrical parameters like the length and location of the internal walls, the size and location of the fans and the welding point position has also been studied. Finally, the effect of the variance of the fluid properties like the Prandtl number due to the temperature changes was determined. Using these results, the best possible configuration amongst the ones considered is chosen for the applied use. This arrangement may result in the best efficiency of the equipment.
2 equations are solved in the left half of the domain. The boundary conditions of the flow are shown in Fig.2. The right boundary (No. 7) is considered to be a symmetry line as stated above. The lower boundary (No. 3) is an impermeable wall with no slip conditions for velocity and zero-flux condition for heat. The outer edges (No. 4) are considered as free boundaries with zero heat flux, based on the assumption that the changes of temperature in the vicinity of these boundaries are negligible relative to the middle point of the domain. The No. 1 and No. 2 boundaries are considered as inlet and outlet boundaries (constant velocity) respectively, with a prescribed temperature in the inlet and an extrapolated temperature for the outlet. For the internal walls (No.6), a no slip condition on velocity and a convective condition for temperature are imposed. Fig 1 - Computational domain of the problem 2) Governing Equations The air flow is assumed to be incompressible, twodimensional and steady. The continuity, x and y momentum and the energy equation for the assumed conditions are + =0 (1) + = 1 + ( + ) (2) + = 1 + ( + ) (3) ( + )= ( + ) (4) It should be mentioned that due to the low velocities in the solution domain, the dissipation in the energy equation (4) has been neglected. To handle the flow problem in a more general fashion, the equations have been non-dimensionalized, by using the following definitions =, =, =, =, = Substituting (5) in (1) to (4) gives (5) + =0 (6) + = + 1 ( + ) (7) + = + 1 ( + ) (8) ( + )=. ( + ) (9) In which the Reynolds and Prandtl numbers are defined as = and =. In the next sections, the star superscripts on the variables are dropped for ready reference. 3) Boundary conditions To reduce the computational cost of the solution, the Fig 2 - Boundary conditions of the domain 4) Solution algorithm Equations (6) to (9) are discretised using the finite volume method, on a collocated grid. Such a grid is based on storing the velocities and pressures in the same locations in the mesh, which is in contrast to the commonly used staggered grid. Incorporating the collocated grid has the advantages of having a lower amount of storage required and also the ease of programming. Its disadvantage is that it may lead to a checker-board pressure distribution in the domain. To overcome such a problem, the Rhie-Chow interpolation scheme [1] has been incorporated to evaluate the volumetric fluxes on the cell surfaces. For the pressure-velocity coupling, the SIMPLE [2] algorithm has been used. To find the quantities on the surfaces of the control volumes, the Upwind method has been used. Using the above mentioned algorithms and techniques, results in a linear system of equations, which have been transformed into a tri-diagonal matrix with the use of the ADI technique and solved with the Tomas Tridiagonal algorithm. For the sake of maximum accuracy, the grid has been refined at the points of larger velocity gradients. A C++ code is developed to implement the mentioned solution process.
3 5) Code verification The developed code has been verified against the wellknown test case of the cavity flow. The geometry of the flow is shown in Fig. 3. In the figure, the side and bottom walls are stationary and the no slip condition is imposed on them. The top boundary on the other hand, is moving with a constant velocity of 1 m/s to the right. The x-velocity profiles on the vertical middle line of the geometry, i.e. =0.5, have been verified with the experimental results of the well-known work of Ghia et. al.[3] and shown in Fig. 5, for the two Reynolds numbers. Fig 3 Geometry of the cavity flow The domain has been discretised using a grid for the Reynolds number of 400 and a grid for the Reynolds number of To show the flow configuration more clearly, the streamlines are shown in Fig. 4, for both of the Reynolds numbers. Fig 5 - x-velocity profiles on = 0.5 line Re=400 Re=1000 Fig 4 - Streamlines Re=400 Re=1000 (with the definition Re=uL x /ν) As can be seen, the results of the simulation are in good agreement with the experimental results. The y-velocity profiles also agree with the empirical data in a similar manner, but to save space the results have not been presented here. 5) Grid independency A structured grid has been used to transform the PDE equations to a set of linear algebraic equations. To achieve an acceptable level of accuracy and also to avoid unnecessarily high computational costs, the five cases of different grid sizes have been considered; , , , and The results for the x-velocity on the = 4.8 line and y- velocity on = 1.2 line (based on the coordinate system origin at the lower left corner of Fig. 1) have been compared for the stated grid sizes, in Fig. 6 and Fig. 6 respectively. As illustrated in the figure, the meshes more refined than the size of , collapse
4 Fig 6 - Grid independency results x-velocity on the = 4.8 line y-velocity on = 1.2 line together. So we choose this mesh size as our basic grid size. To achieve the final maximum accuracy, the grid is refined in the places where high gradients of the velocity are expected, as shown in Fig. 7. 6) Simulation results The velocity variables are nondimensionalized with respect to the inlet velocity so the inlet velocity has been assumed to have a value of 1m/s. Inlet Reynolds number of 200, which is defined, based on the characteristic length of H has been considered. The outlet velocity is calculated based on the ratio of the inlet Reynolds number to that of the outlet (outlet Reynolds number is equal to 80). Flow streamlines and temperature Contours are shown in Fig. 8. and Fig. 9, respectively. As can been seen, two large kidney-shaped circulating flow regions are formed between the internal walls as a result of flow entrapment between them which also prevents effective ventilation. Another point is that according to Fig. 8, the flow has the ability to travel into and out of the domain, due to assumption of the free boundary for the outer edges. Also, since the net outflow of the inlets and outlets together are a negative quantity, the net flow on the outer edges (excluding the inlets and outlets) should be towards the outside of the domain. This forces the streamlines to become denser near the outlets of the Fig 7 - Grid refinement Basic grid locally refined grid refined Fig 8 - Flow streamlines domain. The temperature contours are shown in Fig.9. As can be seen from the figure, the constant temperature lines are perpendicular to the lower wall, in accordance with the insulated B.C. assumption. Also, gradient, far
5 the ambient temperature, i.e. larger distance results in the loss of the flow velocity at the boundaries, without any benefits. Fig 9 - Temperature contour from the welding point. Also there exists a for the side free boundaries show the same behavior, which confirms our prediction of negligible temperature large temperature gradient near the center of the domain, exhibiting effects of the ventilation near the point. 7) Effect of parameters In this section, the effect of various parameters on the flow field is investigated in order to reach at the best possible configuration for the system at hand. 7-1) Boundary conditions Here, three various configurations are considered: The welding stations are positioned in an array, having a large distance apart. This gives us the free boundary condition. The same as previous, but the distances are small so that the symmetry B.C. is applicable. (c) Each welding station is separated from the surrounding using partitions. This results in wall B.C. The flow stream lines and temperature contours are demonstrated in Fig. 10. The differences between the streamlines in parts to (c) is largely due to the large vortices at the top corners of the domain, which result from flow area expansion in the and (c) cases. In case, the flow entrance from the upper boundary and its exhaust from the vertical boundaries, prevents the formation of any large vortices. The temperature contours show a larger temperature for the case of the free boundary. This is due to the outflow of a large amount of the current before getting close enough to the welding point. Amongst the two other configurations the symmetry boundary, which stands for a small-distance array of welding stands in a workshop, gives a lower temperature distribution in the domain. This is in contrast to the no-slip boundary conditions in the case of the wall boundaries (a partition-separated array), that slows down the flow around the welding point, that lowers the effect of the ventilation. These results suggest using an array of stands in the welding workshop, where the stands are positioned in a way that the minimum distance between them is chosen so that the temperature at the boundary just reaches 7-2) Inlet and outlet Reynolds numbers The effects of the changes of the inlet and outlet boundary conditions are shown in Fig. 11. The main case is repeated for ease of comparison. Cases and (c) represent the configurations in which the inlet and outlet Reynolds numbers are the same and are equal to 80 and 20 respectively. In case the equivalence of the inlet and outlet Reynolds numbers, reduces the tendency of the current to cross the boundaries of the domain, which moves the flow towards the formation of circulating vortexes in the domain corners. In case (c) the lowered Reynolds numbers in the inlet and outlet, results in a streamline pattern which interestingly resembles the creeping flow, where the flow lines smoothly turns around the obstacles without separation. The reduction of the Reynolds number from case to and from to (c) results in the shrinkage of the vortexes amongst the internal walls, in such a way that the big vortex in breaks down into two smaller vortices in and finally at (c) it diminishes to a very small size. The velocity reduction in the flow domain results in an increased magnitude of the temperatures in cases and (c). 7-3) Internal wall configuration Regarding the internal wall configuration, one may consider attaching the walls together to prevent harmful effects of the radiation to the surrounding people, but unfortunately this feature may prevent the effective ventilation as well. This conflict motivates us to investigate this parameter. The velocity vectors and temperature contours in the domain, in such a configuration is shown in Fig. 12. As can be seen, the velocity magnitude decreases in the middle region and the current is entrapped between the walls, as was expected, which results in the diminution of the convection cooling and causes in the temperature rise in the vicinity of the welding point from 320 C in Fig. 10 part (aʹ) to 400 C in Fig. 12 part. Since the temperature difference in the two circumstances is too large, the idea of attaching the walls, may not work as a proper solution and we have to look for other methods to circumvent the radiation hazard. 7-4) Dependence on the Prandtl number The analysis here may become subjected to several variations in the parameters. The uncertainty in the parameters can be due to the seasonal temperature changes. Here we investigate the effect of the variation of the Prandtl number in the equations. The temperature contours for two different Prandtl numbers are shown in Fig. 13. As can be seen the lower Prandtl number results in the penetration of the high temperatures in a bigger region in the domain.
6 (c) (aʹ) (bʹ) (cʹ) Fig Effect of different boundary conditions to (c) Streamlines for the free, symmetry and wall B.C.s and (aʹ) to (cʹ) Related temperature contours (c) (aʹ) (bʹ) (cʹ) Fig Effect of inlet and outlet Reynolds numbers to (c) Streamlines with Rein=200, Reout=80 Rein=80, Reout=80 (c) Rein=20, Reout=20 - (aʹ) to (cʹ) Related temperature contours
7 Fig 12 - Effect of the internal wall configuration Velocity contours Temperature contour 8) Conclusions A finite volume CFD code has been developed for the simulation of the ventilation of a welding environment. Since the welding configuration considered in this study is a common installment in many welding stands, its study and optimization may be of a general interest. We have used this code as a means for a virtual laboratory to investigate the effect of various parameters on the function of the system. Decision about the position of welding stands relative to each other, i.e. being separated through partitions or being in near or large distances from each other, can be subjected to several parameters in which the effective ventilation may contribute as a most important parameter. To investigate this issue, the flow has been simulated using three different configurations and it has been found that being apart through partitions can be the best arrangement, which would be hard to identify before the numerical simulation. Upon considering the results of sections 7-2 and 7-3, in order to have a better cooling one has to use more powerful fans and detach the internal walls from each other but has also to account for reverse effects of the more energy consumption and radiation hazards respectively. Fig 13 - Temperature contours for different Prandtls = 0.1 = 0.7 To improve the reliability of the above analysis the following works are recommended: 1- Simulating the flow using the three dimensional geometry. 2- Considering of ventilation from an overhead location in thee dimensional model 3- Considering the effects of variable Prandtl number in solution domain. References [1] Rhie C.M. and Chow W.L., (1983), Numerical Study of Turbulent Flow Past an Airfoil with Trailing Edge Separation, AIAA Journal, 21(11): [2] Versteeg, H., K., Malalasekera, W., H.G., (1995), An Introduction to Computational Fluid Dynamics, Longman Group. [3] Ghia, U., Ghia, K., N., Shin, C., T., (1982), High- Reynolds Solutions for Incompressible Flow Using the Navier-Stokes Equations and a Multigrid Method, J., Comput. Phys., 48(3):
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
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 informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
More informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
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 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 informationA numerical investigation of tip clearance flow in Kaplan water turbines
Published in the proceedings of HYDROPOWER INTO THE NEXT CENTURY - III, 1999. ISBN 9522642 9 A numerical investigation of tip clearance flow in Kaplan water turbines M.Sc. H. Nilsson Chalmers University
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationProject #1 Internal flow with thermal convection
Project #1 Internal flow with thermal convection MAE 494/598, Fall 2017, Project 1 (20 points) Hard copy of report is due at the start of class on the due date. The rules on collaboration will be released
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 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 informationEffect of Periodic Variation of Sol-air Temperature on the Performance of Integrated Solar Collector Storage System
Engineering, 2010, 2, 832-840 doi:10.4236/eng.2010.210106 Published Online October 2010 (http://www.scirp.org/journal/eng) Effect of Periodic Variation of Sol-air Temperature on the Performance of Integrated
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 informationComputation for the Backward Facing Step Test Case with an Open Source Code
Computation for the Backward Facing Step Test Case with an Open Source Code G.B. Deng Equipe de Modélisation Numérique Laboratoire de Mécanique des Fluides Ecole Centrale de Nantes 1 Rue de la Noë, 44321
More informationNumerical studies on natural ventilation flow in an enclosure with both buoyancy and wind effects
Numerical studies on natural ventilation flow in an enclosure with both buoyancy and wind effects Ji, Y Title Authors Type URL Numerical studies on natural ventilation flow in an enclosure with both buoyancy
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 informationDepartment of Energy Science & Engineering, IIT Bombay, Mumbai, India. *Corresponding author: Tel: ,
ICAER 2011 AN EXPERIMENTAL AND COMPUTATIONAL INVESTIGATION OF HEAT LOSSES FROM THE CAVITY RECEIVER USED IN LINEAR FRESNEL REFLECTOR SOLAR THERMAL SYSTEM Sudhansu S. Sahoo* a, Shinu M. Varghese b, Ashwin
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
More informationSTUDY OF A PASSIVE SOLAR WINTER HEATING SYSTEM BASED ON TROMBE WALL
STUDY OF A PASSIVE SOLAR WINTER HEATING SYSTEM BASED ON TROMBE WALL Dr. G.S.V.L.Narasimham Chief Research Scientist, RAC, Dept. of Mechanical Engineering, Indian Institute of Science,Bengaluru- 560012,
More informationEnhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins
Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins Jubin Jose 1, Reji Mathew 2 1Student, Dept. of Mechanical Engineering, M A
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 informationEFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY
THERMAL SCIENCE: Year 2018, Vol. 22, No. 6A, pp. 2413-2424 2413 EFFECT OF THE INLET OPENING ON MIXED CONVECTION INSIDE A 3-D VENTILATED CAVITY by Hicham DOGHMI *, Btissam ABOURIDA, Lahoucin BELARCHE, Mohamed
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationTransient Thermal Flow and Thermal Stress Analysis Coupled NASTRAN and SC/Tetra
Transient Thermal Flow and Thermal Stress Analysis Coupled NASTRAN and SC/Tetra Qin Yin Fan Software CRADLE Co., Ltd. ABSTRACT In SAE paper 2004-01-1345, author focused on how to use a steady state temperature
More informationInternational Journal of Engineering Research and General Science Volume 3, Issue 6, November-December, 2015 ISSN
NUMERICAL AND EXPERIMENTAL INVESTIGATION OF STAGGERED INTERRUPTED FIN ARRANGEMENT IN A NATURAL CONVECTION FIELD Mr.Bhushan S Rane 1, Prof. M D Shende 2 1 (P G Student, Department of Mechanical Engineering,
More informationEntropic Evaluation of Dean Flow Micromixers
COMSOL Conference, Boston, 2013 Brian Vyhnalek, Petru S. Fodor and Miron Kaufman Physics Department Cleveland State University Entropic Evaluation of Dean Flow Micromixers ABSTRACT We study the use of
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 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 informationOn Clean Cooling Systems for Wind Turbine Nacelle operating in Hot Climate
International Conférence en Clean Cooling Technologiesin the ME NA Regions ICT3_MENA 201 Bou Smail, W. Tipaza, 5-6 October 2015 On Clean Cooling Systems for Wind Turbine Nacelle operating in Hot Climate
More informationCOMPUTATIONAL FLOW ANALYSIS THROUGH A DOUBLE-SUCTION IMPELLER OF A CENTRIFUGAL PUMP
Proceedings of the Fortieth National Conference on Fluid Mechanics and Fluid Power December 12-14, 2013, NIT Hamirpur, Himachal Pradesh, India FMFP2013_141 COMPUTATIONAL FLOW ANALYSIS THROUGH A DOUBLE-SUCTION
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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationNumerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions
Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén
More informationApplied CFD Project 1. Christopher Light MAE 598
Applied CFD Project 1 Christopher Light MAE 598 October 5, 2017 Task 1 The hot water tank shown in Fig 1 is used for analysis of cool water flow with the heat from a hot plate at the bottom. For all tasks,
More informationCFD ANALYSIS OF TRIANGULAR ABSORBER TUBE OF A SOLAR FLAT PLATE COLLECTOR
Int. J. Mech. Eng. & Rob. Res. 2013 Basavanna S and K S Shashishekar, 2013 Research Paper ISSN 2278 0149 www.imerr.com Vol. 2, No. 1, January 2013 2013 IJMERR. All Rights Reserved CFD ANALYSIS OF TRIANGULAR
More informationEnhancement of the momentum interpolation method on non-staggered grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 33: 1 22 Enhancement of the momentum interpolation method on non-staggered grids J. Papageorgakopoulos, G. Arampatzis,
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 informationTHERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
More information3D Numerical Simulation of Supercritical Flow in Bends of Channel
3D Numerical Simulation of Supercritical Flow in Bends of Channel Masoud. Montazeri-Namin, Reyhaneh-Sadat. Ghazanfari-Hashemi, and Mahnaz. Ghaeini- Hessaroeyeh Abstract An attempt has been made to simulate
More informationCFD as a Tool for Thermal Comfort Assessment
CFD as a Tool for Thermal Comfort Assessment Dimitrios Koubogiannis dkoubog@teiath.gr G. Tsimperoudis, E. Karvelas Department of Energy Technology Engineering Technological Educational Institute of Athens
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 informationPrediction of Performance Characteristics of Orifice Plate Assembly for Non-Standard Conditions Using CFD
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869, Volume-3, Issue-5, May 2015 Prediction of Performance Characteristics of Orifice Plate Assembly for Non-Standard Conditions
More informationExternal and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationOpen PRAIRIE: Open Public Research Access Institutional Repository and Information Exchange
South Dakota State University Open PRAIRIE: Open Public Research Access Institutional Repository and Information Exchange Theses and Dissertations 2016 The Effects of Diffuser Exit Velocity and Distance
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 informationLattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocks
Lattice Boltzmann Modeling of Wave Propagation and Reflection in the Presence of Walls and Blocs Maysam Saidi, Hassan Basirat Tabrizi, Member, IAENG, and Reza Sepahi Samian Abstract Lattice Boltzmann method
More informationNumerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows
Numerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows T. Do L. Chen J. Tu B. Anderson 7 November 2005 Abstract Flow-induced noise from fully submerged lifting bodies
More informationIJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online):
IJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online): 2321-0613 Performance of Rectangular Baffle Plate Shell and Tube Heat Exchanger using Computational
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 informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationMathematical Modelling of Ceramic Block Heat Transfer Properties
Proceedings of the 3 RD INTERNATIONAL CONFERENCE ADVANCED CONSTRUCTION 18-19 October, 2012, Kaunas, Lithuania Kaunas University of Technology, Faculty of Civil Engineering and Architecture Studentu st.
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 informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationNumerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
More informationNumerical study of the steady state uniform flow past a rotating cylinder
Numerical study of the steady state uniform flow past a rotating cylinder J. C. Padrino and D. D. Joseph December 17, 24 1 Introduction A rapidly rotating circular cylinder immersed in a free stream generates
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 informationComparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
Proceedings of 4 th ICCHMT May 17-0, 005, Paris-Cachan, FRANCE 381 Comparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
More informationNumerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder
Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationA CFD Analysis Of A Solar Air Heater Having Triangular Rib Roughness On The Absorber Plate
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.2, pp 964-971, April-June 2013 ICGSEE-2013[14th 16th March 2013] International Conference on Global Scenario in
More informationVertical Mantle Heat Exchangers for Solar Water Heaters
for Solar Water Heaters Y.C., G.L. Morrison and M. Behnia School of Mechanical and Manufacturing Engineering The University of New South Wales Sydney 2052 AUSTRALIA E-mail: yens@student.unsw.edu.au Abstract
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
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 informationAcoustic Streaming Driven Mixing
Acoustic Streaming Driven Mixing Nitesh Nama, Po-Hsun Huang, Francesco Costanzo, and Tony Jun Huang Department of Engineering Science and Mechanics The Pennsylvania State University, State College, PA,
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 informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationThe Natural Convective Heat Transfer in Rectangular Enclosure Containing Two Inclined Partitions
AASCIT Journal of Energy 2015; 2(3): 36-43 Published online June 10, 2015 (http://www.aascit.org/journal/energy) The Natural Convective Heat Transfer in Rectangular Enclosure Containing Two Inclined Partitions
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationTurbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation Source
International Scientific Colloquium Modelling for Material Processing Riga, September 16-17, 2010 Turbulence Model Affect on Heat Exchange Characteristics Through the Beam Window for European Spallation
More information2-D CFD analysis of passenger compartment for thermal comfort and ventilation
2-D CFD analysis of passenger compartment for thermal comfort and ventilation Ajay Giri 1*, Brajesh Tripathi 2 and H.C. Thakur 3 *1, 3 School of Engineering, Gautam Buddha University, Greater Noida, U.P,
More informationEFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION
EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION A. K. Kansal, P. Suryanarayana, N. K. Maheshwari Reactor Engineering Division, Bhabha Atomic Research Centre,
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 informationGeneral Solution of the Incompressible, Potential Flow Equations
CHAPTER 3 General Solution of the Incompressible, Potential Flow Equations Developing the basic methodology for obtaining the elementary solutions to potential flow problem. Linear nature of the potential
More informationCalculations on a heated cylinder case
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
More informationSTATOR/ROTOR INTERACTION
TASK QUARTERLY 10 No 2, 113 124 CFD MODELLING OF TURBINE STAGE STATOR/ROTOR INTERACTION JERZY ŚWIRYDCZUK Institute of Fluid-Flow Machinery, Polish Academy of Sciences, Fiszera 14, 80-952 Gdansk, Poland
More informationThe effect of Entry Region on Thermal Field
The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving
More informationPVP COMPUTATIONAL FLUID DYNAMICS (CFD) PROVIDES ALTERNATIVE TO CONVENTIONAL HEAT TREATMENT GUIDELINES
Proceedings of the ASME 2017 Pressure Vessels and Piping Conference PVP2017 July 16-20, 2017, Waikoloa, Hawaii, United States PVP2017-65856 COMPUTATIONAL FLUID DYNAMICS (CFD) PROVIDES ALTERNATIVE TO CONVENTIONAL
More informationTHERMAL ANALYSIS OF AN INDUCTION MOTOR BY HYBRID MODELING OF A THERMAL EQUIVALENT CIRCUIT AND CFD
THERMAL ANALYSIS OF AN INDUCTION MOTOR BY HYBRID MODELING OF A THERMAL EQUIVALENT CIRCUIT AND CFD Thiago Voigdlener, voigdlener@gmail.com Daniel Jochelavicius, danieljoche@gmail.com Federal University
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 informationDesign of Test Section for Coriolis Rig and Numerical Simulation of Cooling Flow Through Test Section
University Turbine Systems Research (UTSR) Industrial Fellowship Program 2012 Design of Test Section for Coriolis Rig and Numerical Simulation of Cooling Flow Through Test Section Prepared For: Solar Turbines,
More informationTwo-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet
PHYSICS OF FLUIDS 18, 085103 2006 Two-dimensional model problem to explain counter-rotating vortex pair formation in a transverse jet Suman Muppidi and Krishnan Mahesh Aerospace Engineering & Mechanics,
More informationNumerical Investigation of Conjugate Natural Convection Heat Transfer from Discrete Heat Sources in Rectangular Enclosure
Proceedings of the World Congress on Engineering 4 Vol II, WCE 4, July - 4, 4, London, U.K. Numerical Investigation of Conjugate Natural Convection Heat Transfer from Discrete Heat Sources in Rectangular
More informationCFD Analysis of High Temperature and High Velocity System: Plasma Torch
CFD Analysis of High Temperature and High Velocity System: Plasma Torch Abhishek Pratap Singh Bhadauria 1 1 Department of Mechanical Engineering, K. J. Somaiya College of Engineering, Mumbai, Maharashtra,
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
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 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 informationANALYSIS OF LAMINAR FORCED CONVECTION INSIDE A SQUARE VENTILATED CAVITY USING THE OPENFOAM
ANALYSIS OF LAMINAR FORCED CONVECTION INSIDE A SQUARE VENTILATED CAVITY USING THE OPENFOAM P. C. Mioralli a, V. L. Scalon b, E. Avallone a, and S. A. Verdério Júnior c a,c Instituto Federal de São Paulo
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationMIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER
Suranaree J. Sci. Technol. Vol. 17 No. 2; April - June 2010 139 MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER Md. Mustafizur Rahman 1 *, M. A. Alim 1 and Sumon Saha
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 information