Monte carlo methods for the solution of ODE's and PDE's
|
|
- Baldric McCoy
- 5 years ago
- Views:
Transcription
1 University of Washington ME-535 Computational Techniques in Mechanical Engineering Monte carlo methods for the solution of ODE's and PDE's Authors: Tony Enslow David Jonsson June 9, 2014
2 1 Introduction The Monte Carlo method is an extremely powerful numerical method for high-dimensional and complex problems. It provides the tools necessary to both visualize the process, as well as computing the solution of many engineering problems.although it is not very ecient for lower order dierential equations, it is especially ecient when the computational complexity grows exponentially with increasing dimensions (PDE). The method uses random number generation to approximate numerical solutions. 2 Background The Monte Carlo method, as we know it today, dates back to the 1940s in the development in the atomic bomb. The name came from the casino industry because of the randomness of gambling and its connection to probability. It was therefore named after a popular place in Monaco, renowned for its gambling. The method of using randomness in solving mathematical problems can be traced back to the eighteenth century, to estimate the true value of pi. Another widely used application of the method is also to estimate the area of amorphous shapes. This can be done by enclosing the shape in a unit square and sampling random points. Then by determining the amount of points inside the shape, one can determine the area of the shape by the ratio of points inside and outside of the shape boundaries.this method has also been found to be very useful in stochastic processes. In particular, to estimate investments in nance by evaluating nancial derivatives and investment risks by evaluating the probability of success. This method also has wide applications in computational physics, and physical chemistry due to the randomness of small particles often studied in these elds. It is clear that this method has a wide variety of applications. Brownian motion Brownian motion is the random movement of particles suspended in a uid. It is caused by the collision of larger particles with atoms or molecules in the uid, causing them to alter their movement. This motion is a limit to the stochastic process,the random walk, is often used in Monte Carlo methods. Albert Einstein predicted that the Brownian motion of a particle in uid is characterized by the diusion coecient D = k b T b where b is the Boltzmann constant. The general diusion equation is given by δu δt = D δu δx 2 And from there we can nd the root mean square displacement at time t disp = 2Dt Due to the relationship of Brownian motion and the random walks, this motion is particularly important in many applications of Monte Carlo methods, including the solution of dierential equations 2,3.
3 Random Walks In many cases of engineering, closed form equations can often be hard to derive for a given problem. When this is the case, using a stochastic approach is often an attractive alternative. Using random numbers to calculate the behavior of the problem often gives a good approximation of the solution. A good example of this is the random walk problem of a drunk person walking home from a bar. If we simulate one walk, where a step in any direction can be taken with the same probability, we might end far from the bar or close to the starting point. On the other hand, if we simulate this walk many times, the average distance should give an expected quantitative result to the average displacement of the individual. This is called the method of random walks. The probability of each step will not always be equal and the use of Brownian motion will also come into play where a collision will alter your movement, but as can be seen in the examples below, this method is at the core of many Monte Carlo applications. 3 Method and Precision Any algorithm that uses random numbers to get a quantitative solution to a problem, is a Monte Carlo method. The idea is to initiate a random walk within a domain of a given equation. Probability of moving in each direction is then dened by the dierential equation that is being solved. The walk continues until the particle reaches the boundary of the domain. This can often be a good approximation, given that a sucient number of simulations are run. For low dimensional equations, this is not an ecient method, but it becomes increasingly ecient as the dimensions increase, compared to many other numerical methods. Although there are no general predened steps for a Monte Carlo method, most of them follow a general 4 step trend 1. Estimate the proper range of inputs. For example, if we are trying to estimate the density of air during a process, we would never input a negative number. 2. Generate random numbers in the desired range using either a uniform or non-uniform distribution of randomness. 3. Derive a deterministic function for the given problem, and use the random numbers as an input 4. Iterate and use either a ratio or the average value of all iterations. Continue to iterate until the solution becomes stable. Generally, the precision of the Monte Carlo method is strictly based on the number of iterations. Given a function f(x) with uniform random distribution for the value x 0 x 1, the error can be estimated by error = 1 N N f(x i ) I i=1 Where I is the estimated quantity. Simplifying this equation gives the error approximation error σ f N η(0, 1) where σf 2 = (f(x) I) 2 dx ω
4 and η(0, 1) is the standard random variable. It is worth noting that the Monte Carlo 1 method converges with the rate of N which is rather slow compared to other methods. 4 ODEs The Monte Carlo method is not ecient for standard ODEs with a known analytical solution, as was stated above. It becomes very good for stochastic dierential equations, higher dimension equations, and equations with special conditions or unknown equations. Due to the scope of this paper, this chapter will only cover the last case but the other two cases were briey mentioned above. For an equation with special conditions, we imagine an area with a boundary conditions. We start from some point inside the area and walk towards the boundary. At each step we can go in some direction and the probability of each direction is controlled by the dierential equation at that point. When the boundary is reached, the particle is absorbed or eliminated, depending on the boundary condition. We perform this movement n number of times using uniformly randomized set of numbers in a chosen interval and after those n walks, we will have the solution to our dierential equation inside that perimeter. Example 1 For clarity, let's do an engineering example for a colloid thruster, to demonstrate the application of the Monte Carlo method in the engineering realm. A colloid thruster is a space propulsion system that creates an electrospray of charged liquid droplets to produce propulsion. The problem with solving this system analtically is that the movement of the particles is random in nature, and therefore we cannot assume that they all travel in the same direction during propulsion. Consider the trajectory of a particle inside the electric eld, 0 x L,that can be described using a set of dierential equations dx dt = v x dy dt = v y dv x dt = A x(x) dv r dt = 0 Where v x is the velocity in the x direction and v y is the velocity in the y direction and all particles have dierent, random initial velocity distribution. A x = e dφ(x) is the acceleration induced by the electric eld with the boundary conditions for x m dx { A x0 (1 3x 2 + 2x 3 ), if x L. A x (x) = 0, if x > L. With the initial position of particles x(0) = y(0) = 0 and an initial velocity v x (0) = V 0 cosα 0 and v y (0) = V 0 sinα 0 where α is the angle of angle of the particle. When solving the sets of equation we assume that α is a uniform random variable. We run the code in appendix A and get the average velocities, v x, v y after the particle leaves the electric eld, given in chosen units v x = ± , v y = ±
5 When solving the set of equation we rst generated a random sample of the departure angle with a xed velocity. Then we analyzed every computation for each of the random samples. Finally we found the average velocities and their standard deviation. 5 PDEs Partial dierential equations are very important in many applications of engineering. Monte Carlo methods often give a good solution to these equations. As PDE's can often become very complex and dicult to solve analytically, applying a randomized method such as the Monte Carlo method, can be a good approach to solving them. In order to show the solution process, let's consider two examples. Example 1 We consider the Dirichlet solution to the Laplace equation. Let there be a bounded, connected, region S and a function F(x) which satises Laplace's equation everywhere inside of S where φ is given at the boundary k j=1 = δ2 φ j δx 2 j φ(x) = f(q) = 0 (1) Next we set up the problem by drawing a square around S and splitting it up into a grid. Then we set up our random walks in two dimensions. We start from the center and walk towards the boundary. We solve the dierential equation (1) for two dimensions, by substituting the known dierence equations into the ODEs, and get Which yields 0 = δ2 φ δx + δ2 φ 2 δy φ(p 1) + φ(p 2 ) + φ(p 3 ) + φ(p 4 ) 2 h 2 φ(p ) = φ(p 1) + φ(p 2 ) + φ(p 3 ) + φ(p 4 ) 4φ(P ) 4 From equation 2 we can simulate the solution to the problem. We walk from the center towards the boundary with a choice of going to 4 surrounding nodes (left,right,up,down). (2)
6 After N simulation of random walks, all ending somewhere at the boundary, the solution can be approximated using φ(x, y) N n i f(q i ) i=1 where n is the number of walks terminated at the boundary node Q i. Random walks do not necessarily have an equal probability to walk in each direction. Let's consider a general elliptic dierential equation in the same region S. Which has the dierence equation N δ 2 φ β 11 δx + 2β δ 2 φ 2 12 δxδy + φ β22δ2 δy + 2α δφ 2 1 δx + 2α δφ 2 δy = 0 (4) β 11 xx φ + 2β 12 xy φ + β 22 yy φ + 2α 1 x φ + 2α 2 y φ = 0 We approximate the dierence equation for ve points around (x,y) where each point is a possible destination for the next step. The probabilities for each point are found by deriving the dierence equation and arriving at the function for φ φ(x, y) = φ(x + h, y)p 1 + φ(x, y + h)p 2 + φ(x h, y)p 3 + φ(x, y h)p 4 + φ(x + h, y + h)p 5 Where Where P 1 = β α 1 h 2β 12, P 2 = β α 2 h 2β 12 D D P 3 = β 11 D, P 4 = β 22 D, P 5 = 2β12 D D = 2β11 2β β (α 1 + α 2 )h The solution from these random walks is then nally found using equation 3 above. We have now seen how the Monte Carlo method can be easily applied to two dimensional problems. Let's look at one of the most important partial dierential equation used by engineers; the diusion equation. Since the basic model for the diusion equation of heat uses the idea of random spread of particles in all direction, the Monte Carlo method is good for acquiring the solution as well as giving some insight into the individual movement of each particle. (3) Example 2 Let's consider the one dimensional diusion equation δu δt = δu δx 2 u(x, 0) = u 0 (x) Which we can describe with a the dierence equation u n+1 i = (1 2λ)u n i + λ(u n i+1 + u n i 1 (5) where u n i is the approximation to u(i x, n t). λ = t. The random walk of each x 2 heat particle can then be described by the dierence equation(5) for each time point. Each
7 particle has three possible actions, stay put, go left, and go right. The probability of each option is given respectively P 1 = 1 2λ, P 2 = λ, P 3 = λ Let's now apply this general approach to particles in a uid container. We put particles into a container lled with a dierent substance. We observe their behavior and their relative concentration in the uid and model it using the Monte Carlo method described above. We then compare the solution to the analytical solution of the problem 11. Figure 1: Particles in uid We see that the Monte Carlo method for time steps and random walks gives a good approximation to the diusion equation. The Matlab code for this gure is given in appendix A. 6 Conclusion The Monte Carlo methods demonstrate a good solution for many problems in engineering. In many cases it can be dicult to nd the closed form equation of a process. This was shown in the example of a thrusting power of a colloid thruster, where the particle departure angle is a random process which greatly aects the solution to the problem. We furthermore explored the domain of partial dierential equations where we showed the ecacy of the Monte Carlo method for solving boundary value problems using random walks. We then solved the one dimensional heat equation and demonstrating the usefulness of this method. For higher order equations, this method becomes increasingly attractive compared to other methods, especially when an analytical solution is hard to derive.
8 7 Reference 1. Lecture notes on numerical methods for ODE's, MIT, 2012 http : //web.mit.edu/course/16/16.90/bac 2. Petar KoÄ oviä, Fakultet za obrazovanje rukovodecih kadrova u privredi, Novi Sad,Brownian motion development for Monte Carlo method applied on european style option price forecasting,2011, International journar of economics and law 3. Brownian motion http : //en.wikipedia.org/wiki/brownian_motion 4. Anders à sterling, Diusion equation and Monte Carlo. 5. David Bindel and Jonathan Goodman,2006, Principles of Scientiï c Computing Monte Carlo methods 6. Bernard Lapeyre,Introduction to Monte-Carlo Methods, Monte Carlo methodshttp : //en.wikipedia.org/wiki/monte C arlo m ethod 8. Branislav K. Nikoli Ä, From Random Numbers to Monte Carlo 9. Rubin H Landau, Random Walks 10. Anne GILLE-GENEST,2012,Introduction to the Monte Carlo Methods 11. Moroko, William J., and Russel E. Caisch. "A quasi-monte Carlo approach to particle simulation of the heat equation." SIAM Journal on Numerical Analysis 30.6 (1993): Peter M orters and Yuval Peres, 2008,Brownian Motion 13. Gregory F. Lawler,Random Walk and the Heat Equation 14. Matthew. N. O. Sadiku et.al, 2007, Direct Monte Carlo simulation of time-dependent problems 15. Dr. Eric Ayars, Computational Physics With Python
9 Appendix A Matlab code for colloid thruster 1 function [vx, vr] = thruster(v0, a0, Ax0, L, dt) % initial condition x = 0; r = 0; vx = V0*cos(a0); vr = V0*sin(a0); % time integration using Forward Euler while (x â L) x = x + dt*vx; r = r + dt*vr; vx = vx + dt*ax0*(1-3*x.ë 2 + 2*xx.Ë 3); end % Deterministic (non-random) parameters V0 = 0.1; Ax = 1.0; L = 1.0; dt = 0.001; a0max = 60. * pi / 180.; % 1. Generate a random sample of the input parameters a0mc = rand(10000, 1) * a0max; % Array to store Monte Carlo outputs vxmc = []; vrmc = []; % 2. Analyze (deterministically) each set of inputs in the sample for i = 1:10000 [vx, vr] = thruster(v0, a0mc(i), Ax, L, dt); vxmc = [vxmc; vx]; vrmc = [vrmc; vr]; forward_propulsion(i) = vx; side_propulsion(i) = vr; end % mean muvx = mean(vxmc); muvr = mean(vrmc); % standard deviation sigmavx = std(vxmc); sigmavr = std(vrmc);
10 Matlab code for particles in uid %========================================================================= %******** To do Monte Carlo simulation of particle diffusion ************* %************************ By Dr. Mahesha MG ****************************** % Date: 17/11/2012 %========================================================================= clear clc a1=fopen('particle.dat','w'); %Data will be stored in this file t=0:50000; nop=10000; % Total number of particles x=[nop]; nl=nop; % Asssumed that all particles are on left side of the box at t=0 for i=1:50000 lambda = nl/nop; if rand()<lambda nl=nl-1; else nl=nl+1; end fprintf(a1,'%d\t%d\n',t(i),nl); x=[x;nl]; end y=(nop/2)*(1+exp(-2*t/nop)); %Analytic solution plot(t,x,'b',t,y,'r') legend('monte-carlo','analytic') xlabel('time in s') ylabel('no. of particles on left side') fclose(a1);
Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume Method
Ninth International Conference on Computational Fluid Dynamics (ICCFD9), Istanbul, Turkey, July 11-15, 2016 ICCFD9-0113 Solution of the Two-Dimensional Steady State Heat Conduction using the Finite Volume
More informationFrom Random Numbers to Monte Carlo. Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that
From Random Numbers to Monte Carlo Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that Random Walk Through Life Random Walk Through Life If you flip the coin 5 times you will
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More informationMath 345 Intro to Math Biology Lecture 21: Diffusion
Math 345 Intro to Math Biology Lecture 21: Diffusion Junping Shi College of William and Mary November 12, 2018 Functions of several variables z = f (x, y): two variables, one function value Domain: a subset
More informationPrancing Through Quantum Fields
November 23, 2009 1 Introduction Disclaimer Review of Quantum Mechanics 2 Quantum Theory Of... Fields? Basic Philosophy 3 Field Quantization Classical Fields Field Quantization 4 Intuitive Field Theory
More informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationPROPOSITION of PROJECTS necessary to get credit in the 'Introduction to Parallel Programing in 2014'
PROPOSITION of PROJECTS necessary to get credit in the 'Introduction to Parallel Programing in 2014' GENERALITIES The text may be written in Polish, or English. The text has to contain two parts: 1. Theoretical
More informationt x 0.25
Journal of ELECTRICAL ENGINEERING, VOL. 52, NO. /s, 2, 48{52 COMPARISON OF BROYDEN AND NEWTON METHODS FOR SOLVING NONLINEAR PARABOLIC EQUATIONS Ivan Cimrak If the time discretization of a nonlinear parabolic
More informationCNH3C3 Persamaan Diferensial Parsial (The art of Modeling PDEs) DR. PUTU HARRY GUNAWAN
CNH3C3 Persamaan Diferensial Parsial (The art of Modeling PDEs) DR. PUTU HARRY GUNAWAN Partial Differential Equations Content 1. Part II: Derivation of PDE in Brownian Motion PART II DERIVATION OF PDE
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationMonte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time
Monte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time DESY Summer Student Programme, 214 Ronnie Rodgers University of Oxford, United Kingdom Laura Raes University of
More informationG : Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into
G25.2651: Statistical Mechanics Notes for Lecture 3 I. MICROCANONICAL ENSEMBLE: CONDITIONS FOR THERMAL EQUILIBRIUM Consider bringing two systems into thermal contact. By thermal contact, we mean that the
More informationApplied Mathematics 505b January 22, Today Denitions, survey of applications. Denition A PDE is an equation of the form F x 1 ;x 2 ::::;x n ;~u
Applied Mathematics 505b January 22, 1998 1 Applied Mathematics 505b Partial Dierential Equations January 22, 1998 Text: Sobolev, Partial Dierentail Equations of Mathematical Physics available at bookstore
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationA slow transient diusion in a drifted stable potential
A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive
More informationMarch 16, Abstract. We study the problem of portfolio optimization under the \drawdown constraint" that the
ON PORTFOLIO OPTIMIZATION UNDER \DRAWDOWN" CONSTRAINTS JAKSA CVITANIC IOANNIS KARATZAS y March 6, 994 Abstract We study the problem of portfolio optimization under the \drawdown constraint" that the wealth
More informationPOTENTIAL THEORY AND HEAT CONDUCTION DIRICHLET S PROBLEM
Chapter 6 POTENTIAL THEORY AND HEAT CONDUCTION DIRICHLET S PROBLEM M. Ragheb 9/19/13 6.1 INTRODUCTION The solution of the Dirichlet problem is one of the easiest approaches to grasp using Monte Carlo methodologies.
More informationNovel determination of dierential-equation solutions: universal approximation method
Journal of Computational and Applied Mathematics 146 (2002) 443 457 www.elsevier.com/locate/cam Novel determination of dierential-equation solutions: universal approximation method Thananchai Leephakpreeda
More informationMATH 294???? FINAL # 4 294UFQ4.tex Find the general solution y(x) of each of the following ODE's a) y 0 = cosecy
3.1. 1 ST ORDER ODES 1 3.1 1 st Order ODEs MATH 294???? FINAL # 4 294UFQ4.tex 3.1.1 Find the general solution y(x) of each of the following ODE's a) y 0 = cosecy MATH 294 FALL 1990 PRELIM 2 # 4 294FA90P2Q4.tex
More informationLecture 21: Stochastic Differential Equations
: Stochastic Differential Equations In this lecture, we study stochastic differential equations. See Chapter 9 of [3] for a thorough treatment of the materials in this section. 1. Stochastic differential
More informationMath Camp Notes: Everything Else
Math Camp Notes: Everything Else Systems of Dierential Equations Consider the general two-equation system of dierential equations: Steady States ẋ = f(x, y ẏ = g(x, y Just as before, we can nd the steady
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 information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationSimple Harmonic Motion
Simple Harmonic Motion (FIZ 101E - Summer 2018) July 29, 2018 Contents 1 Introduction 2 2 The Spring-Mass System 2 3 The Energy in SHM 5 4 The Simple Pendulum 6 5 The Physical Pendulum 8 6 The Damped Oscillations
More informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationMat1062: Introductory Numerical Methods for PDE
Mat1062: Introductory Numerical Methods for PDE Mary Pugh January 13, 2009 1 Ownership These notes are the joint property of Rob Almgren and Mary Pugh 2 Boundary Conditions We now want to discuss in detail
More informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More information1 The relation between a second order linear ode and a system of two rst order linear odes
Math 1280 Spring, 2010 1 The relation between a second order linear ode and a system of two rst order linear odes In Chapter 3 of the text you learn to solve some second order linear ode's, such as x 00
More informationOn reaching head-to-tail ratios for balanced and unbalanced coins
Journal of Statistical Planning and Inference 0 (00) 0 0 www.elsevier.com/locate/jspi On reaching head-to-tail ratios for balanced and unbalanced coins Tamas Lengyel Department of Mathematics, Occidental
More informationNotes on Iterated Expectations Stephen Morris February 2002
Notes on Iterated Expectations Stephen Morris February 2002 1. Introduction Consider the following sequence of numbers. Individual 1's expectation of random variable X; individual 2's expectation of individual
More informationf(z)dz = 0. P dx + Qdy = D u dx v dy + i u dy + v dx. dxdy + i x = v
MA525 ON CAUCHY'S THEOREM AND GREEN'S THEOREM DAVID DRASIN (EDITED BY JOSIAH YODER) 1. Introduction No doubt the most important result in this course is Cauchy's theorem. Every critical theorem in the
More informationMonte Carlo Lecture Notes II, Jonathan Goodman. Courant Institute of Mathematical Sciences, NYU. January 29, 1997
Monte Carlo Lecture Notes II, Jonathan Goodman Courant Institute of Mathematical Sciences, NYU January 29, 1997 1 Introduction to Direct Sampling We think of the computer \random number generator" as an
More informationIn this activity, we explore the application of differential equations to the real world as applied to projectile motion.
Applications of Calculus: Projectile Motion ID: XXXX Name Class In this activity, we explore the application of differential equations to the real world as applied to projectile motion. Open the file CalcActXX_Projectile_Motion_EN.tns
More informationA j = 0.1 cm 2 10 cm 10 cm 10 cm. W j Wj. W j W j. W j. 10 cm 10 cm 10 cm. r i
ME 131B Fluid Mechanics Solutions to Week Eight Problem Session: Angular Momentum Principle (3/2/98) 1. In control volume analysis, all governing principles share the same common structure: storage = inow
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 informationREVIEW: Waves on a String
Lecture 14: Solution to the Wave Equation (Chapter 6) and Random Walks (Chapter 7) 1 Description of Wave Motion REVIEW: Waves on a String We are all familiar with the motion of a transverse wave pulse
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More 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 informationThe WENO Method for Non-Equidistant Meshes
The WENO Method for Non-Equidistant Meshes Philip Rupp September 11, 01, Berlin Contents 1 Introduction 1.1 Settings and Conditions...................... The WENO Schemes 4.1 The Interpolation Problem.....................
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationAppendix on complex numbers:
Appendix A Appendix on complex numbers: A. complex numbers We begin with a review of several properties of complex numbers, their representation, and some of their basic properties. The use of complex
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 informationMA 102 Mathematics II Lecture Feb, 2015
MA 102 Mathematics II Lecture 1 20 Feb, 2015 Differential Equations An equation containing derivatives is called a differential equation. The origin of differential equations Many of the laws of nature
More informationX. Hu, R. Shonkwiler, and M.C. Spruill. School of Mathematics. Georgia Institute of Technology. Atlanta, GA 30332
Approximate Speedup by Independent Identical Processing. Hu, R. Shonkwiler, and M.C. Spruill School of Mathematics Georgia Institute of Technology Atlanta, GA 30332 Running head: Parallel iip Methods Mail
More informationThe Definite Integral. Day 6 Motion Problems Strategies for Finding Total Area
The Definite Integral Day 6 Motion Problems Strategies for Finding Total Area ARRIVAL---HW Questions Working in PODS Additional Practice Packet p. 13 and 14 Make good use of your time! Practice makes perfect!
More informationAPM526 CLASSROOM NOTES PART IB. APM526, Spring 2018 Last update: MAR 12
APM526 CLASSROOM NOTES PART IB APM526, Spring 2018 Last update: MAR 12 1 PARTICLE METHODS 1. Particle based methods: Large systems of ODEs for individual particles. 2. Partial differential equations: Equations
More informationCSCI1950V Project 4 : Smoothed Particle Hydrodynamics
CSCI1950V Project 4 : Smoothed Particle Hydrodynamics Due Date : Midnight, Friday March 23 1 Background For this project you will implement a uid simulation using Smoothed Particle Hydrodynamics (SPH).
More informationThe Rocket Equation. Lukas Lundin. 26th January 2016
The Rocket Equation Lukas Lundin 26th January 2016 Abstract In this project we study the basics of rocket propulsion and rocket motion in the vicinity of the Earth. Furthermore we will compare dierent
More informationVorticity and Potential Vorticity
Chapter 4 Vorticity and Potential Vorticity In this chapter we explore a way of solving the shallow water equations for motions with characteristic time scales longer than the rotation period of the earth.
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
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 Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More information14.30 Introduction to Statistical Methods in Economics Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 4.30 Introduction to Statistical Methods in Economics Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationDiscrete particle settling. Flocculent settling. Compression
An introduction to sedimentation theory in wastewater treatment Bengt Carlsson Systems and Control Group Uppsala University Nov 96, rev Okt 98 Abstract This material is made for the course \Wastewater
More informationUNIVERSITY OF MARYLAND Department of Physics College Park, Maryland. PHYSICS Ph.D. QUALIFYING EXAMINATION PART II
UNIVERSITY OF MARYLAND Department of Physics College Park, Maryland PHYSICS Ph.D. QUALIFYING EXAMINATION PART II August 23, 208 9:00 a.m. :00 p.m. Do any four problems. Each problem is worth 25 points.
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 informationPartial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction
Page of 8 Partial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction FILE:Chap 3 Partial Differential Equations-V6. Original: May 7, 05 Revised: Dec 9, 06, Feb 0, 07,
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 informationvestigated. It is well known that the initial transition to convection is accompanied by the appearance of (almost) regular patterns. s the temperatur
Symmetry of ttractors and the Karhunen-Loeve Decomposition Michael Dellnitz Department of Mathematics University of Hamburg D-2046 Hamburg Martin Golubitsky Department of Mathematics University of Houston
More informationMonte Carlo methods for Statistical Inference: Generation of Random Numbers
Monte Carlo methods for Statistical Inference: Generation of Random Numbers Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 25th February 2004 Meet at NS 104 On Wednesday
More informationTime Integration Methods for the Heat Equation
Time Integration Methods for the Heat Equation Tobias Köppl - JASS March 2008 Heat Equation: t u u = 0 Preface This paper is a short summary of my talk about the topic: Time Integration Methods for the
More informationThermodynamics. Basic concepts. Thermal equilibrium and temperature
hermodynamics Basic concepts hermodynamics is a phenomenological description of macroscopic systems with many degrees of freedom. Its microscopic justication is provided by statistical mechanics. Equilibrium
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationMIXED BOUNDARY VALUE PROBLEMS
MIXED BOUNDARY VALUE PROBLEMS M. Ragheb 9/19/013 INTRODUCTION Mixed boundary value problems are practical situations that are met in most potential and other mathematical physics problems. In this case
More informationSimulation of Plume-Spacecraft Interaction
Simulation of Plume-Spacecraft Interaction MATÍAS WARTELSKI Master of Science Thesis in Space Physics Examiner: Prof. Lars Blomberg Space and Plasma Physics School of Electrical Engineering Royal Institute
More informationCarnegie Mellon University Forbes Ave. Pittsburgh, PA 15213, USA. fmunos, leemon, V (x)ln + max. cost functional [3].
Gradient Descent Approaches to Neural-Net-Based Solutions of the Hamilton-Jacobi-Bellman Equation Remi Munos, Leemon C. Baird and Andrew W. Moore Robotics Institute and Computer Science Department, Carnegie
More informationMinimum and maximum values *
OpenStax-CNX module: m17417 1 Minimum and maximum values * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 In general context, a
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 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 informationENGI Partial Differentiation Page y f x
ENGI 3424 4 Partial Differentiation Page 4-01 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
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 informationKinematics and One Dimensional Motion
Kinematics and One Dimensional Motion Kinematics Vocabulary Kinema means movement Mathematical description of motion Position Time Interval Displacement Velocity; absolute value: speed Acceleration Averages
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationPartial Differential Equations and Image Processing. Eshed Ohn-Bar
Partial Differential Equations and Image Processing Eshed Ohn-Bar OBJECTIVES In this presentation you will 1) Learn what partial differential equations are and where do they arise 2) Learn how to discretize
More informationPathwise Construction of Stochastic Integrals
Pathwise Construction of Stochastic Integrals Marcel Nutz First version: August 14, 211. This version: June 12, 212. Abstract We propose a method to construct the stochastic integral simultaneously under
More informationChapter 4: Monte Carlo Methods. Paisan Nakmahachalasint
Chapter 4: Monte Carlo Methods Paisan Nakmahachalasint Introduction Monte Carlo Methods are a class of computational algorithms that rely on repeated random sampling to compute their results. Monte Carlo
More informationEconomics 2010c: Lectures 9-10 Bellman Equation in Continuous Time
Economics 2010c: Lectures 9-10 Bellman Equation in Continuous Time David Laibson 9/30/2014 Outline Lectures 9-10: 9.1 Continuous-time Bellman Equation 9.2 Application: Merton s Problem 9.3 Application:
More informationName: ID.NO: Fall 97. PLEASE, BE NEAT AND SHOW ALL YOUR WORK; CIRCLE YOUR ANSWER. NO NOTES, BOOKS, CALCULATORS, TAPE PLAYERS, or COMPUTERS.
MATH 303-2/6/97 FINAL EXAM - Alternate WILKERSON SECTION Fall 97 Name: ID.NO: PLEASE, BE NEAT AND SHOW ALL YOUR WORK; CIRCLE YOUR ANSWER. NO NOTES, BOOKS, CALCULATORS, TAPE PLAYERS, or COMPUTERS. Problem
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page of 6 Exercise 2 Solve the completely inhomogeneous diffusion problem on the half-line v t kv xx = f(x, t) for < x
More informationMARKOV CHAINS: STATIONARY DISTRIBUTIONS AND FUNCTIONS ON STATE SPACES. Contents
MARKOV CHAINS: STATIONARY DISTRIBUTIONS AND FUNCTIONS ON STATE SPACES JAMES READY Abstract. In this paper, we rst introduce the concepts of Markov Chains and their stationary distributions. We then discuss
More informationThe One Dimensional Heat Equation
The One Dimensional Heat Equation Adam Abrahamsen and David Richards May 22, 2002 Abstract In this document we will study the flow of heat in one dimension through a small thin rod. We will use the derivation
More informationOutline. 1 Why PIC Simulations Make Sense. 2 The Algorithm. 3 Examples Surface High Harmonics Generation. 4 Extensions Of The PIC Algorithm
PIC Simulations an Introduction GRK 1203 Meeting February 12-15, 2008, Oelde Outline 1 Simulations Make Sense 2 3 Surface High Harmonics Generation 4 Of PIC Plasma Physics Is Complex Experiment real thing
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 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 informationREDUNDANT TRINOMIALS FOR FINITE FIELDS OF CHARACTERISTIC 2
REDUNDANT TRINOMIALS FOR FINITE FIELDS OF CHARACTERISTIC 2 CHRISTOPHE DOCHE Abstract. In this paper we introduce so-called redundant trinomials to represent elements of nite elds of characteristic 2. The
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More informationNumerical Schemes from the Perspective of Consensus
Numerical Schemes from the Perspective of Consensus Exploring Connections between Agreement Problems and PDEs Department of Electrical Engineering and Computer Sciences University of California, Berkeley
More informationwhere (E) is the partition function of the uniform ensemble. Recalling that we have (E) = E (E) (E) i = ij x (E) j E = ij ln (E) E = k ij ~ S E = kt i
G25.265: Statistical Mechanics Notes for Lecture 4 I. THE CLASSICAL VIRIAL THEOREM (MICROCANONICAL DERIVATION) Consider a system with Hamiltonian H(x). Let x i and x j be specic components of the phase
More informationProceedings of Eight SIAM Conference on Parallel Processing for Scientific Computing, 10 pages, CD-ROM Format, Minneapolis, Minnesota, March 14-17,
Three Dimensional Monte Carlo Device imulation with Parallel Multigrid olver Can K. andalc 1 C. K. Koc 1. M. Goodnick 2 Abstract We present the results in embedding a multigrid solver for Poisson's equation
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More information