Solution of Differential Equation by Finite Difference Method
|
|
- Adam Phelps
- 5 years ago
- Views:
Transcription
1 NUMERICAL ANALYSIS University of Babylon/ College of Engineering/ Mechanical Engineering Dep. Lecturer : Dr. Rafel Hekmat Class : 3 rd B.Sc Solution of Differential Equation by Finite Difference Method What is a Differential Equation? An equation that consists of derivatives is called a differential equation. Differential equations have applications in all areas of science and engineering. Mathematical formulation of most of the physical and engineering problems leads to differential equations. Types of Differential Equations a- Ordinary Differential Equations (ODE) b- Partial Differential Equations (PDE) An ordinary differential equation is that in which all the derivatives are with respect to a single independent variable. Examples of ordinary differential equations include Ordinary differential equations are classified in terms of order and degree. Order of an ordinary differential equation is the same as the highest derivative and the degree of an ordinary differential equation is the power of highest derivative. Thus the differential equation, is of order 3 and degree 1, whereas the differential equation is of order 1 and degree. 1
2 PDEs are differential equations containing two or more independent variables, so the derivatives are partial derivatives. For example,given a function of two variables (x, y) f(x, y), the partial derivatives with respect to x and y are f x and f y PDE example I: Laplace equation The Laplace equation is a second ordered PDE appearing for example in Fluid Mechanics (potential flows) and Electromagnetics (electromagnetic field in a charge free region): where represents, for example, the velocity potential or the temperature. PDE example II: Heat equation The heat equation is a second order PDE describing how temperature T diffuses through a medium of thermal diffusivity α: PDE example III: Wave equation The wave equation describes how the sound pressure p propagates at speed c through a medium (at rest), and plays a role in acoustics, fluid mechanics, and quantum mechanics: Classification of Partial Differential Equations For analyzing the equations for fluid flow problems, it is convenient to consider the case of a second-order differential equation given in the general form as: In the coefficients A, B, C, D, E and F are either constants or functions of only (x, y) (do not contain φ or its derivatives), it is said to be a linear equation, otherwise it is a non-linear equation. If G =, the aforesaid equation is homogeneous, otherwise it is non-homogeneous.
3 Again for the above mentioned equation if B 4AC =, the equation is parabolic if B 4AC <, the equation is elliptic if B 4AC >, the equation is hyperbolic Example: Classify the following linear second order partial differential equations (PDEs) with solution u(x,y)in the xy-plane. 3
4 Example: Consider the one-dimensional damped wave equation 9uxx= utt+ 6ut Initial and Boundary Conditions Most PDEs have an infinite number of admissible solutions. Thus the PDE alone is not sufficient to get a unique solution. Usually some boundary conditions and initial conditions are required. For the heat equation the simplest boundary conditions are fixed temperatures at both ends: Where l is the length of the rod and h1(t) the temperature at the first end and h(t) the temperature at the second end. The initial conditions specify an arbitrary initial temperature distribution inside the rod: What is the finite difference method? The finite difference method is used to solve ordinary differential equations that have conditions imposed on the boundary rather than at the initial point. These problems are called boundary-value problems. Steps of finite difference solution: Divide the solution region into a grid of nodes, Approximate the given differential equation by finite difference equivalent, 4
5 Solve the differential equations subject to the boundary conditions and/or initial conditions. First Derivative Approximations The derivative of a given function f(x) can be approximated as: df ( x ) f ( x x) f ( x ) dx x df ( x ) f ( x ) f ( x x) dx x (1) () df x f x x f x x ( ) ( ) ( ) (3) dx x Forward Difference Formula Backword Difference Formula Central Difference Formula Second Derivative Approximation f f ( x x) f ( x ) f ( x x) x ( x) Example: Solve f 4 f, x 1. x Subject to the boundary conditions: Assume x=.5 f () f (1) 1 Solution: Discretization of the equation: f ( i 1) f ( i) f ( i 1) 4 ( ) f i ( x) f () f (1) f () 4 f (1) (.5) Since there is one unknown, one equation is enough to find this unknown. 5
6 If x=.5 f (1) f () f () f (1) 1 f () f (1) f () 4 f (1) (.5) f (3) f () f (1) 4 f () (.5) f (4) f (3) f () 4 f (3) (.5) f(1), f() and f(3) are unknowns. Therefore, three equations are written to find three unknowns. Example: Consider the diffusion equation: k t x Where k is a constant. Discretized equation: Where ( i, j 1) ( i, j) ( i 1, j) ( i, j) ( i 1, j) k t ( x) x i x, i,1,,..., n t j t, j,1,,..., n Forward difference for t and central difference for x is used. t r k( x) We can write: ( i, j 1) r ( i 1, j) (1 r) ( i, j) r ( i 1, j) This equation can be calculated in terms of boundary and initial conditions. This explicit formula can be used to find (x,t+ )explicitly in terms of (x,t) Example: The deflection in a simply supported beam with a uniform load q and a tensile axial load T is given by ` where x:location along the beam (in) T:Tension applied (lbs) E:Young s modulus of elasticity of the beam (psi) 6
7 I:Second moment of area (in4) q:uniform loading intensity (lb/in) L:Length of beam (in) Given, Find the deflection of the beam at x=5". Use a step size of Δx=5 and approximate the derivatives by central divided difference approximation. Solution: Substituting the given values, Approximating the derivative y at node i by the central divided difference approximation x We can rewrite the equation as: Since Δx=5, we have 4 nodes as given in Figure below Writing the equation at each node, we get 7
8 Equations above are 4 simultaneous equations with 4 unknowns and can be written in matrix form as Solving the equations we get, Example Take the case of a pressure vessel that is being tested in the laboratory to check its ability to withstand pressure. For a thick pressure vessel of inner radius a and outer radius b, the differential equation for the radial displacement u of a point along the thickness is given by The inner radius a= 5, and the outer radius b= 8, and the material of the pressure vessel is ASTM A36 steel Divide the radial thickness of the pressure vessel into 6 equidistant nodes, and find the radial displacement profile At node i along the radial thickness of the pressure vessel 8
9 Let us break the thickness, b-a, of the pressure vessel into n+1 nodes, that is is node and is node. That r=a is nodes i= and r=b is node i=n. That means we have n+1 unknowns. We can write the above equation for nodes 1,., n-1. This will give us n-1 equations. At the edge nodes, i= and i=n, we use the boundary conditions of We have been asked to do the calculations for n=5, that is a total of 6 nodes. This gives 9
10 Writing the above Equations in matrix form gives After solving Explicit versus implicit Finite Difference Schemes During the lecture we solved the transient (time-dependent) heat equation in 1D In explicit finite difference schemes, the temperature at time n+1 depends explicitly on the temperature at time n. The explicit finite difference discretization of above equation is This can rearranged in the following manner (with all quantities at time n+1 on the left-handside and quantities at time n on the right-hand-side) Since we knowt n n i+1, T i and T n i 1, we can compute T n i. The major advantage of explicit finite difference methods is that they are relatively simple and computationally fast. However, the main drawback is that stable solutions are obtained only when In implicit finite difference schemes, the spatial derivatives T are evaluated (at least partially) x at the new time step. The simplest implicit discretization of heat equation in 1D is This can be rearranged so that unknown terms are on the left and known terms are on the right Where s = K t. The main advantage of implicit finite difference methods is that there are x 1
11 no restrictions on the time step, which is good news if we want to simulate geological processes at high spatial resolution. Taking large time steps, however, may result in an inaccurate solution. Solving an implicit finite difference scheme We solve the transient heat equation on the domain l/ x l/ with the following boundary conditions with the initial condition As usual, the first step is to discretize the spatial domain with nx finite difference points. The implicit finite difference discretization of the temperature equation is and the known right-hand-side vector rhs is 11
Classification 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 informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More information1D Wave PDE. Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) Richard Sear.
Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) November 12, 2018 Wave equation in one dimension This lecture Wave PDE in 1D Method of Separation of
More informationChapter 5 Types of Governing Equations. Chapter 5: Governing Equations
Chapter 5 Types of Governing Equations Types of Governing Equations (1) Physical Classification-1 Equilibrium problems: (1) They are problems in which a solution of a given PDE is desired in a closed domain
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture1/1_1.htm 1 of 1 6/19/2012 4:29 PM The Lecture deals with: Classification of Partial Differential Equations Boundary and Initial Conditions Finite Differences
More informationOrdinary Differential Equations
Chapter 10 Ordinary Differential Equations 10.1 Introduction Relationship between rate of change of variables rather than variables themselves gives rise to differential equations. Mathematical formulation
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationTHE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))
THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the
More informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationVideo 15: PDE Classification: Elliptic, Parabolic and Hyperbolic March Equations 11, / 20
Video 15: : Elliptic, Parabolic and Hyperbolic Equations March 11, 2015 Video 15: : Elliptic, Parabolic and Hyperbolic March Equations 11, 2015 1 / 20 Table of contents 1 Video 15: : Elliptic, Parabolic
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
More informationElements of Matlab and Simulink Lecture 7
Elements of Matlab and Simulink Lecture 7 Emanuele Ruffaldi 12th May 2009 Copyright 2009,Emanuele Ruffaldi. This work is licensed under the Creative Commons Attribution-ShareAlike License. PARTIAL DIFFERENTIAL
More informationAn Introduction to Partial Differential Equations
An Introduction to Partial Differential Equations Ryan C. Trinity University Partial Differential Equations Lecture 1 Ordinary differential equations (ODEs) These are equations of the form where: F(x,y,y,y,y,...)
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
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 informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationSolutions Definition 2: a solution
Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential
More information08.06 Shooting Method for Ordinary Differential Equations
8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method
More informationCIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass
CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional
More informationBEAM DEFLECTION THE ELASTIC CURVE
BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment
More informationStress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy
Stress Analysis Lecture 3 ME 276 Spring 2017-2018 Dr./ Ahmed Mohamed Nagib Elmekawy Axial Stress 2 Beam under the action of two tensile forces 3 Beam under the action of two tensile forces 4 Shear Stress
More informationIntro to Research Computing with Python: Partial Differential Equations
Intro to Research Computing with Python: Partial Differential Equations Erik Spence SciNet HPC Consortium 28 November 2013 Erik Spence (SciNet HPC Consortium) PDEs 28 November 2013 1 / 23 Today s class
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationNumerical Solutions of Partial Differential Equations
Numerical Solutions of Partial Differential Equations Dr. Xiaozhou Li xiaozhouli@uestc.edu.cn School of Mathematical Sciences University of Electronic Science and Technology of China Introduction Overview
More informationAnalysis III for D-BAUG, Fall 2017 Lecture 11
Analysis III for D-BAUG, Fall 2017 Lecture 11 Lecturer: Alex Sisto (sisto@math.ethz.ch) 1 Beams (Balken) Beams are basic mechanical systems that appear over and over in civil engineering. We consider a
More informationSolution: The moment of inertia for the cross-section is: ANS: ANS: Problem 15.6 The material of the beam in Problem
Problem 15.4 The beam consists of material with modulus of elasticity E 14x10 6 psi and is subjected to couples M 150, 000 in lb at its ends. (a) What is the resulting radius of curvature of the neutral
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationOrdinary Differential Equations
Ordinary Differential Equations (MA102 Mathematics II) Shyamashree Upadhyay IIT Guwahati Shyamashree Upadhyay ( IIT Guwahati ) Ordinary Differential Equations 1 / 1 Books Shyamashree Upadhyay ( IIT Guwahati
More informationGEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE
GEOSYNTHETICS ENGINEERING: IN THEORY AND PRACTICE Prof. J. N. Mandal Department of Civil Engineering, IIT Bombay, Powai, Mumbai 400076, India. Tel.022-25767328 email: cejnm@civil.iitb.ac.in Module-13 LECTURE-
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationDeflection profile analysis of beams on two-parameter elastic subgrade
1(213) 263 282 Deflection profile analysis of beams on two-parameter elastic subgrade Abstract A procedure involving spectral Galerkin and integral transformation methods has been developed and applied
More information1.1 Motivation: Why study differential equations?
Chapter 1 Introduction Contents 1.1 Motivation: Why stu differential equations?....... 1 1.2 Basics............................... 2 1.3 Growth and decay........................ 3 1.4 Introduction to Ordinary
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 informationSOLUTION OF EQUATION AND EIGENVALUE PROBLEMS PART A( 2 MARKS)
CHENDU COLLEGE OF ENGINEERING AND TECHNOLOGY (Approved by AICTE New Delhi, Affiliated to Anna University Chennai. Zamin Endathur Village, Madurntakam Taluk, Kancheepuram Dist.-603311.) MA6459 - NUMERICAL
More informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More information6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and
6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa (15.5 10 6 psi) and an original diameter of 3.8 mm (0.15 in.) will experience only elastic deformation when a tensile
More information! 1.1 Definitions and Terminology
! 1.1 Definitions and Terminology 1. Introduction: At times, mathematics aims to describe a physical phenomenon (take the population of bacteria in a petri dish for example). We want to find a function
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationMathematical Computing
IMT2b2β Department of Mathematics University of Ruhuna A.W.L. Pubudu Thilan Differential Equations Types of Differential Equations Differential equations can basically be classified as ordinary differential
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
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 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 informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More information4.MECHANICAL PROPERTIES OF MATERIALS
4.MECHANICAL PROPERTIES OF MATERIALS The diagram representing the relation between stress and strain in a given material is an important characteristic of the material. To obtain the stress-strain diagram
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 13 Saarland University 5. Januar 2015 c Daria Apushkinskaya (UdS) PDE and BVP lecture 13 5. Januar 2015 1 / 35 Purpose of Lesson To interpretate
More informationDefinition of differential equations and their classification. Methods of solution of first-order differential equations
Introduction to differential equations: overview Definition of differential equations and their classification Solutions of differential equations Initial value problems Existence and uniqueness Mathematical
More informationPartial Differential Equations (PDEs)
C H A P T E R Partial Differential Equations (PDEs) 5 A PDE is an equation that contains one or more partial derivatives of an unknown function that depends on at least two variables. Usually one of these
More information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationSolutions to Laplace s Equations- II
Solutions to Laplace s Equations- II Lecture 15: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Laplace s Equation in Spherical Coordinates : In spherical coordinates
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion
EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 3 Torsion Introduction Stress and strain in components subjected to torque T Circular Cross-section shape Material Shaft design Non-circular
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
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 informationComputational Neuroscience. Session 1-2
Computational Neuroscience. Session 1-2 Dr. Marco A Roque Sol 05/29/2018 Definitions Differential Equations A differential equation is any equation which contains derivatives, either ordinary or partial
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationBefore we look at numerical methods, it is important to understand the types of equations we will be dealing with.
Chapter 1. Partial Differential Equations (PDEs) Required Readings: Chapter of Tannehill et al (text book) Chapter 1 of Lapidus and Pinder (Numerical Solution of Partial Differential Equations in Science
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
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 informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
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 informationLecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018
Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationModule 4 : Deflection of Structures Lecture 4 : Strain Energy Method
Module 4 : Deflection of Structures Lecture 4 : Strain Energy Method Objectives In this course you will learn the following Deflection by strain energy method. Evaluation of strain energy in member under
More informationLECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD
ECTURE 19: SEPARATION OF VARIABES, HEAT CONDUCTION IN A ROD The idea of separation of variables is simple: in order to solve a partial differential equation in u(x, t), we ask, is it possible to find a
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationtwenty one concrete construction: shear & deflection ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2014 lecture
ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2014 lecture twenty one concrete construction: Copyright Kirk Martini shear & deflection Concrete Shear 1 Shear in Concrete
More informationPartial Differential Equations
Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but
More informationLecture 2. Introduction to Differential Equations. Roman Kitsela. October 1, Roman Kitsela Lecture 2 October 1, / 25
Lecture 2 Introduction to Differential Equations Roman Kitsela October 1, 2018 Roman Kitsela Lecture 2 October 1, 2018 1 / 25 Quick announcements URL for the class website: http://www.math.ucsd.edu/~rkitsela/20d/
More informationMathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 2 The wave equation Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.0 28/09/2018 1 Learning objectives of this lecture Understand the fundamental properties of the wave equation
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
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 informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationSTRENGTH OF MATERIALS-I. Unit-1. Simple stresses and strains
STRENGTH OF MATERIALS-I Unit-1 Simple stresses and strains 1. What is the Principle of surveying 2. Define Magnetic, True & Arbitrary Meridians. 3. Mention different types of chains 4. Differentiate between
More information