Lesson 9: Diffusion of Heat (discrete rod) and ode45
|
|
- Bathsheba Norris
- 6 years ago
- Views:
Transcription
1 Lesson 9: Diffusion of Heat (discrete rod) and ode Applied Problem. Consider the temperature in a thin rod such as a wire with electrical current. Assume the ends have a fixed temperature and the current generates heat in the rod. Then the temperature will not be uniform and must depend both on space and time, u(x,t). The problem is to approximate the temperature given the initial temperature and the temperature at the ends of the rod. Partitioning the rod into a finite number of segments can form a system of ordinary differential equations, whose solution can be approximated using the Matlab command ode Differential Equation Model. The model for the diffusion of heat in space is based on empirical observations. Discrete Fourier Heat Law: (a). heat flows from hot to cold, (b). the change in heat is proportional to the cross sectional area, change in time and (change in temperature)/(change space normal to the cross section). The last term is a good approximation provided the change in space is small. The proportionality constant, K, is called the thermal conductivity. The K varies with the particular material and with the temperature. Here we will assume the temperature varies over a smaller range so that K is approximately a constant.
2 2 Consider a thin wire so that there is diffusion in only one direction, x. The wire will have a current going through it with a source of heat, f, which is from the electrical resistance of the wire. The f has units of (heat)/(volume*time). We will assume the ends of the wire are kept at zero temperature, and the initial temperature is given. The goal is to be able to predict the temperature inside the wire for any future time and space location. In order to develop a model to predict the temperature = u(x,t), we will discretize space and let u(ih,t) be approximated by u i where h = L/n and L is the length of the wire. The model will have the general form change in heat content (heat from the source) + (heat from diffusion through the right end) + (heat from diffusion through the left end). This is depicted in the figure below where the time step has not been indicated. For time we can choose either t or t+ t. u i-1 * u u i i+1 * * h A h is the change in length A is the cross section area u i is the approximate temperature. Figure: Heat Diffusion in a Thin Wire The heat diffusing in the right face (when (u i+1 - u i ) /h > 0) is A t K (u i+1 - u i ) /h.
3 3 The heat diffusing out the left face (when (u i - u i-1 ) /h > 0) is A t K (u i - u i-1 ) /h. Therefore, the heat from diffusion is A t K (u i+1 - u i ) /h - A t K (ui - u i-1 ) /h. The heat from the source is Ah t f. The heat content of the volume Ah at time t is ρc u i Ah where ρ is the density and c is the specific heat. By combining these we have the following approximation of the change in the heat for the small volume Ah: ρc u i (t+ t) Ah - ρc u i (t) Ah = Ah t f + A t K (u i+1 - u i )/h - A t K (u i - u i-1 )/h. Now, divide by ρcah, define α = (K/ρc) /h 2 and let t go the zero to the get the system of ordinary differential equations where the left side gives the time derivative, u i '. ODE Model for Heat Diffusion. u i ' = 1/ρc f + α(u i+1 + ui-1 ) - 2α u i where (1) i = 1,...,n-1, u i (0) = 0 for i = 1,...,n-1 and (2) u 0 (t) = u n (t) = 0 for t > 0. (3) Equation (2) is the initial temperature set equal to zero, and (3) is the temperature at the left and right ends set equal to zero.
4 4 Equation (1) may be put into the matrix version of a system of ODEs. For example, if the wire is divided into four equal parts, then n = 4 and (1) may be written as three scalar equations u 1 ' = 1/ρc f + α(u 2 + 0) - 2α u1 u 2 ' = 1ρc f + α(u 1 + u3 ) - 2α u 2 u 3 ' = 1/ρc f + α(0 + u 2 ) - 2α u 3 Or, using a vector equation and a matrix u' + A u = b where ' u 1 u ' u' = u2, u u2, b (1/ ρc) f 1 and A α = = =. ' u 3 u A steady state solution is given by setting the derivative vector equal to the zero vector and solving the algebraic system Au = bwhere ' u 1 0 u 1 ' u' = u 2 0 = and u = u2 is constant vector. ' u 3 0 u 3 Thus, the steady state solution is inv(a)*b or A\b. u -1 = A b. This can be done in Matlab by either Another model uses partial derivatives with respect to time and space. PDE Model for Heat Diffusion. ρcu t - (Ku x ) x = f u(x,0) = given and u(0,t), u(1,t) = given, a partial differential equation, an initial condition, boundary conditions.
5 5 The steady state solution is given by setting u t = 0 so that the temperature depends only on space and - (Ku x ) x = f and u(0), u(1) are given. For example, if K = 0.001, f = 1, u(0) = 0 and u(1) = 0, then u(x) must be a quadratic function of x u(x) = x 2 /2 + a x + b. The boundary conditions imply the choice of a and b so that u(x) = -500x(1 x). 9.3 Method of Solution. In this lesson we will use the Matlab command ode45 to solve our systems of differential equation. This command is a robust implementation for systems of differential equations, which uses a variable step size method and the fourth and fifth order Runge-Kutta method. 9.4 Matlab Implementation. Since there are three unknowns and three differential equations and we wish to use Matlab's ode45 scheme, the file, ypheatode.m, must contain the three right sides of the differential equations where y(1) = temperature in the left segment, y(2) = temperature in the center segment and y(3) = temperature in the right segment. The heatode.m file contains the call to ode45, and the Matlab command plot generates the graphs of the three population groups. The initial temperatures are 1.0 and stored in the uo array in the heatode.m file.
6 6 m-files ypheatode.m function ypheatode = ypheatode(t,u) K =.001; rho = 1.; c = 1.; L = 1.; dx = L/4; const = (K/(rho*c))/(dx*dx); f = [1 1 1]/(rho*c); ypheatode(1) = -const*(2*u(1) - u(2)) + f(1); ypheatode(2) = -const*(-u(1) + 2*u(2) - u(3)) + f(2); ypheatode(3) = -const*(2*u(3) - u(2)) + f(3); ypheatode = [ypheatode(1) ypheatode(2) ypheatode(3)]'; m-files heatode.m. clear; clf K =.001; rho = 1.; c = 1.; L = 1.; dx = L/4; const = (K/(rho*c))/(dx*dx); f = [1 1 1]/(rho*c); uo = [1 1 1]; to = 0; tf = 400; [t u] = ode45('ypheatode',[to tf],uo); maxk = size(t,1); x = 0:.25:1; figure(1) for k = 1:5:maxk plot(x,[0 u(k,:) 0]) title('temperatures at various times') xlabel('position in space') axis([ ]) hold on pause end t(1:5:maxk) A = const*[2-1 0;-1 2-1;0-1 2]; d = [1 1 1]'/(rho*c); sssol = A\d plot(x, [0 sssol' 0]) figure(2) plot(t,u(:,2)) title('temperature in the center') xlabel('time')
7 7 9.5 Numerical Experiments. The first graphical output has multiple curves each associated with a particular time as given in the numerical output times = sssol = Note the last curve approximates the steady state solution found by the command A\d. The second graph is for the temperature at the center segment versus time. This graph indicates the approximation of 125.0, which is the steady state solution at the center node.
8 8 temperatures at various times position in space 140 temperature in the center time
9 9 9.6 Additional Calculations. Consider the ODE heat diffusion model with K = 0.001, rho = 1.0, c = 1.0, L = 1.0 and variable heat source f. (a). (b). (c). (d). State the system of differential equations. Modify the ypheatode.m and heatode.m files. For f = [ ]' execute the heatode.m file. Use the hold on command and repeat (c) using f = [ ]' and f = [ ]' (e). Examine the three sets of three curves and numerical outputs. What happens as the heat source f increases?
Lesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationLesson 4: Population, Taylor and Runge-Kutta Methods
Lesson 4: Population, Taylor and Runge-Kutta Methods 4.1 Applied Problem. Consider a single fish population whose size is given by x(t). The change in the size of the fish population over a given time
More informationFalling Mass and Euler Methods
Lesson 2: Falling Mass and Euler Methods 2.1 Applied Problem. Consider a mass, which is falling through a viscous medium. One example is a falling rock, and a second example is a person in parachute. These
More informationLesson 14: Van der Pol Circuit and ode23s
Lesson 4: Van der Pol Circuit and ode3s 4. Applied Problem. A series LRC circuit when coupled via mutual inductance with a triode circuit can generate a sequence of pulsing currents that have very rapid
More informationLesson 13: Rapid Cooling and Euler-trapezoid
Lesson 13: Rapid Cooling and Euler-trapezoid 13.1 Applied Problem. If an object is being cooled very rapidly, and the model is Newton's law of cooling, then the constant c in u t = c(usur - u) will be
More informationDEPARTMENT OF MATHEMATICS. The University of Queensland. MATH3203 Scientific Computation: Algorithmic Design and Implementation
DEPARTMENT OF MATHEMATICS The University of Queensland MATH303 Scientific Computation: Algorithmic Design and Implementation Section A Boundary Value Problem in dimension Semester, 00 A simplified problem
More informationLesson 13: Rapid Cooling and Euler-trapezoid
Lesson 3: Rapid Cooling and Euler-trapezoid 3. Applied Problem. If an object is being cooled very rapidly, and the model is Newton's law of cooling, then the constant c in u t = c(usur - u) will be large.
More informationThe One Dimensional Heat Equation
The One Dimensional Heat Equation Adam Abrahamsen and David Richards May 22, 2002 Abstract In this document we will study the flow of heat in one dimension through a small thin rod. We will use the derivation
More informationFalling Mass and Improved Euler
Lesson 3: Falling Mass and Improved Euler 3.1 Applied Problem. Consider a mass which is falling through a viscous medium. One example is a falling rock, and a second example is a person in parachute. These
More informationDiffusion Processes. Lectures INF2320 p. 1/72
Diffusion Processes Lectures INF2320 p. 1/72 Lectures INF2320 p. 2/72 Diffusion processes Examples of diffusion processes Heat conduction Heat moves from hot to cold places Diffusive (molecular) transport
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationChapter 10: Partial Differential Equations
1.1: Introduction Chapter 1: Partial Differential Equations Definition: A differential equations whose dependent variable varies with respect to more than one independent variable is called a partial differential
More informationFinite Difference Method
Finite Difference Method for BVP ODEs Dec 3, 2014 1 Recall An ordinary differential equation is accompanied by auxiliary conditions. In the analytical method, these conditions are used to evaluate the
More informationDiffusion - The Heat Equation
Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationThe Heat Equation John K. Hunter February 15, The heat equation on a circle
The Heat Equation John K. Hunter February 15, 007 The heat equation on a circle We consider the diffusion of heat in an insulated circular ring. We let t [0, ) denote time and x T a spatial coordinate
More informationNumerical modeling of rock deformation: 06 FEM Introduction
Numerical modeling of rock deformation: 6 FEM Introduction Stefan Schmalholz schmalholz@erdw.ethz.ch NO E 6 AS 9, Thursday -, NO D Introduction We have learned that geodynamic processes can be described
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationA First Course on Kinetics and Reaction Engineering Example S5.1
Example S5.1 Problem Purpose This example illustrates the use of the MATLAB template file SolvBVDif.m to solve a second order boundary value ordinary differential equation. Problem Statement Solve the
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (35 points) Leslie Leroy Irvin bails out of an airplane at the altitude of 16,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s
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 informationIYGB Mathematical Methods 1
IYGB Mathematical Methods Practice Paper A Time: 3 hours Candidates may use any non programmable, non graphical calculator which does not have the capability of storing data or manipulating algebraic expressions
More informationOver Determined Algebraic Systems: Curve Fitting to Data
Over Determined Algebraic Systems: Curve Fitting to Data Introduction. In this section we investigate algebraic systems that have more equations than unknowns. In general, these will not have a solution
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 informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE201: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationLecture 8: Calculus and Differential Equations
Lecture 8: Calculus and Differential Equations Dr. Mohammed Hawa Electrical Engineering Department University of Jordan EE21: Computer Applications. See Textbook Chapter 9. Numerical Methods MATLAB provides
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationHOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, Maximum score 7.0 pts.
HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, 2010. Maximum score 7.0 pts. Three problems are to be solved in this homework assignment. The
More informationMA/CS 615 Spring 2019 Homework #2
MA/CS 615 Spring 019 Homework # Due before class starts on Feb 1. Late homework will not be given any credit. Collaboration is OK but not encouraged. Indicate on your report whether you have collaborated
More informationJanuary 18, 2008 Steve Gu. Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB,
Introduction to MATLAB January 18, 2008 Steve Gu Reference: Eta Kappa Nu, UCLA Iota Gamma Chapter, Introduction to MATLAB, Part I: Basics MATLAB Environment Getting Help Variables Vectors, Matrices, and
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationPDE: The Method of Characteristics Page 1
PDE: The Method of Characteristics Page y u x (x, y) + x u y (x, y) =, () u(, y) = cos y 2. Solution The partial differential equation given can be rewritten as follows: u(x, y) y, x =, (2) where = / x,
More informationOrdinary Differential Equations (ODE)
Ordinary Differential Equations (ODE) Why study Differential Equations? Many physical phenomena are best formulated mathematically in terms of their rate of change. Motion of a swinging pendulum Bungee-jumper
More informationA First Course on Kinetics and Reaction Engineering Example S6.2
Example S6.2 Problem Purpose This example illustrates the use of the MATLAB template file SolvBVDifS.m to solve a second order boundary value ordinary differential equation with a singularity at the boundary
More informationLecture 1: Newton s Cooling Law, Fourier s Law and Heat Conduction
Lecture 1: Newton s Cooling Law, Fourier s Law and Heat Conduction V.Vuorinen Aalto University School of Engineering Heat and Mass Transfer Course, Autumn 2016 October 31 st 2017, Otaniemi ville.vuorinen@aalto.fi
More informationNumerical methods. Examples with solution
Numerical methods Examples with solution CONTENTS Contents. Nonlinear Equations 3 The bisection method............................ 4 Newton s method.............................. 8. Linear Systems LU-factorization..............................
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationIntegration by Parts Logarithms and More Riemann Sums!
Integration by Parts Logarithms and More Riemann Sums! James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 16, 2013 Outline 1 IbyP with
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More informationProject Two. Outline. James K. Peterson. March 27, Cooling Models. Estimating the Cooling Rate k. Typical Cooling Project Matlab Session
Project Two James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 27, 2018 Outline Cooling Models Estimating the Cooling Rate k Typical Cooling
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs 2e: Section 3 - Exercise Page of 6 Exercise Carefully derive the equation of a string in a medium in which the resistance is proportional to the velocity Solution There are two ways (among
More informationMath 5587 Lecture 2. Jeff Calder. August 31, Initial/boundary conditions and well-posedness
Math 5587 Lecture 2 Jeff Calder August 31, 2016 1 Initial/boundary conditions and well-posedness 1.1 ODE vs PDE Recall that the general solutions of ODEs involve a number of arbitrary constants. Example
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More information2 Scalar Ordinary Differential Equation (ODE s)
Dr S. Roberts/Dr L. Stals 2006 Department of Mathematics Australian National University APAC-ANU Teaching Module Scilab ODE s Scilab Numerical Methods Tutorials 1 Introduction This tutorial is designed
More informationThe One-Dimensional Heat Equation
The One-Dimensional Heat Equation R. C. Trinity University Partial Differential Equations February 24, 2015 Introduction The heat equation Goal: Model heat (thermal energy) flow in a one-dimensional object
More informationBoundary conditions. Diffusion 2: Boundary conditions, long time behavior
Boundary conditions In a domain Ω one has to add boundary conditions to the heat (or diffusion) equation: 1. u(x, t) = φ for x Ω. Temperature given at the boundary. Also density given at the boundary.
More informationChapter 12 Partial Differential Equations
Chapter 12 Partial Differential Equations Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 12.1 Basic Concepts of PDEs Partial Differential Equation A
More informationLesson 11: Mass-Spring, Resonance and ode45
Lesson 11: Mass-Spring, Resonance and ode45 11.1 Applied Problem. Trucks and cars have springs and shock absorbers to make a comfortable and safe ride. Without good shock absorbers, the truck or car will
More information1.5 Application Linear Equations and Temperature Oscillations
1.5 Application Linear Equations and Temperature Oscillations Section 1.5 of the text describes the following 4-step algorithm for solving the linear first-order differential equation 1. Begin by calculating
More informationA First Course on Kinetics and Reaction Engineering Example 23.1
Example 23.1 parameter. Problem Purpose This problem illustrates the transient analysis of a CSTR following a change in an operating Problem Statement Recall the isothermal 4430 cm 3 steady-state chemostat
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationINTRODUCTION TO COMPUTER METHODS FOR O.D.E.
INTRODUCTION TO COMPUTER METHODS FOR O.D.E. 0. Introduction. The goal of this handout is to introduce some of the ideas behind the basic computer algorithms to approximate solutions to differential equations.
More informationCHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN
CHAPTER 6 STATE SPACE: FREQUENCY RESPONSE, TIME DOMAIN 6. Introduction Frequency Response This chapter will begin with the state space form of the equations of motion. We will use Laplace transforms to
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationLecture 8: Differential Equations. Philip Moriarty,
Lecture 8: Differential Equations Philip Moriarty, philip.moriarty@nottingham.ac.uk NB Notes based heavily on lecture slides prepared by DE Rourke for the F32SMS module, 2006 8.1 Overview In this final
More informationFinite Volume Method for Scalar Advection in 2D
Chapter 9 Finite Volume Method for Scalar Advection in D 9. Introduction The purpose of this exercise is to code a program to integrate the scalar advection equation in two-dimensional flows. 9. The D
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 informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationA First Course on Kinetics and Reaction Engineering Example 26.3
Example 26.3 unit. Problem Purpose This problem will help you determine whether you have mastered the learning objectives for this Problem Statement A perfectly insulated tubular reactor with a diameter
More informationLaboratory 1. Solving differential equations with nonzero initial conditions
Laboratory 1 Solving differential equations with nonzero initial conditions 1. Purpose of the exercise: - learning symbolic and numerical methods of differential equations solving with MATLAB - using Simulink
More informationNumerically Solving Partial Differential Equations
Numerically Solving Partial Differential Equations Michael Lavell Department of Applied Mathematics and Statistics Abstract The physics describing the fundamental principles of fluid dynamics can be written
More informationBake, shake or break - and other applications for the FEM. 5: Do real-life experimentation using your FEM code
Bake, shake or break - and other applications for the FEM Programming project in TMA4220 - part 2 by Kjetil André Johannessen TMA4220 - Numerical solution of partial differential equations using the finite
More informationSolution. Your sketch should show both x and y axes marked from 5 to 5 and circles of radius 1, 3, and 5 centered at the origin.
Solutions of the Sample Problems for the First In-Class Exam Math 246, Fall 208, Professor David Levermore () (a) Sketch the graph that would be produced by the following Matlab command. fplot(@(t) 2/t,
More informationUse of Differential Equations In Modeling and Simulation of CSTR
Use of Differential Equations In Modeling and Simulation of CSTR JIRI VOJTESEK, PETR DOSTAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nám. T. G. Masaryka
More information1. Differential Equations (ODE and PDE)
1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable
More information1 Introduction to MATLAB
L3 - December 015 Solving PDEs numerically (Reports due Thursday Dec 3rd, carolinemuller13@gmail.com) In this project, we will see various methods for solving Partial Differential Equations (PDEs) using
More informationPH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract
PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed (110108461) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary
More informationEGR 111 Heat Transfer
EGR 111 Heat Transfer The purpose of this lab is to use MATLAB to determine the heat transfer in a 1-dimensional system. New MATLAB commands: (none) 1. Heat Transfer 101 Heat Transfer is thermal energy
More informationMarching on the BL equations
Marching on the BL equations Harvey S. H. Lam February 10, 2004 Abstract The assumption is that you are competent in Matlab or Mathematica. White s 4-7 starting on page 275 shows us that all generic boundary
More informationA Hypothesis about Infinite Series, of the minimally centered variety. By Sidharth Ghoshal May 17, 2012
A Hypothesis about Infinite Series, of the minimally centered variety By Sidharth Ghoshal May 17, 2012 Contents Taylor s Theorem and Maclaurin s Special Case: A Brief Introduction... 3 The Curious Case
More informationSolutions of separable equations often are in an implicit form; that is, not solved for the dependent variable.
Section 2.2 Separable DEs Key Terms: Separable DEs; the DE can be rearranged so that all terms containing the independent variable are on one side of the equal sign and all terms containing the dependent
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 informationThe Heat Equation with partial differential equations
Applications of Differential Equations 1/100 The Heat Equation with partial differential equations David Richards and Adam Abrahamsen Introduction The purpose of our presentation is to derive the heat
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationMA 102 Mathematics II Lecture Feb, 2015
MA 102 Mathematics II Lecture 1 20 Feb, 2015 Differential Equations An equation containing derivatives is called a differential equation. The origin of differential equations Many of the laws of nature
More informationHOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, Maximum score 6.0 pts.
HOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, 2012. Maximum score 6.0 pts. Read Strang s (old) book, Sections 3.1-3.4 and 6.4-6.5.
More informationLecture 7. Heat Conduction. BENG 221: Mathematical Methods in Bioengineering. References
BENG 221: Mathematical Methods in Bioengineering Lecture 7 Heat Conduction References Haberman APDE, Ch. 1. http://en.wikipedia.org/wiki/heat_equation Lecture 5 BENG 221 M. Intaglietta Heat conduction
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationPartial Differential Equations - part of EM Waves module (PHY2065)
Partial Differential Equations - part of EM Waves module (PHY2065) Richard Sear February 7, 2013 Recommended textbooks 1. Mathematical methods in the physical sciences, Mary Boas. 2. Essential mathematical
More informationAnalysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017
Analysis III (BAUG Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017 Question 1 et a 0,..., a n be constants. Consider the function. Show that a 0 = 1 0 φ(xdx. φ(x = a 0 + Since the integral
More informationAbsorbing Boundary Conditions for Nonlinear Wave Equations
Absorbing Boundary Conditions for Nonlinear Wave Equations R. B. Gibson December, 2010 1 1 Linear Wave Equations 1.1 Boundary Condition 1 To begin, we wish to solve where with the boundary condition u
More information12.7 Heat Equation: Modeling Very Long Bars.
568 CHAP. Partial Differential Equations (PDEs).7 Heat Equation: Modeling Very Long Bars. Solution by Fourier Integrals and Transforms Our discussion of the heat equation () u t c u x in the last section
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationCOMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE
COMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE Rafael ÅMAN*, Heikki HANDROOS*, Pasi KORKEALAAKSO** and Asko ROUVINEN** * Laboratory
More informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep Copyright
More informationAPPENDIX 1 MATLAB AND ANSYS PROGRAMS
APPENDIX 1 MATLAB AND ANSYS PROGRAMS This appendix lists all the MATLAB and ANSYS codes used in each chapter, along with a short description of the purpose of each. MATLAB codes have the suffix.m and the
More informationLecture 5: Travelling Waves
Computational Biology Group (CoBI), D-BSSE, ETHZ Lecture 5: Travelling Waves Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015 26. Oktober 2016 2 / 68 Contents 1 Introduction to Travelling Waves
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 informationEffect of Variable Viscosity on Hydro Magnetic Flow and Heat Transfer Over a Stretching Surface with Variable Temperature
37 Effect of Variable Viscosity on Hydro Magnetic Flow and Heat Transfer Over a Stretching Surface with Variable Temperature M. Y. Akl Department of Basic Science, Faculty of Engineering (Shopra Branch),
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationProject Two. James K. Peterson. March 26, Department of Biological Sciences and Department of Mathematical Sciences Clemson University
Project Two James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University March 26, 2019 Outline 1 Cooling Models 2 Estimating the Cooling Rate k 3 Typical
More informationSolve the set of equations using the decomposition method and AX=B:
PROBLEM 1.1 STATEMENT Solve the set of equations using the decomposition method and AX=B: 5x 1 + 3x 2 + 4x 3 = 12 6x 1 + 3x 2 + 4x 3 = 15 7x 1 + 9x 2 + 2x 3 = 10 PROBLEM 1.2 STATEMENT 1 4 5 3 4 5 If A
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 information