Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
|
|
- Madlyn Bryant
- 5 years ago
- Views:
Transcription
1 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 this project is to develop and apply a computer program to solve the twodimensional unsteady Navier-Stokes equations for an elliptic flow in a square cavity with top wall sliding at a given velocity. We assume that the fluid is incompressible and has constant viscosity. The two-dimensional unsteady Navier-Stokes equations for flow with constant density and viscosity are listed as follows: u x + v y = 0 (1) ρ( u t + u u x + v u y ) = P x + µ( 2 u x u y 2 ) (2) ρ( v t + u v x + v v y ) = P y + µ( 2 v x v y 2 ) (3) where u and v are the two Cartesian velocities, x and y are spatial coordinates, t is the time, ρ is the density, P is the static pressure, and µ is the dynamic viscosity. The side length of the square cavity is set to 1 meter and a grid is used in the numerical computation. The grid layout is shown in Figure 1. There are ghost nodes located outside the domain of interest. The boundary conditions are that v velocity at all the domain boundaries is zero, while u velocity of the flow near the top wall equals the top wall velocity and it remains zero at other boundaries. Numerical results for different Reynolds numbers of 50, 100 and 200 will be shown. The Reynolds number is based on the side length of the square cavity, top wall velocity and the kinematic viscosity 1
2 Figure 1: Grid layout of the fluid, as is shown in Equation 4: Re = ρ V L µ (4) where V is the top wall velocity and L is the side length. The time step size t yields to the following stability limit: 1 t u x + v y + (5) 2ν + 2ν x 2 y 2 In our project, a time step of 0.002s is used in all cases which always yields above stability limit given that flow density ρ = 1kg/m 3, dynamic viscosity µ = 0.001kg/m s, x = y = 0.025m. A collocated scheme is used for velocity discretizations and pressure-poisson equations are formulated to avoid checker-board splitting. 2 Differencing Scheme 2.1 Discretization First order time and second order space discretizations are used in this project. An explicit time discretization scheme has been used for both convection and diffusion, and pressure-poisson equations are formulated for determining the pressure field at the new time step. Relevant discretizations are shown as follows: ( u x )n un i+1/2,j un i 1/2,j x (6) 2
3 ( v y )n vn i,j+1/2 un i,j 1/2 y ( u t )n un+1 i,j u n i,j t ( v t )n vn+1 i,j vi,j n t (u u x )n u n u n i+1,j un i 1,j i,j 2 x (v u y )n vi,j n u n i,j+1 un i,j 1 2 y (u v x )n u n vi+1,j n vn i 1,j i,j 2 x (v v y )n vi,j n vi,j+1 n vn i,j 1 2 y ( P x )n P i+1,j n P i 1,j n 2 x ( P y )n P i,j+1 n P i,j 1 n 2 y (7) (8) (9) (10) (11) (12) (13) (14) (15) 2.2 Collocated Fractional Step Algorithm ( 2 u x 2 )n = un i+1,j 2un i,j + un i 1,j x 2 (16) ( 2 u y 2 )n = un i,j+1 2un i,j + un i,j 1 y 2 (17) ( 2 v x 2 )n = vn i+1,j 2vn i,j + vn i 1,j x 2 (18) ( 2 v y 2 )n = vn i,j+1 2vn i,j + vn i,j 1 y 2 (19) A collocated fractional step algorithm is used in the simulation. Main procedures of this algorithm can be shown as follows: Prescribe domain size Prescribe solution and flow parameters 3
4 Initialize flow variables Start time loop The following steps will be executed inside each time loop: Solve for û, ˆv at each node (i, j) inside the domain ρûi,j u n i,j t = C n x + D n x = = ρ un i,j 2 x (un i+1,j u n i 1,j) ρ vn i,j 2 y (un i,j+1 u n i,j 1) + µ x 2 (un i+1,j 2u n i,j + u n i 1,j) + µ y 2 (un i,j+1 2u n i,j + u n i,j 1) (20) ρ ˆv i,j v n i,j t = C n y + D n y = Evaluate û, ˆv at midpoints = ρ un i,j 2 x (vn i+1,j v n i 1,j) ρ vn i,j 2 y (vn i,j+1 v n i,j 1) + µ x 2 (vn i+1,j 2v n i,j + v n i 1,j) + µ y 2 (vn i,j+1 2v n i,j + v n i,j 1) (21) û i+1/2,j 1 2 (û i,j + û i+1,j ) (22) û i 1/2,j 1 2 (û i,j + û i 1,j ) (23) ˆv i,j+1/2 1 2 (ˆv i,j + û i,j+1 ) (24) ˆv i,j 1/2 1 2 (ˆv i,j + û i,j 1 ) (25) Discrete Pressure-Poisson Equations ρ t U = ρ û i 1/2,j t (ûi+1/2,j x + ˆv i,j+1/2 ˆv i,j 1/2 ) y = ( P i+1,j n+1 n+1 2Pi,j + Pi 1,j n+1 x 2 ) + ( P i,j+1 n+1 n+1 2Pi,j + Pi,j 1 n+1 y 2 ) (26) Note that on the domain boundary since u 1.5,j and u nx+1/2,j are fixed for all j, v i,1.5 and u i,ny+1/2 are fixed for all i, small adjustment needs to be made to the original pressure- 4
5 Poisson equation. For example at (nx, j): {u nx+1/2,j [û nx 1/2,j t n+1 nx,j (P P n+1 ρ x nx 1,j )]}/ x +[ˆv nx,j+1/2 ˆv nx,j 1/2 t n+1 (P ρ nx,j+1 P nx,j n+1 y ) + t n+1 (P ρ nx,j P n+1 y nx,j 1 )]/ y = 0 (27) The discrete pressure-poisson equations will be solved by SOR or other iterative solvers. Correct grid node velocities u n+1 i,j, v n+1 i,j with û i,j, ˆv i,j and P n+1 i 1,j, P n+1 i+1,j, P n+1 i,j+1, P n+1 i,j 1 : u n+1 i,j v n+1 i,j = û i,j t ρ = ˆv i,j t ρ n+1 i+1,j (P P i 1,j n+1 2 x n+1 i,j+1 (P P i,j 1 n+1 2 y ) (28) ) (29) Modify the velocities at the ghost nodes outside the domain so that they can be used at the beginning of next time loop. (End of the time loop) In this project, we use SOR with an over-relaxation factor of 1.4 to solve the pressure-poisson equations. Convergence is assumed to be reached when error residual is less than Error residual is defined as: E = D (LHS RHS)2 nx ny where LHS and RHS are the left-hand side and the right-hand side of the pressure-poisson equation at each node and D denotes the sum of all the square of difference between LHS and RHS at each node within the domain. nx ny is the number of nodes within the domain. Similarly, we assume that steady state is reached when convergence residuals for both u and v velocities are less than 1e-6. Convergence residual is defined as: C = D (OLD NEW )2 nx ny where OLD and NEW are respectively the old value from last loop and the new value calculated in current loop and D denotes the sum of all the square of difference between OLD and NEW at each node within the domain. nx ny is the number of nodes within the domain. 5
6 3 Results and Discussion Figure 2 shows the contours of u velocity for Re=50 at 4 time intervals (time steps = 500, 1000, 2000, and 5000; time step size equals 0.002s). Figure 3 shows the contours of v velocity for Re=50 at 4 time intervals (time steps = 500, 1000, 2000, and 5000; time step size equals 0.002s). Figure 4 shows the streamlines of velocity vector for Re=50 at 4 time intervals (time steps = 500, 1000, 2000, and 5000; time step size equals 0.002s). It can be clearly observed that there is a vortex formed inside the square cavity and its center moves downward until the flow field becomes steady. The rotational direction of the vortex is clockwise, which is in accord with the moving direction of the top wall. (a) t=1s Figure 2: (a)(b)(c)(d): Contours of u velocity for Re=50 at 4 time intervals 6
7 (a) t=1s Figure 3: (a)(b)(c)(d): Contours of v velocity for Re=50 at 4 time intervals (a) t=1s Figure 4: (a)(b)(c)(d): Streamlines of velocity vector for Re=50 at 4 time intervals 7
8 Figure 5, 6 and 7 show respectively the contours of u velocity, contours of v velocity and streamlines of velocity vector for Re=100 at 4 time intervals (time steps = 500, 1000, 2000, and 5000; time step size equals 0.002s). (a) t=1s Figure 5: (a)(b)(c)(d): Contours of u velocity for Re=100 at 4 time intervals 8
9 (a) t=1s Figure 6: (a)(b)(c)(d): Contours of v velocity for Re=100 at 4 time intervals (a) t=1s Figure 7: (a)(b)(c)(d): Streamlines of velocity vector for Re=100 at 4 time intervals 9
10 Figure 8, 9 and 10 show respectively the contours of u velocity, contours of v velocity and streamlines of velocity vector for Re=200 at 4 time intervals (time steps = 500, 1000, 2000, and 5000; time step size equals 0.002s). (a) t=1s Figure 8: (a)(b)(c)(d): Contours of u velocity for Re=200 at 4 time intervals 10
11 (a) t=1s Figure 9: (a)(b)(c)(d): Contours of v velocity for Re=200 at 4 time intervals (a) t=1s Figure 10: (a)(b)(c)(d): Streamlines of velocity vector for Re=200 at 4 time intervals 11
12 Numbers of iterations it takes to reach steady state for different Reynolds numbers and the corresponding time costed are shown in Table 1. Reynolds Number Iterations Time (s) Residual of u velocity Residual of v velocity e e e e e e-7 Table 1: Important parameters when steady state is reached for different Reynolds numbers As can be observed from Table 1, for Re=50 and 100, at t=10s (5000 time steps) the flow inside the cavity is already or very close to steady. But for Re=200, it takes as long as 24.9s (12431 time steps) to reach steady state. Figure 11 displays a comparison of contours of u and v velocity at steady state for different Reynolds numbers. Figure 12 shows the u velocity profiles at the vertical centerline and v velocity profiles at the horizontal centerline at steady state for different Reynolds numbers. Figure 13 displays the streamlines of velocity vector at steady state for different Reynolds numbers. When Re=200, at steady state one major vortex whose rotational direction is clockwise is at the middle of the cavity. Two vortices whose rotational direction is counterclockwise are observed at two bottom corners of the cavity (Figure 13 (c1)(c2)). No obvious vortices can be observed at the corners for lower Reynolds numbers, even at steady state. In Figure 14 the profiles of u velocity at the vertical centerline and v velocity at the horizontal centerline we generated for Re=100 at steady state are compared with those generated by Omari[2]. They are in good agreement with each other. In Figure 15 the streamlines of velocity vector we generated for Re=200 at steady state are compared with those generated by Baloch et al.[3]. Good agreement can be clearly observed. 12
13 (a1) Re=50, u velocity (a2) Re=50, v velocity (b1) Re=100, u velocity (b2) Re=100, v velocity (c1) Re=200, u velocity (c2) Re=200, v velocity Figure 11: (a)(b)(c): Contours of u and v velocity at steady state for different Reynolds numbers 13
14 (a1) Re=50, vertical centerline u velocity (a2) Re=50, horizontal centerline v velocity (b1) Re=100, vertical centerline u velocity (b2) Re=100, horizontal centerline v velocity (c1) Re=200, vertical centerline u velocity (c2) Re=200, horizontal centerline v velocity Figure 12: (a)(b)(c): Profiles of u and v velocity at the centerline at steady state for different Reynolds numbers 14
15 (a) Re=50 (b) Re=100 (c1) Re=200 (c2) Re=200, detail of left bottom corner Figure 13: (a)(b)(c): Streamlines of velocity vector at steady state for different Reynolds numbers 15
16 (a1) Our u velocity profile (a2) Omari[2] (b1) Our v velocity profile (b2) Omari[2] Figure 14: (a)(b): Comparison of profiles of u velocity at the vertical centerline and v velocity at the horizontal centerline generated by our codes and by Omari[2] for Re=100 at steady state (a) Our result (b) Baloch et al.[3] Figure 15: (a)(b): Comparison of streamlines of velocity vector generated by our codes and by Baloch et al.[3] for Re=200 at steady state 16
17 4 Summary In this project, we successfully developed a program to solve the two-dimensional unsteady Navier- Stokes equations for an elliptic flow in a square cavity with top wall sliding at a given velocity. Numerical results for flow of different Reynolds numbers at different time intervals and steady state have been shown. Satisfactory agreement with results from other literatures has been achieved. References [1] Joel H. Ferziger and Milovan Peric, Computational Methods for Fluid Dynamics, 3rd edition. [2] Omari R., CFD Simulations of Lid Driven Cavity Flow at Moderate Reynolds Number, European Scientific Journal, 2013, 9(15). [3] Baloch Z. A. K., Baloch M. E., Shaikh H., Computer Simulation of Lid Driven shallow Cavity Flows: Effects of Fluid Inertia and Aspect Ratios. 17
Parabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
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 informationHandout 8 MATLAB Code for the Projection Method and BC Details for the Lid-Driven Cavity Problem
Handout 8 MATLAB Code for the Projection Method and BC Details for the Lid-Driven Cavity Problem Let s simplify the equations of the last page of Handout 7 for the case h = Δx = Δy because the code given
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 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 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 informationNOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL
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 informationFine Grid Numerical Solutions of Triangular Cavity Flow
Published in : The European Physical Journal - Applied Physics (2007) Eur. Phys. J. Appl. Phys. 2007; Vol 38: pp 97-105 Fine Grid Numerical Solutions of Triangular Cavity Flow Ercan Erturk 1 and Orhan
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
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 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 informationFine grid numerical solutions of triangular cavity flow
Eur. Phys. J. Appl. Phys. 38, 97 105 (2007) DOI: 10.1051/epjap:2007057 THE EUROPEAN PHYSICAL JOURNAL APPLIED PHYSICS Fine grid numerical solutions of triangular cavity flow E. Erturk 1,a and O. Gokcol
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
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 informationThomas Pierro, Donald Slinn, Kraig Winters
Thomas Pierro, Donald Slinn, Kraig Winters Department of Ocean Engineering, Florida Atlantic University, Boca Raton, Florida Applied Physics Laboratory, University of Washington, Seattle, Washington Supported
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 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 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 informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
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 informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationLattice Boltzmann Method Solver Documentation
Lattice Boltzmann Method Solver Documentation Release 0.0.1 Sayop Kim May 03, 2016 Contents 1 Code Instruction 3 1.1 Quick instruction for running the simulation.............................. 3 2 Results
More informationNumerical solutions of 2-D incompressible driven cavity flow with wavy bottom surface
American Journal of Applied Mathematics 015; (1-1): 0-4 Published online December 5, 014 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.s.01500101.14 ISSN: 0-004 (Print); ISSN: 0-006X
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More informationDIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER
Conference on Turbulence and Interactions TI26, May 29 - June 2, 26, Porquerolles, France DIRECT NUMERICAL SIMULATION IN A LID-DRIVEN CAVITY AT HIGH REYNOLDS NUMBER E. Leriche, Laboratoire d Ingénierie
More informationDevelopment of a Parallel, 3D, Lattice Boltzmann Method CFD Solver for Simulation of Turbulent Reactor Flow
DOE-FIU SCIENCE & TECHNOLOGY WORKFORCE DEVELOPMENT PROGRAM STUDENT SUMMER INTERNSHIP TECHNICAL REPORT June 4, 2012 to August 10, 2012 Development of a Parallel, 3D, Lattice Boltzmann Method CFD Solver
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 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 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 informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationDUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationComputational Fluid Dynamics (CFD, CHD)*
1 / 1 Computational Fluid Dnamics (CFD, CHD)* PDE (Shocks 1st); Part I: Basics, Part II: Vorticit Fields Rubin H Landau Sall Haerer, Producer-Director Based on A Surve of Computational Phsics b Landau,
More informationNumerical investigation of natural convection inside a wavy enclosure
International Journal of Thermal ciences 42 (2003) 397 406 www.elsevier.com/locate/ijts Numerical investigation of natural convection inside a wavy enclosure Prodip Kumar Das 1, hohel Mahmud,2 Department
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 informationMETHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION
TASKQUARTERLY9No1,17 35 METHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION ZBIGNIEW KOSMA Institute of Applied Mechanics, Radom University of Technology,
More informationSIMULATION OF MIXED CONVECTIVE HEAT TRANSFER USING LATTICE BOLTZMANN METHOD
International Journal of Automotive and Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 2, pp. 130-143, July-December 2010 Universiti Malaysia Pahang DOI: http://dx.doi.org/10.15282/ijame.2.2010.3.0011
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 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 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 informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationEULER AND SECOND-ORDER RUNGE-KUTTA METHODS FOR COMPUTATION OF FLOW AROUND A CYLINDER
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,
More informationRotating-surface-driven non-newtonian flow in a cylindrical enclosure
Korea-Australia Rheology Journal Vol. 22, No. 4, December 2010 pp. 265-272 Rotating-surface-driven non-newtonian flow in a cylindrical enclosure Yalin Kaptan 1, *, Ali Ecder 2 and Kunt Atalik 2 1 Hansung
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 informationA STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES
http://doi.org/10.4867/jpe-017-01-11 JPE (017) Vol.0 (1) Mohapatra, C. R. Preliminary Note A STUDY ON NATURAL CONVECTION HEAT TRANSFER IN COMPLEX BOUNDARIES Received: 3 February 017 / Accepted: 01 April
More informationNumerical Study of the Moving Boundary Problem During Melting Process in a Rectangular Cavity Heated from Below
American Journal of Applied Sciences 4 (4): 25-256, 2007 ISSN 546-9239 2007 Science Publications Corresponding Author: Numerical Study of the Moving Boundary Problem During Melting Process in a Rectangular
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 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 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 informationLATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY
LATTICE BOLTZMANN SIMULATION OF FLUID FLOW IN A LID DRIVEN CAVITY M. Y. Gokhale, Ignatius Fernandes Maharashtra Institute of Technology, Pune 4 38, India University of Pune, India Email : mukundyg@yahoo.co.in,
More informationFINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS
7. ANKARA INTERNATIONAL AEROSPACE CONFERENCE AIAC-2013-143 11-13 September 2013 - METU, Ankara TURKEY FINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS Fatih ÇEVİK
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationAli J. Chamkha* Ahmed Ali Shaker
Progress in Computational Fluid Dynamics, Vol. 12, No. 5, 2012 309 Conduction-combined forced and natural convection in a lid-driven parallelogram-shaped enclosure divided by a solid partition Ali J. Chamkha*
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 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 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 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 informationSingle Curved Fiber Sedimentation Under Gravity. Xiaoying Rong, Dewei Qi Western Michigan University
Single Curved Fiber Sedimentation Under Gravity Xiaoying Rong, Dewei Qi Western Michigan University JunYong Zhu, Tim Scott USDA Forest Products Laboratory ABSTRACT Dynamics of single curved fiber sedimentation
More informationSENSITIVITY ANALYSIS OF THE FACTORS AFFECTING FORCE GENERATION BY WING FLAPPING MOTION
Proceedings of the ASME 2013 International Mechanical Engineering Congress and Exposition IMECE2013 November 15-21, 2013, San Diego, California, USA IMECE2013-65472 SENSITIVITY ANALYSIS OF THE FACTORS
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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationNATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN ENCLOSURES CONTAINING AN INTERNAL LOCAL HEAT SOURCE
Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 24, No. 03, pp. 375-388, July - September, 2007 NATURAL CONVECTION HEAT TRANSFER IN PARTIALLY OPEN
More informationNumerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall
Heat Transfer Research, 2011, Vol. 42, No. 3 Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall SAEID JANI, 1* MEYSAM AMINI, 2 and MOSTAFA MAHMOODI
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationFlow Structure Investigations in a "Tornado" Combustor
Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas
More informationCFD Modeling of Reciprocating Flow around a Bend in Pulse Tube Cryocoolers
CFD Modeling of Reciprocating Flow around a Bend in Pulse Tube Cryocoolers I.Nachman 1, N. Pundak 1, and G. Grossman 2 1 Ricor Cryogenic and Vacuum Systems En Harod Ihud 18960, Israel 2 Faculty of Mechanical
More informationCOMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS
HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 214 Orlando, Florida COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationThe investigation on SIMPLE and SIMPLER algorithm through lid driven cavity
29, Issue 1 (2017) 10-23 Journal of Advanced Research in Fluid Mechanics and Thermal Sciences Journal homepage: www.akademiabaru.com/arfmts.html ISSN: 2289-7879 The investigation on SIMPLE and SIMPLER
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
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 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 informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
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 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 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 informationJMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman
JMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman The exercises will be carried out on PC s in the practicum rooms. Several (Matlab and Fortran) files are required. How these can be obtained
More informationA Cartesian-Grid Discretisation Scheme Based on Local Integrated RBFNs for Two-Dimensional Elliptic Problems
Copyright 2009 Tech Science Press CMES, vol.51, no.3, pp.213-238, 2009 A Cartesian-Grid Discretisation Scheme Based on Local Integrated RBFNs for Two-Dimensional Elliptic Problems N. Mai-Duy 1 and T. Tran-Cong
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 informationApplication of Viscous Vortex Domains Method for Solving Flow-Structure Problems
Application of Viscous Vortex Domains Method for Solving Flow-Structure Problems Yaroslav Dynnikov 1, Galina Dynnikova 1 1 Institute of Mechanics of Lomonosov Moscow State University, Michurinskiy pr.
More informationNumerical study of 2D heat transfer in a scraped surface heat exchanger
Computers & Fluids 33 (2004) 869 880 www.elsevier.com/locate/compfluid Numerical study of 2D heat transfer in a scraped surface heat exchanger K.-H. Sun a, *, D.L. Pyle a, A.D. Fitt b, C.P. Please b, M.J.
More informationResearch Article Direct Simulation of Low-Re Flow around a Square Cylinder by Numerical Manifold Method for Navier-Stokes Equations
Applied Mathematics Volume 22, Article ID 465972, 4 pages doi:.55/22/465972 Research Article Direct Simulation of Low-Re Flow around a Square Cylinder by Numerical Manifold Method for Navier-Stokes Equations
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (24 ) 55 556 th International Conference on Mechanical Engineering, ICME 23 Analysis of heat transfer and flow due to natural
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 informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
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 informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationUNSTEADY MIXED CONVECTION IN A POROUS MEDIA FILLED LID-DRIVEN CAVITY HEATED BY A SEMI-CIRCULAR HEATERS
THERMAL SCIENCE, Year 2015, Vol. 19, No. 5, pp. 1761-1768 1761 UNSTEADY MIXED CONVECTION IN A POROUS MEDIA FILLED LID-DRIVEN CAVITY HEATED BY A SEMI-CIRCULAR HEATERS by Md Mustafizur RAHMAN a,c, Hakan
More informationAlternative and Explicit Derivation of the Lattice Boltzmann Equation for the Unsteady Incompressible Navier-Stokes Equation
ISSN (e): 2250 3005 Volume, 06 Issue, 12 December 2016 International Journal of Computational Engineering Research (IJCER) Alternative and Explicit Derivation of the Lattice Boltzmann Equation for the
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 information