Numerical Astrophysics: hydrodynamics
|
|
- Gavin Denis Singleton
- 5 years ago
- Views:
Transcription
1 Numerical Astrophysics: hydrodyamics
2 Part 1: Numerical solutios to the Euler Equatios
3 Outlie Numerical eperimets are a valuable tool to study astrophysical objects (where we ca rarely do direct eperimets). Numerical simulatios i astrophysics cover a wide rage of situatios, with may techiques ad various degrees of compleity. The most widely used i astrophysics ca be roughly divided ito: N-Body dyamics; Fluid dyamics (hydrodyamics, magetohydrodyamics, SPH...). This first series of lectures will provide a hads o approach to umerical hydrodyamics. The goal is to write a series of codes to solve the equatios of gas dyamics i oe dimesio, or i multi-dimesios. Material (templates ad sample codes) will be made available i fortra 90, as well as basic plottig scripts i guplot or pytho (matplotlib).
4 Hydrodyamics I order to describe a system of particles as a fluid several coditios must be met The mea free path has to be smaller tha the scale-legth of the fluctuatios of the macroscopic variables (such as desity, pressure...). mfp L The time betwee collisios must be small compared with the time-scale of chages i the flow. t coll t flow The mea distace betwee particles has to be smaller tha the scale-legth of the macroscopic variables. l = 1/3 L
5 The Euler equatios We cosider a series of fied volume elemets (Euleria mesh). The equatios that describe the fluid dyamics i a fied volume elemet (disregardig viscous effects) are: Mass coservatio + r ( u) Mometum (Newto s 2d u) + r ( uu)+rp = f Eergy + r [u (E + P )] = L + f et u Alog with a equatio of state of 55 PDE E = 1 2 u2 + P 1 they form a system
6 Cartesia coordiates, writte i matri @z = S U = u v w E 3 7 5, F = u u 2 + P uv uw u(e + P ) 3 7 5, G = v vu v 2 + P vw v(e + P ) 3 7 5, H = w wu wv w 2 + P w(e + P ) 3 7 5, S = G 0 f f y f z L + f u ad u=(u, v, w) coserved variables,y, ad z flues source terms
7 Fiite differeces Discretize the coordiate i N poits i = 0 + i, i =1, 2,...N cosider a fuctio f()! f( i )=f i the derivative of f() ca be approimated by fiite differeces f 0 fwd() = f 0 back() = f( + ) f() f() f( ) = f i+1 f i = f i f i 1 f f 0 cet( i ) f i+1 f i 1 2 f i f i+1 f() f i 1 f 0 cet() = f( + ) f( ) 2 = f i+1 f i 1 2 Remember the formal defiitio of the derivative f 0 () = lim h!0 f( + h) f() h i 1 i i+1
8 Applicatio of fiite differeces: FTCS Discretize space as: i = 0 + i ; i =1, 2... Ad time as: t = t 0 + t;,=1, 2,... With this otatio U(, t) =U i The we apply forward differece i time (FT), ad cetered differecig i U+1 i t F i+1 F i @ = S ) U+1 i = U i t 2 F i+1 F i 1 + ts However, it turs out to be ucoditioally ustable for hyperbolic equatios (i.e. useless for our case)
9 The La-Friedrichs method La fied the problem of stability of the FTCS by replacig U i! 1 2 U i+1 + U i 1 Thus the La Method reads: U +1 i = 1 2 U i+1 + U i 1 t 2 F i+1 F i 1 + ts However, we ca ot take arbitrarily large values of restricted by a stability criterium. t, they are
10 The advectio equatio Cosider the liear + =0 with iitial coditios has ) = S U = u, F = au, ad S =0 u(, t = 0) = f() u(, t) =f( at) f a<0 f(, t = 0) a>0
11 Vo Neuma aalysis Cosider the advectio equatio If oe assumes that oe has: Thus the error follows Apply La Method to oe Fourier mode ) +1 i = 1 2 ˆ +1 i+1 + i 1 + a @t + =0 u =ū + i+1 i 1 true solutio ˆ = 1 2 ejk +e jk + a t 2 ejk +e jk error(e.g. trucatio) (, t) = i = ˆ (t)e jk = ˆ e jk ; j = p 1 = cos(k) j a t si(k) The error grows if ˆ +1 ˆ > 1 The coditio for stability becomes a t = a t = Co apple 1 This is the famous Courat-Friedrichs-Lewy stability coditio
12 Geeral solutio strategy (1D): 1. Discretize domai i N (+2) poits i =0, 1, 2...N +(N + 1) 2. Iitial coditios i the etire domai: 3. Determie the timestep allowed by the CFL coditio t = C,withC<1 a 4. Calculate the Flues F i = F(U i ) 5. Advace the solutio, e.g.: U +1 i = 1 2 U i+1 + U i 1 t 2 6. Apply boudary coditios U +1 0, U +1 N+1 F i+1 F i At this poit = +1, (t = t + t) ad we iterate ts
13 Eample of the advectio equatio solved with the La method Advectio of a regular pulse Periodic Boudary coditios: U +1 0 = U +1 N, U +1 N +1 = U+1 1 i =0 i =1 i =2... i = Ni =(N+ 1) v Advectio eq. w/la method t= 0 t t= 2 t t= 4 t t= 6 t Eample code i fortra: advectio_la.f90 Guplot plottig script plot_advectio.gpl pytho (matpotlib) plottig script plot_advectio.py Dowload from: Click o the right o the tab: Hydrodyamics ICTP-SAIFR School
14 Pgas Task 1: Write a program (suggestio modify the advectio_la.f90) to solve the Euler equatios. Iitial coditios: classic Sod tube test Euler equatios w/la s Method, t= L =1.0 P L = u R =0.1 u L = u R =0 0.0 P R = You will eed to chage u(0:+1)to u(3,0:+1) u(1,i):desity at _i u(2,i): mometum at _i u(3,i): total eergy desity (same with F, UP) Set the iitial coditios (see fig) Chage output to write desity, velocity, ad Pressure The ew CFL coditio is Co t =mi u + cs Boudary coditios: ouflow Itegrate oly to all domai U +1 0 = U +1 1, U +1 N +1 = U+1 N t ma =0.2
15 Here s how the Sod tube test solutio looks like after t=0.2 Euler equatios w/la s Method, t=0.2 Rarefactio wave u Pgas cotact discotiuity Shock
Lecture 2: Finite Difference Methods in Heat Transfer
Lecture 2: Fiite Differece Methods i Heat Trasfer V.Vuorie Aalto Uiversity School of Egieerig Heat ad Mass Trasfer Course, Autum 2016 November 1 st 2017, Otaiemi ville.vuorie@aalto.fi Overview Part 1 (
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationComputational Fluid Dynamics. Lecture 3
Computatioal Fluid Dyamics Lecture 3 Discretizatio Cotiued. A fourth order approximatio to f x ca be foud usig Taylor Series. ( + ) + ( + ) + + ( ) + ( ) = a f x x b f x x c f x d f x x e f x x f x 0 0
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More informationHydro code II. Simone Recchi
Hydro code II Simoe Recchi 16.06.2011 Literature E.F. Toro: Riema Solvers ad Numerical Methods for Fluid Dyamics, Spriger R.J. LeVeque: Fiite Volume Methods for Hyperbolic Problems, Cambridge Uiversity
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationAnalysis of a Numerical Scheme An Example
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 3 Jauary 5, 7 Aalysis of a Numerical Scheme A Example Grétar Tryggvaso Numerical Aalysis Example Use the leap-frog method (cetered
More informationSolutions to Final Exam Review Problems
. Let f(x) 4+x. Solutios to Fial Exam Review Problems Math 5C, Witer 2007 (a) Fid the Maclauri series for f(x), ad compute its radius of covergece. Solutio. f(x) 4( ( x/4)) ( x/4) ( ) 4 4 + x. Sice the
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationFree Surface Hydrodynamics
Water Sciece ad Egieerig Free Surface Hydrodyamics y A part of Module : Hydraulics ad Hydrology Water Sciece ad Egieerig Dr. Shreedhar Maskey Seior Lecturer UNESCO-IHE Istitute for Water Educatio S. Maskey
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationChapter 10 Partial Differential Equations and Fourier Series
Math-33 Chapter Partial Differetial Equatios November 6, 7 Chapter Partial Differetial Equatios ad Fourier Series Math-33 Chapter Partial Differetial Equatios November 6, 7. Boudary Value Problems for
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationCHAPTER 8 SYSTEMS OF PARTICLES
CHAPTER 8 SYSTES OF PARTICLES CHAPTER 8 COLLISIONS 45 8. CENTER OF ASS The ceter of mass of a system of particles or a rigid body is the poit at which all of the mass are cosidered to be cocetrated there
More informationAP Calculus BC Review Applications of Derivatives (Chapter 4) and f,
AP alculus B Review Applicatios of Derivatives (hapter ) Thigs to Kow ad Be Able to Do Defiitios of the followig i terms of derivatives, ad how to fid them: critical poit, global miima/maima, local (relative)
More informationFinite Dierence Schemes
MATH-459 Numerical Methods for Coservatio Laws by Prof. Ja S. Hesthave Solutio set 2: Fiite Dierece Schemes Exercise 2. Cosistecy A method is cosistet if its local trucatio error T k satises T k (x, t)
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differetial Equatios Eric de Sturler Uiversity of Illiois at Urbaa-Champaig Cosider the liear first order hyperbolic equatio Øu Øu Øt + a(x, t) Øx = 0 I spite of its simple
More informationPartial Differential Equations
Partial Differetial Equatios Part 2 Massimo Ricotti ricotti@astro.umd.edu Uiversity of Marylad Partial Differetial Equatios p.1/15 Upwid differecig I additio to amplitude errors (istability or dampig),
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More information1. Hydrogen Atom: 3p State
7633A QUANTUM MECHANICS I - solutio set - autum. Hydroge Atom: 3p State Let us assume that a hydroge atom is i a 3p state. Show that the radial part of its wave fuctio is r u 3(r) = 4 8 6 e r 3 r(6 r).
More informationMatsubara-Green s Functions
Matsubara-Gree s Fuctios Time Orderig : Cosider the followig operator If H = H the we ca trivially factorise this as, E(s = e s(h+ E(s = e sh e s I geeral this is ot true. However for practical applicatio
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationNotes on iteration and Newton s method. Iteration
Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f
More informationSalmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations
3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)
More informationThe Advection-Diffusion equation!
ttp://www.d.edu/~gtryggva/cf-course/! Te Advectio-iffusio equatio! Grétar Tryggvaso! Sprig 3! Navier-Stokes equatios! Summary! u t + u u x + v u y = P ρ x + µ u + u ρ y Hyperbolic part! u x + v y = Elliptic
More informationExample 2. Find the upper bound for the remainder for the approximation from Example 1.
Lesso 8- Error Approimatios 0 Alteratig Series Remaider: For a coverget alteratig series whe approimatig the sum of a series by usig oly the first terms, the error will be less tha or equal to the absolute
More informationElectromagnetic wave propagation in Particle-In-Cell codes
Electromagetic wave propagatio i Particle-I-Cell codes Remi Lehe Lawrece Berkeley Natioal Laboratory (LBNL) US Particle Accelerator School (USPAS) Summer Sessio Self-Cosistet Simulatios of Beam ad Plasma
More informationFor example suppose we divide the interval [0,2] into 5 equal subintervals of length
Math 1206 Calculus Sec 1: Estimatig with Fiite Sums Abbreviatios: wrt with respect to! for all! there exists! therefore Def defiitio Th m Theorem sol solutio! perpedicular iff or! if ad oly if pt poit
More informationAn efficient time integration method for extra-large eddy simulations
A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler Departmet of Mathematics Master s Thesis A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationLesson 03 Heat Equation with Different BCs
PDE & Complex Variables P3- esso 3 Heat Equatio with Differet BCs ( ) Physical meaig (SJF ) et u(x, represet the temperature of a thi rod govered by the (coductio) heat equatio: u t =α u xx (3.) where
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationFluid Physics 8.292J/12.330J % (1)
Fluid Physics 89J/133J Problem Set 5 Solutios 1 Cosider the flow of a Euler fluid i the x directio give by for y > d U = U y 1 d for y d U + y 1 d for y < This flow does ot vary i x or i z Determie the
More informationFinite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 39-764 Ide Copericus Value (3): 64 Impact Factor (3): 4438 Fiite Differece Approimatio for Trasport Equatio with Shifts Arisig i Neuroal Variability
More informationNumerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance
Iteratioal oural of Sciece ad Research (ISR) ISSN (Olie): 39-74 Ide Copericus Value (3): 4 Impact Factor (3): 4438 Numerical Solutio of the First-Order Hyperbolic Partial Differetial Equatio with Poit-Wise
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Hyperbolic PDEs: Coupled system/noliear coservatio laws/a oliear Lax-Wedroff scheme (Lecture 18, Week 6 Markus Schmuck Departmet of Mathematics ad Maxwell Istitute for Mathematical
More informationMaximum and Minimum Values
Sec 4.1 Maimum ad Miimum Values A. Absolute Maimum or Miimum / Etreme Values A fuctio Similarly, f has a Absolute Maimum at c if c f f has a Absolute Miimum at c if c f f for every poit i the domai. f
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More informationLecture 25 (Dec. 6, 2017)
Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed
More informationLinear Elliptic PDE s Elliptic partial differential equations frequently arise out of conservation statements of the form
Liear Elliptic PDE s Elliptic partial differetial equatios frequetly arise out of coservatio statemets of the form B F d B Sdx B cotaied i bouded ope set U R. Here F, S deote respectively, the flux desity
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
More informationCOMM 602: Digital Signal Processing
COMM 60: Digital Sigal Processig Lecture 4 -Properties of LTIS Usig Z-Trasform -Iverse Z-Trasform Properties of LTIS Usig Z-Trasform Properties of LTIS Usig Z-Trasform -ve +ve Properties of LTIS Usig Z-Trasform
More informationP.3 Polynomials and Special products
Precalc Fall 2016 Sectios P.3, 1.2, 1.3, P.4, 1.4, P.2 (radicals/ratioal expoets), 1.5, 1.6, 1.7, 1.8, 1.1, 2.1, 2.2 I Polyomial defiitio (p. 28) a x + a x +... + a x + a x 1 1 0 1 1 0 a x + a x +... +
More informationAE/ME 339 Computational Fluid Dynamics (CFD)
AE/ME 339 Computatioal Fluid Dyamics (CFD 0//004 Topic0_PresCorr_ Computatioal Fluid Dyamics (AE/ME 339 Pressure Correctio Method The pressure correctio formula (6.8.4 Calculatio of p. Coservatio form
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationFinite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts
Iteratioal Joural of Scietific Egieerig ad Research (IJSER) wwwiseri ISSN (Olie): 347-3878, Impact Factor (4): 35 Fiite Differece Approimatio for First- Order Hyperbolic Partial Differetial Equatio Arisig
More information1 Adiabatic and diabatic representations
1 Adiabatic ad diabatic represetatios 1.1 Bor-Oppeheimer approximatio The time-idepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationLecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables
Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio
More informationNonequilibrium Excess Carriers in Semiconductors
Lecture 8 Semicoductor Physics VI Noequilibrium Excess Carriers i Semicoductors Noequilibrium coditios. Excess electros i the coductio bad ad excess holes i the valece bad Ambiolar trasort : Excess electros
More informationMore Elementary Aspects of Numerical Solutions of PDEs!
ttp://www.d.edu/~gtryggva/cfd-course/ Outlie More Elemetary Aspects o Numerical Solutios o PDEs I tis lecture we cotiue to examie te elemetary aspects o umerical solutios o partial dieretial equatios.
More informationTHEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS. Marcel Migdalovici 1 and Daniela Baran 2
ICSV4 Cairs Australia 9- July, 007 THEORETICAL RESEARCH REGARDING ANY STABILITY THEOREMS WITH APPLICATIONS Marcel Migdalovici ad Daiela Bara Istitute of Solid Mechaics, INCAS Elie Carafoli, 5 C-ti Mille
More information2. Fourier Series, Fourier Integrals and Fourier Transforms
Mathematics IV -. Fourier Series, Fourier Itegrals ad Fourier Trasforms The Fourier series are used for the aalysis of the periodic pheomea, which ofte appear i physics ad egieerig. The Fourier itegrals
More informationMathematical Series (You Should Know)
Mathematical Series You Should Kow Mathematical series represetatios are very useful tools for describig images or for solvig/approimatig the solutios to imagig problems. The may be used to epad a fuctio
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationmx bx kx F t. dt IR I LI V t, Q LQ RQ V t,
Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio
More information1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations
. Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a
More informationAccuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics
http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture Jauary 3, 7 Grétar Tryggvaso It is clear that although the umerical solutio is qualitatively similar to the aalytical solutio,
More informationHonors Calculus Homework 13 Solutions, due 12/8/5
Hoors Calculus Homework Solutios, due /8/5 Questio Let a regio R i the plae be bouded by the curves y = 5 ad = 5y y. Sketch the regio R. The two curves meet where both equatios hold at oce, so where: y
More informationFrequency Domain Filtering
Frequecy Domai Filterig Raga Rodrigo October 19, 2010 Outlie Cotets 1 Itroductio 1 2 Fourier Represetatio of Fiite-Duratio Sequeces: The Discrete Fourier Trasform 1 3 The 2-D Discrete Fourier Trasform
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationPrecalculus MATH Sections 3.1, 3.2, 3.3. Exponential, Logistic and Logarithmic Functions
Precalculus MATH 2412 Sectios 3.1, 3.2, 3.3 Epoetial, Logistic ad Logarithmic Fuctios Epoetial fuctios are used i umerous applicatios coverig may fields of study. They are probably the most importat group
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More information1. pn junction under bias 2. I-Vcharacteristics
Lecture 10 The p Juctio (II) 1 Cotets 1. p juctio uder bias 2. I-Vcharacteristics 2 Key questios Why does the p juctio diode exhibit curret rectificatio? Why does the juctio curret i forward bias icrease
More informationMATH 31B: MIDTERM 2 REVIEW
MATH 3B: MIDTERM REVIEW JOE HUGHES. Evaluate x (x ) (x 3).. Partial Fractios Solutio: The umerator has degree less tha the deomiator, so we ca use partial fractios. Write x (x ) (x 3) = A x + A (x ) +
More informationFFTs in Graphics and Vision. The Fast Fourier Transform
FFTs i Graphics ad Visio The Fast Fourier Trasform 1 Outlie The FFT Algorithm Applicatios i 1D Multi-Dimesioal FFTs More Applicatios Real FFTs 2 Computatioal Complexity To compute the movig dot-product
More informationPC5215 Numerical Recipes with Applications - Review Problems
PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationIntroduction to Astrophysics Tutorial 2: Polytropic Models
Itroductio to Astrophysics Tutorial : Polytropic Models Iair Arcavi 1 Summary of the Equatios of Stellar Structure We have arrived at a set of dieretial equatios which ca be used to describe the structure
More informationTHE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION
Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 S93 THE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION by
More informationCHAPTER 1 SEQUENCES AND INFINITE SERIES
CHAPTER SEQUENCES AND INFINITE SERIES SEQUENCES AND INFINITE SERIES (0 meetigs) Sequeces ad limit of a sequece Mootoic ad bouded sequece Ifiite series of costat terms Ifiite series of positive terms Alteratig
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationA large class of initial value (time-evolution) PDEs in one space dimension can be cast into the form of a flux-conservative equation, u (19.1.
834 Chapter 19. Partial Differetial Equatios egieerig; these methods allow cosiderable freedom i puttig computatioal elemets where you wat them, importat whe dealig with highly irregular geometries. Spectral
More informationJournal of Computational Physics 149, (1999) Article ID jcph , available online at
Joural of Computatioal Physics 149, 418 422 (1999) Article ID jcph.1998.6131, available olie at http://www.idealibrary.com o NOTE Defiig Wave Amplitude i Characteristic Boudary Coditios Key Words: Euler
More informationChapter 4 : Laplace Transform
4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic
More informationNumerical 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 informationFundamental Concepts: Surfaces and Curves
UNDAMENTAL CONCEPTS: SURACES AND CURVES CHAPTER udametal Cocepts: Surfaces ad Curves. INTRODUCTION This chapter describes two geometrical objects, vi., surfaces ad curves because the pla a ver importat
More informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More informationBeurling Integers: Part 2
Beurlig Itegers: Part 2 Isomorphisms Devi Platt July 11, 2015 1 Prime Factorizatio Sequeces I the last article we itroduced the Beurlig geeralized itegers, which ca be represeted as a sequece of real umbers
More informationLecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box
561 Fall 013 Lecture #5 page 1 Last time: Lecture #5: Begi Quatum Mechaics: Free Particle ad Particle i a 1D Box u 1 u 1-D Wave equatio = x v t * u(x,t): displacemets as fuctio of x,t * d -order: solutio
More informationENGI 9420 Engineering Analysis Assignment 3 Solutions
ENGI 9 Egieerig Aalysis Assigmet Solutios Fall [Series solutio of ODEs, matri algebra; umerical methods; Chapters, ad ]. Fid a power series solutio about =, as far as the term i 7, to the ordiary differetial
More informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationHOMEWORK #10 SOLUTIONS
Math 33 - Aalysis I Sprig 29 HOMEWORK # SOLUTIONS () Prove that the fuctio f(x) = x 3 is (Riema) itegrable o [, ] ad show that x 3 dx = 4. (Without usig formulae for itegratio that you leart i previous
More informationLecture 6: Integration and the Mean Value Theorem. slope =
Math 8 Istructor: Padraic Bartlett Lecture 6: Itegratio ad the Mea Value Theorem Week 6 Caltech 202 The Mea Value Theorem The Mea Value Theorem abbreviated MVT is the followig result: Theorem. Suppose
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationALGEBRAIC GEOMETRY COURSE NOTES, LECTURE 5: SINGULARITIES.
ALGEBRAIC GEOMETRY COURSE NOTES, LECTURE 5: SINGULARITIES. ANDREW SALCH 1. The Jacobia criterio for osigularity. You have probably oticed by ow that some poits o varieties are smooth i a sese somethig
More informationOffice: JILA A709; Phone ;
Office: JILA A709; Phoe 303-49-7841; email: weberjm@jila.colorado.edu Problem Set 5 To be retured before the ed of class o Wedesday, September 3, 015 (give to me i perso or slide uder office door). 1.
More informationNumerical Methods for Partial Differential Equations
Numerical Methods for Partial Differetial Equatios CAAM 45 Sprig 005 Lecture 4 -step time-steppig methods: stability, accuracy Ruge-Kutta Methods, Istructor: Tim Warburto Today Recall AB stability regios
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationQuantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.
Quatum Mechaics I 1 April, 14 Assigmet 5: Solutio 1 For a particle icidet o a potetial step with E < V, show that the magitudes of the amplitudes of the icidet ad reflected waves fuctios are the same Fid
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationShock-Turbulence Interaction
Shock-Turbulece Iteractio A.Sakurai ad M.Tsukamoto Tokyo Deki Uiversity, Nishikicho -, Kada, Chiyoda-ku, Tokyo, Japa Abstract. For the geeral purpose of ivestigatig pheomeo of shock-turbulece iteractio,
More information17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)
7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.
More information