FINITE VOLUME METHOD: BASIC PRINCIPLES AND EXAMPLES
|
|
- Hortense Golden
- 6 years ago
- Views:
Transcription
1 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
2 ACKNOWLEDGMENTS I would like to thank my mentor, Professor G. Biswas without whose help, this presentation would not have been possible. My presentation is a result of his lecture notes and his consistent guidance. Participating in the winter academy has been a learning experience and I would like to thank all the members of the academy for giving me this opportunity. It has been a pleasure working on this presentation.
3 FINITE VOLUME METHOD: AN OUTLINE Introduction to Finite Volume Method 2 approaches -Weighted Residual Approach -Physical Approach Boundary Conditions in Finite Volume Method Equations with First Order Derivatives Equations with Second Order Derivatives
4 FINITE VOLUME METHOD Finite Volume Method is a sub domain method with piecewise definition of the field variable in the neighborhood of chosen control volumes. The total solution domain is divided into many small control volumes which are usually rectangular in shape. Nodal points are used within these control volumes for interpolating the field variable and usually, a single node at the centre of the control volume is used for each control volume. This method was developed by Patankar and Spalding and they proposed the use of physical approach for deriving the nodal equations.
5 FINITE VOLUME METHOD Consider a 2-D, steady heat conduction in rectangular geometry (Figure 1).The 2-D heat conduction equation is (1) Where T(x,y) is the temperature field, k is the thermal conductivity and Q is the heat generation per unit volume. The two alternative ways of settings up the nodal equations are the weighted residual approach and the physical approach
6 WEIGHTED RESIDUAL APPROACH Weighted residual approach for 2-D heat conduction equation: (2) Where the weight W i =1 within the i th control volume W i =0 outside the i th control volume Thus, we get, for each i= 1, n (3) Integrating the above equation by parts, we get:
7 WEIGHTED RESIDUAL APPROACH Typical control volume Figure 1: Grid Arrangement for the Finite Volume Method Where the Gauss divergence theorem has been used to convert the volume integral to a surface integral. (4)
8 WEIGHTED RESIDUAL APPROACH The meaning of Equ. (4) is that the net heat generation rate in the control volume is equal to the net sum of the rate of heat energy going out of the control volume ) where c i is the boundary of the control volume
9 PHYSICAL APPROACH Figure 2: Balance of Heat Flux in a Control Volume
10 PHYSICAL APPROACH Considering the balance of heat flux in Figure 2 for the control volume (for unit depth in z direction): Net rate of heat leaving the control volume through the boundary = Rate of heat generation within the control volume (CV) at steady state: (5) which is the same equation as obtained before.
11 PHYSICAL APPROACH In the implementation of the finite volume method, the heat fluxes are represented in terms of nodal temperatures (e.g. T E etc. at the CV centers). Thus, assuming temperature to have linear variation between points E and P, the heat flux q e can be evaluated as follows: (6) While deriving equation (6) it has been assumed that the cell size is, constant in x-direction. Similarly, q w is given by (7)
12 PHYSICAL APPROACH Using similar expressions for q n and q s also, the nodal equation for point P becomes: (8) This equation can be rewritten in the familiar form used in finite difference as: Where, (9) During numerical implementation, the subscripts E, W, etc. will be changed to numerical indices of i, j and solved in the some way.
13 BOUNDARY CONDITIONS: When heat flux at a boundary is specified Figure 3: Prescribed Heat Flux at the Boundary The control volumes adjacent to the x = 0 boundary as shown Figure 3, the term q w will be substituted by q 0 in equation (5).Thus, (10)
14 BOUNDARY CONDITIONS: When boundary temperature is specified When boundary temperature is specified, the control volume shapes near the boundary can be changed to facilitate the implementation of the boundary conditions. For example: consider the BC, T=T L at boundary x=l. For nodes on boundary an imaginary extension of CV s outside the actual domain can be considered. We consider the boundary to be the center of the boundary cells and their width to be x/2, thus reducing the width of the adjacent cells to 3 x/4. Figure 4: Boundary Condition, at x=l, T= T L
15 BOUNDARY CONDITIONS The nodal equation for the adjacent cell P will be written considering a shortened control volume (11) Where Note that T L has been substituted instead of T W is the above equation, so the boundary condition is being directly applied.
16 EQUATION WITH FIRST DERIVATIVES General first-order equation: (12) For the above equation is the inviscid momentum equation in the x-directions. In a similar manner, for x direction viscous momentum equation, (13)
17 EQUATION WITH FIRST DERIVATIVES Figure-5: Finite Volume Mesh System
18 EQUATION WITH FIRST DERIVATIVES We consider the area integral of general first order equation over Ω p (14) Ω p is the finite volume (quadrilateral) ABCD shown in figure 5 Green s theorem: (15)
19 EQUATION WITH FIRST DERIVATIVES Using Green s theorem the equation (14) becomes: (16) where H = (F,G) and n is the outward unit normal of segment ds. On a counter clockwise contour dyi dx j 2 2, where ds = dx dy n Hence H.n ds = F dy G dx. ds The finite volume method is a discretization of the governing equation in integral form, in contrast to the finite difference method, which is unusually applied to the governing equation in differential form. (17)
20 EQUATION WITH FIRST DERIVATIVES Where A is the area of the quadrilateral ABCD in Figure-5, and the average value of E over the quadrilateral is represented by E i,j and the remaining terms are approximations for the line integral over segments AB, BC, CD and DA respectively Figure-5: Finite Volume Mesh System
21 EQUATION WITH FIRST DERIVATIVES With similar expressions for time, then:, etc. If A is not a function of (18) If the grid-mesh is uniform and coincides with lines of constant x and y, the above equation (18) becomes: (19) which coincides with a central difference representation for the spatial terms in equation (12)
22 EQUATIONS WITH SECOND DERIVATIVES Let us consider Laplace s equation: (20) The finite volume method follows the application of the subdomain method to equation (20). where Green s theorem has again been used and now: (21)
23 EQUATIONS WITH SECOND DERIVATIVES The above equation can be evaluated approximately over the segments AB, BC,CD and DA by: (22)
24 EQUATIONS WITH SECOND DERIVATIVES Figure 4: Finite Volume for a Curvilinear Grid Here is evaluated as the mean value over the area. in Figure 4.
25 EQUATIONS WITH SECOND DERIVATIVES Thus we can write: (23) (24) Using Green s theorem: A similar expression can be written for dx A' B' C' D' i, j 1 x A' B' x B B' C' i, j x C' D' x A D' A'
26 EQUATIONS WITH SECOND DERIVATIVES If the mesh is not too distorted, Hence, y y x x S A' B' B' C' A' B' B' C' A, B y y x x S C' D' D' A' C' D' D' A' A' B' C' D' y y x x AB j 1, j AB j 1, j A' B' B' C' (25)
27 EQUATIONS WITH SECOND DERIVATIVES Therefore the equation (23) & (24) becomes: (26) (27) Now, the first line of equation (22) can be evaluated as: (28)
28 EQUATIONS WITH SECOND DERIVATIVES In a similar manner, one can evaluate the following: (29) (30)
29 EQUATIONS WITH SECOND DERIVATIVES If the mesh is not too distorted, The equivalent expressions for similar terms in equation (22), finally yields: and other (31)
30 EQUATIONS WITH SECOND DERIVATIVES Where the geometrical parameters are: In equation (31) are evaluated as the average of the four surrounding nodal values. Thus
31 EQUATIONS WITH SECOND DERIVATIVES Substitution into equation (31) generates the following nine-points discretization of equation (20): (32) We need to calculate the geometric quantities like M AB, N AB etc. only once for a given grid and can use the obtained values for all subsequent calculations.
32 EQUATIONS WITH SECOND DERIVATIVES Equation (32) is solved conveniently using a Successive Over- Relaxation (SOR) technique. (33) and the improved better value is: where is the relaxation parameter. (34) The discretised equation (32) reduces to the centre finite difference scheme on a uniform rectangular grid (35)
33 SUMMARY Finite Volume Method is a sub domain method with piecewise definition of the field variable in the neighborhood of chosen control volumes. The total solution domain is divided into many small control volumes which are usually rectangular in shape. Boundary Conditions may be applied in finite volume equations by substituting for heat fluxes at boundaries or by changing the control volume spacing (while substituting for boundary temperatures). The finite volume method is a discretization of the governing equation in integral form, in contrast to the finite difference method, which is unusually applied to the governing equation in differential form. Finite volume method can be applied in first and second order equations and the discretized equation finally reduces to the central finite difference scheme on a uniform rectangular grid. One attractive feature of the finite volume method is that it can handle Neumann boundary condition as readily as the Dirichlet boundary condition.
34 REFERENCES Lecture notes of Professor G. Biswas Fletcher, 'Computational Techniques for Fluid Dynamics Peyret and Taylor, Computational Methods for Fluid Flow'
35 Thank You!
Module 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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
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 informationChapter 4. Two-Dimensional Finite Element Analysis
Chapter 4. Two-Dimensional Finite Element Analysis general two-dimensional boundary-value problems 4.1 The Boundary-Value Problem 2nd-order differential equation to consider α α β φ Ω (4.1) Laplace, Poisson
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 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 informationDivergence Theorem and Its Application in Characterizing
Divergence Theorem and Its Application in Characterizing Fluid Flow Let v be the velocity of flow of a fluid element and ρ(x, y, z, t) be the mass density of fluid at a point (x, y, z) at time t. Thus,
More informationDiscretization 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 informationFinite volume method for CFD
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:
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 informationSolving the Generalized Poisson Equation Using the Finite-Difference Method (FDM)
Solving the Generalized Poisson Equation Using the Finite-Difference Method (FDM) James R. Nagel September 30, 2009 1 Introduction Numerical simulation is an extremely valuable tool for those who wish
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
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 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 informationApplications in Fluid Mechanics
CHAPTER 8 Applications in Fluid 8.1 INTRODUCTION The general topic of fluid mechanics encompasses a wide range of problems of interest in engineering applications. The most basic definition of a fluid
More information8. Introduction to Computational Fluid Dynamics
8. Introduction to Computational Fluid Dynamics We have been using the idea of distributions of singularities on surfaces to study the aerodynamics of airfoils and wings. This approach was very powerful,
More informationBrief Review of Vector Algebra
APPENDIX Brief Review of Vector Algebra A.0 Introduction Vector algebra is used extensively in computational mechanics. The student must thus understand the concepts associated with this subject. The current
More informationModule 2 : Convection. Lecture 12 : Derivation of conservation of energy
Module 2 : Convection Lecture 12 : Derivation of conservation of energy Objectives In this class: Start the derivation of conservation of energy. Utilize earlier derived mass and momentum equations for
More informationElements of Vector Calculus : Line and Surface Integrals
Elements of Vector Calculus : Line and Surface Integrals Lecture 2: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay In this lecture we will talk about special functions
More informationNumerical Solution for Two Dimensional Laplace Equation with Dirichlet Boundary Conditions
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, Volume 6, Issue 4 (May. - Jun. 2013), PP 66-75 Numerical Solution for Two Dimensional Laplace Equation with Dirichlet Boundary
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More informationChapter 3. Formulation of FEM for Two-Dimensional Problems
Chapter Formulation of FEM for Two-Dimensional Problems.1 Two-Dimensional FEM Formulation Many details of 1D and 2D formulations are the same. To demonstrate how a 2D formulation works we ll use the following
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationdifferential problems
A numerical study of 2D integrated RBFNs incorporating Cartesian grids for solving 2D elliptic differential problems N. Mai-Duy and T. Tran-Cong Computational Engineering and Science Research Centre Faculty
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationexample consider flow of water in a pipe. At each point in the pipe, the water molecule has a velocity
Module 1: A Crash Course in Vectors Lecture 1: Scalar and Vector Fields Objectives In this lecture you will learn the following Learn about the concept of field Know the difference between a scalar field
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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
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 informationSolution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method
Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July 11-15, 2016 ICCFD9-0113 Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
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 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 informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationMath 212-Lecture 20. P dx + Qdy = (Q x P y )da. C
15. Green s theorem Math 212-Lecture 2 A simple closed curve in plane is one curve, r(t) : t [a, b] such that r(a) = r(b), and there are no other intersections. The positive orientation is counterclockwise.
More informationLecture 27: Structural Dynamics - Beams.
Chapter #16: Structural Dynamics and Time Dependent Heat Transfer. Lectures #1-6 have discussed only steady systems. There has been no time dependence in any problems. We will investigate beam dynamics
More informationIndustrial Applications of Computational Mechanics Extended Civil Engineering Problems. Prof. Dr.-Ing. Casimir Katz SOFiSTiK AG
Industrial Applications of Computational Mechanics Extended Civil Engineering Problems Prof. Dr.-Ing. Casimir Katz SOFiSTiK AG Beyond Structural Analysis Multiphysics» Heat Conduction, Radiation» Potential
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 informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
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 informationNumerical Solution of Laplace Equation By Gilberto E. Urroz, October 2004
Numerical Solution of Laplace Equation By Gilberto E. Urroz, October 004 Laplace equation governs a variety of equilibrium physical phenomena such as temperature distribution in solids, electrostatics,
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationFluid Thermal Interaction of High Speed Compressible Viscous Flow Past Uncooled and Cooled Structures by Adaptive Mesh
The 20th Conference of echanical Engineering Network of Thailand 18-20 October 2006, Nakhon Ratchasima, Thailand luid Thermal nteraction of High peed Compressible Viscous low Past Uncooled and Cooled tructures
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationSection 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow
Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Anderson: Chapter 2 pp. 41-54 1 Equation of State: Section 1 Review p = R g T " > R g = R u M w - R u = 8314.4126
More informationPreamble. Circuit Analysis II. Mesh Analysis. When circuits get really complex methods learned so far will still work,
Preamble Circuit Analysis II Physics, 8 th Edition Custom Edition Cutnell & Johnson When circuits get really complex methods learned so far will still work, but they can take a long time to do. A particularly
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 informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationIntroduction to Computational Fluid Dynamics
AML2506 Biomechanics and Flow Simulation Day Introduction to Computational Fluid Dynamics Session Speaker Dr. M. D. Deshpande M.S. Ramaiah School of Advanced Studies - Bangalore 1 Session Objectives At
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
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 informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationMultivariable Calculus Lecture #12 Notes
Multivariable alculus Lecture #12 Notes In this lecture, we will develop a list of statements equivalent to a vector field being conservative and state and prove Green s heorem to help connect these facts.
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
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 informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More informationConservation and dissipation principles for PDEs
Conservation and dissipation principles for PDEs Modeling through conservation laws The notion of conservation - of number, energy, mass, momentum - is a fundamental principle that can be used to derive
More informationWave Theory II (6) Finite Element Method
Wave Theory II (6) Finite Method Jun-ichi Takada (takada@ide.titech.ac.jp) In this lecture, the finite element method (FEM) is described. The Helmholtz equation and the boundary condition are transformed
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 informationAE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Finite Element Methods Michael Bader Technical University of Munich Summer 2016 Part I Looking Back: Discrete Models for Heat Transfer and the Poisson Equation Modelling
More informationComputational Fluid Dynamics Overview, Scope and Relevance
Computational Fluid Dynamics Overview, Scope and Relevance P. R. Naren School of Chemical & Biotechnology SASTRA University Thanjavur 613401 E-mail: prnaren@scbt.sastra.edu at Faculty Development Program
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
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 information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationME19b. FINAL REVIEW SOLUTIONS. Mar. 11, 2010.
ME19b. FINAL REVIEW SOLTIONS. Mar. 11, 21. EXAMPLE PROBLEM 1 A laboratory wind tunnel has a square test section with side length L. Boundary-layer velocity profiles are measured at two cross-sections and
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 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 informationarxiv: v2 [physics.comp-ph] 4 Feb 2014
Fast and accurate solution of the Poisson equation in an immersed setting arxiv:1401.8084v2 [physics.comp-ph] 4 Feb 2014 Alexandre Noll Marques a, Jean-Christophe Nave b, Rodolfo Ruben Rosales c Abstract
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture 1-20 September, 2017 Institute of Structural Engineering Method of Finite Elements II 1 Course
More information1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM)
1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM) 1.1 Introduction In this work, performances of two most widely
More informationY. M. TEMIS AND V. V. KARABAN 1
J. KSIAM Vol.5, No.2, 39-51, 2001 BOUNDARY ELEMENT TECHNIQUE IN TORSION PROBLEMS OF BEAMS WITH MULTIPLY CONNECTED CROSS-SECTIONS Y. M. TEMIS AND V. V. KARABAN 1 Abstract. This paper shows how boundary
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 informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More information26.4. Basic Complex Integration. Introduction. Prerequisites. Learning Outcomes
Basic omplex Integration 6.4 Introduction omplex variable techniques have been used in a wide variety of areas of engineering. This has been particularly true in areas such as electromagnetic field theory,
More informationPhysics 9 WS M5 (rev. 1.0) Page 1
Physics 9 WS M5 (rev. 1.0) Page 1 M-3. Faraday s Law Questions for discussion 1. In the figure below, there is a non-uniform magnetic field pointing into the page. a) If you move the metal loop to the
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. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationLecture 12: Finite Elements
Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature
More informationi.e. the conservation of mass, the conservation of linear momentum, the conservation of energy.
04/04/2017 LECTURE 33 Geometric Interpretation of Stream Function: In the last class, you came to know about the different types of boundary conditions that needs to be applied to solve the governing equations
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
More information