PowerPoints organized by Dr. Michael R. Gustafson II, Duke University
|
|
- Ethelbert Griffin
- 6 years ago
- Views:
Transcription
1 Part 6 Chapter Boundary-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
2 Chapter Objectives Understanding the difference between initial-value and boundary-value problems. Knowing how to express an n th order ODE as a system of n first-order ODEs. Knowing how to implement the shooting method for linear ODEs by using linear interpolation to generate accurate shots. Understanding how derivative boundary conditions are incorporated into the shooting method.
3 Objectives (cont) Knowing how to solve nonlinear ODEs with the shooting method by using root location to generate accurate shots. Knowing how to implement the finite-difference method. Understanding how derivative boundary conditions are incorporated into the finite-difference method. Knowing how to solve nonlinear ODEs with the finite-difference method by using root location methods for systems of nonlinear algebraic equations.
4 Boundary-Value Problems Boundary-value problems are those where conditions are not known at a single point but rather are given at different values of the independent variable. Boundary conditions may include values for the variable or values for derivatives of the variable. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
5 Higher Order Systems MATLAB s ODE solvers are based on solving first-order differential equations only. To solve an n th order system (n>1), the system must be written as n first-order equations: dt z dx dz ht dx dt h T T 0 dx T Each first-order equation needs an initial value or boundary value to solve.
6 The Shooting Method One method for solving boundary-value problems - the shooting method - is based on converting the boundary-value problem into an equivalent initial-value problem. Generally, the equivalent system will not have sufficient initial conditions and so a guess is made for any undefined values. The guesses are changed until the final solution satisfies all the boundary conditions. For linear ODEs, only two shots are required - the proper initial condition can be obtained as a linear interpolation of the two guesses. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
7 Example Steady-state temperature distribution for a long, thin rod positioned between two constant temperature walls. Solve T a dt h T T 0 dx x=0 L with L=10 m, h =0.05 m -, T =00 K, T(0) = 300 K, and T(10) = 400 K. Assume no radial temperature gradient. First - break into two equations: convection conduction T T b dt z dx dz T dx dt h T T 0 dx
8 Example Code Code for derivatives: function dy=ex30(x,y) %linear ODE boundary-value problem heat conduction/convection p. 547 Chapra dy=[y();-0.05*(00-y(1))]; Code for 1 st shot try z(0)=z a1 =-5 [x,y]=ode45(@ex30, [0 10], [300-5]); Tb1=y(length(y)) >>ans Tb1 = Code for nd shot try z(0)=z a =-0 [x,y]=ode45(@ex30, [0 10], [300-0]); Tb=y(length(y)) >>ans Tb =
9 Interpolation Because it s a linear ODE, use linear interpolation to solve for the exact z(0)=z a ( z z ) ( z z ) ( T T ) z z ( z z ) ( T T ) ( T T ) ( T T ) z a a a1 a a1 b b1 a a1 a a1 b b1 b b1 b b1 ( ) 5 ( 0 ( 5)) K/m ( ) Now solve for the T(x) distribution: [t,y]=ode45(@ex30,[0 10],[ ]);
10 Boundary Conditions Dirichlet boundary conditions are those where a fixed value of a variable is known at a particular location. Neumann boundary conditions are those where a derivative is known at a particular location. Shooting methods can be used for either kind of boundary condition.
11 Example Shooting Method with derivative boundary conditions (p. 549) Same as previous, except convection is on the left side. Solve: dt h T T 0 dx T convection conduction h=1 J/(m K s), k = 00 J/ (s m K) First - break into two equations: dt h T T 0 dx convection T = 00 K x=0 L a b dt z dx dz T dx T b
12 Boundary condition at a: Solve for gradient =z a : Guess T a1 = 300K: function dy=ex30(x,y) %linear ODE boundary-value problem heat conduction/convection %p. 547 Chapra dy=[y();-0.05*(00-y(1))]; end dt hac( T T (0)) kac (0) dx dt h (0) ( T (0) T ) dx k dt 1 1 K za 1 (0) ( T (0) T ) (300 00) 0.5 dx m First guess >> [x,y]=ode45(@ex30,[0 10],[300,0.5]); >> Tb1=y(length(y)) Tb1 = Result is higher than known T b = 400 K, so make a new guess.
13 nd guess T a = 150 K, z a = -0.5 [x,y]=ode45(@ex30,[0 10],[150,-0.5]); >> Tb=y(length(y)) Tb = First guess : Ta1 = 300K, za1 =0.5K/m Tb1= K Second guess: Ta= 150K, za = -0.5 Tb= K So we know that Ta is in between 300 and 150K. Using linear interpolation, ( T T ) ( T T ) ( T T ) T T ( T T ) ( T T ) ( T T ) ( T T ) T a a a1 a a1 a a1 a a1 b b1 b b1 b b1 b b ( ) K And z a = K/m.
14 Difference between Dirichlet and Nuemann BC problems Difference in the last two problems was that when the derivative B.C. is used, we don t know Ta. Therefore, must guess Ta and solve for za, and then run. In the first problem, Ta is known, so we guess za only.
15 The Shooting Method for Nonlinear ODEs For nonlinear ODEs, interpolation between two guesses will not necessarily result in an accurate estimate of the required boundary condition. Instead, the boundary condition can be used to write a roots problem with the estimate as a variable.
16 Example dt dx Solve with =.7x10-9 K -3 m -, L=10 m, h =0.05 m-, T =00 K, T(0) = 300 K, and T(10) = 400 K. d T dx h T T T 4 T 4 0 First - break into two equations: dt z 4 4 dx ht T T T 0 dz dx T T
17 Example Code Code for derivatives: function dy=dydxn(x,y) dy=[y(); -0.05*(00-y(1))-.7e-9*(1.6e9-y(1)^4)]; Code for residual: function r=res(za) [0 10], [300 za]); r=y(length(x),1)-400; Code for finding root of residual: -50)%50 is the guess Code for solving system: [0 10], [ ) ]);
18
19 Finite-Difference Methods The most common alternatives to the shooting method are finite-difference approaches. In these techniques, finite differences are substituted for the derivatives in the original equation, transforming a linear differential equation into a set of simultaneous algebraic equations.
20 Finite-Difference Example Convert: d T dx h T T 0 into n-1 simultaneous equations at each interior point using centered difference equations: dt Ti 1T i Ti 1 dx x Ti 1T i Ti 1 h T Ti 0 x T hx T T hx T i1 i i1 (.16)
21 Finite-Difference Example (cont) Since T 0 and T n are known, they will be on the right-hand-side of the linear algebra system (in this case, in the first and last entries, respectively): h x 1 1 h x 1 1 h x T 1 T T n1 h x T T 0 h x T h x T T n
22 Example of finite-difference method for linear ODE Use the finite-difference approach to solve the same problem as before. Use 4 interior nodes with segment length of x = m. For node 1: -T 0 +.T 1 T = 40 Substituting T 0 = 300 gives Assemble all the equations:.t 1 T = T T T T4 440
23 >> A=[ ; ; ; ]; >> b=[ ]'; >> T=A\b T =
24 Derivative Boundary Conditions Neumann boundary conditions are resolved by solving the centered difference equation at the point and rewriting the system equation accordingly. Use Chapra Tables p For example, if there is a Neumann condition at the T 0 point, dt T1 T 1 dt T 1 T1 x dx 0 x dx 0 T hx T T hx T dt T1 x hx T0 T1 hx T dx 0 dt hx T0 T1 hx T x dx 0 Substitute T -1 into.16 (.18)
25 Example Incorporating Derivative Boundary Conditions Problem: Generate the finite-difference solution for a 10-m rod with x = m, h = 0.05 m -, T = 00 K, and the boundary conditions: T a = 0 and T b = 400 K. Note that the first condition means that the slope of the solution should approach zero at the rod s left end. Aside from this case, also generate the solution for dt/dx = -0 at x =0. Solution. Eq (.18) can be used to represent node 0 as.t 0 - T 1 = 0.05* *00-0=40
26 Write Eq.(.16) for the interior nodes. For example, for node 1, -T 0 +.T 1 T = 40 A similar approach can be used for the remaining interior nodes. The final system of equations can be assembled in matrix form as. T T T T T 440 4
27 >> A=[ ; ; ; ; ]; >> b=[ ] '; T=A\b T =
28 For derivative at x = 0 is set to -0 K/m, the simultaneous equations are Solution: >>T = T T T T T 440 4
29 Plot of Soln to Example.6 (derivative boundary conditions)
30 Finite-Difference Method for Nonlinear ODEs Root location methods for systems of equations may be used to solve nonlinear ODEs. Another method is to adapt a successive substitution algorithm to calculate the values of the interior points. dt h T T T T dx Ti 1T i Ti 1 0 h T T T T x Collecting terms, 4 4 T ( hx ) T T hx T x T T 4 4 i1 i i1 i
31 Solve for the left-side T i in terms of the nonlinear T i on the right side: T i hx T x T T T T 4 4 i i1 i1 hx Assume initial values of T i all = 0. Solve for T 1. Substitute this new value of T 1, and solve for T. Continue solving for all T n. This is one iteration. Now use these new values to complete a second iteration. Continue these iterations until tolerance is met for all T i.
32 Example of Finite-Difference Approach for Nonlinear ODEs Solve the previous convection/radiation problem using finiteelement and successive substitution. with =.7x10-9 K -3 m -, L=10 m, h =0.05 m-, T =00 K, T(0) = 300 K, and T(10) = 400 K. Use four interior nodes with a segment length of x = m. Recall that we solved the same problem with the shooting method in Example.4.
33 T T T T () 00.7 x 10 () (00 0 ) () 00.7 x 10 () (00 0 ) () 00.7 x 10 () (00 0 ) () 0.05() 0.05() () 00.7 x 10 () (00 0 ) () Continue until converging to final result: T T T T T T
Ordinary 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 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 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 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 informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 17 Numerical Integration Formulas PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 4 Chapter 14 General Linear Least-Squares ares and Nonlinear Regression PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationGeneral Linear Least-Squares and Nonlinear Regression
Part 4 Chapter 14 General Linear Least-Squares and Nonlinear Regression PowerPoints organized by Dr. Michael R. Gustafson II, Duke University Revised by Prof. Jang, CAU All images copyright The McGraw-Hill
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 19 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 3 Chapter 10 LU Factorization PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 1 Chapter 4 Roundoff and Truncation Errors PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationInternal Forced Convection. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Internal Forced Convection Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Introduction Pipe circular cross section. Duct noncircular cross section. Tubes small-diameter
More informationParallel Circuits. Chapter
Chapter 5 Parallel Circuits Topics Covered in Chapter 5 5-1: The Applied Voltage V A Is the Same Across Parallel Branches 5-2: Each Branch I Equals V A / R 5-3: Kirchhoff s Current Law (KCL) 5-4: Resistance
More informationIntroduction to Heat and Mass Transfer. Week 8
Introduction to Heat and Mass Transfer Week 8 Next Topic Transient Conduction» Analytical Method Plane Wall Radial Systems Semi-infinite Solid Multidimensional Effects Analytical Method Lumped system analysis
More information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationFault Tolerant Computing CS 530DL
Fault Tolerant Computing CS 530DL Additional Lecture Notes Modeling Yashwant K. Malaiya Colorado State University March 8, 2017 1 Quantitative models Derived from first principles: Arguments are actual
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationHEAT TRANSFER FROM FINNED SURFACES
Fundamentals of Thermal-Fluid Sciences, 3rd Edition Yunus A. Cengel, Robert H. Turner, John M. Cimbala McGraw-Hill, 2008 HEAT TRANSFER FROM FINNED SURFACES Mehmet Kanoglu Copyright The McGraw-Hill Companies,
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 informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationMTF071 Computational Fluid Dynamics of Turbulent
CHALMERS TEKNISKA HÖGSKOLA Termo- och Fluiddynamik 42 96 Göteborg MTF07 Computational Fluid Dynamics of Turbulent Flow http://wwwtfdchalmersse/gr-kurs/mtf07 Task K2 Lars Davidson 2003-02-03 In Task K you
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 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 informationComputational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #11 Fundamentals of Discretization: Finite Difference
More informationAppendix F. + 1 Ma 1. 2 Ma Ma Ma ln + K = 0 (4-173)
5:39p.m. Page:949 Trimsize:8.5in 11in Appendix F F.1 MICROSOFT EXCEL SOLVER FOR NON-LINEAR EQUATIONS The Solver is an optimization package that finds a maximum, minimum, or specified value of a target
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 informationChapter 2 HEAT CONDUCTION EQUATION
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 2 HEAT CONDUCTION EQUATION Mehmet Kanoglu University of Gaziantep
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 informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More 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 informationOrdinary Differential Equations (ode)
Ordinary Differential Equations (ode) Numerical Methods for Solving Initial condition (ic) problems and Boundary value problems (bvp) What is an ODE? =,,...,, yx, dx dx dx dx n n 1 n d y d y d y In general,
More information2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13
2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13 REVIEW Lecture 12: Classification of Partial Differential Equations (PDEs) and eamples with finite difference discretizations Parabolic PDEs Elliptic
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 1 INTRODUCTION AND BASIC CONCEPTS
Heat and Mass Transfer: Fundamentals & Applications 5th Edition in SI Units Yunus A. Çengel, Afshin J. Ghajar McGraw-Hill, 2015 Chapter 1 INTRODUCTION AND BASIC CONCEPTS Mehmet Kanoglu University of Gaziantep
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 informationChapter 5 Time-Dependent Conduction
Chapter 5 Time-Dependent Conduction 5.1 The Lumped Capacitance Method This method assumes spatially uniform solid temperature at any instant during the transient process. It is valid if the temperature
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 informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationCHAPTER 7. An Introduction to Numerical Methods for. Linear and Nonlinear ODE s
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 1 A COLLECTION OF HANDOUTS ON FIRST ORDER ORDINARY DIFFERENTIAL
More informationQuasi-linear first order equations. Consider the nonlinear transport equation
Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +
More informationOrdinary Differential Equations
Ordinary Differential Equations Part 4 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Ordinary Differential Equations p. 1/23 Two-point Boundary Value Problems NRiC 17. BCs specified at two
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 informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
More informationMotion Along a Straight Line
Chapter 2 Motion Along a Straight Line PowerPoint Lectures for University Physics, Twelfth Edition Hugh D. Young and Roger A. Freedman Lectures by James Pazun Copyright 2008 Pearson Education Inc., publishing
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationCS545 Contents XVI. Adaptive Control. Reading Assignment for Next Class. u Model Reference Adaptive Control. u Self-Tuning Regulators
CS545 Contents XVI Adaptive Control u Model Reference Adaptive Control u Self-Tuning Regulators u Linear Regression u Recursive Least Squares u Gradient Descent u Feedback-Error Learning Reading Assignment
More informationMinimizing & Maximizing Functions
Minimizing & Maximizing Functions Nonlinear functions may have zero to many minima and maxima. Example: find the minimum of y = 3x 2 2x + 1 Minima & maxima occur in functions where the slope changes sign
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationLecture 16 : Fully Developed Pipe flow with Constant Wall temperature and Heat Flux
Module 2 : Convection Lecture 16 : Fully Developed Pipe flow with Constant Wall temperature and Heat Flux Objectives In this class: The fully developed temperature profile for uniform circumferential heating
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 informationMethods for Computing Periodic Steady-State Jacob White
Introduction to Simulation - Lecture 15 Methods for Computing Periodic Steady-State Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Periodic Steady-state problems Application
More informationChapter. Triangles. Copyright Cengage Learning. All rights reserved.
Chapter 3 Triangles Copyright Cengage Learning. All rights reserved. 3.5 Inequalities in a Triangle Copyright Cengage Learning. All rights reserved. Inequalities in a Triangle Important inequality relationships
More informationq x = k T 1 T 2 Q = k T 1 T / 12
Conductive oss through a Window Pane q T T 1 Examine the simple one-dimensional conduction problem as heat flow through a windowpane. The window glass thickness,, is 1/8 in. If this is the only window
More information7.5 Operations with Matrices. Copyright Cengage Learning. All rights reserved.
7.5 Operations with Matrices Copyright Cengage Learning. All rights reserved. What You Should Learn Decide whether two matrices are equal. Add and subtract matrices and multiply matrices by scalars. Multiply
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationCharacteristic finite-difference solution Stability of C C (CDS in time/space, explicit): Example: Effective numerical wave numbers and dispersion
Spring 015 Lecture 14 REVIEW Lecture 13: Stability: Von Neumann Ex.: 1st order linear convection/wave eqn., F-B scheme Hyperbolic PDEs and Stability nd order wave equation and waves on a string Characteristic
More informationWhat are Numerical Methods? (1/3)
What are Numerical Methods? (1/3) Numerical methods are techniques by which mathematical problems are formulated so that they can be solved by arithmetic and logic operations Because computers excel at
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationDifferential equations. Background mathematics review
Differential equations Background mathematics review David Miller Differential equations First-order differential equations Background mathematics review David Miller Differential equations in one variable
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 information1D heat conduction problems
Chapter 1D heat conduction problems.1 1D heat conduction equation When we consider one-dimensional heat conduction problems of a homogeneous isotropic solid, the Fourier equation simplifies to the form:
More informationPhysics 432, Physics of fluids: Assignment #1 Solutions
Physics 432, Physics of fluids: Assignment #1 Solutions 1. This is an investigation of the two-dimensional draining bathtub problem described in lecture, where the tub has fluid level h, width 2L and drain
More informationBOUNDARY VALUE PROBLEMS
BOUNDARY VALUE PROBLEMS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 BOUNDARY VALUE PROBLEMS 3 NEWTONS SHOOTING METHOD 4 SUMMARY OUTLINE 1 REVIEW 2 BOUNDARY VALUE PROBLEMS 3 NEWTONS SHOOTING
More informationqxbxg. That is, the heat rate within the object is everywhere constant. From Fourier s
PROBLEM.1 KNOWN: Steady-state, one-dimensional heat conduction through an axisymmetric shape. FIND: Sketch temperature distribution and explain shape of curve. ASSUMPTIONS: (1) Steady-state, one-dimensional
More informationConstructing Taylor Series
Constructing Taylor Series 8-8-200 The Taylor series for fx at x = c is fc + f cx c + f c 2! x c 2 + f c x c 3 + = 3! f n c x c n. By convention, f 0 = f. When c = 0, the series is called a Maclaurin series.
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationPhysics of Everyday Phenomena. Chapter 2
Physics of Everyday Phenomena W. Thomas Griffith Juliet W. Brosing Chapter 2 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Question 2.1 Ben leaves his home
More information(A) Opening Problem Newton s Law of Cooling
Lesson 55 Numerical Solutions to Differential Equations Euler's Method IBHL - 1 (A) Opening Problem Newton s Law of Cooling! Newton s Law of Cooling states that the temperature of a body changes at a rate
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 1.3 MODEL ANALYSIS
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 informationPhysics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =...
Physics 606, Quantum Mechanics, Final Exam NAME Please show all your work. (You are graded on your work, with partial credit where it is deserved.) All problems are, of course, nonrelativistic. 1. Consider
More informationQuiescent Steady State (DC) Analysis The Newton-Raphson Method
Quiescent Steady State (DC) Analysis The Newton-Raphson Method J. Roychowdhury, University of California at Berkeley Slide 1 Solving the System's DAEs DAEs: many types of solutions useful DC steady state:
More informationFundamentals of Heat Transfer Muhammad Rashid Usman
Fundamentals of Heat Transfer Muhammad Rashid Usman Institute of Chemical Engineering and Technology University of the Punjab, Lahore. Figure taken from: http://heatexchanger-design.com/2011/10/06/heat-exchangers-6/
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationASSUMPTIONS: (1) Homogeneous medium with constant properties, (2) Negligible radiation effects.
PROBEM 5.88 KNOWN: Initial temperature of fire clay bric which is cooled by convection. FIND: Center and corner temperatures after 50 minutes of cooling. ASSUMPTIONS: () Homogeneous medium with constant
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationPROBLEM 3.10 KNOWN: Dimensions and surface conditions of a plate thermally joined at its ends to heat sinks at different temperatures. FIND: (a) Differential equation which determines temperature distribution
More informationMon Jan Improved acceleration models: linear and quadratic drag forces. Announcements: Warm-up Exercise:
Math 2250-004 Week 4 notes We will not necessarily finish the material from a given day's notes on that day. We may also add or subtract some material as the week progresses, but these notes represent
More informationSECTION 7: CURVE FITTING. MAE 4020/5020 Numerical Methods with MATLAB
SECTION 7: CURVE FITTING MAE 4020/5020 Numerical Methods with MATLAB 2 Introduction Curve Fitting 3 Often have data,, that is a function of some independent variable,, but the underlying relationship is
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationThe Shooting Method for Boundary Value Problems
1 The Shooting Method for Boundary Value Problems Consider a boundary value problem of the form y = f(x, y, y ), a x b, y(a) = α, y(b) = β. (1.1) One natural way to approach this problem is to study the
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationENERGY ANALYSIS OF CLOSED SYSTEMS
Thermodynamics: An Engineering Approach Seventh Edition in SI Units Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 4 ENERGY ANALYSIS OF CLOSED SYSTEMS Mehmet Kanoglu University of Gaziantep
More informationFIND: (a) Sketch temperature distribution, T(x,t), (b) Sketch the heat flux at the outer surface, q L,t as a function of time.
PROBLEM 5.1 NOWN: Electrical heater attached to backside of plate while front surface is exposed to convection process (T,h); initially plate is at a uniform temperature of the ambient air and suddenly
More informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationNovel determination of dierential-equation solutions: universal approximation method
Journal of Computational and Applied Mathematics 146 (2002) 443 457 www.elsevier.com/locate/cam Novel determination of dierential-equation solutions: universal approximation method Thananchai Leephakpreeda
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Current Technical Advisor: Dr. Kelly O. Homan Original Technical Advisor: Dr. D.C. Look
ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT Current Technical Advisor: Dr. Kelly O. Homan Original Technical Advisor: Dr. D.C. Look (Updated on 8 th December 2009, Original version: 3 rd November 2000) 1
More information