Finite volume method for CFD
|
|
- Ursula Patrick
- 6 years ago
- Views:
Transcription
1 Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor: Dr. Suman Chakraborty, IIT Kharagpur Dr. Vivek V. Buwa, IIT Delhi
2 Overview Conservation Laws Discretisation Concept Concept of Finite Volume Method FVM for 1D Steady state diffusion FVM for Convection and Diffusion Calculation of Flow Field 1. Vorticity Approach 2. Preliminary Variable approach SIMPLE Algorithm SIMPLER Algorithm 1D Unsteady State Heat Conduction
3 Conservation Laws Conservation of Mass Conservation of Momentum Conservation of Energy
4 Conservation of Chemical Species ( ρm t l ) + div ( ρum l + J l ) = R t δz δy Rate of change of mass Convection Flux Diffusion flux δx *m l : mass fraction of chemical species u : velocity field R t : rate of generation of the chemical species per unit volume. By the use of Flick's Law J 1 = Γ l (grad m l ). Γ l : diffusion coefficient.
5 Conservation of Momentum ( ρ t u) + div ( ρ uu ) = div ( μ gradu ) p x + + Bx Vx Transient term Advective term Diffusive term μ: : Coefficient of viscosity p: pressure B x : x-direction x body force per unit volume. Vx: Viscous terms in addition to those expressed by div(μ grad u).
6 Conservation of Energy For a steady low-velocity flow with negligible viscous dissipation is Energy equation is div (ρuh) = div( kgradt) + S h h: specific enthalpy k: thermal conductivity T : Temperature S h : Volumetric Rate of heat generation For ideal gases and for liquids and solids, h= ct, for a constant c. Equation reduces to: div(ρut ) k = div( gradt ) c + S c h
7 General Differential Equation ( ρφ ) t + div( ρuφ ) = div( Γgrad φ ) + S Unsteady term Convection term Diffusion term Source term Γ : diffusion coefficient Φ : Dependent variable, stand for different quantities, like mass fraction of chemical species, enthalpy or temperature, velocity component
8 Nature of Coordinates Φ = Φ( ( x,y,z,t) With respect to x: 1D, 2D, 3D With respect to t: unsteady and steady Computational Saving is achieved by decreasing independent variable, eg: flow around an airplane, axisymmetric flow in circular pipe etc.. One way (parabolic) and two way (elliptic) coordinates.
9 Inference from Conservation laws Conservation equations are governed by partial differential equation. These equations can be represented as a general equation for the variable φ. Analytical Solutions exist for few cases and many a times complex. Numerical methods have to applied. Numerical method should employ the values of Φ as the primary unknown.
10 Discretization Concept Transferring continuous models and equations into discreet counterparts.
11 Concept of Finite Volume Method Subdivide the problem domain into non overlapping control volumes Integrate the governing equation over each of these control volumes. Representative value of variable = value of the variable at the geometric centre of CV Piecewise profiles expressing the variation of φ between the grid points are used to evaluate the integrals. Discretized equation obtained in this manner, represents the conservation principle for the finite control volume.
12 Finite Volume Grid δ x w δ x e N Δ y W P E δ y n δ y s S Δ x
13 FVM for 1D steady State diffusion Governing equation of diffusion is d dx d φ ( Γ ) + S = dx 0 Γ : Diffusion coefficient S: Source term Boundary values of A and B are prescribed
14 Contd. Step 1: Grid generation w Δx e W Step 2: Discretization P E Δv d ( Γ dx dφ ) dv dx + Δv SdV dφ = ( ΓA ) dx e dφ ( ΓA ) dx w + S ΔV = 0 A: cross-sectional area of control volume ΔV: Volume of Control Volume S : Average value of Source over control volume.
15 Contd. Interpolating linearly Γ w = ( ΓW + Γ P ) ( Γ P + Γ E, Γ e = 2 2 Diffusive Flux terms are as follows ( ΓA ( ΓA dφ ( φ E φ P ) ) e = Γe Ae dx δx PE dφ ( φw φ P ) ) w = Γw Aw dx δx WP ) Using Central differencing scheme
16 Contd. Finally the equation is: a p Φ p = a w Φ w + a e Φ E + S u a W a E a P ΓW A δ x WP W Γ E A δ x Now, Discretized equations must be setup at each of the nodal points PE E a W + a E S P
17 FVM for Convection and Diffusion problem Governing equation for the process is: d dx ( ρuφ ) The Continuity equation : = d dx ( Γ dφ ) dx d dx ( ρ u ) = 0 There are various approaches to find the numerical solution of a convection diffusion problem.
18 Discretization Scheme Central Differencing Scheme Upwind differencing scheme Hybrid differencing scheme Power Law Scheme Quick Scheme
19 Calculation of Flow Field Flow field needs to be calculated from appropriate governing equations Velocity component is governed by momentum equations Difficulties: 1. Pressure gradient forms a part of the source term for the momentum equation. 2. There is no obvious equation for obtaining Pressure. Once the correct pressure field is substituted in momentum equation, resultant Velocity field satisfies Continuity equation
20 Vorticity Based Methods Pressure gradient term is eliminated from the momentum equation by cross differentiating the 2 components of the momentum equation (in a 2-D flow) and subtracting them. This gives rise to the following equation involving Vorticity vector ζ. The equation for the stream function ψ(x,y) is given as ] [ y x y v x u t + = + + ζ ζ ρ μ ζ ζ ζ ζ ψ ψ = y x
21 Assessment of Vorticity based Merits: method 1. Pressure makes no appearance 2. Need to solve only two equations to obtain ζ and and ψ. 3. Poisson s equation can be used to derive pressure. Demerits: 1. Value of vorticity at a wall is difficult to specify and causes trouble in getting converged solution 2. Pressure is required for calculation of density and other fluid properties. 3. Extraction of pressure from vorticity offsets computational saving 4. Cannot be extended to 3-D situations for which a streamfunction does not exist.
22 Primitive Variables Method To overcome the shortcomings of the vorticity method primitive variable method is used Dependent variables are the velocity components and pressure. Need to solve the velocity equations and the continuity equation simultaneously. Main task is to convert the indirect information in the continuity equation into direct algorithm for the calculation of pressure.
23 Representation of Pressure Gradient term W w P e E Integration of pressure gradient over control volume gives P w - P e Assuming linear profile: P w -P e = (P W + P P )/2 (P P + P E )/2 = (P W P E )/2 Momentum equation contains the pressure difference between alternate grid points.
24 Contd It is a zig-zag zag pressure field Formula predicts it as uniform pressure field. Similar is the case with 2D and 3D fields.
25 Representation of the Continuity Equation Continuity equation for 1D can be written as u x = 0 Integrating over control volume u E -u W = 0 Again the continuity equation wants the equality of velocities at alternate grid points Where is the Flaw? The flaw is in the fact that we are calculating all the variables for the same grid points.
26 REMEDY: The Staggered Grid Arrange the velocity components on a different grid. This grid is called the staggered grid. Velocity components are calculated at CV faces which form the staggered grid Pressure is evaluated at the main grid points.
27 Staggered Grid Dotted line represents the Staggered Grid. Velocity Components lie on the faces of CV. u point v point P point
28 Advantages of Staggered grid Discretized Continuity equation will contain difference of adjacent velocity components. Pressure difference between two adjacent grid points becomes the driving force for the velocity component located between them.
29 Momentum Equation Discretised x momentum equation for the staggered CV shown in the fig is: a e u e = a nb u nb + b + (P P -P E )A e subscript nb represents neighboring term. Discretised y momentum equation for the staggered CV shown in the fig is a n u n = a nb u nb + b + (P P -P N )A n Staggered grid for Velocity
30 Contd. Momentum equation can be solved only if pressure field is given or somehow estimated. Unless the correct pressure field is known the resulting velocity field will not satisfy the continuity equation. Any inaccurate velocity field resulting from a guessed value p*,will be denoted by u*,v*,w*. They are related by the following momentum equations : a e u e * = a nb u nb * + b + (P P *- P E *)A e a n u n * = a nb u nb * + b + (P P *- P N *)A n a t w t * = a nb w nb * + b + (P P *- P t *)A t Aim: find a way of improving these guessed values p* such that the resulting value of velocities satisfy the continuity equation. Hence we employ a pressure correction and the corresponding velocity correction.
31 Velocity Corrections Suppose we employ a correction p` to the guessed value p* of the pressure and the corresponding corrections in velocity be u',v',w'.then using the momentum equation the following relation can be derived. If effect of neighbors is neglected u e = u e * + d e ( P P` - P E`) v n = v n * + d n ( P P` - P N`) w n = w n * + d t ( P P` - P t`) where d e = A e /a e a e u e` = a nb u nb` + (P P`- P E`)A e
32 Pressure-Correction Equation We get the correction for pressure from the continuity equation. ρ ( ρu) ( ρv) ( ρw) = 0 t x y z Substitution of corrected velocities in the discretised form of pressure equation gives the equation of the following form a p p P`= a E p`e + a W p`w p` + a S p`s + a T p`t + a B p B` ` + b where a $ = ρ $ d $ A #, A # = Area perpendicular to line p # a P = a E + a W + a N + a S + a T + a B b = f( u*, v*, w*)
33 Keep the corrected pressure p as p* SIMPLE ALGORITHM Start Guess pressure field Solve Discretised momentum equation to obtain u*, v*, w* Solve the p` equation Calculate p by adding p` to previous value of p* Calculate u,v,w from the velocity correction relation Solve all other discretized transport equations No Convergence? YES STOP
34 Discussion of Pressure-correction Equation Σa nb u`nb is neglected in the velocity correction equation, which enables to cast the pressure equation in a general conservative form. The algorithm is called semi implicit because, the term Σa nb u`nb which represents the implicit influence of the pressure correction on velocities is dropped. On convergence all the corrections tend to zero and there is no error induced on dropping Σa nb u`nb For high pressure flows the pressure correction formula should accommodate density correction term. Pressure obtained is a relative variable as an outcome of this algorithm and not a absolute quantity.
35 Drawbacks of SIMPLE Approximation of the pressure correction by neglecting the term Σa nb u nb in the SIMPLE algorithm leads to a rather exaggerated pressure correction. Due to the omission of the neighboring velocity corrections, pressure correction has to take the entire burden of correcting the velocities. This does a pretty good job of correcting velocity it does a poor job of correcting pressure. Hence we have to do many iterations before convergence could result.
36 SIMPLER: Revised SIMPLE Separate pressure equation is formulated to calculate the exact pressure field Velocity field is still calculated using pressure correction formula The Pressure equation: The momentum equation can be written as u e = Σa nb u a nb e + b + d e( Pp PE ) where d e = A e /a e
37 Contd. Define û e ( pseudo velocity ) = ( a nb u nb + b )/a e Hence the velocity is u e = û e + d e ( p P p E ) Now if we substitute this velocity field in the discretised continuity equation we get the following pressure equation. a p p P = a E p E + a W p W + a N p N + a S p S + a T p T + a B p B + b where all a s are same as before.
38 SIMPLER Algorithm START Guess Velocity Field Calculate coefficients of discretised momentum equation Calculate û,v and w by substituting u nb Calculate coefficients of pressure equation and solve it to get pressure field P Use P as P* and solve for u*,v*,w* Solve pressure correction equation Correct velocity field But not pressure and solve other dicretised equations No Convergence Yes STOP
39 SIMPLER v/s SIMPLE Although SIMPLER has been found to give faster convergence than SIMPLE, it should be recognized that one iteration of SIMPLER involves much more computational work. It involves solving the pressure equation. Calculation of û, v and w which is not there in SIMPLE. However since SIMPLER requires fewer iterations for convergence, the additional effort required is justified.
40 1D Unsteady State Heat Conduction Governing Equation is T T ρc = ( k ) + t x x S W w P e E W W δx we
41 Contd. Integrating over control volume and over a time interval from t to t +Δt+ t + Δt t + Δt t CV T T ρc dvdt = ( k ) dvdt t x x This can be written as t CV + t +Δt t CV sdvdt e w t+ Δt t+ Δt T [ ρc t t dt] dv = t T [( ka x ) e T ( ka x ) w ] dt + t+δt t _ S ΔVdt
42 Contd. Temperature at node is assumed to prevail over whole control volume CV t+δt T o [ ρc dt] dv = ρc( Tp Tp ) ΔV t t ρc( T P T o P ) ΔV = t+δt t [( k e TE T A δx PE P ) ( k w TP T A δx WP W )] dt + t+δt t SΔVdt (applying central differencing scheme)
43 Contd. Integration wrt time can be carried out in the following way t + Δt Tdt θ = 0, Explicit Scheme θ = 1, Implicit Scheme t = [ θt (1 θ ) T ] Δt θ = 0.5, Crank-Nicolson scheme P o P
44 Contd. The Final Discretized equation is a P T P + [ a o P = aw[ θt W (1 θ) a W + (1 θ) T o W (1 θ) a E ] + ] T a o P E [ θt + b E + (1 θ) T o E ] a P = θ(a W +a E )+a o P a W = K w /δx WP, a E = K e /δx PE, b = SΔx
45 Conclusion The basic conservation equations can be written in the form of a differential equation for a general variable φ. The concept of discretization can be used as a tool to solving these differential equations over a domain of interest. Finite Volume Method can be used to find the discreet solution of flow fields as well as convection diffusion problems. We discussed the various schemes under the FVM to find out the solutions of a given flow field. Finally we conclude that Finite Volume Method is a useful tool to find the solution of a given flow field.
46 Thank You
Discretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
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 informationChapter 5. The Finite Volume Method for Convection-Diffusion Problems. The steady convection-diffusion equation is
Chapter 5 The Finite Volume Method for Convection-Diffusion Problems Prepared by: Prof. Dr. I. Sezai Eastern Mediterranean University Mechanical Engineering Department Introduction The steady convection-diffusion
More informationFINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS
IJRET: International Journal of Research in Engineering and Technology eissn: 39-63 pissn: 3-738 FINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS Asish Mitra Reviewer: Heat
More informationFINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES SHRUTI JAIN B.Tech III Year, Electronics and Communication IIT Roorkee Tutors: Professor G. Biswas Professor S. Chakraborty ACKNOWLEDGMENTS I would like
More informationSolution Methods. Steady convection-diffusion equation. Lecture 05
Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationOn A Comparison of Numerical Solution Methods for General Transport Equation on Cylindrical Coordinates
Appl. Math. Inf. Sci. 11 No. 2 433-439 (2017) 433 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.18576/amis/110211 On A Comparison of Numerical Solution Methods
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationDraft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
Draft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Instructor: Jayathi Y. Murthy School of Mechanical Engineering Purdue University Spring 00 c 1998 J.Y. Murthy and S.R. Mathur.
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 informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
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 informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationVALIDATION OF ACCURACY AND STABILITY OF NUMERICAL SIMULATION FOR 2-D HEAT TRANSFER SYSTEM BY AN ENTROPY PRODUCTION APPROACH
Brohi, A. A., et al.: Validation of Accuracy and Stability of Numerical Simulation for... THERMAL SCIENCE: Year 017, Vol. 1, Suppl. 1, pp. S97-S104 S97 VALIDATION OF ACCURACY AND STABILITY OF NUMERICAL
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 informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationUnsteady MHD Couette Flow with Heat Transfer in the Presence of Uniform Suction and Injection
Mechanics and Mechanical Engineering Vol. 12, No. 2 (2008) 165 176 c Technical University of Lodz Unsteady MHD Couette Flow with Heat Transfer in the Presence of Uniform Suction and Injection Hazem A.
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
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 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 informationIntroduction to Fluid Dynamics
Introduction to Fluid Dynamics Roger K. Smith Skript - auf englisch! Umsonst im Internet http://www.meteo.physik.uni-muenchen.de Wählen: Lehre Manuskripte Download User Name: meteo Password: download Aim
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More 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 informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationModule 2: Introduction to Finite Volume Method Lecture 14: The Lecture deals with: The Basic Technique. Objectives_template
The Lecture deals with: The Basic Technique file:///d /chitra/nptel_phase2/mechanical/cfd/lecture14/14_1.htm[6/20/2012 4:40:30 PM] The Basic Technique We have introduced the finite difference method. In
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 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 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 informationFluid Flow Modelling with Modelica
NE T E C H N I C A L N O T E Marco Bonvini 1*, Mirza Popovac 2 1 Dipartimento di Elettronica e Informazione, Politecnico di Milano, Via Ponzio 34/5, 20133 Milano, Italia; 1 * bonvini@elet.polimi.it 2 Austrian
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #12 Fundamentals of Discretization: Finite Volume Method
More informationCFD Analysis On Thermal Energy Storage In Phase Change Materials Using High Temperature Solution
CFD Analysis On Thermal Energy Storage In Phase Change Materials Using High Temperature Solution Santosh Chavan 1, M. R. Nagaraj 2 1 PG Student, Thermal Power Engineering, PDA College of Engineering, Gulbarga-585102,
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 informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More 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 informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationOpenFOAM selected solver
OpenFOAM selected solver Roberto Pieri - SCS Italy 16-18 June 2014 Introduction to Navier-Stokes equations and RANS Turbulence modelling Numeric discretization Navier-Stokes equations Convective term {}}{
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
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 informationAnswers to Exercises Computational Fluid Dynamics
Answers to Exercises Computational Fluid Dynamics Exercise - Artificial diffusion upwind computations.9 k=. exact.8 k=..7 k=..6 k=..5.4.3.2...2.3.4.5.6.7.8.9 x For k =.,. and., and for N = 2, the discrete
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationDevelopment of unsteady algorithms for pressurebased unstructured solver for two-dimensional incompressible flows
Graduate Theses and Dissertations Graduate College 2009 Development of unsteady algorithms for pressurebased unstructured solver for two-dimensional incompressible flows Angela Dwi Lestari Iowa State University
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 information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationBERNOULLI EQUATION. The motion of a fluid is usually extremely complex.
BERNOULLI EQUATION The motion of a fluid is usually extremely complex. The study of a fluid at rest, or in relative equilibrium, was simplified by the absence of shear stress, but when a fluid flows over
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationTotal energy in volume
General Heat Transfer Equations (Set #3) ChE 1B Fundamental Energy Postulate time rate of change of internal +kinetic energy = rate of heat transfer + surface work transfer (viscous & other deformations)
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 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 informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationNumerical modelling of the dynamic contact angle
GR Numerical modelling of the dynamic contact angle Simon van Mourik height(x,t) 0.055 0.05 0.045 0.04 0.035 0.03 4 0.025 0 0.01 0.02 0.03 0 1 2 time 3 space Department of Mathematics GR Master s Thesis
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationME 608 Numerical Methods in Heat, Mass, and Momentum Transfer. q 0 = εσ ( T 4 T 4) Figure 1: Computational Domain for Problem 1.
ME 68 Numerical Methods in Heat, Mass, and Momentum Transfer Mid-Term Examination Solution Date: March 1, 21 6: 8: PM Instructor: J. Murthy Open Book, Open Notes Total: 1 points 1. Consider steady 1D conduction
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
More informationINDIAN INSTITUTE OF TECHNOOGY, KHARAGPUR Date: -- AN No. of Students: 5 Sub. No.: ME64/ME64 Time: Hours Full Marks: 6 Mid Autumn Semester Examination Sub. Name: Convective Heat and Mass Transfer Instructions:
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationIntroduction to Heat and Mass Transfer. Week 10
Introduction to Heat and Mass Transfer Week 10 Concentration Boundary Layer No concentration jump condition requires species adjacent to surface to have same concentration as at the surface Owing to concentration
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 informationDevelopment of an explicit pressure-based unstructured solver for three-dimensional incompressible flows with graphics hardware acceleration
Graduate Theses and Dissertations Graduate College 2012 Development of an explicit pressure-based unstructured solver for three-dimensional incompressible flows with graphics hardware acceleration Jordan
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 informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More 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 informationINTEGRAL ANALYSIS OF LAMINAR INDIRECT FREE CONVECTION BOUNDARY LAYERS WITH WEAK BLOWING FOR SCHMIDT NO. 1
INTEGRA ANAYSIS OF AMINAR INDIRECT FREE CONVECTION BOUNDARY AYERS WITH WEAK BOWING FOR SCHMIDT NO. Baburaj A.Puthenveettil and Jaywant H.Arakeri Department of Mechanical Engineering, Indian Institute of
More information1. INTRODUCTION TO CFD SPRING 2018
1. INTRODUCTION TO CFD SPRING 018 1.1 What is computational fluid dynamics? 1. Basic principles of CFD 1.3 Stages in a CFD simulation 1.4 Fluid-flow equations 1.5 The main discretisation methods Appendices
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 informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
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 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 informationV/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0
UNIT III FLOW THROUGH PIPES 1. List the types of fluid flow. Steady and unsteady flow Uniform and non-uniform flow Laminar and Turbulent flow Compressible and incompressible flow Rotational and ir-rotational
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More information6.2 Modeling of Systems and Components
Chapter 6 Modelling of Equipment, Processes, and Systems 61 Introduction Modeling is one of the most important elements of thermal system design Most systems are analyzed by considering equations which
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 informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationThe Advanced Research WRF (ARW) Dynamics Solver
Dynamics: Introduction The Advanced Research WRF (ARW) Dynamics Solver 1. What is a dynamics solver? 2. Variables and coordinates 3. Equations 4. Time integration scheme 5. Grid staggering 6. Advection
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationLocal Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf**
Controlled CO 2 Diversified fuels Fuel-efficient vehicles Clean refining Extended reserves Local Time Step for a Finite Volume Scheme I.Faille F.Nataf*, F.Willien, S.Wolf** *: Laboratoire J.L.Lions **:Université
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 informationTHREE-DIMENSIONAL FINITE DIFFERENCE MODEL FOR TRANSPORT OF CONSERVATIVE POLLUTANTS
Pergamon Ocean Engng, Vol. 25, No. 6, pp. 425 442, 1998 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0029 8018/98 $19.00 + 0.00 PII: S0029 8018(97)00008 5 THREE-DIMENSIONAL FINITE
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More information