EULER AND SECOND-ORDER RUNGE-KUTTA METHODS FOR COMPUTATION OF FLOW AROUND A CYLINDER
|
|
- Alison Greene
- 6 years ago
- Views:
Transcription
1
2 EULER AND SEOND-ORDER RUNGE-KUTTA METHODS FOR OMPUTATION OF FLOW AROUND A YLINDER László Daróczy, László Baranyi MSc student, Professor University of Miskolc, Hungary Department of Fluid and Heat Engineering, University of Miskolc, Hungary ABSTRAT The two-dimensional low-reynolds number flow around a circular cylinder is investigated numerically using two different time discretisation methods ( st order Euler and nd order Runge-Kutta) and two different types of hardware (PU- and GPGPU-based) using the finite difference method. omputations were done at Reynolds number and 5 for different dimensionless time steps (.;.;.4;.5) and different mesh sizes (545; 366). For all results flow properties such as drag, lift and base pressure coefficients are analysed and compared with each other. omputational results obtained using the four different types of code agree well, and are also very similar to values in the literature. INTRODUTION When computing flow around a bluff body such as a cylinder a proper balance needs to be found between computational time and accuracy. Higher order computational methods are more accurate than lower order methods when the temporal and spatial disretisations are the same. If a higher order method is chosen then it is possible to obtain the same accuracy as for a lower order method even at a larger time step and or with a coarser mesh. The Euler ( st order) method, although considered old-fashioned by many, is still in use, even for problems involving highly unsteady phenomenon such as wake flows of bluff bodies. Among the nd order methods the most frequently used are the Runge- Kutta and the Adams-Bashforth methods []; several versions of each eist. There are also 4 th order Runge Kutta methods but they are not so frequently used in the study of wake flow. For this sudy we tried the four most often recommended versions of nd order Runge-Kutta (RK) method: the Heun, Midpoint, Ralston and Taylor series methods [, 3]. st and nd order results are compared at different time steps, grid sizes, computational domains and hardware. The D computations are for a low-reynolds number flow around a stationary circular cylinder placed in a uniform stream. EULER AND RUNGE-KUTTA METHODS We briefly introduce the computational procedure developed by the second author [4, 5] incorporating the Euler method. A non-inertial system fied to the accelerating cylinder is used to compute two-dimensional (D) low Reynolds
3 number flow around a circular cylinder placed perpendicular to an otherwise uniform flow. The governing equations are the Navier-Stokes equations, the equations of continuity and a Poisson equation for pressure. On the cylinder surface no-slip boundary condition is used for the velocity and a Neumann type boundary condition for the pressure. At the far region potential flow is assumed. Boundaryfitted coordinates are used to impose boundary conditions accurately. The physical domain, which consists of two concentric circles (with radii R, R ), is mapped into a rectangular computational domain. The governing equations and boundary conditions are also transformed to the computational domain (not shown here) and solved by finite difference method. Space derivatives are approimated by fourthorder central differences, ecept for the convective terms, for which a third-order modified upwind scheme is used. The Poisson equation for pressure is solved by the successive over-relaation (SOR) method. The Navier-Stokes equations are integrated using the Euler method to yield the velocity field. The continuity equation is satisfied at every time step. For further details see [4] and [5]. We tested all four Runge-Kutta methods mentioned, but finally Heun's method was chosen for implementation as a compromise between speed and accuracy. The dimensionless Navier-Stokes equation for a non-inertial system fied to the accelerating cylinder is v = ( v ) v p+ v a, () t Re where v is the velocity vector, p is the pressure, Re = Ud / ν is the Reynolds number based on the free stream velocity U, cylinder diameter d and kinematic viscosity of the fluid ν, a is the cylinder acceleration vector ( here), and t is the dimensionless time. The pressure Poisson equation can be written as [4] D u v u v = p t y y. () With superscript n denoting a time level and t a time step, the sequence of temporal integration of equation () using Heun s RK method is as follows. Step. Determine RK coefficient k and intermediate velocity v and pressure p n n n n k= ( v ) v + p v Re. (3) Using k and equation () the intermediate velocity v can be obtained v = v + k t v + v, (4) n n n + where t is the time step and v is the velocity of the cylinder. Now we can determine an intermediate pressure p using v and equation ()
4 n D u v u v = p t y y, (5) n n where D = div( u ). In this way while solving the Poisson equation (5) we also enforce the continuity condition of D = div( v ) =. n+ n Step. Determine RK coefficient k and final velocity u and pressure p + k = ( v ) v + p v Re, (6) n n n n v + = v + ( k+ k ) t v + + v, (7) n n+ n+ n+ n+ D u v u v n+ = p t y y. (8) As can be seen from this sequence, the SOR method, which requires the large bulk of the computational time, has to be applied twice in a single time step when using the Runge-Kutta method, while with the Euler method it is applied only once. However, tests showed that in most cases the computational time does not double. DIFFERENT TYPES OF ODES In the current study FlowFD software, developed by the first author, provided the framework for Euler and Runge-Kutta methods, running on either PU or GPU. The software package includes many post-processing options, including visualization of flow and limit cycles in D and 3D, analysis of different coefficients, and management for scheduling simulations. GPU codes use a GPGPU (General-purpose computing on graphics processing units) programming method. The code itself is ++ based, while the computation kernels are written in NVIDIA UDA. Instead of the original SOR method [4], the GPU version uses a red-black Gauss-Seidel iteration with successive overrelaation (SOR) [6]. The results presented in the present article were all run on Tesla 5 using 448 cores with double precision and E on. The study uses four different variations of FD codes. Types and are the recreation in ++ of the original FORTRAN code developed by the second author [4], using st order Euler-type time discretisation; Type uses the PU, and Type runs on the GPU. Types 3 and 4 are based on the same code, slightly modified for the Runge-Kutta method; Type 3 runs on the PU and Type 4 on the GPU. TESTED ASES & EVALUATION METHODS and 5 cases were tested for a stationary cylinder. For the GPU codes a 545 (peripheral radial) grid was used at both Reynolds numbers, and for PU codes two grids were used: for 5 a 545 grid with R / R = 6, while the same was used for, in addition to a smaller grid of 366 ( R / R = 6).
5 For all cases 4 4 t = was used with ( ) t=,,4,5 non-dimensional time step. 4 For Reynolds number time step t 6 failed to converge. As Fast Fourier Transform (FFT) analysis revealed only one dominant frequency in the spectrum, for evaluation of the results the built-in Maut algorithm was used, which calculates the period time of the oscillations using the local maima of the coefficients. The vorticity contours were the same for all cases belonging to the same Reynolds number (see Figure ). RESULTS Figure Vorticity contours () The results are presented in Tables and for the different settings. Drag coefficient ( D, mean), lift coefficients (amplitude of L ) and base pressure coefficients (, mean) are given. The results are compared to the values pb presented in the work of Posdziech and Grundmann, who used a systematic approach to determine force coefficients and Strouhal number of two-dimensional flow around a circular cylinder using the spectral element method [7]. We can see that the results correlate well: for eample, values are presented in [7] as.33 5 Hardware D Table Results - PU Time step Euler (st order) Runge-Kutta (nd order) Δt (dimensionless) PU D pb L D pb L D pb L
6 (for their largest domain size, at ), which is similar to the range of values presented in Tables and ( for both PU and GPU). They also gave a value of.99 for 5 (around.3 for us), and values are D presented as.69 (for the largest domain size, ) and.87 for 5. Hardware Table Results - GPU GPU Time step Euler (st order) Runge-Kutta (nd order) Δt (dimensionless) D pb L D pb L The non-dimensional vorte shedding frequency St = fd / U is also compared with the results in [7]. In case of, for comparison we will use the value of.633 ([7], Table 4, largest domain). The values were in the range of , with relative errors in the range of.576%.636%. The Runge- Kutta codes resulted in smaller relative error. However, for the smaller grid (366) we obtained results in the range of , corresponding to the relative error of.69%.9%, slightly higher than in the case of the larger grid. For 5, [7] suggests several formulas to evaluate St for different Reynolds numbers. Using the non-integer polynomial (lowest error) we get omparing these values with the tested cases, we see that results are between.8356 and.8365, which is a very low error (between.6% and.58%). For L amplitudes the presented values in [7] are.36 (, largest domain size) and.53 (5, largest domain). The amplitude of the lift coefficient compared well with the values in [7]. In the case of, for the smaller mesh values were in the range of.35.34, while for the larger mesh they were For 5, our simulations gave values in the range of ONLUSIONS The second order results agree very well with the results of the first order temporal discretisation for the time steps investigated. This indicates that the first order results are accurate enough for practical cases under the conditions investigated. Both sets of results agree well with the spectral element results given in [7]. In addition, performing the computations using PU and GPGPU revealed no pb
7 significant differences in the results, although the computational time required is roughly 8- times lower for the GPU. While we were able to implement the Runge-Kutta method into our eisting and proven code we found that we were unable to take advantage of the benefits of higher order methods, i.e., being able to use larger time steps and still obtain the same accuracy. When the time step was increased the SOR method failed to converge. It appears that the SOR method and the Runge-Kutta method are somewhat incompatible. We are currently investigating ways to replace SOR. The PU Euler code is the identical to the original FORTRAN version, which has been used many times in publications and articles, and whose results correlate well with the literature [5]. From these results it can be clearly concluded that the GPGPU and PU-based code yield identical results up to the defined precision 5 ( ). omparing the results of GPU codes we can see the same tendency. AKNOWLEDGEMENTS The authors would like to say special thanks to NVIDIA for donating an NVIDIA Tesla 5 within the Academic Partnership Program. The support provided by the Hungarian Research Foundation (OTKA project K 768) is gratefully acknowledged. The work was carried out as part of the TÁMOP-4...B- //KONV-- project in the framework of the New Hungarian Development Plan. The realization of this project is supported by the European Union, co-financed by the European Social Fund. REFERENES [] FERZIGER, J.H., PERIĆ, M.: omputational Methods for Fluid Dynamics (3rd edition). Springer-Verlag, Berlin (). [] HEEVER, E.: The nd order Runge-Kutta method. [3] KAW, A.K., KALU, E.E.: Numerical Methods with Applications. Ed. World Wide Web, 9. [4] BARANYI, L.: omputation of unsteady momentum and heat transfer from a fied circular cylinder in laminar flow. Journal of omputational and Applied Mechanics 4() (3), 3-5. [5] BARANYI, L.: Numerical simulation of flow around an orbiting cylinder at different ellipticity values. Journal of Fluids and Structures 4 (8), [6] SHINN, A.F.: omputational Fluid Dynamics (FD) using Graphics Processing Units. Accelerators for Science and Engineering Applications: GPUs and Multicores [7] POSDZIEH, O., GRUNDMANN, R.: A systematic approach to the numerical calculation of fundamental quantities of the two-dimensional flow over a circular cylinder. Journal of Fluids and Structures 3 (7),
Numerical Simulation of Mechanical Energy Transfer between Fluid and a Circular Cylinder Forced to Follow an Elliptical Path
OS16-3-3 The 12th International Symposium on Fluid Control, Measurement and Visualization November183, 213, Nara, Japan FLUCOM213 Numerical Simulation of Mechanical nergy Transfer between Fluid and a Circular
More informationffiloo [g PROCEEDINGS OF THE SEVENTH CONFERENCE ON MECHANICAL ENGINEERING Made by IIIVTCD .liiiiillilillllillllliillillillllim' I
PROCEEDINGS OF THE SEVENTH CONFERENCE ON MECHANICAL ENGINEERING May 25-26. 2010 Budapest. Hungary ) ffiloo [g Made by IIIVTCD ISBN 978-963-31 3-007-0 @: t o t '.liiiiillilillllillllliillillillllim' I Budapest
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
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 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 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 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 informationNumerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement
Numerical Investigation of Vortex Induced Vibration of Two Cylinders in Side by Side Arrangement Sourav Kumar Kar a, 1,, Harshit Mishra a, 2, Rishitosh Ranjan b, 3 Undergraduate Student a, Assitant Proffessor
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 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 informationNUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER ON A CARTESIAN GRID
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2000 Barcelona, 11-14 September 2000 c ECCOMAS NUMERICAL COMPUTATION OF A VISCOUS FLOW AROUND A CIRCULAR CYLINDER
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 informationVORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS
THERMAL SCIENCE, Year 2012, Vol. 16, No. 5, pp. 1395-1399 1395 VORTEX SHEDDING PATTERNS IN FLOW PAST INLINE OSCILLATING ELLIPTICAL CYLINDERS by Li-Zhong HUANG a* and De-Ming NIE b a State Key Laboratory
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationPredicting vortex-induced vibration from driven oscillation results
Applied Mathematical Modelling 3 (26) 196 112 www.elsevier.com/locate/apm Predicting vortex-induced vibration from driven oscillation results J.S. Leontini *, B.E. Stewart, M.C. Thompson, K. Hourigan Department
More informationAn example of the Rvachev function method
arxiv:1603.00320v1 [physics.flu-dyn] 1 Mar 2016 An example of the Rvachev function method Alexander V. Proskurin Altai State University, Altai State Technical University, k210@list.ru Anatoly M. Sagalakov
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 informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 1, No 3, 2010
CFD analysis of 2D unsteady flow around a square cylinder Gera.B, Pavan K. Sharma, Singh R.K Reactor Safety Division, Bhabha Atomic Research Centre, Trombay, Mumbai, India 400085 pa1.sharma@gmail.com ABSTRACT
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 informationNumerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow
Fluid Structure Interaction VII 147 Numerical investigation on vortex-induced motion of a pivoted cylindrical body in uniform flow H. G. Sung 1, H. Baek 2, S. Hong 1 & J.-S. Choi 1 1 Maritime and Ocean
More informationStatistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders
Statistical Analysis of the Effect of Small Fluctuations on Final Modes Found in Flows between Rotating Cylinders Toshiki Morita 1, Takashi Watanabe 2 and Yorinobu Toya 3 1. Graduate School of Information
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20
2.29 Numerical Fluid Mechanics Fall 2011 Lecture 20 REVIEW Lecture 19: Finite Volume Methods Review: Basic elements of a FV scheme and steps to step-up a FV scheme One Dimensional examples d x j x j 1/2
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 informationApplication of a Virtual-Boundary Method for the Numerical Study of Oscillations Developing Behind a Cylinder Near A Plane Wall
Fluid Dynamics, Vol. 39, No. 1, 2004, pp. 61 68. Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2004, pp. 69 77. Original Russian Text Copyright 2004 by Kit, Nikitin,
More informationVortex Induced Vibrations
Vortex Induced Vibrations By: Abhiroop Jayanthi Indian Institute of Technology, Delhi Some Questions! What is VIV? What are the details of a steady approach flow past a stationary cylinder? How and why
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More informationModule 3: Velocity Measurement Lecture 15: Processing velocity vectors. The Lecture Contains: Data Analysis from Velocity Vectors
The Lecture Contains: Data Analysis from Velocity Vectors Velocity Differentials Vorticity and Circulation RMS Velocity Drag Coefficient Streamlines Turbulent Kinetic Energy Budget file:///g /optical_measurement/lecture15/15_1.htm[5/7/2012
More informationSHEAR LAYER REATTACHMENT ON A SQUARE CYLINDER WITH INCIDENCE ANGLE VARIATION
Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 9- December 9 SHEAR LAYER REATTACHMENT ON A SQUARE CYLINDER WITH INCIDENCE ANGLE VARIATION Priyanka
More informationDay 24: Flow around objects
Day 24: Flow around objects case 1) fluid flowing around a fixed object (e.g. bridge pier) case 2) object travelling within a fluid (cars, ships planes) two forces are exerted between the fluid and the
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationON PARTITIONED AND MONOLITHIC COUPLING STRATEGIES IN LAGRANGIAN VORTEX METHODS FOR 2D FSI PROBLEMS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK ON PARTITIONED AND MONOLITHIC COUPLING STRATEGIES
More informationVortex structures in the wake of a buoyant tethered cylinder at moderate to high reduced velocities
European Journal of Mechanics B/Fluids 23 (2004) 127 135 Vortex structures in the wake of a buoyant tethered cylinder at moderate to high reduced velocities K. Ryan, M.C. Thompson, K. Hourigan Fluids Laboratory
More informationA STRONG COUPLING SCHEME FOR FLUID-STRUCTURE INTERACTION PROBLEMS IN VISCOUS INCOMPRESSIBLE FLOWS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2005 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2005 A STRONG COUPLING SCHEME
More informationSimulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls
Fluid Structure Interaction V 85 Simulation analysis using CFD on vibration behaviors of circular cylinders subjected to free jets through narrow gaps in the vicinity of walls K. Fujita Osaka City University,
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 information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 4
2.29 Spring 2015 Lecture 4 Review Lecture 3 Truncation Errors, Taylor Series and Error Analysis Taylor series: 2 3 n n i1 i i i i i n f( ) f( ) f '( ) f ''( ) f '''( )... f ( ) R 2! 3! n! n1 ( n1) Rn f
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 informationINTRODUCTION OBJECTIVES
INTRODUCTION The transport of particles in laminar and turbulent flows has numerous applications in engineering, biological and environmental systems. The deposition of aerosol particles in channels and
More informationSuppression of the unsteady vortex wakes of a circular cylinder pair by a doublet-like counter-rotation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2009) Published online in Wiley InterScience (www.interscience.wiley.com)..2075 Suppression of the unsteady vortex wakes
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 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 informationA. Bottaro, D. Venkataraman & F. Negrello Università di Genova, Italy
A. Bottaro, D. Venkataraman & F. Negrello Università di Genova, Italy A. Bottaro, D. Venkataraman & F. Negrello & G. Tadmor Università di Genova, Italy Focus on passive actuators, what works, why it does,
More informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
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 informationExperimental and Numerical Investigation of Flow over a Cylinder at Reynolds Number 10 5
Journal of Modern Science and Technology Vol. 1. No. 1. May 2013 Issue. Pp.52-60 Experimental and Numerical Investigation of Flow over a Cylinder at Reynolds Number 10 5 Toukir Islam and S.M. Rakibul Hassan
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 informationNavier-Stokes equations
1 Navier-Stokes equations Introduction to spectral methods for the CSC Lunchbytes Seminar Series. Incompressible, hydrodynamic turbulence is described completely by the Navier-Stokes equations where t
More informationCFD DESIGN OF A GENERIC CONTROLLER FOR VORTEX-INDUCED RESONANCE
Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 9-11 December 2009 CFD DESIGN OF A GENERIC CONTROLLER FOR VORTEX-INDUCED RESONANCE Andrew A. ANTIOHOS,
More informationLength Learning Objectives Learning Objectives Assessment
Universidade Federal Fluminense PGMEC Course: Advanced Computational Fluid Dynamics Coordinator: Vassilis Theofilis Academic Year: 2018, 2 nd Semester Length: 60hrs (48hrs classroom and 12hrs tutorials)
More informationUNIT IV BOUNDARY LAYER AND FLOW THROUGH PIPES Definition of boundary layer Thickness and classification Displacement and momentum thickness Development of laminar and turbulent flows in circular pipes
More informationOPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK
OPEN CHANNEL FLOW Numerical Methods and Computer Applications Roland Jeppson TECHNISCHE INFORMATlONSBiBUOTHEK UNIVERSITATSB'BUOTHEK HANNOVER Si. i. CRC Press Taylor &.Francis Group Boca Raton London New
More informationCFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2
CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration G.S.T.A. Bangga 1*, W.A. Widodo 2 1,2 Department of mechanical engineering Field of study energy conversion Institut
More informationExperimental Aerodynamics. Experimental Aerodynamics
Lecture 3: Vortex shedding and buffeting G. Dimitriadis Buffeting! All structures exposed to a wind have the tendency to vibrate.! These vibrations are normally of small amplitude and have stochastic character!
More informationDYNAMIC ANALYSIS OF WIND EFFECTS BY USING AN ARTIFICAL WIND FUNCTION
2008/3 PAGES 21 33 RECEIVED 15. 3. 2008 ACCEPTED 10. 7. 2008 J. GYÖRGYI, G. SZABÓ DYNAMIC ANALYSIS OF WIND EFFECTS BY USING AN ARTIFICAL WIND FUNCTION J. Györgyi, Prof Budapest University of Technology
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationNumerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation
, pp.49-58 http://dx.doi.org/10.1457/ijast.016.9.06 Numerical Investigation of the Fluid Flow around and Past a Circular Cylinder by Ansys Simulation Mojtaba Daneshi Department of Mechanical Engineering,
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
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 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 informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationStudy on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration Time
World Journal of Mechanics, 2018, 8, 301-310 http://www.scirp.org/journal/wjm ISSN Online: 2160-0503 ISSN Print: 2160-049X Study on Non-Uniqueness of Taylor Vortex Flow Changing Inner Cylinder Acceleration
More informationDNS of the Taylor-Green vortex at Re=1600
DNS of the Taylor-Green vortex at Re=1600 Koen Hillewaert, Cenaero Corentin Carton de Wiart, NASA Ames koen.hillewaert@cenaero.be, corentin.carton@cenaero.be Introduction This problem is aimed at testing
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
More informationMA3232 Numerical Analysis Week 9. James Cooley (1926-)
MA umerical Analysis Week 9 James Cooley (96-) James Cooley is an American mathematician. His most significant contribution to the world of mathematics and digital signal processing is the Fast Fourier
More informationSuppression of 3D flow instabilities in tightly packed tube bundles
Suppression of 3D flow instabilities in tightly packed tube bundles Nicholas Kevlahan kevlahan@mcmaster.ca Department of Mathematics & Statistics CSFD, June 13 15 2004 p.1/33 Collaborators CSFD, June 13
More informationNumerical simulations of the edge tone
Numerical simulations of the edge tone I. Vaik, G. Paál Department of Hydrodynamic Systems, Budapest University of Technology and Economics, P.O. Box 91., 1521 Budapest, Hungary, {vaik, paal}@vizgep.bme.hu
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationVortex-induced vibration of a slender single-span cylinder
Vortex-induced vibration of a slender single-span cylinder N. Oikou Delft University of Technology, the Netherlands The goal of this paper is to study the vortex-induced vibration of slender cylindrical
More informationChapter 4. MAC Scheme. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes. u = 0
Chapter 4. MAC Scheme 4.1. MAC Scheme and the Staggered Grid. The MAC scheme is a numerical method used to solve the incompressible Navier-Stokes equation in the velocity-pressure formulation: (4.1.1)
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationINVESTIGATION OF THE FLOW OVER AN OSCILLATING CYLINDER WITH THE VERY LARGE EDDY SIMULATION MODEL
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationDedalus: An Open-Source Spectral Magnetohydrodynamics Code
Dedalus: An Open-Source Spectral Magnetohydrodynamics Code Keaton Burns University of California Berkeley Jeff Oishi SLAC National Accelerator Laboratory Received ; accepted ABSTRACT We developed the magnetohydrodynamic
More information1. Fluid Dynamics Around Airfoils
1. Fluid Dynamics Around Airfoils Two-dimensional flow around a streamlined shape Foces on an airfoil Distribution of pressue coefficient over an airfoil The variation of the lift coefficient with the
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More informationSimulation of Cross Flow Induced Vibration
Simulation of Cross Flow Induced Vibration Eric Williams, P.Eng Graduate Student, University of New Brunswic, Canada Andrew Gerber, PhD, P.Eng Associate Professor, University of New Brunswic, Canada Marwan
More informationL.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133, ORSAY CEDEX, FRANCE fax number:
Large Eddy Simulations of a spatially developing incompressible 3D mixing layer using the v-ω formulation. C. TENAUD, S. PELLERIN, A. DULIEU and L. TA PHUOC L.I.M.S.I. - U.P.R. C.N.R.S. 3251, B.P. 133,
More informationChalmers Publication Library
Chalmers Publication Library Quantitative numerical analysis of flow past a circular cylinder at Reynolds number between 5 and 2 This document has been downloaded from Chalmers Publication Library (CPL).
More informationNumerical Investigation of Laminar Flow over a Rotating Circular Cylinder
International Journal of Mechanical & Mechatronics Engineering IJMME-IJENS Vol:13 No:3 32 Numerical Investigation of Laminar Flow over a Rotating Circular Cylinder Ressan Faris Al-Maliky Department of
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationA numerical study on effect of corner radius and Reynolds number on fluid flow over a square cylinder
Sādhanā Vol. 42, No. 7, July 2017, pp. 1155 1165 DOI 10.1007/s12046-017-0680-2 Ó Indian Academy of Sciences A numerical study on effect of corner radius and Reynolds number on fluid flow over a square
More informationNumerical Simulation and Its Possibility of Application to Green Environment
Numerical Simulation and Its Possibility of Application to Green Environment Anna Kuwana 1,a, Tetuya Kawamura 2,b 1 Faculty of Science and Technology, Gunma University, 1-5-1 Tenjin-cho, Kiryu City, Gunma
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationLaminar flow over a steadily rotating circular cylinder under the influence of uniform shear
Laminar flow over a steadily rotating circular cylinder under the influence of uniform shear Sangmo Kang Citation: Phys. Fluids 18, 047106 (2006); doi: 10.1063/1.2189293 View online: http://dx.doi.org/10.1063/1.2189293
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationAerodynamics. Professor: Luís Eça
Professor: Luís Eça 1. Introduction Aerodynamical forces. Flow description. Dependent variables and physical principles that govern the flow 2. Incompressible, Viscous Flow Analytical solutions of the
More informationNumAn2014 Conference Proceedings
OpenAccess Proceedings of the 6th International Conference on Numerical Analysis, pp 198-03 Contents lists available at AMCL s Digital Library. NumAn014 Conference Proceedings Digital Library Triton :
More informationSide-View Mirror Vibrations Induced Aerodynamically by Separating Vortices
Open Journal of Fluid Dynamics, 2016, 6, 42-56 Published Online March 2016 in SciRes. http://www.scirp.org/journal/ojfd http://dx.doi.org/10.4236/ojfd.2016.61004 Side-View Mirror Vibrations Induced Aerodynamically
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationNUMERICAL SIMULATION OF FLUID FLOW BEHAVIOUR ON SCALE UP OF OSCILLATORY BAFFLED COLUMN
Journal of Engineering Science and Technology Vol. 7, No. 1 (2012) 119-130 School of Engineering, Taylor s University NUMERICAL SIMULATION OF FLUID FLOW BEHAVIOUR ON SCALE UP OF OSCILLATORY BAFFLED COLUMN
More information2011 Christopher William Olenek
2011 Christopher William Olenek STUDY OF REDUCED ORDER MODELS FOR VORTEX-INDUCED VIBRATION AND COMPARISON WITH CFD RESULTS BY CHRISTOPHER WILLIAM OLENEK THESIS Submitted in partial fulfillment of the requirements
More informationABSTRACT INTRODUCTION
Numerical simulation of pulsating flow around a cube C. Dargent, D. Dartus, J. George Institut de Mecanique des Fluides de Toulouse, Avenue du Professeur Camille Soula, 31400 Toulouse, France ABSTRACT
More informationSpectral Vanishing Viscosity Method for Large-Eddy Simulation of Turbulent Flows
Journal of Scientific Computing ( 2006) DOI: 10.1007/s10915-005-9029-9 Spectral Vanishing Viscosity Method for Large-Eddy Simulation of Turbulent Flows Richard Pasquetti 1 Received October 6, 2004; accepted
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationNumerical simulation of the flow behind a rotary oscillating circular cylinder
PHYSICS OF FLUIDS VOLUME 10, NUMBER 4 APRIL 1998 Numerical simulation of the flow behind a rotary oscillating circular cylinder Seung-Jin Baek and Hyung Jin Sung a) Department of Mechanical Engineering;
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationBasic Concepts: Drag. Education Community
Basic Concepts: Drag 011 Autodesk Objectives Page Introduce the drag force that acts on a body moving through a fluid Discuss the velocity and pressure distributions acting on the body Introduce the drag
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More information