Linear Advection Equation
|
|
- Gregory Cannon
- 6 years ago
- Views:
Transcription
1 Linear Advection Equation Numerical Methods Laboratory Modelling the Climate System ARC Centre of Excellence for Climate System Science 2nd Annual Winter School 18 June Analytical Solution The linear advection equation is given by u t + σ u x = 0, (1) where u = u(x, t) and σ = constant. To solve the equation we choose cyclic boundary conditions: u(0, t) = u(l, t) and the initial condition: u(x, 0) = f(x) = Ae ikx for 0 x L, f(x + L) = f(x). This means the initial condition is just a plane wave. Using separation of variables, we have u(x, t) = X(x)T (t). meaning Eqn. (1) becomes X(x) T t + σt (t) X x = 0, 1
2 and dividing by X(x)T(t) 1 dt (t) T (t) = σ 1 X(x) dx(x). This can only have a solution if both sides are constant, as the left hand side (LHS) is a function of t only and the right hand side (RHS) is a function of x only. Let so that Therefore, 1 dx X dx = λ 1 dt T = σλ. dx dx = λx and dt = σλt. Which have solutions: so X(x) = X 0 e λx and T (t) = T 0 e σλt u(x, t) = X(x)T (t) = X 0 T 0 e λx σλt = u 0 e λ(x σt). The initial condition gives so u(x, 0) = u 0 e λ(x 0) = Ae ikx A = u 0 and λ = ik, u(x, t) = u 0 e ik(x σt) This solution is that of a plane wave that propagates forward in the x direction with velocity σ. This is why the equation is called the linear advection equation the initial condition is transported forward by the bulk flow with velocity σ. 2
3 2 Numerical Solution Using Centered Differences We now analyze how the centered differences (in space) method changes the solution to the advection equation. Write the linear advection equation from (1) as where u j t = σ u j+1 u j 1 u j = U(t)e ijk. Here we have introduced a spatial discretization into the equation. Then { } du(t) e ijk = σ U(t)e ik(j+1) U(t)e ik(j 1) du(t) Euler s formula tells us so = σu(t)e ijk { { } e ik e ik = σu(t). e ±iζ = cos ζ ± i sin ζ, du(t) Now define 2i sin k = σu(t) 1 du(t) U(t) σ sin k = σ k, = d[ln(u(t))] } e ik e ik sin k = ikσu(t) k sin k = ikσ k. and this gives d[ln(u(t))] = ikσ, ln(u(t)) ln(u 0 ) = ikσ (t t 0 ). 3
4 Therefore ( ) U(t) ln = ikσ t, since u(0) = u 0. u 0 U(t) = u 0 e ikσ t Substituting back into u j u j = u 0 e ikσ t e ijk = u 0 e ik(j σ t). But from the analytical solution we have u(j, t) = u 0 e ik(j σt) and since σ σ there is a discretization error. Since σ = f(k), the error is a function of the wavenumber k. This means each wavenumber will behave differently and may not propagate at the same speed σ as in the analytical solution. This is called numerical dispersion. The result derived here implies that centered differences is a dispersive numerical scheme. Let s have a closer look at the properties of σ. The smallest wavelength that can be resolved in a grid is λ s = (Figure 1). Since k = 2π/λ, and k s = 2π/ = π/, σ s = σ sin(k s) k s = σ sin π π = 0. No matter what the true phase speed, the shortest waves will have a numerical phase speed of 0. This implies large error are present for the shortest resolvable waves! What about longer waves? Say λ = 4, λ = 4 k = π σ = σ sin(π/2) π/2 = σ. 4
5 Figure 1: Representation of the minimum resolvable wavelength (λ, the distance from peak to peak) on a numerical grid. Or say λ = 8, λ = 8 k = π 4 σ = σ sin(π/4) π/4 = 0.9σ. This means the longer the wave, the smaller the error in propagation speed. 3 More on Phase Speed General Case Let λ = n σ = σ sin(2π/n) 2π/n as n σ n = lim n σ sin(2π/n) 2π/n. Use L Hôpital s rule: f(x) lim x a g(x) = f (x) g (x). 5
6 So σ n = σ 2π n 2 σ n = σ. cos(2π/n) n= 2π n 2 So the larger n is (ie. the wavelength) the smaller the error! Group Velocity In a wave packet (an overlay of many different waves with varying wave numbers) energy is transported with the group velocity σ g = ω k where ω = frequency. The analytical solution is where u = u 0 e ik(x σt) = u 0 e i(kx ωt) ω = kσ. σ g = (kσ) = σ. k This tells us that the group velocity equals the phase speed, so the wave is nondispersive! Numerical Solution Lets apply group velocity to the numerical solution: so ω = kσ = kσ sin(k) k = σ sin(k) σ g = ω k = σ cos(k) 6
7 Example: λ = k = π/ σ g = σ cos π = σ. Therefore, the numerical group velocity is not equal to the phase velocity as in the analytical solution. In fact it has the opposite sign! 4 Discretization Error of Backwards Differencing We now analyse the backward in space differencing scheme. We start by writing the advection equation from (1) as u t = σ u x Let u j = U(t)e ijk du(t) e ijk = σ u j t = σ u j u j 1. ( 1 du(t) = σ U(t) = σ = ikσ ( U(t)e ijk U(t)e i(j 1)k ) 1 e ik ( 1 (cos k i sin k) ( 1 (cos k i sin k) ik = ikσ where ( ) 1 cos k + i sin k σ = σ. ik Therefore and U(t) = u 0 e ikσ t u j = U(t)e ijk = u 0 e ik(j σ t) 7 ) ) )
8 Analytic Solution: u = u 0 e ik(x σt) at x = j u = u 0 e ik(j σt). Discretization error since σ σ. Investigate σ a bit more: u j = u 0 e = u 0 e 1 cos k i sin k ik[j σt( ik ) σt( ik )] 1 cos k σt( sin k ) e ik[j σt( k )]. Term containing sin k/k numerical dispersion. Term containing (1 cos k)/ damping of solution for σ > 0, or growth and hence instability for σ < 0. 8
Understand the basic principles of spectroscopy using selection rules and the energy levels. Derive Hund s Rule from the symmetrization postulate.
CHEM 5314: Advanced Physical Chemistry Overall Goals: Use quantum mechanics to understand that molecules have quantized translational, rotational, vibrational, and electronic energy levels. In a large
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationPartial differential equations (ACM30220)
(ACM3. A pot on a stove has a handle of length that can be modelled as a rod with diffusion constant D. The equation for the temperature in the rod is u t Du xx < x
More informationu tt = a 2 u xx u tt = a 2 (u xx + u yy )
10.7 The wave equation 10.7 The wave equation O. Costin: 10.7 1 This equation describes the propagation of waves through a medium: in one dimension, such as a vibrating string u tt = a 2 u xx 1 This equation
More informationNumerical dispersion and Linearized Saint-Venant Equations
Numerical dispersion and Linearized Saint-Venant Equations M. Ersoy Basque Center for Applied Mathematics 11 November 2010 Outline of the talk Outline of the talk 1 Introduction 2 The Saint-Venant equations
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationWaves Primer II Basic wave form. Waves are: 1. Important 2. Easy 3. Neat
Waves Primer II Waves are: 1. Important 2. Easy 3. Neat 23.1 Basic wave form The wave form that you will almost always use is a = Re[A c e i(kx+ly ωt) ] (23.1) Where a variable A c waveamplitude(complexinthiscase)
More informationEXAM 3 MAT 167 Calculus I Spring is a composite function of two functions y = e u and u = 4 x + x 2. By the. dy dx = dy du = e u x + 2x.
EXAM MAT 67 Calculus I Spring 20 Name: Section: I Each answer must include either supporting work or an explanation of your reasoning. These elements are considered to be the main part of each answer and
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
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 informationMARN 5898 Fourier Analysis.
MAR 5898 Fourier Analysis. Dmitriy Leykekhman Spring 2010 Goals Fourier Series. Discrete Fourier Transforms. D. Leykekhman - MAR 5898 Parameter estimation in marine sciences Linear Least Squares 1 Complex
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 5 Fourier Analysis of Finite Difference Methods In this lecture, we determine the stability of PDE discretizations using Fourier analysis. First, we begin with Fourier analysis of PDE s, and then
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationProblem 1. Part a. Part b. Wayne Witzke ProblemSet #4 PHY ! iθ7 7! Complex numbers, circular, and hyperbolic functions. = r (cos θ + i sin θ)
Problem Part a Complex numbers, circular, and hyperbolic functions. Use the power series of e x, sin (θ), cos (θ) to prove Euler s identity e iθ cos θ + i sin θ. The power series of e x, sin (θ), cos (θ)
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
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 informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationFunctions of Several Variables
Functions of Severl Vribles Sketching Level Curves Sections Prtil Derivtives on every open set on which f nd the prtils, 2 f y = 2 f y re continuous. Norml Vector x, y, 2 f y, 2 f y n = ± (x 0,y 0) (x
More informationDispersion relations, linearization and linearized dynamics in PDE models
Dispersion relations, linearization and linearized dynamics in PDE models 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationMATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:
MATH 251 Final Examination December 19, 2012 FORM A Name: Student Number: Section: This exam has 17 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all
More informationCalculus of Variations
ECE 68 Midterm Exam Solution April 1, 8 1 Calculus of Variations This exam is open book and open notes You may consult additional references You may even discuss the problems (with anyone), but you must
More informationThe dynamics of a simple troposphere-stratosphere model
The dynamics of a simple troposphere-stratosphere model by Jessica Duin under supervision of Prof. Dr. A. Doelman, UvA Dr. W.T.M. Verkley, KNMI August 31, 25 Universiteit van Amsterdam Korteweg-de vries
More informationAPPM 1350 Final Exam Fall 2017
APPM 350 Final Exam Fall 207. (26 pts) Evaluate the following. (a) Let g(x) cos 3 (π 2x). Find g (π/3). (b) Let y ( x) x. Find y (4). (c) lim r 0 e /r ln(r) + (a) (9 pt) g (x) 3 cos 2 (π 2x)( sin(π 2x))(
More informationFundamental Solution
Fundamental Solution onsider the following generic equation: Lu(X) = f(x). (1) Here X = (r, t) is the space-time coordinate (if either space or time coordinate is absent, then X t, or X r, respectively);
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationLECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD
ECTURE 19: SEPARATION OF VARIABES, HEAT CONDUCTION IN A ROD The idea of separation of variables is simple: in order to solve a partial differential equation in u(x, t), we ask, is it possible to find a
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationTraveling Harmonic Waves
Traveling Harmonic Waves 6 January 2016 PHYC 1290 Department of Physics and Atmospheric Science Functional Form for Traveling Waves We can show that traveling waves whose shape does not change with time
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationAdvection, Conservation, Conserved Physical Quantities, Wave Equations
EP711 Supplementary Material Thursday, September 4, 2014 Advection, Conservation, Conserved Physical Quantities, Wave Equations Jonathan B. Snively!Embry-Riddle Aeronautical University Contents EP711 Supplementary
More informationDynamics of beams. Modes and waves. D. Clouteau. September 16, Department of Mechanical and Civil Engineering Ecole Centrale Paris, France
Dynamics of and waves Department of Mechanical and Civil Engineering Ecole Centrale Paris, France September 16, 2008 Outline 1 2 Pre-stressed Bernoulli-Euler in moving frames Kinematical assumption Resultant
More informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationChapter 1: Limits and Continuity
Chapter 1: Limits and Continuity Winter 2015 Department of Mathematics Hong Kong Baptist University 1/69 1.1 Examples where limits arise Calculus has two basic procedures: differentiation and integration.
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationSelected HW Solutions
Selected HW Solutions HW1 1 & See web page notes Derivative Approximations. For example: df f i+1 f i 1 = dx h i 1 f i + hf i + h h f i + h3 6 f i + f i + h 6 f i + 3 a realmax 17 1.7014 10 38 b realmin
More informationSound Waves Sound Waves:
3//18 Sound Waves Sound Waves: 1 3//18 Sound Waves Linear Waves compression rarefaction Inference of Sound Wave equation: Sound Waves We look at small disturbances in a compressible medium (note: compressible
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More informationTransmission Lines, Stacked Insulated Washers Lines, Tesla Coils and the Telegrapher s Equation
Transmission Lines, Stacked Insulated Washers Lines, Tesla oils and the Telegrapher s Equation Kurt Nalty August, 2008 Abstract Tesla coils have more in common with transmission lines than transformers.
More informationMth Review Problems for Test 2 Stewart 8e Chapter 3. For Test #2 study these problems, the examples in your notes, and the homework.
For Test # study these problems, the examples in your notes, and the homework. Derivative Rules D [u n ] = nu n 1 du D [ln u] = du u D [log b u] = du u ln b D [e u ] = e u du D [a u ] = a u ln a du D [sin
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationProblem 1 In each of the following problems find the general solution of the given differential
VI Problem 1 dt + 2dy 3y = 0; dt 9dy + 9y = 0. Problem 2 dt + dy 2y = 0, y(0) = 1, y (0) = 1; dt 2 y = 0, y( 2) = 1, y ( 2) = Problem 3 Find the solution of the initial value problem 2 d2 y dt 2 3dy dt
More informationDavid J. Starling Penn State Hazleton PHYS 214
All the fifty years of conscious brooding have brought me no closer to answer the question, What are light quanta? Of course today every rascal thinks he knows the answer, but he is deluding himself. -Albert
More informationWAVE PACKETS & SUPERPOSITION
1 WAVE PACKETS & SUPERPOSITION Reading: Main 9. PH41 Fourier notes 1 x k Non dispersive wave equation x ψ (x,t) = 1 v t ψ (x,t) So far, we know that this equation results from application of Newton s law
More informationMath 121 Winter 2010 Review Sheet
Math 121 Winter 2010 Review Sheet March 14, 2010 This review sheet contains a number of problems covering the material that we went over after the third midterm exam. These problems (in conjunction with
More informationNanomaterials for Photovoltaics (v11) 11 Quantum Confinement. Quantum confinement (I) The bandgap of a nanoparticle increases as the size decreases.
Quantum Confinement Quantum confinement (I) The bandgap of a nanoparticle increases as the size decreases. A rule of thumb is E g d where d is the particle size. Quantum confinement (II) The Schrodinger
More informationComplex Variables & Integral Transforms
Complex Variables & Integral Transforms Notes taken by J.Pearson, from a S4 course at the U.Manchester. Lecture delivered by Dr.W.Parnell July 9, 007 Contents 1 Complex Variables 3 1.1 General Relations
More informationLaminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations
Laminar Boundary Layers Answers to problem sheet 1: Navier-Stokes equations The Navier Stokes equations for d, incompressible flow are + v ρ t + u + v v ρ t + u v + v v = 1 = p + µ u + u = p ρg + µ v +
More information2003 Mathematics. Advanced Higher. Finalised Marking Instructions
2003 Mathematics Advanced Higher Finalised Marking Instructions 2003 Mathematics Advanced Higher Section A Finalised Marking Instructions Advanced Higher 2003: Section A Solutions and marks A. (a) Given
More informationSolitons : An Introduction
Solitons : An Introduction p. 1/2 Solitons : An Introduction Amit Goyal Department of Physics Panjab University Chandigarh Solitons : An Introduction p. 2/2 Contents Introduction History Applications Solitons
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationMath 180, Exam 2, Practice Fall 2009 Problem 1 Solution. f(x) = arcsin(2x + 1) = sin 1 (3x + 1), lnx
Math 80, Exam, Practice Fall 009 Problem Solution. Differentiate the functions: (do not simplify) f(x) = x ln(x + ), f(x) = xe x f(x) = arcsin(x + ) = sin (3x + ), f(x) = e3x lnx Solution: For the first
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationECE 301 Fall 2011 Division 1. Homework 1 Solutions.
ECE 3 Fall 2 Division. Homework Solutions. Reading: Course information handout on the course website; textbook sections.,.,.2,.3,.4; online review notes on complex numbers. Problem. For each discrete-time
More informationWaves, the Wave Equation, and Phase Velocity. We ll start with optics. The one-dimensional wave equation. What is a wave? Optional optics texts: f(x)
We ll start with optics Optional optics texts: Waves, the Wave Equation, and Phase Velocity What is a wave? f(x) f(x-) f(x-) f(x-3) Eugene Hecht, Optics, 4th ed. J.F. James, A Student's Guide to Fourier
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationFinite difference approximation
Finite difference approximation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen January
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationSOLUTIONS TO PROBLEMS FROM ASSIGNMENT 5. Problems 3.1:6bd
SOLUTIONS TO PROBLEMS FROM ASSIGNMENT 5 Statement. Solve the problem Problems 3.1:6bd u t = u xx, (t ), with the initial condition u(x, ) = f(x), where the functions u(x, t) and f(x) are assumed to be
More informationNumerics: The analysis and implementation of numerical methods for solving differential equations
Numerics: The analysis and implementation of numerical methods for solving differential equations Dr Hilary Weller, Lecturer version November 6, 217 Timetable Week Chapters to
More informationSeparation of Variables. A. Three Famous PDE s
Separation of Variables c 14, Philip D. Loewen A. Three Famous PDE s 1. Wave Equation. Displacement u depends on position and time: u = u(x, t. Concavity drives acceleration: u tt = c u xx.. Heat Equation.
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationLecture 1: sine-gordon equation and solutions
Lecture 1: sine-gordon equation and solutions Equivalent circuit Derivation of sine-gordon equation The most important solutions plasma waves a soliton! chain of solitons resistive state breather and friends
More informationLES of Turbulent Flows: Lecture 3
LES of Turbulent Flows: Lecture 3 Dr. Jeremy A. Gibbs Department of Mechanical Engineering University of Utah Fall 2016 1 / 53 Overview 1 Website for those auditing 2 Turbulence Scales 3 Fourier transforms
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 34 Outline 1 Infinite Order Finite Differences
More informationMath 201 Assignment #11
Math 21 Assignment #11 Problem 1 (1.5 2) Find a formal solution to the given initial-boundary value problem. = 2 u x, < x < π, t > 2 u(, t) = u(π, t) =, t > u(x, ) = x 2, < x < π Problem 2 (1.5 5) Find
More informationBASIC WAVE CONCEPTS. Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, Giancoli?
1 BASIC WAVE CONCEPTS Reading: Main 9.0, 9.1, 9.3 GEM 9.1.1, 9.1.2 Giancoli? REVIEW SINGLE OSCILLATOR: The oscillation functions you re used to describe how one quantity (position, charge, electric field,
More information= k, (2) p = h λ. x o = f1/2 o a. +vt (4)
Traveling Functions, Traveling Waves, and the Uncertainty Principle R.M. Suter Department of Physics, Carnegie Mellon University Experimental observations have indicated that all quanta have a wave-like
More informationChapter 4. Nonlinear Hyperbolic Problems
Chapter 4. Nonlinear Hyperbolic Problems 4.1. Introduction Reading: Durran sections 3.5-3.6. Mesinger and Arakawa (1976) Chapter 3 sections 6-7. Supplementary reading: Tannehill et al sections 4.4 and
More informationInverse Kinematics. Mike Bailey. Oregon State University. Inverse Kinematics
Inverse Kinematics Mike Bailey mjb@cs.oregonstate.edu inversekinematics.pptx Inverse Kinematics Forward Kinematics solves the problem if I know the link transformation parameters, where are the links?.
More informationLecture 24. Scott Pauls 5/21/07
Lecture 24 Department of Mathematics Dartmouth College 5/21/07 Material from last class The heat equation α 2 u xx = u t with conditions u(x, 0) = f (x), u(0, t) = u(l, t) = 0. 1. Separate variables to
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationIntroduction to Condensed Matter Physics
Introduction to Condensed Matter Physics The Reciprocal Lattice M.P. Vaughan Overview Overview of the reciprocal lattice Periodic functions Reciprocal lattice vectors Bloch functions k-space Dispersion
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationMath 113 Winter 2005 Departmental Final Exam
Name Student Number Section Number Instructor Math Winter 2005 Departmental Final Exam Instructions: The time limit is hours. Problem consists of short answer questions. Problems 2 through are multiple
More informationMath 266 Midterm Exam 2
Math 266 Midterm Exam 2 March 2st 26 Name: Ground Rules. Calculator is NOT allowed. 2. Show your work for every problem unless otherwise stated (partial credits are available). 3. You may use one 4-by-6
More informationUniversity of Toronto at Scarborough Department of Computer and Mathematical Sciences, Mathematics MAT C46S 2013/14.
University of Toronto at Scarborough Department of Computer and Mathematical Sciences, Mathematics MAT C46S 213/14 Problem Set #4 Due date: Thursday, March 2, 214 at the beginning of class Part A: 1. 1.
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationMath Assignment 14
Math 2280 - Assignment 14 Dylan Zwick Spring 2014 Section 9.5-1, 3, 5, 7, 9 Section 9.6-1, 3, 5, 7, 14 Section 9.7-1, 2, 3, 4 1 Section 9.5 - Heat Conduction and Separation of Variables 9.5.1 - Solve the
More informationENGIN 211, Engineering Math. Laplace Transforms
ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving
More informationWave Equation With Homogeneous Boundary Conditions
Wave Equation With Homogeneous Boundary Conditions MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 018 Objectives In this lesson we will learn: how to solve the
More informationDefinitions of the Laplace Transform (1A) Young Won Lim 1/31/15
Definitions of the Laplace Transform (A) Copyright (c) 24 Young W. Lim. Permission is granted to copy, distriute and/or modify this document under the terms of the GNU Free Documentation License, Version.2
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 informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationM412 Assignment 9 Solutions
M4 Assignment 9 Solutions. [ pts] Use separation of variables to show that solutions to the quarter-plane problem can be written in the form u t u xx ; t>,
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More informationIntroduction to Waves in Structures. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil
Introduction to Waves in Structures Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Waves in Structures Characteristics of wave motion Structural waves String Rod Beam Phase speed, group velocity Low
More information