Courant and all that. Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy
|
|
- Estella Bailey
- 5 years ago
- Views:
Transcription
1 Consistency, Convergence Stability Numerical Dispersion Computational grids and numerical anisotropy c dx ε 1 he goal of this lecture is to understand how to find suitable parameters (i.e., grid spacing dx and time increment ) for a numerical simulation and knowing what can go wrong. 1
2 A simple example: Newtonian Cooling Numerical solution to first order ordinary differential equation d = f (, t) We can not simply integrate this equation. We have to solve it numerically! First we need to discretise time: t = t 0 + and for emperature = ( t )
3 A finite-difference approximation Let us try a forward difference: d t= t 1 = + + O( )... which leads to the following explicit scheme : + f (, t + 1 ) his allows us to calculate the emperature as a function of time and the forcing inhomogeneity f(,t). Note that there will be an error O() which will accumulate over time.
4 Coffee? Let s try to apply this to the Newtonian cooling problem: Air Cappucino How does the temperature of the liquid evolve as a function of time and temperature difference to the air?
5 Newtonian Cooling he rate of cooling (d/) will depend on the temperature difference ( cap - air ) and some constant (thermal conductivity). his is called Newtonian Cooling. With = cap - air being the temperature difference and τ the time scale of cooling then f(,t)=-/ τ and the differential equation describing the system is d = /τ with initial condition = i at t=0 and τ>0.
6 Analytical solution his equation has a simple analytical solution: ( t) = i exp( t / τ ) How good is our finite-difference appoximation? For what choices of will we obtain a stable solution? Our FD approximation is: + 1 = + 1 τ = = (1 (1 ) τ ) τ
7 FD Algorithm + 1 = (1 ) τ 1. Does this equation approximation converge for -> 0? 2. Does it behave like the analytical solution? With the initial condition = 0 at t=0: 2 1 = 0(1 ) τ = 1( 1 ) = 0 (1 )(1 τ τ = 0(1 ) τ leading to : ) τ
8 Understanding the numerical approximation ) 0(1 τ = Let us use =t / where t is the total time up to time step : t + = τ 1 0 his can be expanded using the binomial theorem = t t τ τ
9 ! )! (! r r r =... where we are interested in the case that -> 0 which is equivalent to -> r r r + = 1) 2)...( 1)( ( )! (! as a result r! r r
10 substituted into the series for we obtain: 2 2 t t ! τ 2! τ which leads to 2 t 1 t τ 2! τ... which is the aylor expansion for = 0 exp( t / τ )
11 So we conclude: For the Newtonian Cooling problem, the numerical solution converges to the exact solution when the time step gets smaller. How does the numerical solution behave? = 0 exp( t / τ ) he analytical solution decays monotonically! + 1 = (1 ) τ What are the conditions so that +1 <?
12 Convergence Convergence means that if we decrease time and/or space increment in our numerical algorithm we get closer to the true solution.
13 Stability + 1 = (1 ) τ +1 < requires 0 1 < 1 τ or 0 < τ he numerical solution decays only montonically for a limited range of values for! Again we seem to have a conditional stability.
14 + 1 = (1 ) τ if τ < < 2τ then ( 1 ) < 0 τ the solution oscillates but converges as 1-/τ <1 if > 2τ then / τ > 2 1-/τ<-1 and the solution oscillates and diverges... now let us see how the solution looks like...
15 Matlab solution % Matlab Program - Newtonian Cooling % initialise values nt=10; t0=1. tau=.7; =1. % initial condition =t0; % time extrapolation for i=1:nt, (i+1)=(i)-/tau*(i); end % plotting plot()
16 Example 1 Blue numerical Red - analytical
17 Example 2 Blue numerical Red - analytical
18 Example 3 Blue numerical Red - analytical
19 Example 4 Blue numerical Red - analytical
20 Stability Stability of a numerical algorithm means that the numerical solution does not tend to infinite (or zero) while the true solution is finite. In many case we obtain algoriths with conditional stability. hat means that the solution is stable for values in well-defined intervals.
21 Grid anisotropy In which direction is the solution most accurate and why?
22 Grid anisotropy with ac2d.m FD FD
23 and more a Green s function FD
24 Grids Cubed sphere Hexahedral grid Spectral element implementation etrahedral grid Discontinuous Galerkin method
25 Grid anisotropy
26 Grid anisotropy Grid anisotropy is the directional dependence of the accuracy of your numerical solution if you do not use enough points per wavelength. Grid anisotropy depends on the actual grid you are using. Cubic grids display anisotropy, hexagonal grids in 2D do not. In 3D there are no grid with isotropic properties! Numerical solutions on unstructured grids usually do not have this problem because of averaging effects, but they need more points per wavelength!
27 Real problems Alluvial Basin V S = 300m/s f max = 3Hz λ min = V S /f max = 100m 36 km 30 km Bedrock V S = 3200m/s f max = 3Hz λ min = V S /f max = m
28 Hexahedral Grid Generation
29 Courant ime step v P dx ε ε dx v P Largest velocity Smallest grid size
30 - Summary Any numerical solution has to be checked if it converges to the correct solution (as we have seen there are different options when using FD and not all do converge!) he number of grid points per wavelength is a central concept to all numerical solutions to wave like problems. his desired frequency for a simulation imposes the necessary space increment. he Courant criterion, the smallest grid increment and the largest velocity determine the (global or local) time step of the simulation Examples on the exercise sheet.
Numerical Methods in Geophysics: Implicit Methods
Numerical Methods i Geophysics: What is a implicit scheme? Explicit vs. implicit scheme for Newtoia oolig rak-nicholso Scheme (mixed explicit-implicit Explicit vs. implicit for the diffusio equatio Relaxatio
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More information3.4 Introduction to power series
3.4 Introduction to power series Definition 3.4.. A polynomial in the variable x is an expression of the form n a i x i = a 0 + a x + a 2 x 2 + + a n x n + a n x n i=0 or a n x n + a n x n + + a 2 x 2
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationChecking the Radioactive Decay Euler Algorithm
Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More information(Refer Slide Time: 03: 09)
Computational Electromagnetics and Applications Professor Krish Sankaran Indian Institute of Technology Bombay Lecture No 26 Finite Volume Time Domain Method-I Welcome back in the precious lectures we
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationMATH3283W LECTURE NOTES: WEEK 6 = 5 13, = 2 5, 1 13
MATH383W LECTURE NOTES: WEEK 6 //00 Recursive sequences (cont.) Examples: () a =, a n+ = 3 a n. The first few terms are,,, 5 = 5, 3 5 = 5 3, Since 5
More informationFINITE DIFFERENCES. Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations.
FINITE DIFFERENCES Lecture 1: (a) Operators (b) Forward Differences and their calculations. (c) Backward Differences and their calculations. 1. Introduction When a function is known explicitly, it is easy
More informationDOING PHYSICS WITH MATLAB WAVE MOTION
DOING PHYSICS WITH MATLAB WAVE MOTION THE [1D] SCALAR WAVE EQUATION THE FINITE DIFFERENCE TIME DOMAIN METHOD Ian Cooper School of Physics, University of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY
More information16.4. Power Series. Introduction. Prerequisites. Learning Outcomes
Power Series 6.4 Introduction In this Section we consider power series. These are examples of infinite series where each term contains a variable, x, raised to a positive integer power. We use the ratio
More informationECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM. Prof. Peter Bermel February 6, 2017
ECE 695 Numerical Simulations Lecture 12: Thermomechanical FEM Prof. Peter Bermel February 6, 2017 Outline Static Equilibrium Dynamic Equilibrium Thermal transport mechanisms Thermal transport modeling
More informationseismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model
seismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model properties inversion of seismic observables for 3D Earth structure
More informationReplacing the a in the definition of the derivative of the function f at a with a variable x, gives the derivative function f (x).
Definition of The Derivative Function Definition (The Derivative Function) Replacing the a in the definition of the derivative of the function f at a with a variable x, gives the derivative function f
More informationStep Response of First-Order Systems
INTRODUCTION This tutorial discusses the response of a first-order system to a unit step function input. In particular, it addresses the time constant and how that affects the speed of the system s response.
More informationModule 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series
Module 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series Lecture 27 : Series of functions [Section 271] Objectives In this section you will learn
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationTaylor Series and Numerical Approximations
Taylor Series and Numerical Approximations Hilary Weller h.weller@reading.ac.uk August 7, 05 An introduction to the concept of a Taylor series and how these are used in numerical analysis to find numerical
More informationThe Derivative Function. Differentiation
The Derivative Function If we replace a in the in the definition of the derivative the function f at the point x = a with a variable x, we get the derivative function f (x). Using Formula 2 gives f (x)
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationGeneralized Fourier Analyses of Semi-Discretizations of the Advection-Diffusion Equation
Generalized Fourier Analyses of Semi-Discretizations of the Advection-Diffusion Equation Sponsored by: ASCI Algorithms - Advanced Spatial Discretization Project Mark A. Christon, homas E. Voth Computational
More informationThe Finite Difference Method
The Finite Difference Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 32 Outline 1 Introduction Motivation
More informationMath 1b Sequences and series summary
Math b Sequences and series summary December 22, 2005 Sequences (Stewart p. 557) Notations for a sequence: or a, a 2, a 3,..., a n,... {a n }. The numbers a n are called the terms of the sequence.. Limit
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationMATH 117 LECTURE NOTES
MATH 117 LECTURE NOTES XIN ZHOU Abstract. This is the set of lecture notes for Math 117 during Fall quarter of 2017 at UC Santa Barbara. The lectures follow closely the textbook [1]. Contents 1. The set
More information* = 2 = Probability distribution function. probability of finding a particle near a given point x,y,z at a time t
Quantum Mechanics Wave functions and the Schrodinger equation Particles behave like waves, so they can be described with a wave function (x,y,z,t) A stationary state has a definite energy, and can be written
More information9.3. Total number of phonon modes, total energy and heat capacity
Phys50.nb 6 E = n = n = exp - (9.9) 9... History of the Planck distribution or the Bose-Einstein distribution. his distribution was firstly discovered by Planck in the study of black-body radiation. here,
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationAn explicit time-domain finite-element method for room acoustics simulation
An explicit time-domain finite-element method for room acoustics simulation Takeshi OKUZONO 1 ; Toru OTSURU 2 ; Kimihiro SAKAGAMI 3 1 Kobe University, JAPAN 2 Oita University, JAPAN 3 Kobe University,
More informationFLUID MECHANICS 3 - LECTURE 4 ONE-DIMENSIONAL UNSTEADY GAS
FLUID MECHANICS 3 - LECTURE 4 ONE-DIMENSIONAL UNSTEADY GAS Consider an unsteady 1-dimensional ideal gas flow. We assume that this flow is spatially continuous and thermally isolated, hence, it is isentropic.
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationSeries. 1 Convergence and Divergence of Series. S. F. Ellermeyer. October 23, 2003
Series S. F. Ellermeyer October 23, 2003 Convergence and Divergence of Series An infinite series (also simply called a series) is a sum of infinitely many terms a k = a + a 2 + a 3 + () The sequence a
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
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 informationMITOCW MIT18_02SCF10Rec_61_300k
MITOCW MIT18_02SCF10Rec_61_300k JOEL LEWIS: Hi. Welcome back to recitation. In lecture, you've been learning about the divergence theorem, also known as Gauss's theorem, and flux, and all that good stuff.
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
More information1. Transition dipole moment
1. Transition dipole moment You have measured absorption spectra of aqueous (n=1.33) solutions of two different chromophores (A and B). The concentrations of the solutions were the same. The absorption
More informationGregory's quadrature method
Gregory's quadrature method Gregory's method is among the very first quadrature formulas ever described in the literature, dating back to James Gregory (638-675). It seems to have been highly regarded
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week 2 November 6 November Deadline to hand in the homeworks: your exercise class on week 9 November 13 November Exercises (1) Let X be the following space of piecewise
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 informationCosmology. Lecture Topics. Colliding Galaxy Simulations The expanding universe Cosmology. The Age of the Universe The Big Bang
Cosmology Lecture 27: Cosmology 27 1 APOD: GOODS project Lecture Topics Colliding Galaxy Simulations The expanding universe Cosmology The cosmological principle The Age of the Universe The Big Bang 2 27-1
More informationLECTURE 4-1 INTRODUCTION TO ORDINARY DIFFERENTIAL EQUATIONS
130 LECTURE 4-1 INTRODUCTION TO ORDINARY DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS: A differential equation (DE) is an equation involving an unknown function and one or more of its derivatives. A differential
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationG12MAN Mathematical Analysis Fully-justified answers to Revision Quiz questions
G12MAN Mathematical Analysis Fully-justified answers to Revision Quiz questions Remember that, unless otherwise specified or not applicable for some reason), you MUST justify your answers to questions
More information1. The growth of a cancerous tumor can be modeled by the Gompertz Law: dn. = an ln ( )
1. The growth of a cancerous tumor can be modeled by the Gompertz Law: ( ) dn b = an ln, (1) dt N where N measures the size of the tumor. (a) Interpret the parameters a and b (both non-negative) biologically.
More information35. RESISTIVE INSTABILITIES: CLOSING REMARKS
35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More information31.1.1Partial derivatives
Module 11 : Partial derivatives, Chain rules, Implicit differentiation, Gradient, Directional derivatives Lecture 31 : Partial derivatives [Section 31.1] Objectives In this section you will learn the following
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationA Brief Introduction to Numerical Methods for Differential Equations
A Brief Introduction to Numerical Methods for Differential Equations January 10, 2011 This tutorial introduces some basic numerical computation techniques that are useful for the simulation and analysis
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More informationChapter 18. Remarks on partial differential equations
Chapter 8. Remarks on partial differential equations If we try to analyze heat flow or vibration in a continuous system such as a building or an airplane, we arrive at a kind of infinite system of ordinary
More informationMATH20411 PDEs and Vector Calculus B
MATH2411 PDEs and Vector Calculus B Dr Stefan Güttel Acknowledgement The lecture notes and other course materials are based on notes provided by Dr Catherine Powell. SECTION 1: Introctory Material MATH2411
More informationEuler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.
Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be
More informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationEntropy Changes & Processes
Entropy Changes & Processes Chapter 4 of Atkins: he Second Law: he Concepts Section 4.3 Entropy of Phase ransition at the ransition emperature Expansion of the Perfect Gas Variation of Entropy with emperature
More informationMethods of solution. 2.1 Introduction
.1 Introduction The numerical solution of non-linear problems is quite different from linear problems. The solution can not be obtained with a single step but by a superposition of linear solutions. Because
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationREVIEW: The Matching Method Algorithm
Lecture 26: Numerov Algorithm for Solving the Time-Independent Schrödinger Equation 1 REVIEW: The Matching Method Algorithm Need for a more general method The shooting method for solving the time-independent
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 informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
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 informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationEntropy Changes & Processes
Entropy Changes & Processes Chapter 4 of Atkins: he Second Law: he Concepts Section 4.3, 7th edition; 3.3, 8th and 9th editions Entropy of Phase ransition at the ransition emperature Expansion of the Perfect
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationNew Scaling Factors of 2-D Isotropic-Dispersion Finite Difference Time Domain (ID-FDTD) Algorithm for Lossy Media
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 56, NO., FEBRUARY 8 63 is shown by squares. Each point is computed by averaging reflection coefficient values calculated for each component of the regular
More informationModule 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series
Module 9 : Infinite Series, Tests of Convergence, Absolute and Conditional Convergence, Taylor and Maclaurin Series Lecture 26 : Absolute convergence [Section 261] Objectives In this section you will learn
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationSection 3.2 : Sequences
Section 3.2 : Sequences Note: Chapter 11 of Stewart s Calculus is a good reference for this chapter of our lecture notes. Definition 52 A sequence is an infinite ordered list a 1, a 2, a 3,... The items
More informationChem 452 Mega Practice Exam 1
Last Name: First Name: PSU ID #: Chem 45 Mega Practice Exam 1 Cover Sheet Closed Book, Notes, and NO Calculator The exam will consist of approximately 5 similar questions worth 4 points each. This mega-exam
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationINTRODUCTION TO MICROWAVE REMOTE SENSING - II. Dr. A. Bhattacharya
1 INTRODUCTION TO MICROWAVE REMOTE SENSING - II Dr. A. Bhattacharya The Radiation Framework The information about features on the Earth s surface using RS depends on measuring energy emanating from the
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More information8.1 Sequences. Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = 1.
8. Sequences Example: A sequence is a function f(n) whose domain is a subset of the integers. Notation: *Note: n = 0 vs. n = Examples: 6. Find a formula for the general term a n of the sequence, assuming
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More informationSupplementary material for: Self-similar slip instability on interfaces with rate- and state-dependent friction
rspa.royalsocietypublishing.org Supplementary material for: Self-similar slip instability on interfaces with rate- and state-dependent friction Robert C. Viesca 1 1 Tufts University, Medford, MA 0155 c
More information1D Heat equation and a finite-difference solver
1D Heat equation and a finite-difference solver Guillaume Riflet MARETEC IST 1 The advection-diffusion equation The original concept applied to a property within a control volume from which is derived
More informationA Few Examples. A Few Examples
Section 3.2 : Sequences Note: Chapter of Stewart s Calculus is a good reference for this chapter of our lecture notes. Definition 52 A sequence is an infinite ordered list A Few Examples (( ) n + ) n=
More informationLecture 2 Blackbody radiation
Lecture 2 Blackbody radiation Absorption and emission of radiation What is the blackbody spectrum? Properties of the blackbody spectrum Classical approach to the problem Plancks suggestion energy quantisation
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 6: January 30, 2018 Inverse z-transform Lecture Outline! z-transform " Tie up loose ends " Regions of convergence properties! Inverse z-transform " Inspection " Partial
More informationComputational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. SumanChakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #11 Fundamentals of Discretization: Finite Difference
More informationGradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice
1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationRelaxation. Ravinder Reddy
Relaxation Ravinder Reddy Relaxation What is nuclear spin relaxation? What causes it? Effect on spectral line width Field dependence Mechanisms Thermal equilibrium ~10-6 spins leads to NMR signal! T1 Spin-lattice
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationNumerical Methods for Partial Differential Equations CAAM 452. Spring 2005
Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks
More informationAbstract. Introduction
Stability Criterion for Explicit Schemes (Finite-Difference Method) on the Solution of the Advection-Diffusion Equation L.F. León, P.M. Austria Mexican Institute of Water Technology, Paseo Cuauhnáhuac
More information