PHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison
|
|
- Carol Price
- 5 years ago
- Views:
Transcription
1 INTRODUCTION PHY37F/47F - Coputational Physics Background Material for Expt 3 - Heat Equation David Harrison In the Pendulu Experient, we studied the Runge-Kutta algorith for solving ordinary differential equations (ODEs) Here, we study techniques for solving partial differential equations (PDEs) The general proble of solving PDEs is huge; Press et al in Nuerical Recipes clai it requires an entire book However, there are soe siple cases Separation of Variables: Consider the Schrödinger equation: h / Ψ(x,t) x +V(x)Ψ(x,t)ih Ψ(x,t) t () If we assue that the wavefunction can be written as a function of position ties a function of tie: then: Ψ(x,t) ψ(x)φ(t) h d / ψ(x) ψ(x) dx +V(x)ψ(x) ih φ(t) dφ c () dt where c is a constant Thus we have turned the proble into two ordinary differential equations First order PDEs: If the equations are linear, then a closed solution can easily be found If the equations are non-linear, a faily of paraetric solutions can be found The standard Matheatica package Calculus PDSolve handles this latter case Second order PDEs: There are two general classes of these: Poisson, soeties called "elliptic" 7
2 -- Cauchy, soeties called "initial value" The Poisson equation is, of course: u u + x y ρ(x,t) (3) This type of equation turns out to be very easy to solve using a technique called a relaxation ethod The Cauchy class of second order PDEs is subdivided further: i The wave equation (a hyperbolic for): ii u v u t x The heat equation (a parabolic for): (4) u u c (5) t x x For both of these, you will discover that stability is a proble: any reasonable looking algoriths just don t work unless care is taken In this Experient, we will concentrate on the heat equation Usually, u is the teperature We will assue that we are solving the equation for a one diensional slab of width L We will usually assue that c is a constant so the heat equation becoes: u(x,t) u(x,t) c t x We will adopt units where x/l x and tc/l t, so the heat equation is now: u(x,t) u(x,t) t x x These notes are organised as follows: I Explicit algorith II III Iplicit algorith Crank-Nicholson algorith (6)
3 -3- IV Schrödinger s equation V References Note that the "usual" code listing is not included in this docuent That is because you will be odifying soe of the code to produce a new procedure, and cutting and pasting the code within the notebook is the siplest way to do this Thus, the code listing appears in the notebook instead of these suppleentary notes for the Experient I EXPLICIT ALGORITHM We iagine that we divide our one diensional slab into n interior points, and will call the distance between each point x We will step through tiesteps We will be using the notation that: u j u(x j j x,t ) j n+ (n+) x,,, (I) (tie ) n- n- n n+ j (position j x) We will write the initial condition as: u j g (x) (I) and the boundary conditions as: u α β u n+ It is reasonable to write the left hand side of the heat equation, Equation (6), as: (I3a) (I3b)
4 -4- u u + j u j t We write the right hand side of Equation (6) as: u x xj u x xj + u x xj (I4) / x (I5) Note that we write u/ x evaluated at two points that are not on our grid This doesn t atter because: Thus, Equation (I5) becoes: u x xj + u x xj u(x,t) x xj u j+ x u j u j u j x x u j+ u j +u j (I6) The above way of solving the second-order partial derivative is called the ethod of finite differences so: The heat equation states that Equation (I4) equals Equation (I6): u j + u j x u j+ u j +u j (I7) u + j u j + x u j+ u j +u j (I8) Note that the right hand side of Equation (I8) contains values of u at tie t, and these ters are cobined to give the value of u at the later tie t + (+) Once the values of u at tie t + (+) are known we can then find the values at t + (+), and so on Thus we can solve the heat equation for all values of t using this explicit algorith The deterination of u j + depends on knowing the values of u at three positions at the earlier tie:
5 -5- + j- j j+ Note that this schee does not allow one to deterine u + because we have no values for u But that does not atter: u + is the value of u at the left side of the slab, and is given by the boundary condition Equation (I3a) Siilarly the boundary condition Equation (I3b) gives the value for u + n+ So Equation (I8) is used only to evaluate the interior values of u + The above way of solving the heat equation is pretty siple Of the three algoriths you will investigate to solve the heat equation, this one is also the fastest and also can give the ost accurate result However, the result will be accurate only if you choose tiesteps and a space grid size x such that: ( x) µ Solutions produced in violation of Equation (I9) will be unstable, often producing ridiculous results Note that this eans that for a given tiestep the space grid size x can not be too sall For a given space grid size, the tiestep can not be too big We will discuss later in this docuent the fact that even if the tiestep and space grid satisfy Equation (I9), the explicit algorith can still produce wrong results for soe physical systes The reasons for the instability of the explicit algorith is soewhat beyond the level of this course; consult the references for further inforation However, we can justify Equation (I9) (I9) Golub and Ortega, 8
6 -6- The heat equation also governs the diffusion of, say, a sall quantity of perfue in the air You probably already know that diffusion is a for of rando walk so after a tie t we expect the perfue has diffused a distance x t One solution to the heat equation gives the density of the gas as a function of position and tie: u (x,t) ρ(x,t) e x σ σ (I) where: σ ct and c is the heat constant defined in Equation (5), for which we usually have been choosing units so that it equals In a diffusion context c is often called the diffusion constant Equation (I) shows that the density of the diffusing gas is a Gaussian, and that the standard deviation describing the width of the Gaussian increases as the square root of t This eans that in a tie, the olecules travel a distance x on the order of: x c Thus, if we are solving the equations using an explicit algorith, for a given : x (I) since otherwise the distribution is not spreading as fast as we know it should A trivial rearrangeent of Equation (I) gives: which is just Equation (I9) II IMPLICIT ALGORITHM ( x) Consider the following equation and copare it to Equation (I7) fro the previous section: u + u j ( x) (u + j+ u + j +u + j ) (II) Here we are evaluating u/ x using finite differences at tie t + instead of t A oent s reflection should convince you that Equation (II) is just as reasonable as Equation (I7) We rearrange Equation (II):
7 -7- (+µ)u + j µ(u + j+ +u + j ) u j (II) where recall that: µ ( x) Thus we are trying to find u at ties + at three different positions j-, j and j+ fro a single value of u at the earlier tie + j- j j+ At first glance, this is ipossible However, consider an (n+) (n+) identity atrix Id: Id an (n+) (n+) (,-) triangular atrix A: A and an n+ diensional colun atrix b :
8 -8- b µα µβ Then the entire set of equations for all values of j given by Equation (II) can be written as: (Id + µa)u + u +b (II3) where: u u u u u n u n+ If we know u, the values of u at tie, Equation (II3) ay be solved for u + using, for exaple, the LU decoposition that you studied in the Exercise The above technique for solving the heat equation is called an iplicit algorith You will discover in the Experient that this algorith is stable for all values of µ That is the good news There is also soe bad news: the iplicit technique is uch slower than the explicit one and is also uch less accurate for a given value of tiestep and space grid size x III CRANK-NICHOLSON ALGORITHM The explicit algorith begins with: u + u j ( x) (u j+ u j +u j ) while the iplicit one begins with: u + u j ( x) (u + j+ u + j +u + j ) Average these two equations:
9 -9- u + u j ( x) (u j+ u j +u j +u + j+ u + j +u + j )) (III) In ters of the atrices defined in the previous section, we can write Equation (III) as: µ (Id + A)u + µ (Id A)u +b (III) This equation is, of course, solvable using LU decoposition We are finding the values of u + j, u + j and u + j+ fro the values of u j,u j and u j+ + j- j j+ The above is the Crank-Nicholson algorith You will discover that it aintains the absolute stability of the iplicit algorith while recovering soe of the accuracy of the explicit one IV SCHRÖDINGER S EQUATION Assuing h and /, Schrödinger s equation is: i Ψ(x,t) HΨ(x,t) t H +V(x) x (IV) In the late 96 s, Goldberg, Schey and Schwartz had a good idea for a lecture deonstration in quantu echanics 3 They decided to solve this equation using thennew coputer technology to produce a ovie The initial conditions for Equation (IV) are 3 Goldberg, Schey and Schwartz, Aer Jour Phys 35, (967) 77
10 -- Ψ(x,)g(x) and the boundary conditions are: Ψ(,t) Ψ(,t) The solution to Equation (IV) can be written as: Ψ(x,t) e iht Ψ(x, ) (IV) Expand Equation (IV) and throw out higher order ters: Ψ(x,t + ) ( iht)ψ(x,t) (IV3) Thus if we know Ψ(x,t) we can produce a tie series of values for Ψ at ties t+, t+, Using finite differences to evaluate the / x ter in the Hailtonian H, the right hand side of Equation (IV3) will give a ter involving (Ψ j Ψ j +Ψ j+ ) This should be recognisable to you as an explicit algorith We can ultiply Equation (IV) by e iht and expand to get: ( + iht)ψ(x,t + ) Ψ(x,t) (IV4) Now evaluating the Hailtonian using finite differences will give a ter involving (Ψ + j Ψ + j +Ψ + j+ ) on the left hand side This is, of course, just an iplicit algorith Everything we have said about explicit and iplicit algoriths applies to these cases too However, the physics of Quantu Mechanics puts an additional constraint on Ψ: x Ψ * (x,t)ψ(x,t)dx for all ties t This is called unitarity and physically eans that the probability that the object being described by the wave function is soewhere between x and x is one It has been known since long before the invention of the coputer that Equations (IV3) and (IV4) violate unitarity This eans that neither the explicit or iplicit algorith can be used to solve Schrödinger s equation because they both violate a physical principle A long-known equation that does not violate unitarity is called Cayley s for Multiply Equation (IV) by e iht/ and expand both sides: (IV5)
11 -- ( + i Ht)Ψ(x,t +)( i Ht)Ψ(x,t) (IV6) Using finite differences to evaluate the / x ters in the Hailtonian on both sides of the equation will give us a Crank-Nicholson algorith The lesson to be learned here is that just knowing the nuerical ethods is soeties not sufficient Just as for the pendulu, where the physics pointed us to the syplectic algorith, here the physics points us to a Crank-Nicholson algorith when solving Schrödinger s equation V REFERENCES Gene H Golub and Jaes M Ortega, Scientific Coputing and Differential Equations (Acadeic Press, 99), Chapter 7 Willia H Press, Brian P Flannery, Saul A Teukolsky, and Willia T Vetterling, Nuerical Recipes: The Art of Scientific Coputing or Nuerical Recipes in C: The Art of Scientific Coputing (Cabridge Univ Press), Chapter 7 Exaples of solving Poisson s equation using relaxation ethods, and a visualisation of Schrödinger s equation solved using a Crank-Nicholson algorith ay be found at: ringland/instruction/pde/ Copyright 998 David M Harrison This is version, date (/d/y) /8/98
Donald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationQuantum algorithms (CO 781, Winter 2008) Prof. Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search
Quantu algoriths (CO 781, Winter 2008) Prof Andrew Childs, University of Waterloo LECTURE 15: Unstructured search and spatial search ow we begin to discuss applications of quantu walks to search algoriths
More informationlecture 36: Linear Multistep Mehods: Zero Stability
95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,
More informationNow multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:
Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationPhysically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationSome Perspective. Forces and Newton s Laws
Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More informationKinematics and dynamics, a computational approach
Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationAN APPLICATION OF CUBIC B-SPLINE FINITE ELEMENT METHOD FOR THE BURGERS EQUATION
Aksan, E..: An Applıcatıon of Cubıc B-Splıne Fınıte Eleent Method for... THERMAL SCIECE: Year 8, Vol., Suppl., pp. S95-S S95 A APPLICATIO OF CBIC B-SPLIE FIITE ELEMET METHOD FOR THE BRGERS EQATIO by Eine
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationFinite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields
Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for
More informationExplicit Approximate Solution for Finding the. Natural Frequency of the Motion of Pendulum. by Using the HAM
Applied Matheatical Sciences Vol. 3 9 no. 1 13-13 Explicit Approxiate Solution for Finding the Natural Frequency of the Motion of Pendulu by Using the HAM Ahad Doosthoseini * Mechanical Engineering Departent
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More informationCombining Classifiers
Cobining Classifiers Generic ethods of generating and cobining ultiple classifiers Bagging Boosting References: Duda, Hart & Stork, pg 475-480. Hastie, Tibsharini, Friedan, pg 246-256 and Chapter 10. http://www.boosting.org/
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationPHY 171. Lecture 14. (February 16, 2012)
PHY 171 Lecture 14 (February 16, 212) In the last lecture, we looked at a quantitative connection between acroscopic and icroscopic quantities by deriving an expression for pressure based on the assuptions
More informationVector Spaces in Physics 8/6/2015. Chapter 4. Practical Examples.
Vector Spaces in Physics 8/6/15 Chapter 4. Practical Exaples. In this chapter we will discuss solutions to two physics probles where we ae use of techniques discussed in this boo. In both cases there are
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationPhysics 139B Solutions to Homework Set 3 Fall 2009
Physics 139B Solutions to Hoework Set 3 Fall 009 1. Consider a particle of ass attached to a rigid assless rod of fixed length R whose other end is fixed at the origin. The rod is free to rotate about
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lessons 7 20 Dec 2017 Outline Artificial Neural networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationlecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationSupplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators
Suppleentary Inforation for Design of Bending Multi-Layer Electroactive Polyer Actuators Bavani Balakrisnan, Alek Nacev, and Elisabeth Sela University of Maryland, College Park, Maryland 074 1 Analytical
More informationOcean 420 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers
Ocean 40 Physical Processes in the Ocean Project 1: Hydrostatic Balance, Advection and Diffusion Answers 1. Hydrostatic Balance a) Set all of the levels on one of the coluns to the lowest possible density.
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationPhysics 221A: HW3 solutions
Physics 22A: HW3 solutions October 22, 202. a) It will help to start things off by doing soe gaussian integrals. Let x be a real vector of length, and let s copute dxe 2 xt Ax, where A is soe real atrix.
More informationMeasuring Temperature with a Silicon Diode
Measuring Teperature with a Silicon Diode Due to the high sensitivity, nearly linear response, and easy availability, we will use a 1N4148 diode for the teperature transducer in our easureents 10 Analysis
More informationThe path integral approach in the frame work of causal interpretation
Annales de la Fondation Louis de Broglie, Volue 28 no 1, 2003 1 The path integral approach in the frae work of causal interpretation M. Abolhasani 1,2 and M. Golshani 1,2 1 Institute for Studies in Theoretical
More informationNumerical Solution of the MRLW Equation Using Finite Difference Method. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.1401 No.3,pp.355-361 Nuerical Solution of the MRLW Equation Using Finite Difference Method Pınar Keskin, Dursun Irk
More informationThis model assumes that the probability of a gap has size i is proportional to 1/i. i.e., i log m e. j=1. E[gap size] = i P r(i) = N f t.
CS 493: Algoriths for Massive Data Sets Feb 2, 2002 Local Models, Bloo Filter Scribe: Qin Lv Local Models In global odels, every inverted file entry is copressed with the sae odel. This work wells when
More informationBALLISTIC PENDULUM. EXPERIMENT: Measuring the Projectile Speed Consider a steel ball of mass
BALLISTIC PENDULUM INTRODUCTION: In this experient you will use the principles of conservation of oentu and energy to deterine the speed of a horizontally projected ball and use this speed to predict the
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2016 Lessons 7 14 Dec 2016 Outline Artificial Neural networks Notation...2 1. Introduction...3... 3 The Artificial
More informationQuantum Chemistry Exam 2 Take-home Solutions
Cheistry 60 Fall 07 Dr Jean M Standard Nae KEY Quantu Cheistry Exa Take-hoe Solutions 5) (0 points) In this proble, the nonlinear variation ethod will be used to deterine an approxiate solution for the
More informationUCSD Spring School lecture notes: Continuous-time quantum computing
UCSD Spring School lecture notes: Continuous-tie quantu coputing David Gosset 1 Efficient siulation of quantu dynaics Quantu echanics is described atheatically using linear algebra, so at soe level is
More informationOn Lotka-Volterra Evolution Law
Advanced Studies in Biology, Vol. 3, 0, no. 4, 6 67 On Lota-Volterra Evolution Law Farruh Muhaedov Faculty of Science, International Islaic University Malaysia P.O. Box, 4, 570, Kuantan, Pahang, Malaysia
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More information10.3 Spectral Stability Criterion for Finite-Difference Cauchy Problems
Finite-Difference Schees for Partial Differential Equations 349 1.3 Spectral Stability Criterion for Finite-Difference Cauchy Probles Perhaps the ost widely used approach to the analysis of stability for
More information1 Brownian motion and the Langevin equation
Figure 1: The robust appearance of Robert Brown (1773 1858) 1 Brownian otion and the Langevin equation In 1827, while exaining pollen grains and the spores of osses suspended in water under a icroscope,
More informationNumerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda
Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji
More informationIntelligent Systems: Reasoning and Recognition. Artificial Neural Networks
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley MOSIG M1 Winter Seester 2018 Lesson 7 1 March 2018 Outline Artificial Neural Networks Notation...2 Introduction...3 Key Equations... 3 Artificial
More information4 = (0.02) 3 13, = 0.25 because = 25. Simi-
Theore. Let b and be integers greater than. If = (. a a 2 a i ) b,then for any t N, in base (b + t), the fraction has the digital representation = (. a a 2 a i ) b+t, where a i = a i + tk i with k i =
More informationOSCILLATIONS AND WAVES
OSCILLATIONS AND WAVES OSCILLATION IS AN EXAMPLE OF PERIODIC MOTION No stories this tie, we are going to get straight to the topic. We say that an event is Periodic in nature when it repeats itself in
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationGeneralized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,
Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University
More informationTutorial Exercises: Incorporating constraints
Tutorial Exercises: Incorporating constraints 1. A siple pendulu of length l ass is suspended fro a pivot of ass M that is free to slide on a frictionless wire frae in the shape of a parabola y = ax. The
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationPhysics 215 Winter The Density Matrix
Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it
More informationBipartite subgraphs and the smallest eigenvalue
Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.
More informationThe Solution of One-Phase Inverse Stefan Problem. by Homotopy Analysis Method
Applied Matheatical Sciences, Vol. 8, 214, no. 53, 2635-2644 HIKARI Ltd, www.-hikari.co http://dx.doi.org/1.12988/as.214.43152 The Solution of One-Phase Inverse Stefan Proble by Hootopy Analysis Method
More informationThe Weierstrass Approximation Theorem
36 The Weierstrass Approxiation Theore Recall that the fundaental idea underlying the construction of the real nubers is approxiation by the sipler rational nubers. Firstly, nubers are often deterined
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationDRAFT. Memo. Contents. To whom it may concern SVN: Jan Mooiman +31 (0) nl
Meo To To who it ay concern Date Reference Nuber of pages 219-1-16 SVN: 5744 22 Fro Direct line E-ail Jan Mooian +31 )88 335 8568 jan.ooian@deltares nl +31 6 4691 4571 Subject PID controller ass-spring-daper
More informationSeismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.
Seisic Analysis of Structures by K Dutta, Civil Departent, II Delhi, New Delhi. Module 5: Response Spectru Method of Analysis Exercise Probles : 5.8. or the stick odel of a building shear frae shown in
More informationLecture 20 November 7, 2013
CS 229r: Algoriths for Big Data Fall 2013 Prof. Jelani Nelson Lecture 20 Noveber 7, 2013 Scribe: Yun Willia Yu 1 Introduction Today we re going to go through the analysis of atrix copletion. First though,
More informationi ij j ( ) sin cos x y z x x x interchangeably.)
Tensor Operators Michael Fowler,2/3/12 Introduction: Cartesian Vectors and Tensors Physics is full of vectors: x, L, S and so on Classically, a (three-diensional) vector is defined by its properties under
More informationDispersion. February 12, 2014
Dispersion February 1, 014 In aterials, the dielectric constant and pereability are actually frequency dependent. This does not affect our results for single frequency odes, but when we have a superposition
More informationA note on the realignment criterion
A note on the realignent criterion Chi-Kwong Li 1, Yiu-Tung Poon and Nung-Sing Sze 3 1 Departent of Matheatics, College of Willia & Mary, Williasburg, VA 3185, USA Departent of Matheatics, Iowa State University,
More informationHee = ~ dxdy\jj+ (x) 'IJ+ (y) u (x- y) \jj (y) \jj (x), V, = ~ dx 'IJ+ (x) \jj (x) V (x), Hii = Z 2 ~ dx dy cp+ (x) cp+ (y) u (x- y) cp (y) cp (x),
SOVIET PHYSICS JETP VOLUME 14, NUMBER 4 APRIL, 1962 SHIFT OF ATOMIC ENERGY LEVELS IN A PLASMA L. E. PARGAMANIK Khar'kov State University Subitted to JETP editor February 16, 1961; resubitted June 19, 1961
More informationProblem Set 2. Chapter 1 Numerical:
Chapter 1 Nuerical: roble Set 16. The atoic radius of xenon is 18 p. Is that consistent with its b paraeter of 5.15 1 - L/ol? Hint: what is the volue of a ole of xenon atos and how does that copare to
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More informationOptimizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2
Research Paper 223 Optiizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2 Background: Success in solving the protein structure prediction
More informationLecture 8 Symmetries, conserved quantities, and the labeling of states Angular Momentum
Lecture 8 Syetries, conserved quantities, and the labeling of states Angular Moentu Today s Progra: 1. Syetries and conserved quantities labeling of states. hrenfest Theore the greatest theore of all ties
More informationStability Analysis of the Matrix-Free Linearly Implicit 2 Euler Method 3 UNCORRECTED PROOF
1 Stability Analysis of the Matrix-Free Linearly Iplicit 2 Euler Method 3 Adrian Sandu 1 andaikst-cyr 2 4 1 Coputational Science Laboratory, Departent of Coputer Science, Virginia 5 Polytechnic Institute,
More informationASSUME a source over an alphabet size m, from which a sequence of n independent samples are drawn. The classical
IEEE TRANSACTIONS ON INFORMATION THEORY Large Alphabet Source Coding using Independent Coponent Analysis Aichai Painsky, Meber, IEEE, Saharon Rosset and Meir Feder, Fellow, IEEE arxiv:67.7v [cs.it] Jul
More information1 Bounding the Margin
COS 511: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #12 Scribe: Jian Min Si March 14, 2013 1 Bounding the Margin We are continuing the proof of a bound on the generalization error of AdaBoost
More informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 4
Massachusetts Institute of Technology Quantu Mechanics I (8.04) Spring 2005 Solutions to Proble Set 4 By Kit Matan 1. X-ray production. (5 points) Calculate the short-wavelength liit for X-rays produced
More informationEstimating Parameters for a Gaussian pdf
Pattern Recognition and achine Learning Jaes L. Crowley ENSIAG 3 IS First Seester 00/0 Lesson 5 7 Noveber 00 Contents Estiating Paraeters for a Gaussian pdf Notation... The Pattern Recognition Proble...3
More informationBernoulli Wavelet Based Numerical Method for Solving Fredholm Integral Equations of the Second Kind
ISSN 746-7659, England, UK Journal of Inforation and Coputing Science Vol., No., 6, pp.-9 Bernoulli Wavelet Based Nuerical Method for Solving Fredhol Integral Equations of the Second Kind S. C. Shiralashetti*,
More informationMulti-Scale/Multi-Resolution: Wavelet Transform
Multi-Scale/Multi-Resolution: Wavelet Transfor Proble with Fourier Fourier analysis -- breaks down a signal into constituent sinusoids of different frequencies. A serious drawback in transforing to the
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.it.edu 16.323 Principles of Optial Control Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 16.323 Lecture 10 Singular
More informationI. Understand get a conceptual grasp of the problem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departent o Physics Physics 81T Fall Ter 4 Class Proble 1: Solution Proble 1 A car is driving at a constant but unknown velocity,, on a straightaway A otorcycle is
More informationComparison of Charged Particle Tracking Methods for Non-Uniform Magnetic Fields. Hann-Shin Mao and Richard E. Wirz
42nd AIAA Plasadynaics and Lasers Conferencein conjunction with the8th Internati 27-30 June 20, Honolulu, Hawaii AIAA 20-3739 Coparison of Charged Particle Tracking Methods for Non-Unifor Magnetic
More informationRESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS
BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di
More informationFirst of all, because the base kets evolve according to the "wrong sign" Schrödinger equation (see pp ),
HW7.nb HW #7. Free particle path integral a) Propagator To siplify the notation, we write t t t, x x x and work in D. Since x i, p j i i j, we can just construct the 3D solution. First of all, because
More informationDepartment of Physics, Sri Venkateswara University, Tirupati Range Operations, Satish Dhawan Space Centre SHAR, ISRO, Sriharikota
Trajectory Estiation of a Satellite Launch Vehicle Using Unscented Kalan Filter fro Noisy Radar Measureents R.Varaprasad S.V. Bhaskara Rao D.Narayana Rao V. Seshagiri Rao Range Operations, Satish Dhawan
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationExplicit Analytic Solution for an. Axisymmetric Stagnation Flow and. Heat Transfer on a Moving Plate
Int. J. Contep. Math. Sciences, Vol. 5,, no. 5, 699-7 Explicit Analytic Solution for an Axisyetric Stagnation Flow and Heat Transfer on a Moving Plate Haed Shahohaadi Mechanical Engineering Departent,
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationProblem Set 8 Solutions
Physics 57 Proble Set 8 Solutions Proble The decays in question will be given by soe Hadronic atric eleent: Γ i V f where i is the initial state, V is an interaction ter, f is the final state. The strong
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationModeling Chemical Reactions with Single Reactant Specie
Modeling Cheical Reactions with Single Reactant Specie Abhyudai Singh and João edro Hespanha Abstract A procedure for constructing approxiate stochastic odels for cheical reactions involving a single reactant
More informationReed-Muller Codes. m r inductive definition. Later, we shall explain how to construct Reed-Muller codes using the Kronecker product.
Coding Theory Massoud Malek Reed-Muller Codes An iportant class of linear block codes rich in algebraic and geoetric structure is the class of Reed-Muller codes, which includes the Extended Haing code.
More informationAn earlier article in this column considered the problem
--- CALC CORNER Estiating nternal Air Cooling Teperature Reduction in a Closed Box Utilizing Theroelectrically Enhanced Heat Rejection Previously published in February, 2013 Bob Sions BM Retired The following
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informationA Simple Regression Problem
A Siple Regression Proble R. M. Castro March 23, 2 In this brief note a siple regression proble will be introduced, illustrating clearly the bias-variance tradeoff. Let Y i f(x i ) + W i, i,..., n, where
More informationFinding Rightmost Eigenvalues of Large Sparse. Non-symmetric Parameterized Eigenvalue Problems. Abstract. Introduction
Finding Rightost Eigenvalues of Large Sparse Non-syetric Paraeterized Eigenvalue Probles Applied Matheatics and Scientific Coputation Progra Departent of Matheatics University of Maryland, College Par,
More informationarxiv: v1 [stat.ml] 31 Jan 2018
Increental kernel PCA and the Nyströ ethod arxiv:802.00043v [stat.ml] 3 Jan 208 Fredrik Hallgren Departent of Statistical Science University College London London WCE 6BT, United Kingdo fredrik.hallgren@ucl.ac.uk
More informationMidterm 1 Sample Solution
Midter 1 Saple Solution NOTE: Throughout the exa a siple graph is an undirected, unweighted graph with no ultiple edges (i.e., no exact repeats of the sae edge) and no self-loops (i.e., no edges fro a
More information