Theory of! Partial Differential Equations-I!
|
|
- Lorraine Nicholson
- 5 years ago
- Views:
Transcription
1 hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and Nonlinear Advecion Equaions! Quasi-linear Second Order Equaions!!- Classiicaion: hyperbolic, parabolic, ellipic! - Domain o Dependence/Inluence! Ill-Posed Problems! Conservaive orm! Navier-Sokes equaions! Eamples o equaions! + U = 0 D = 0 c = 0 Advecion! Diusion! Wave propagaion! Evoluion in ime! Navier-Sokes equaions! Diusion par! + u + v = 1 P ρ + µ u ρ + u Advecion par! + = 0 Laplace equaion! + v = 0 Laplace par! Deiniions! The order o PDE is deermined by he highes derivaives! Linear i no powers or producs o he unknown uncions or is parial derivaives are presen.! + y =, + + = 0 Quasi-linear i i is rue or he parial derivaives o highes order.!! + = =, + y Quasi-linear irs order! parial dierenial equaions!
2 Consider he quasi-linear irs order equaion! a + b = c where he coeiciens are uncions o,y, and, bu no he derivaives o : a = a(, y, ) b = b(, y, ) c = c(, y, ) The soluion o his equaions deines a single valued surace (,y) in hree-dimensional space:! y =(,y) An arbirary change in is given by:! d = d + dy dy y d d The original equaion and he condiions or a small change can be rewrien as:! a + b = c d = d + dy,, 1 (a,b,c) = 0,, 1 (d,dy, d ) = 0 Normal o he surace! =(,y) Boh (a,b,c) and (d,dy,d) are in he surace!! The normal vecor o he curve =(,y) 1! Picking he displacemen in he direcion o (a,b,c)! (d,dy,d ) = ds (a,b, c) -1! Same argumens in he y-direcion. Thus =(,y) n =,, 1 Separaing he componens! d ds = a; dy ds = b; d dy = a b d ds = c;
3 The hree equaions speciy lines in he -y plane! d ds = a; Given he iniial condiions:! dy ds = b; Characerisics! d ds = c; = (s, o ); y = y(s, o ); = (s, o ); he equaions can be inegraed in ime! y d ds = a; dy ds = b; d ds = c; (s); y(s); (s) (0); y(0); (0) slope! d dy = a b Consider he linear advecion equaion! + U = 0 The characerisics are given by:! or! ds = 1; d ds = U; d = U; d = 0; d ds = 0; Which shows ha he soluion moves along sraigh characerisics wihou changing is value! Graphically:! Noice ha hese resuls are speciic or:! + U = 0 (,) = = 0 ( U) The soluion is hereore:! (,) = g U ( ) g() = (, = 0) where! This can be veriied by direc subsiuion:! Se! η(,) = U = g η Then! η = g η ( U ) and! Subsiue ino he original equaion! + U = g η U ( ) + U g = g η η = g η ( 1) η = 0 Disconinuous iniial daa! + U = 0 Since he soluion propagaes along characerisics compleely independenly o he soluion a he ne spaial poin, here is no requiremen ha i is diereniable or even coninuous! (,) = g( U)
4 Add a source! +U = The characerisics are given by:! ds =1; d ds = U; d ds = ; or! d = U; d = = (0)e Moving wave wih! decaying ampliude! Consider a nonlinear (quasi-linear) advecion equaion! + = 0 The characerisics are given by:! or! d =1; ds ds = ; d ds = 0; d = ; d = 0; The slope o he characerisics depends on he value o (,). 3 Muli-valued soluion Why unphysical soluions?! - Because mahemaical equaion neglecs! some physical process (dissipaion)! + ε = 0 Burgers Equaion! 1 Addiional condiion is required o pick ou he physically! Relevan soluion! Correc soluion is epeced rom Burgers equaion wih ε 0 Enropy Condiion hp://users.wpi.edu/~grear/me61.hml! Consrucing physical soluions! + = 0 In mos cases he soluion is no allowed o be muliple valued and he physical soluion mus be reconsruced using conservaion o! The disconinuous soluion propagaes wih a shock speed ha is dieren rom he slope o he characerisics on eiher side! Theory o! Parial Dierenial Equaions-II! Gréar Tryggvason! Spring 010!
5 Ouline! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and Nonlinear Advecion Equaions! Quasi-linear Second Order Equaions!!- Classiicaion: hyperbolic, parabolic, ellipic! - Domain o Dependence/Inluence! Ill-Posed Problems! Conservaive orm! Navier-Sokes equaions! Quasi-linear Second order parial dierenial equaions! Consider! where! a + b + c = d a = a(, y,,, y ) b = b(, y,,, y ) c = c(, y,,, y ) d = d(, y,,, y ) Firs wrie he second order PDE as a sysem o irs order equaions! Deine! hen! v = a + b + c = d = and! w = The second equaion is obained rom! a v + b v + c w = d w v = 0 Thus! Is equivalen o:! where! a + b + c = d a v + b v + c w = d w v = 0 v = w = Any high-order PDE can be rewrien as a sysem o irs order equaions!! a v + b v + c w = d In mari orm! or! v + b a c a v = d a w 1 0 w 0 u + Au y = s w v = 0 u = v w Are here lines in he -y plane, along which he soluion is deermined by an ordinary dierenial equaion?!
6 The oal derivaive is! dv d = v + v = v +α v where! Rae o change o v wih, along he line y=y()! I here are lines (deermined by α) where he soluion is governed by ODEʼs, hen i mus be possible o rewrie he equaions in such a way ha he resul conains only α and he oal derivaives.! α = Add he original equaions:! v + b v a + c a Is his ever equal o! w + l w v v + α v w + l + α w d = a For some lʼs and α = d a Compare he erms:! v + b v a + c a w + l w v v + α v + l w w +α Thereore, we mus have:! b a l = α c a = l α = d a = d a Characerisic lines eis i:! Or, in mari orm:! b a α c a b a l = α c a = l α 1 l 1 = 0 α 0 l The original equaion is:! v + b a c a v = d a w 1 0 w 0 A! b a α 1 = 0 Rewrie! c a α 0 as! l l b a 1 α 0 = 0 c a 0 0 α 0 A T! l ( A T αi )l = 0 The equaion has a soluion only i he deerminan is zero! b a α 1 c a α l The deerminan is:! A T αi = 0 α b a α + c a = 0 Or, solving or α! α = 1 ( a b ± b 4ac) = 0 0
7 α = 1 ( a b ± b 4ac) b 4ac > 0 Two real characerisics! Eamples! b 4ac = 0 One real characerisics! b 4ac < 0 No real characerisics! Eamples! Eamples! c = 0 Comparing wih he sandard orm! a + b + c = d shows ha! a = 1; b = 0; c = c ; d = 0; b 4ac = c = 4c > 0 Hyperbolic! + = 0 Comparing wih he sandard orm! a + b + c = d shows ha! a = 1; b = 0; c = 1; d = 0; b 4ac = = 4 < 0 Ellipic! Eamples! Eamples! = D Comparing wih he sandard orm! a + b + c = d shows ha! a = 0; b = 0; c = D; d = 0; b 4ac = D = 0 Parabolic! Diusion equaion! Laplace equaion! c = 0 Hyperbolic! = D Parabolic! + = 0 Ellipic!
8 α = α = 1 ( a b ± b 4ac) Wave Equaion! Hyperbolic! Parabolic! Ellipic! c = 0 Firs wrie he equaion as a sysem o irs order equaions! Inroduce! yielding! u = ; v = ; c v = 0 v = 0 since! rom he pde! = To ind he characerisics! c v = 0 v = 0 v ( ) = de de A T αi We can also use! α = 1 a b ± b 4ac α c α = ±c + 0 c α = α c v ( ) = 0 ( ) wih! b = 0; a =1; c = c = 0 Two characerisic lines! P d = + 1 c ; d = 1 c d = + 1 c d = 1 c To ind he soluion we need α = ±c o ind he eigenvecors! v c = 0 α 1 v +l c α = 0 Take =1! For! α = +c For! α = c c l = 0 l = c +c l = 0 l = +c l = 0
9 For! α = +c =1 l = c 1 c v = 0 c v = 0 + c c v + c v = 0 du c dv = 0 on d = +c Similarly:! For! α = c =1 l = +c du + c dv = 0 on d = c Add he equaions! Relaion beween he oal derivaive on he characerisic! du c dv = 0 on d = +c du + c dv = 0 on d = c For consan c! dr 1 = 0 on d = +c where r = u cv 1 dr = 0 on d = c where r = u + cv r 1 and r are called he Rieman invarians! Wave equaion-general! The general soluion can hereore be wrien as:! (,) = r 1 ( c) + r ( + c) Two characerisic lines! d = + 1 c ; d = 1 c where! r 1 ( ) = c r ( ) = + c = 0 = 0 P d = + 1 c d = 1 c Can also be veriied by direc subsiuion! hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-III! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and Nonlinear Advecion Equaions! Quasi-linear Second Order Equaions!!- Classiicaion: hyperbolic, parabolic, ellipic! - Domain o Dependence/Inluence! Ill-Posed Problems! Navier-Sokes equaions!
10 Ill-posed problems! Consider he iniial value problem:! = This is simply Laplaceʼs equaion! + = 0 which has a soluion i / or are given on he boundaries! Here, however, his equaion appeared as an iniial value problem, where he only boundary condiions available are a = 0. Since his is a second order equaion we will need wo condiions, which we may assume are ha and / are speciied a =0.! The general soluion can be wrien as:! ( ) = a k ()e ik, Look or soluions o he ype:! = a k ()e ik k where he aʼs depend on he iniial condiions! Subsiue ino:! + = 0 o ge:! d a k = k a k d a k = k a k General soluion! a k () = Ae k + Be k A, B deermined by iniial condiions! a( 0); ika(0) Long wave wih shor wave perurbaions! Generally, boh A and B are non-zero! Thereore:! a() as! Ill-posed Problem!
11 Ill-posed problems generally appear when he iniial or boundary daa and he equaion ype do no mach.! Frequenly arise because small bu imporan higher order eecs have been negleced! Summary! Ill-posedness generally maniess isel in he eponenial growh o small perurbaions so ha he soluion does no depend coninuously on he iniial daa! Summary! Navier-Sokes equaions! Why is he classiicaion Imporan?! 1. Iniial and boundary condiions!. Dieren physics! 3. Dieren numerical mehod apply! Parabolic par! + u + v = 1 P ρ + µ u ρ + u + v = 0 Hyperbolic par! Ellipic equaion! Summary! The Navier-Sokes equaions conain hree equaion ypes ha have heir own characerisic behavior! Depending on he governing parameers, one behavior can be dominan! The dieren equaion ypes require dieren soluion echniques! For inviscid compressible lows, only he hyperbolic par survives! Summary! In he ne several lecures we will discuss numerical soluions echniques or each class:! Hyperbolic equaions, including! soluions o he Euler equaions! Parabolic equaions! Ellipic equaions! Then we will consider advecion/diusion equaions and he special consideraions needed here.! Finally we will reurn o he ull Navier-Sokes equaions!
Theory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationA Note on Fractional Electrodynamics. Abstract
Commun Nonlinear Sci Numer Simula 8 (3 589 593 A Noe on Fracional lecrodynamics Hosein Nasrolahpour Absrac We invesigae he ime evoluion o he racional elecromagneic waves by using he ime racional Maxwell's
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationENGI 9420 Engineering Analysis Assignment 2 Solutions
ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion
More informationScientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS
Scieniic Research o he Insiue o Mahemaics and Compuer Science DIERENT VARIANTS O THE BOUNDARY ELEMENT METHOD OR PARABOLIC EQUATIONS Ewa Majchrzak,, Ewa Ładyga Jerzy Mendakiewicz, Alicja Piasecka Belkhaya
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More information8. Basic RL and RC Circuits
8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationINDEX. Transient analysis 1 Initial Conditions 1
INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera
More informationdy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationRobotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.
Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationMath 2214 Solution Test 1A Spring 2016
Mah 14 Soluion Tes 1A Spring 016 sec Problem 1: Wha is he larges -inerval for which ( 4) = has a guaraneed + unique soluion for iniial value (-1) = 3 according o he Exisence Uniqueness Theorem? Soluion
More informationLecture #6: Continuous-Time Signals
EEL5: Discree-Time Signals and Sysems Lecure #6: Coninuous-Time Signals Lecure #6: Coninuous-Time Signals. Inroducion In his lecure, we discussed he ollowing opics:. Mahemaical represenaion and ransormaions
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationMath 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.
1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be
More informationu(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x
. 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih
More information4.5 Constant Acceleration
4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),
More information1 Solutions to selected problems
1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen
More informationMath 221: Mathematical Notation
Mah 221: Mahemaical Noaion Purpose: One goal in any course is o properly use he language o ha subjec. These noaions summarize some o he major conceps and more diicul opics o he uni. Typing hem helps you
More information1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.
. Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.
More informationHamilton Jacobi equations
Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
More informationR.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#
.#W.#Erickson# Deparmen#of#Elecrical,#Compuer,#and#Energy#Engineering# Universiy#of#Colorado,#Boulder# Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance,
More informationLinear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation
Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationConvolution. Lecture #6 2CT.3 8. BME 333 Biomedical Signals and Systems - J.Schesser
Convoluion Lecure #6 C.3 8 Deiniion When we compue he ollowing inegral or τ and τ we say ha he we are convoluing wih g d his says: ae τ, lip i convolve in ime -τ, hen displace i in ime by seconds -τ, and
More informationSMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.
SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationLinear Control System EE 711. Design. Lecture 8 Dr. Mostafa Abdel-geliel
Linear Conrol Sysem EE 7 MIMO Sae Space Analysis and Design Lecure 8 Dr. Mosafa Abdel-geliel Course Conens Review Sae Space SS modeling and analysis Sae feed back design Oupu feedback design Observer design
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationChallenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k
Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationA New Perturbative Approach in Nonlinear Singularity Analysis
Journal of Mahemaics and Saisics 7 (: 49-54, ISSN 549-644 Science Publicaions A New Perurbaive Approach in Nonlinear Singulariy Analysis Ta-Leung Yee Deparmen of Mahemaics and Informaion Technology, The
More informationMath 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm
Mah Spring 6 Differenial Equaions Final Exam Due Dae: Tuesday, June 6, 5pm Your name (please prin): Insrucions: This is an open book, open noes exam. You are free o use a calculaor or compuer o check your
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationIn this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should
Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationChapter 2: Principles of steady-state converter analysis
Chaper 2 Principles of Seady-Sae Converer Analysis 2.1. Inroducion 2.2. Inducor vol-second balance, capacior charge balance, and he small ripple approximaion 2.3. Boos converer example 2.4. Cuk converer
More informationdi Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.
di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x
More informationTHE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI
THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach
More informationin Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology
Risk and Saey in Engineering Pro. Dr. Michael Havbro Faber ETH Zurich, Swizerland Conens o Today's Lecure Inroducion o ime varian reliabiliy analysis The Poisson process The ormal process Assessmen o he
More informationASTR415: Problem Set #5
ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua
More informationExam 1 Solutions. 1 Question 1. February 10, Part (A) 1.2 Part (B) To find equilibrium solutions, set P (t) = C = dp
Exam Soluions Februar 0, 05 Quesion. Par (A) To find equilibrium soluions, se P () = C = = 0. This implies: = P ( P ) P = P P P = P P = P ( + P ) = 0 The equilibrium soluion are hus P () = 0 and P () =..
More informationMath 2214 Solution Test 1B Fall 2017
Mah 14 Soluion Tes 1B Fall 017 Problem 1: A ank has a capaci for 500 gallons and conains 0 gallons of waer wih lbs of sal iniiall. A soluion conaining of 8 lbsgal of sal is pumped ino he ank a 10 galsmin.
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationSolutions for homework 12
y Soluions for homework Secion Nonlinear sysems: The linearizaion of a nonlinear sysem Consider he sysem y y y y y (i) Skech he nullclines Use a disincive marking for each nullcline so hey can be disinguished
More informationMEI STRUCTURED MATHEMATICS 4758
OXFORD CAMBRIDGE AND RSA EXAMINATIONS Advanced Subsidiary General Cerificae of Educaion Advanced General Cerificae of Educaion MEI STRUCTURED MATHEMATICS 4758 Differenial Equaions Thursday 5 JUNE 006 Afernoon
More informationFitting Parameters in a Differential Equation With a Non-Analytical Solution
NDSolve u.nb 1 Fiing Parameers in a Dierenial Equaion Wih a Non-Analyical Soluion Mahew Smih 9-Jan-005 In[]:= Remove "Global` " O General::"spell" O General::"spell1" Graphics` Oen dierenial equaions do
More information1 st order ODE Initial Condition
Mah-33 Chapers 1-1 s Order ODE Sepember 1, 17 1 1 s order ODE Iniial Condiion f, sandard form LINEAR NON-LINEAR,, p g differenial form M x dx N x d differenial form is equivalen o a pair of differenial
More informationKinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.
Kinemaics Vocabulary Kinemaics and One Dimensional Moion 8.1 WD1 Kinema means movemen Mahemaical descripion of moion Posiion Time Inerval Displacemen Velociy; absolue value: speed Acceleraion Averages
More informationCHAPTER 12 DIRECT CURRENT CIRCUITS
CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As
More informationLecture 10: Wave equation, solution by spherical means
Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationPlasma Astrophysics Chapter 3: Kinetic Theory. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Asrophysics Chaper 3: Kineic Theory Yosuke Mizuno Insiue o Asronomy Naional Tsing-Hua Universiy Kineic Theory Single paricle descripion: enuous plasma wih srong exernal ields, imporan or gaining
More informationUndetermined coefficients for local fractional differential equations
Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani
More information( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.
Andrei Tokmakoff, MIT Deparmen of Chemisry, 3/14/007-6.4 PERTURBATION THEORY Given a Hamilonian H = H 0 + V where we know he eigenkes for H 0 : H 0 n = E n n, we can calculae he evoluion of he wavefuncion
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More informationENV 6015 Solution to Mixing Problem Set
EN 65 Soluion o ixing Problem Se. A slug of dye ( ) is injeced ino a single ank wih coninuous mixing. The flow in and ou of he ank is.5 gpm. The ank volume is 5 gallons. When will he dye concenraion equal
More informationCHBE320 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang
CHBE320 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Spring 208 Dep. of Chemical and Biological Engineering CHBE320 Process Dynamics and Conrol 4- Road Map of he Lecure
More informationPhysics 1402: Lecture 22 Today s Agenda
Physics 142: ecure 22 Today s Agenda Announcemens: R - RV - R circuis Homework 6: due nex Wednesday Inducion / A curren Inducion Self-Inducance, R ircuis X X X X X X X X X long solenoid Energy and energy
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationNon-uniform circular motion *
OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform
More informationL p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity
ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.
More informationSolution of Integro-Differential Equations by Using ELzaki Transform
Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More information