APM526 CLASSROOM NOTES PART IB. APM526, Spring 2018 Last update: MAR 12
|
|
- Bryan Webster
- 5 years ago
- Views:
Transcription
1 APM526 CLASSROOM NOTES PART IB APM526, Spring 2018 Last update: MAR 12 1
2 PARTICLE METHODS 1. Particle based methods: Large systems of ODEs for individual particles. 2. Partial differential equations: Equations for the density of particles. Example: N particles moving with a speed, dependent on their location. ξ n = a(ξ n ), n = 1 : N The corresponding density function: u(x, t) = N n=1 δ(x ξ n (t)) 2
3 DISTRIBUTIONS g(x) distribution, defined by ψ(x)g(x) dx = α(ψ), ψ C 0 ψ : test function. Example: g(x) = 5δ (x) + cos(x). Example: Discretization of δ (x). Integration against a test function corresponds to physically taking a measurement! 3
4 THE CONCEPT OF A WEAK SOLUTION FOR CON- SERVATION LAWS t u(x, t) + x F (u, x, t) = 0 (1) Integrate against a test function ψ ( the microscope ) ψ(x) t u(x, t) dx F (u, x, t) x ψ(x) dx = 0, ψ C 0 (2) Definition: We call the u(x, t) a weak solution of (1) if (2) holds ψ C 0. This allows for finding solutions if u discontinuous and x u does not exist in the classical sense! 4
5 t u(x, t) = N n=1 δ (x ξ n (t))a(ξ n (t)) (3) u(x, t) = N n=1 δ(x ξ n (t)) t u = N n=1 δ (x ξ n (t))a(ξ n (t)), x u = integrated against any test function ψ: ψ(x) t u(x, t) dx = N n=1 δ (x ξ n (t)) ψ (x)a(x)u(x, t), ψ Satisfies the density equation weakly: t u(x, t) + x [a(x)u(x, t)] = 0 5
6 HOW TO COMPUTE DENSITIES FROM PARTICLES Define a grid and grid cells C n = [x n, x n+1 ]. Define indicator functions χ m (x) = ( ) 1 x Cm 0 else Use χ m as test functions χ m u(x, t) dx = N n=1 χ m (ξ n (t)) Counts the number of particles in C m. ( u(x, t) dx = N). 6
7 THE MEAN FIELD PROBLEM IID case: Independent and Identically Distributed particles t ξ n = a(ξ n ) Assume: particles interact t ξ n = a n ( ξ), ξ = (ξ1,.., ξ n ) Density: u(x, t) = N n=1 δ(x ξ n (t)) ψ(x) t u(x, t) dx = ψ (x) N n=1 a n ( ξ(t))δ(x ξ n (t)) dx 7
8 Mean field: N n=1 a n ( ξ(t))δ(x ξ n (t)) N n=1 a mf (u)δ(x ξ n (t)), u = for N. a mf is called the mean field! N n=1 δ(x ξ n (t)) Given a mean field a mf (u)(???) and N, we have (weakly) t u + x [a mf (u)u] = 0 Example: Gravity, Electrodynamics 8
9 PARTICLES AND DIFFUSION (Einstein theory of Brownian motion) Consider the pure transport (hyperbolic) case d ξ(t) = a(ξ(t)), ξ(t + t) = ξ(t) + ta(ξ(t)) (4) dt Alternative (probabilistic) interpretation: Consider only one single particle. u(x, t) dx = dp[ξ(t) = x]. This gives for the density u(x, t) = δ(x ξ(t)) (in weak form) the equation t u(x, t) + x [a(x)u] = 0 (5) Now consider N >> 1 identical particles, all evolving according to (4), and u(x, t) = N 1 n δ(x ξ n (t). Equation (5) still holds. So, for a system of N >> 1 identical particles, (5) describes the evolution of the density u((x, t). 9
10 Brownian motion: Replace a(ξ) with a random variable q(ξ) ξ(t + t) = ξ(t) + tq(ξ(t)), dp[q(ξ) = x] = φ(x, ξ)dx q(ξ) has a mean a(ξ) and variance σ 2 (ξ) φ(x, ξ) dx = 1, xφ(x, ξ) dx = a(ξ), (x a(ξ)) 2 φ(x, ξ) dx = σ 2 (ξ) Normalize φ: φ(x, ξ) = σ 1 φ 0( x a(ξ) σ(ξ), ξ) with φ 0 normalized distribution with mean 0 and variance 1. φ 0 (x, ξ) dx = 1, xφ 0 (x, ξ) dx = 0, x 2 φ 0 (x, ξ) dx = 1, ξ 10
11 Let u(x, t)dx = dp[ξ(t) = x] Sum up over all possibilities: u(x, t + t) = δ(y + tq x)φ(q, y)u(y, t) dydq = δ(y + tq x) 1 σ(y) φ 0( q a(y), y)u(y, t) dydq σ(y) weak formulation: ψ(x)u(x, t+ t) dx = for all test functions ψ C 0. ψ(x)u(x, t+ t) dx = Taylor expand the test function: ψ(x)u(x, t + t) dx = ψ(y+ tq) 1 σ(y) φ 0( q a(y), y)u(y, t) dydq σ(y) ψ(y+ ta+ tσz)φ 0 (z, y)u(y, t) dzdy [ψ(y)+( ta+ tσz)ψ (y)+ 1 2 ( ta+ tσz)2 ψ (y)]φ 0 (z, y)u(y, t) dzdy+h 11
12 [ψ(y)+ ta(y)ψ (y)+ 1 2 ( t2 a(y) 2 + t 2 σ(y) 2 )ψ (y)]u(y, t) dy+h.o.t. strong form: u(x, t+ t) = u(x, t) t x [au(x, t)] x [( t2 a(x) 2 + t 2 σ(x) 2 )u(x, t)] set σ 2 t σ2. so q(ξ) has a mean a(ξ) and variance σ2 (ξ) t Limiting equation for t 0: t u(x, t) = x [au(x, t)] x[σ(x) 2 )u(x, t)] 12
13 Theorem: A stochastic process of the form ξ(t + t) = ξ(t) + tq(ξ(t)) where q(ξ) is a random variable with dp[q(ξ) = x] = φ(x, ξ)dx and E(q(ξ)) = a(ξ), var(q(ξ)) = σ(ξ)2 t i.e. φ(x, ξ) dx = 1, xφ(x, ξ) dx = a(ξ), results in the limit t 0 in the equation t u(x, t) = x [au(x, t)] x [σ(x)2 u(x, t)] for the probability density u(x, t)dx = dp[ξ(t) = x] (x a(ξ)) 2 φ(x, ξ) dx = σ(ξ)2 t 13
14 THINK IN REVERSE: Given a PDE t u(x, t) = x [au(x, t)] x[σ(x) 2 )u(x, t)] Define a Monte Carlo particle method in the following way: Use N particles ξ n (t), n = 1 : N solve ξ n (t + t) = ξ n (t) + tq n (ξ n (t)) with dp[q n (ξ) = y] = φ(y, ξ)dx (1, x, x 2 )φ(x, ξ) dx = (1, a(ξ), a(ξ) 2 + σ(ξ)2 t ) 14
15 JOB: Given a probability distribution φ in [a, b] (a, b can be ±!) φ(x), φ(x) 0, b a φ dx = 1 and a random number generator, producing numbers uniformly distributed in (0, 1), compute random numbers q with dp[q = x] = φ(x)dx SOLUTION: set q = g(y), y (0, 1). g monotone (0, 1) (a, b) P[q < x] = P[g(y) < x] = P[y < g 1 (x)] = g 1 (x) Φ (x) = φ(x), Φ(x) = P[q < x] = g 1 (x) g(y) = Φ 1 (y) Algorithm: Take φ(x), compute the anti - derivative Φ(x) = x a φ(y) dy, Φ(a) = 0, Φ(b) = 1. set g(x) = Φ 1 (x), g(x) (a, b), compute y (0, 1) (with c.f. Matlab rand) set q = g(y) 15
16 EXAMPLE 1: exponential distribution a = 0, b =, φ(x) = e x anti - derivative: Φ(x) = H(x)[1 e x ] computing the inverse: g(y) = ln(1 y). compute y with rand and set q = ln(1 y) EXAMPLE 2: Gaussian: a =, b =, φ(x) = 1 2π exp( x2 2 ) anti - derivative: Φ(x) = 1 x exp( y2 ) dy = erf(x). 2π compute y with rand and set q = erf 1 (y) 2 16
17 2 dimensional distributions: f(x, y) JOB: compute a random pair (q 1, q 2 ) distributed according to f(x, y) conditional probabilities: f(x, y) = φ(y x)φ 1 (x), φ 1 (x) = f(x, y) dy φ 1 (x) is the probability density for x. φ(y x) is the probability density for y when x is already known. 1. compute a random number q 1 according to φ given q 1, compute q according to φ(y q 1 ). 17
18 % compute an exponential distribution with N particles clear format N=1000 % number of particles K=10 % number of grid points X=3% length of interval dx=x/k xi=-log(rand(1,n)) x=x*(0:k)/k % hh(k,n)=(x(k) xi(n)) 18
19 h=(x *ones(1,n) ones(k+1,1)*xi) hh=sum(h,2) u=diff(hh) u=u/n/dx figure(1) plot(x(1:k),u)
20 19
21 RELATION TO DIFFERENCE SCHEMES Given a discretization of t u = x (au)+ x (b x u) of the form u n (t+ t) = u n + ta 2 x (u n 1 u n+1 )+ tb x 2(u n+1 2u n +u n 1 )(t) u n (t + t) = αu n 1 + βu n + γu n 1 (t) α = ta 2 x + tb x2, β = 1 2 tb x 2, γ = ta 2 x + tb x 2, α + β + γ = 1 rewrite as a stochastic process: u n (t + t) = u n 1 (t) u n (t) u n+1 (t) P = α P = β P = γ CFL condition: α, β, γ 0 Modeling option: Add interaction between cells α = α( u(t))... Works for first order in time schemes (like LaxFriedrichs) and not for any higher order in time schemes (like LaxWendroff) 20
22 CONVERGENCE ISSUES Q: What are good approximations to t, x and good values for t, x for difference schemes as t, x 0? Q: What are good choices for φ(x, ξ) and the random variable q(ξ) for the random motion of the particle as N? 21
23 CONVERGENCE OF PARTICLE METHODS: Start with a single particle at position ξ ξ(t + t) = ξ(t) + tq(ξ), dp[q(ξ) = x] = φ(x, ξ)dx Probability: u(x, t)dx = dp[ξ(t) = x] Observable: Mesh: x j, indicator function χ j (x) = 1 if x {x j, x j+1 }, χ j (x) = 0, else Take a Counter: g j (t) = 1 x χ j(ξ(t)) Expectation: E[g j ] = 1 x x j+1 x j χ j (x)u(x, t) dx u(x j, t) 22
24 Compute N realizations of the process: Blackbox: given ξ = (ξ 1,.., ξ N ), compute N identically distributed random variables q n with dp[q = x] = φ(x, ξ(t)) dx Compute averages: G j (t) = 1 N n g n j (t) = 1 N 1 x G(x, t) should approx- Since all realizations are independent, imate E[g(x, t)] u(x, t) n χ j (ξ n (t)) 23
25 LAW OF LARGE NUMBERS AND CENTRAL LIMIT THE- OREM: Let ξ n, n = 1 : N be identically distributed random variables with mean = 0 and variance = σ 2 dp[ξ n = x] = f(x)dx, n, Compute the average q = 1 N Nn=1 ξ n. (1, x, x 2 )f(x) dx = (1, 0, σ 2 ) Then F (y)dy = dp[q = y] is given for N approximately by the Gaussian by N F (y) σ exp( Ny2 2σ 2 ) independent of the structure of f(x)! 24
26 COROLLARY: Let ξ n, n = 1 : N be identically distributed random variables with mean = a and variance = σ 2 dp[ξ n = x] = f(x)dx, n, Compute the average q = 1 N Nn=1 ξ n. (1, x, x 2 )f(x) dx = (1, a, a 2 + σ 2 ) Then F (y)dy = dp[q = y] is given for N approximately by the Gaussian by N a)2 F (y) exp( N(y σ 2σ 2 ) independent of the structure of f(x)! 25
27 THE 1 N PROBLEM: Compute N realizations of a stochastic process. Compute an observable of the process by averaging. The observable will converge to the correct mean a with a rate 1. N 26
28 MARKOV PROCESSES: ω: frequency (time scale). t: time step ξ(t + t) = (1 κ)ξ(t) + κ[ξ(t) + q], κ {0, 1}, P[κ = 1] = ω t, P[q = y] = φ(y, ξ(t)) Define u(x, t) dx = dp[ξ(t) = x] t u(x, t) = ωφ(x, t) 27
29 STABILITY IN THE - NORM U = max n U n Matrix norms: A 1 = max n A = max x 0 m Ax x A mn, A = max m Ax A x n A mn, A 2 = λ max (A T A) 28
30 Difference scheme (linear): u(x, t + t) = n α(x)t n u(x, t) A(m, m + n) = α(x m ) u(x m, t) := U m (t) U m (t + t) = n α m U m+n (t) U(t + t) (max x α(x) ) U(t) Consistency: m α(x m ) = 1 + O( t) If α(x m ) 0 we have α(x) = α(x) and stability U(t + t) (1 + O( t)) U(t). α m 0 gives the CFL condition! 29
31 MONOTONE OPERATORS U(t + t) = Φ(U(t)), U = (u 1,.., u N ), Φ = (φ 1,.., φ N ) Definition: Φ is monotone : consistency: φ j u k 0, j, k, t φ j (C 1) = C, C max j min j u j (t+ t) = max j u j (t+ t) = min j φ j (U(t)) max j φ j (U(t)) min j φ j ( 1 max k φ j ( 1 min k U(t + t) U(t) {u k(t)}) = max {u k(t)} {u k (t)}) = min{u k (t)} Example: Viscous Burger s equation t u+ x ( 1 2 u2 s x u) = 0 k k Remark: Explicit monotone schemes retain positivity! 30
32 Remark: Explicit monotone schemes can be interpreted as the solution of a particle method where particles jump only with multiples of x! Remark: Monotonicity works also for nonlinear scalar equations but not in general for systems!
33 INVERSE MONOTONE OPERATORS AND IMPLICIT SCHEMES Definition: B is inverse monotone (an M-matrix) if ( B jk 0 k B jk K > 0 j for j k (strictly diagonally dominant) ) Theorem: B is an M matrix B
34 Generalization: B is inverse monotone (an M-matrix) if B jk 0 k B jk 0 j α : k B αk > 0 for j k diagonally dominant) Remark: α usually given by the boundary conditions. 32
35 IMPLICIT SCHEMES (LINEAR) BU(t + t) = AU(t) B M-matrix and A > 0. B jk 0 k B jk 1 + tk B, K B > 0 j A jk 0 k A jk 1 tk A, K A 0, j for j k (strictly diagonally dominant) B 1 A 0 and U(t + t 1 tk A 1 + tk B U(t + t 33
36 Example: (damped heat equation) t u = 2 xu cu u(x, t+ t) = u(x, t)+ t x 2(T 2+T 1 )u(x, t + t) c tu(x, t + t) involves solving a system of equations at each time step write this as k k B jk U k (t + t) = k A jk U k (t) B jk = 1 + c t, B jj 0, B jk 0 for j k, A = I, U(t + t) c t U(t) This scheme requires the solution of equations at each time step, but is unconditionally stable (no CFL condition) 34
37 IMPLICIT SCHEMES (NONLINEAR) φ U φ(u(t + t)) = ψ(u(t)) M-matrix and ψ U > 0. φ j U 0 for j k k φ(c 1) (1 + tk φ )c, K φ > 0 c ψ j U 0 k ψ(c 1) c(1 tk ψ ), K ψ 0, c U(t + t 1 tk φ 1 + tk ψ U(t + t Example: Cranck - Nicholson (trapezoidal rule) t u = xu 2 u(t+ t) t 2 x 2(T 2+T 1 )u(t+ t) = u(t)+ t 2 x 2(T 2+T 1 )u(t) 35
38 ARTIFICIAL DISPERSION Consider t u + a x u = 0. Plane wave ansatz: u(x, t) = exp[iξ(x vt) rt] Exact: v = a, r = 0 Difference scheme: v = v(ξ) (incorrect wave speeds for high frequency components!) Not a problem in the presence of artificial diffusion (r > 0)! Example: Cranck Nicholson for t u + a x u = 0 will give the wrong wave speeds due the lack of artificial diffusion. 36
Numerical 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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationMATH 220: Problem Set 3 Solutions
MATH 220: Problem Set 3 Solutions Problem 1. Let ψ C() be given by: 0, x < 1, 1 + x, 1 < x < 0, ψ(x) = 1 x, 0 < x < 1, 0, x > 1, so that it verifies ψ 0, ψ(x) = 0 if x 1 and ψ(x)dx = 1. Consider (ψ j )
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationAPPLIED MATHEMATICS. Part 1: Ordinary Differential Equations. Wu-ting Tsai
APPLIED MATHEMATICS Part 1: Ordinary Differential Equations Contents 1 First Order Differential Equations 3 1.1 Basic Concepts and Ideas................... 4 1.2 Separable Differential Equations................
More informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationMath 2a Prac Lectures on Differential Equations
Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments
More informatione (x y)2 /4kt φ(y) dy, for t > 0. (4)
Math 24A October 26, 2 Viktor Grigoryan Heat equation: interpretation of the solution Last time we considered the IVP for the heat equation on the whole line { ut ku xx = ( < x
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationx ct x + t , and the characteristics for the associated transport equation would be given by the solution of the ode dx dt = 1 4. ξ = x + t 4.
. The solution is ( 2 e x+ct + e x ct) + 2c x+ct x ct sin(s)dx ( e x+ct + e x ct) + ( cos(x + ct) + cos(x ct)) 2 2c 2. To solve the PDE u xx 3u xt 4u tt =, you can first fact the differential operat to
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationNumerical discretization of tangent vectors of hyperbolic conservation laws.
Numerical discretization of tangent vectors of hyperbolic conservation laws. Michael Herty IGPM, RWTH Aachen www.sites.google.com/michaelherty joint work with Benedetto Piccoli MNCFF 2014, Bejing, 22.5.2014
More informationA New Wavelet-based Expansion of a Random Process
Columbia International Publishing Journal of Applied Mathematics and Statistics doi:10.7726/jams.2016.1011 esearch Article A New Wavelet-based Expansion of a andom Process Ievgen Turchyn 1* eceived: 28
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationA NEW NONLINEAR FILTER
COMMUNICATIONS IN INFORMATION AND SYSTEMS c 006 International Press Vol 6, No 3, pp 03-0, 006 004 A NEW NONLINEAR FILTER ROBERT J ELLIOTT AND SIMON HAYKIN Abstract A discrete time filter is constructed
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 informationModeling Rarefaction and Shock waves
April 30, 2013 Inroduction Inroduction Rarefaction and shock waves are combinations of two wave fronts created from the initial disturbance of the medium. Inroduction Rarefaction and shock waves are combinations
More informationRandom Variables. Random variables. A numerically valued map X of an outcome ω from a sample space Ω to the real line R
In probabilistic models, a random variable is a variable whose possible values are numerical outcomes of a random phenomenon. As a function or a map, it maps from an element (or an outcome) of a sample
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationNumerical integration and importance sampling
2 Numerical integration and importance sampling 2.1 Quadrature Consider the numerical evaluation of the integral I(a, b) = b a dx f(x) Rectangle rule: on small interval, construct interpolating function
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 2: Monte Carlo simulation methods Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
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 informationContinuous Random Variables and Continuous Distributions
Continuous Random Variables and Continuous Distributions Continuous Random Variables and Continuous Distributions Expectation & Variance of Continuous Random Variables ( 5.2) The Uniform Random Variable
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationA Lévy-Fokker-Planck equation: entropies and convergence to equilibrium
1/ 22 A Lévy-Fokker-Planck equation: entropies and convergence to equilibrium I. Gentil CEREMADE, Université Paris-Dauphine International Conference on stochastic Analysis and Applications Hammamet, Tunisia,
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
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 informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationWEYL S LEMMA, ONE OF MANY. Daniel W. Stroock
WEYL S LEMMA, ONE OF MANY Daniel W Stroock Abstract This note is a brief, and somewhat biased, account of the evolution of what people working in PDE s call Weyl s Lemma about the regularity of solutions
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 informationInverse Problems in the Bayesian Framework
Inverse Problems in the Bayesian Framework Daniela Calvetti Case Western Reserve University Cleveland, Ohio Raleigh, NC, July 2016 Bayes Formula Stochastic model: Two random variables X R n, B R m, where
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 informationLecture 2. Introduction to Differential Equations. Roman Kitsela. October 1, Roman Kitsela Lecture 2 October 1, / 25
Lecture 2 Introduction to Differential Equations Roman Kitsela October 1, 2018 Roman Kitsela Lecture 2 October 1, 2018 1 / 25 Quick announcements URL for the class website: http://www.math.ucsd.edu/~rkitsela/20d/
More 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 informationFirst order wave equations. Transport equation is conservation law with J = cu, u t + cu x = 0, < x <.
First order wave equations Transport equation is conservation law with J = cu, u t + cu x = 0, < x
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
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 informationPARTIAL DIFFERENTIAL EQUATIONS. Lecturer: D.M.A. Stuart MT 2007
PARTIAL DIFFERENTIAL EQUATIONS Lecturer: D.M.A. Stuart MT 2007 In addition to the sets of lecture notes written by previous lecturers ([1, 2]) the books [4, 7] are very good for the PDE topics in the course.
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More informationCLASSROOM NOTES PART II: SPECIAL TOPICS. APM526, Spring 2018 Last update: Apr 11
CLASSROOM NOTES PART II: SPECIAL TOPICS APM526, Spring 2018 Last update: Apr 11 1 Function Space Methods General Setting: Projection into finite dimensional subspaces t u = F (u), u(t = 0) = u I, F : B
More informationLecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University
Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationPartial Differential Equations, Winter 2015
Partial Differential Equations, Winter 215 Homework #2 Due: Thursday, February 12th, 215 1. (Chapter 2.1) Solve u xx + u xt 2u tt =, u(x, ) = φ(x), u t (x, ) = ψ(x). Hint: Factor the operator as we did
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More information2) Let X be a compact space. Prove that the space C(X) of continuous real-valued functions is a complete metric space.
University of Bergen General Functional Analysis Problems with solutions 6 ) Prove that is unique in any normed space. Solution of ) Let us suppose that there are 2 zeros and 2. Then = + 2 = 2 + = 2. 2)
More information2 Functions of random variables
2 Functions of random variables A basic statistical model for sample data is a collection of random variables X 1,..., X n. The data are summarised in terms of certain sample statistics, calculated as
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationFirst order Partial Differential equations
First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationMonte Carlo Solution of Integral Equations
1 Monte Carlo Solution of Integral Equations of the second kind Adam M. Johansen a.m.johansen@warwick.ac.uk www2.warwick.ac.uk/fac/sci/statistics/staff/academic/johansen/talks/ March 7th, 2011 MIRaW Monte
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationPartial Differential Equations (PDEs)
10.34 Numerical Methods Applied to Chemical Engineering Fall 2015 Final Exam Review Partial Differential Equations (PDEs) 1. Classification (a) Many PDEs encountered by chemical engineers are second order
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationIEOR 4703: Homework 2 Solutions
IEOR 4703: Homework 2 Solutions Exercises for which no programming is required Let U be uniformly distributed on the interval (0, 1); P (U x) = x, x (0, 1). We assume that your computer can sequentially
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationLecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225
Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225 Dr. Asmaa Al Themairi Assistant Professor a a Department of Mathematical sciences, University of Princess Nourah bint Abdulrahman,
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationS6880 #13. Variance Reduction Methods
S6880 #13 Variance Reduction Methods 1 Variance Reduction Methods Variance Reduction Methods 2 Importance Sampling Importance Sampling 3 Control Variates Control Variates Cauchy Example Revisited 4 Antithetic
More informationLINEAR SECOND-ORDER EQUATIONS
LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationBrownian Motion and Stochastic Calculus
ETHZ, Spring 17 D-MATH Prof Dr Martin Larsson Coordinator A Sepúlveda Brownian Motion and Stochastic Calculus Exercise sheet 6 Please hand in your solutions during exercise class or in your assistant s
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationThird part: finite difference schemes and numerical dispersion
Tird part: finite difference scemes and numerical dispersion BCAM - Basque Center for Applied Matematics Bilbao, Basque Country, Spain BCAM and UPV/EHU courses 2011-2012: Advanced aspects in applied matematics
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationLecture 1 Measure concentration
CSE 29: Learning Theory Fall 2006 Lecture Measure concentration Lecturer: Sanjoy Dasgupta Scribe: Nakul Verma, Aaron Arvey, and Paul Ruvolo. Concentration of measure: examples We start with some examples
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More information