Hydro code II. Simone Recchi
|
|
- Lily Davidson
- 5 years ago
- Views:
Transcription
1 Hydro code II Simoe Recchi
2 Literature E.F. Toro: Riema Solvers ad Numerical Methods for Fluid Dyamics, Spriger R.J. LeVeque: Fiite Volume Methods for Hyperbolic Problems, Cambridge Uiversity Press flash.uchicago.edu/site/flashcode irfu.cea.fr/projects/site_ramses/ramses.html
3 Visual outlie set up (grid based) Hydro code evolutio grid geeratio iitial values iput output storage solutio advacemet boudary values time stepper hydro solver source terms solutio update
4 set up (grid based) Hydro code evolutio grid geeratio iitial values iput output storage solutio advacemet boudary values time stepper hydro solver source terms solutio update
5 1 D grid geeratio Remider of uiform grid geeratio x 1/2 x=a x N+1/2 x=b C 1 C 0 C 1 C 2 C N 1 C N C N+1 C N+2 u 1 u 0 u 1 u 2 u N 1 u N+1 u N u N+2 The domai [a,b] has bee divided ito N equal cells C i. The scalar value u i is defied i the midpoit of C i (or it represets a average of u(x) i the cell C i ). C 1, C 0, C N+1, C N+2 : ghost (guard) cells
6 Variable grids x 1/2 x=a x i 3/2 x i 1/2 x i+1/2 x i+3/2 x i 1 x i x i+1 C 1 C 0 C 1 C 2 C i 1 u 1 u 0 u 1 u 2 u i 1 C i u i C i+1 u i+1 x i = x i 1/ 2 x i 1/ 2 All expressios ecoutered last week could be adapted by substitutig the cell size x with x i!
7 Variable grids All expressios ecoutered last week could be adapted by substitutig x with x i! For istace: average over cell u i 1 x i 1/2 x x u x, t dx i 1/2 i stadard discretizatio scheme for coserved laws time step determiatio ad so o u i 1 =u i t x i [ f i 1 /2 f i 1/ 2 ] t=b x i a v
8 Variable grids Notice that last week i the determiatio of the evolutio of desity with slope limiters =[ f v i 1/ 2 i 1/ 2 i 1 = i 3 t f i 1/ 2 f i 1 /2 i 1 v i 1/ r i 1/ 2 r i 1/ 2 i 1 2 v i 1/ 2 1 c ] S r 2 i 1/ 2 i 1 / 2 we have ot assumed that the grids are uiform!
9 Variable grids x i 3/2 x i 1/2 x i+1/2 x i+3/2 x i 1 x i x i+1 C 1 C 0 C 1 C 2 C i 1 u 1 u 0 u 1 u 2 u i 1 C i u i C i+1 u i+1 x i = x i 1/ 2 x i 1/ 2 Thigs to otice: 1. I reflectig boudaries mirrorig cells must be equal, e.g. x 0 = x 1, x 1 = x 2 2. Half idex cell sizes might be eeded (i.e. for the mometum advectio; see also below). They are give by: x i 1 /2 =x i x i 1 =0.5 x i 1/ 2 x i 3 / 2
10 set up (grid based) Hydro code evolutio grid geeratio iitial values iput output storage solutio advacemet boudary values time stepper hydro solver source terms solutio update
11 set up (grid based) Hydro code evolutio grid geeratio iitial values iput output storage solutio advacemet boudary values time stepper hydro solver source terms solutio update
12 Hydro solver outlie More o the Euler equatios Global solutios ad Riema solvers
13 Euler equatios t v =0 m m v P =0 t E t I these equatios we have eglected Exteral forces (e.g. gravity!) Viscosity! Possible source terms [ E P v ]=0
14 Navier Stokes equatios with sources t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė : Viscous stress tesor g: Gravity q: heat flux
15 Navier Stokes equatios with sources t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė : Viscous stress tesor g: Gravity q: heat flux It depeds o the divergece of v ad o the divergece of spatial derivatives of velocities. I 1 D it is proportioal to v xx!
16 Navier Stokes equatios with sources : Viscous stress tesor g: Gravity q: heat flux It depeds o the divergece of v ad o the divergece of spatial derivatives of velocities. I 1 D it is proportioal to v xx! Remember the viscous advectio equatio u t au x = u xx to prevet wave breakig. Itroducig such terms is a cheap way to cosider viscous effects (remember the artificial viscosity of the exercise!) = mi 0, C 1 v r v r
17 Navier Stokes equatios with sources Remember the viscous advectio equatio u t au x = u xx to prevet wave breakig. Itroducig such terms is a cheap way to cosider viscous effects (remember the artificial viscosity of the exercise!) Remember also the modified equatio of the first order upwid scheme u t au x = 1 2 a x 1 c u xx= um u xx it is agai a viscous term! This shows that some form of artificial viscosity is always preset i umerical schemes (i particular firstorder oes). Sometimes this is ot eough ad some extra artificial viscosity must be added. If the built i artificial viscosity is eough to prevet oscillatios, a extra artificial viscosity term is ot eeded.
18 Navier Stokes equatios with sources t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė : Viscous stress tesor g: Gravity q: heat flux g might deped o some (fixed) exteral field (e.g. a fixed dark matter halo) > costat vector! It might be also self gravity of the gas. I this case the (elliptic!) Poisso equatio must be solved! > Complicated 2 =4 G
19 Navier Stokes equatios with sources g might deped o some (fixed) exteral field (e.g. a fixed dark matter halo) > costat vector! It might be also self gravity of the gas. I this case the (elliptic!) Poisso equatio must be solved! > Complicated 2 =4 G i 1 D the calculatio is much simpler sice oe has: g=g M R R 2 g is a vector quatity (defied at the itercell boudary i a staggered grid); M is the mass withi R The self gravity ca be thus very easily implemeted i the 1 D Superova explosio code. However it will appear evidet that gravity is (almost) egligible i this problem.
20 Navier Stokes equatios with sources t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė : Viscous stress tesor g: Gravity q: heat flux Fourier's law relates the heat flux to the temperature gradiet, i.e. q= T (parabolic problem!) The thermal coductivity ca be costat or variable (see exercise) > Solved through Crak Nicolso method
21 Serial vs. global solutio t v =0 m m v P =0 t t [ P v ]=0 This system ca be solved i two ways: serially: three scalar advectio equatios are solved separately globally; the sigle vector advectio equatio W t [ F W ] x = 0, W t A W W x = 0 where F represets the fluxes of the ivolved variables, is solved (bars are omitted hereafter).
22 Serial vs. global solutio W t [ F W ] x =0, W = m, F W = m m u P u P To obtai the global solutio of this vector equatio, we eed some tool to uderstad the global behaviour of the thermodyamic variables. This tool is offered by the solutio of the Riema Problem (RP): W t [ F W ] x =0 W x, 0 =W 0 x = W L W R x 0 x 0 Two (costat!) thermodyamic states i cotact.
23 Riema Problem Replace t=0 with t ; x=0 with x i 1/2 (chage of referece frame), W L with W i 1, W R with W i ad you will immediately realize why is the Riema problem so relevat to the solutio of the Euler equatios. t t 1 W t [ F W ] x =0 W x, t =W 0 x = W i 1 x x i 1 / 2 W i x x i 1 /2 1 /2 t W i 1 W i Costat states i cotact! t x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x
24 Riema Problem W t [ F W ] x =0 W x, 0 =W 0 x = W L W R x 0 x 0 Remember Sod shock problem (Nigel's lecture), m, = 1, 0, 2.5 x , 0, 0.25 x 0.5 Simple RP (iitial velocities are zero).
25 Riema Problem, m, = 1, 0, 2.5 x , 0, 0.25 x 0.5 Simple RP (iitial velocities are zero). Nedwich's solutio.
26 Riema Problem New itermediate states have bee produced! The solutio of the RP is self similar! I.e. it develops alog rays x x itf t t itf =costat
27 RP ad Goduov method How ca we exploit the solutio of the RP? We solve the RP at each iterface ad average the ew state at each cell! t t 1 W i 1 1 /2 t t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x (remember that each state ca be defied as cell average, e.g.) W 1 i 1 x i 1/ 2 x x W x, t 1 dx i 1/ 2 i
28 RP ad Goduov method t t 1 W i 1 1 /2 t t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x W 1 i 1 x i 1/ 2 x x W x, t 1 dx i 1/ 2 i Problems: This itegral is ofte difficult to calculate We do ot kow what to do if rays itersect, i.e. Courat coditio is: a v t x 1 2
29 RP ad Goduov method t t 1 R i 1 /2 R i 1 /2 1 /2 t t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x To simplify we ca cosider oly the solutio of the RP at the iterfaces! (We idicated it with R i 1/2 ). At this poit we cosider the (Goduov) flux with state W equal to these itermediate states, i.e.: W i 1 =W i t God x i [ F i 1 /2 F i 1 / 2 ], F i 1 /2 =F i 1 / 2 = F R i 1/ 2 Courat umbers must be (as usual) smaller tha 1 ad ot of 0.5!
30 RP ad Goduov method t t 1 R i 1 /2 R i 1 /2 1 /2 t t u i 1 u i u i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x Example: advectio equatio u t au x =0 with a>0 (remember: advected quatities remai ualtered!). Clearly: R i 1 /2 =u i 1 u 1 i =u i t > Upwid scheme!, R i 1 / 2 =u i x i [ F u i 1 F u i ]=u i c u i 1 u i
31 RP ad Goduov method t t 1 R i 1 /2 R i 1 /2 1 /2 t t u i 1 u i u i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x Example: advectio equatio u t au x =0 with a<0. R i 1 /2 =u i, R i 1/ 2 =u i 1 u 1 i =u i t x i [ F u i F u i 1 u 1 i =u i c u i u i 1 ]=u i c u i u i 1 c u i 1 u i
32 RP ad Goduov method t t 1 R i 1 /2 R i 1 /2 1 /2 t t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x Example: cotiuity equatio t v x =0 states W i = i, v i T with piecewise costat This is ot a simple extesio of the previous example: R i 1/2 ad R i+1/2 represet ew states! Two waves are emaated from the iterfaces, oe with velocity v i (v i+1 ) ad oe with velocity 0!
33 RP ad Goduov method Two waves are emaated from the iterfaces, oe with velocity v i (v i+1 ) ad oe with velocity 0! Thik at coveyor belts with differet speeds t v i v i+1 v i v i+1 v i v i+1 v i v i+1
34 RP ad Goduov method Two waves are emaated from the iterface, oe with velocity v i (v i+1 ) (shock!) ad oe with velocity 0 (cotact!) A ew state has bee created! t cotact shock v i v i+1
35 RP ad Goduov method t R i 1 /2 W i 1/ 2 R i 1 /2 W i 1/ 2 t 1 1 /2 t t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x How ca we solve the RP (i.e. calculate the itermediate states W*) the? > Coservatio of desity flux at the iterface! i v i = i 1 /2 v i 1, i 1/ 2 = v i i v i 1 i 1 = i t x i [ f i 1/ 2 f i 1/ 2 ]= i t x i i 1 v i 1 i v i
36 RP ad WAF methods Secod way of exploitig the solutio of the RP: we average the fluxes at half timestep! (Weighted average flux WAF methods) t t 1 1 /2 t WAF F i 1 / 2 WAF F i 1 / 2 t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x
37 RP ad WAF methods t t 1 1 /2 t WAF F i 1 / 2 WAF F i 1 / 2 t W i 1 W i W i 1 x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x i 1 x i 3 / 2 x (remember that we wished to defie fluxes at half iteger time idices staggered grid!) The WAF flux ca be defied as: F WAF i 1 / 2 1 x i x x i 1 i 1 /2 F [W x, t 1 / 2 ] dx where W(x,t +1/2 ) is the solutio of the RP evaluated after half a temporal timestep.
38 RP ad WAF methods F i 1 / 2 WAF t 1 /2 t 2 1 F F 3 i 1 / 2 i 1 / 2 F i 1 / 2 4 F i 1 / 2 t W i 1 W i x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x If N waves emaate from a RP, the WAF flux ca be approximated as: F WAF i 1 / 2 1 x i x x F [W x, t 1 / 2 N 1 k i 1 ] dx = k k =1 F i 1/ 2 i 1 /2 F (k) : itermediate fluxes (fluxes of itermediate states); k : weights
39 RP ad WAF methods W i 1 /2 WAF t 1 /2 t 2 1 W W 3 i 1 /2 i 1 /2 W i 1 /2 4 W i 1 /2 t W i 1 W i x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x Equivalet approach: recostruct a average WAF state W WAF ad the calculate its flux W WAF i 1 /2 1 x i x xi W x, t 1/ 2 dx 1 = k=1 i 1/ 2 F WAF i 1/ 2 = F W i 1/ 2 WAF N 1 k k W i 1 / 2
40 RP ad WAF methods t 1 /2 t 2 1 W W 3 i 1 /2 i 1 /2 W i 1 /2 4 W i 1 /2 A 0 A 1 A 2 A 3 A t W i 1 W i x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x Weights k are easy to calculate if the slopes k (propagatio speeds) are kow. They are proportioal to the legth of the segmets A k 1 A k, i.e.: k = A k 1 A k x i 1/ 2
41 RP ad WAF methods Weights k are easy to calculate if the slopes k (propagatio speeds) are kow (removed idex from x for simplicity). It is: k = A k 1 A k x k = 1 2 c i 1/ 2 Oe ca thus obtai: k 1 k ; c i 1 / 2 = t k x, c0 = 1, c N 1 =1 c i 1 /2 F WAF N 1 k i 1 / 2 = k k =1 F i 1 / 2 = 1 2 F i 1 F i 1 2 N k=1 k k c i 1/ 2 F i 1/ 2 k k F i =F W i, F i 1 = F W i 1, F i 1 /2 =F 1 k i 1 / 2 F i 1/ 2 (otice that F i 1/2 (1) =F i 1, F i 1/2 (N+1) =F i ) Alteratively: W WAF i 1 /2 = 1 2 W i 1 W i 1 2 k =1 k k 1 W i 1 /2 =W k i 1/ 2 W i 1/ 2 N k k c i 1 / 2 W i 1 / 2, F WAF WAF i 1/ 2 = F W i 1 /2
42 RP ad WAF methods t 1 /2 t 1 u i 1/ 2 2 u i 1/ 2 A 0 A 1 A 2 1 =a t u i 1 u i x i 3 / 2 x i 1 x i 1 /2 x i x i 1 /2 x Example: advectio equatio k = A k 1 A k x, 1 = c, 2 = c f WAF i 1/ 2 = 1 2 [ au i 1 1 c au i 1 c ] This is the Lax Wedroff flux!
43 Flux limited WAF methods k = A k 1 A k x, 1 = c, 2 = c f WAF i 1/ 2 = 1 2 [ au i 1 1 c au i 1 c ] This is the Lax Wedroff flux! How ca we costruct flux limited (i.e. o oscillatig) WAF methods? We remid the costructio of scalar flux limited methods: we realized that the Lax Wedroff method ca be obtaied from a piecewise liear data recostructio with slopes give by the dowwid slopes. We the modified the data jumps S through a flux limiter fuctio i order to make the iitial data mootoe. f i 1/ 2 =a u i 1 a u i 1 2 a 1 c S i 1 / 2 S i 1 /2 = i 1 / 2 u i 1 /2
44 Flux limited WAF methods f i 1/ 2 =a u i 1 a u i 1 2 a 1 c S i 1 / 2 S i 1 /2 = i 1 / 2 u i 1 /2 sice: F WAF i 1 / 2 = 1 2 F i 1 F i 1 2 k=1 N k c i 1/ 2 k F i 1 /2 by aalogy with the scalar case we ca assume that the (mootoicity preservig) WAF flux is: WAF = 1 2 F i 1 F i 1 2 k=1 F i 1 / 2 S i 1/2 (k) : limited versio of the k th jump, i.e. k S i 1 /2 N k k [1 1 c ] k i 1 / 2 S i 1/ 2 k k = i 1 / 2 W i 1/ 2
45 Flux limited WAF methods WAF = 1 2 F i 1 F i 1 2 k=1 F i 1 / 2 N k k [1 1 c ] k i 1 / 2 S i 1/ 2 Alteratively: W WAF i 1 /2 = 1 2 W i 1 W i 1 2 k =1 WAF WAF = F W i 1/ 2 F i 1/ 2 N k sig c i 1/ 2 k k [1 1 c i 1 /2 ]S i 1/ 2 S i 1/2 (k) : still limited versio of the k th jump, i.e. k S i 1 /2 k k = i 1 / 2 W i 1/ 2
46 Flux limited WAF methods S i 1/2 (k) : still limited versio of the k th jump, i.e. k S i 1 /2 k k = i 1 / 2 W i 1/ 2 For scalar problems we had: i 1/ 2 = u I 1/ 2, I = u i 1/ 2 i 1 a 0 i 1 a 0 : ratio of the upwid jump to the local jump. I vector problems is the ratio of wave jumps, but wave jumps are vectors! How ca we defie their divisio? (Geerally they are ot parallel!)
47 Flux limited WAF methods : ratio of the upwid jump to the local jump. I vector problems is the ratio of wave jumps, but wave jumps are vectors! How ca we defie their divisio? 1. Oe selected sigle quatity which is kow to chage across waves (for istace desity) is chose to calculate the jumps, i.e. k i 1/ 2 = k I 1 / 2 i 1 /2, I= i 1 c i 1/ 2 0 k k i 1 c i 1/ The vectors are projected oto a preferetial directio (for istace the oe of the local jump), therefore oe has: k i 1/ 2 = W k I 1/ 2 k W i 1/ 2 k k W i 1/ 2, I = i 1 c k i 1/ 2 0 k k W i 1 / 2 i 1 c i 1/ 2 0
48 Flux limited WAF methods k i 1/ 2 k i 1/ 2 = k I 1 / 2 i 1 /2 = W k I 1/ 2 k W i 1/ 2, I= i 1 c i 1/ 2 0 k k i 1 c i 1/ 2 0 Limiter fuctios are the same as i the scalar problem, for istace: k k W i 1/ 2, I = i 1 c k i 1/ 2 0 k k W i 1 / 2 i 1 c i 1/ 2 0 Va Leer limiter: mimod = 1 =max [0, mi 1, ] superbee =max [0, mi 1, 2, mi 2, ]
49 Flux limited WAF methods Thigs to otice: 1. W ca be the vector of coserved quatities W=(,m, ) T, but it ca be also the vector of primitive quatities W=(,v,P) T! 2. To calculate the timestep to advace the solutio as usual the Courat umber must be smaller tha 1. We must however take ito accout the speeds of all the waves, emaatig from all the iterfaces, take the fastest at the time level (S max () ) ad calculate the timestep accordig to geerally S max () is well approximated by S max t=b x S max =max i [ v i a i ] but fastest waves are always possible. Commo practice: start with a very low b, icrease it i the first few timesteps util reachig 0.4 or 0.5
50 Flux limited WAF methods: summary Calculate the timestep x t=b max i [ v i a i ] assumig a small value of b durig the first timesteps For each pair of data (W i 1,W i ) solve the RP. Calculate itermediate (k) (k) (k) states W i 1/2, fluxes F i 1/2, jumps W i 1/2 ad wave speeds k Calculate Courat umbers related to each wave: c i 1 / 2 k = t k x
51 Flux limited WAF methods: summary Calculate the wave jump ratios accordig to oe of the followig formulae: k i 1/ 2 k i 1/ 2 Calculate limited jumps Va Leer limiter: mimod = k I 1 / 2 i 1 /2 = W k I 1/ 2 k W i 1/ 2 k S i 1 /2 k, I= i 1 c i 1/ 2 0 k k i 1 c i 1/ 2 0 k W i 1/ 2, I = i 1 c k i 1/ 2 0 k k W i 1 / 2 i 1 c i 1/ 2 0 k k = i 1 / 2 W i 1/ 2 = 1 =max [0, mi 1, ] superbee =max [0, mi 1, 2, mi 2, ]
52 Flux limited WAF methods: summary Calculate the itercell flux accordig to or: WAF = 1 2 F i 1 F i 1 2 k=1 F i 1 / 2 W WAF i 1 /2 = 1 2 W i 1 W i 1 2 k =1 WAF WAF = F W i 1/ 2 F i 1/ 2 Update the solutio Update the boudaries Go to the ext time level N N k k [1 1 c ] k i 1 / 2 S i 1/ 2 k sig c i 1/ 2 W i 1 =W i t x i [ F i 1 /2 F i 1 / 2 ] k k [1 1 c i 1 /2 ]S i 1/ 2
53 Solutio of the RP i pills The RP solutio is self similar! (It develops alog rays of costat x/t) The Euler equatio produces three waves, propagatig at speeds 1 = v a, 2 =v, 3 =v a The cetral wave is always a cotact! The other two waves ca be shocks or rarefactio waves
54 Solutio of the RP i pills Across a cotact oly the desity chages! Velocity ad pressure remai costat (remember the coveyor belts). P ad v are thus also called Geeralized Riema ivariats of the cotact wave. t cotact shock v i v i+1
55 Solutio of the RP i pills Shocks are zero legth discotiuities. Across them, the fluxes of mass, mometum ad eergy are coserved. This leads to the Rakie Hugoiot jump shock coditios! (M: Mach umber) P 2 = 2 P 1 1 M v 2 = 1 v M = 1 1 2/ M 2 Rarefactio waves are regios of liear variatios of the thermodyamic quatities. The Geeralized Riema ivariats of rarefactio waves are: v 2a / 1, v 2a / 1, across 1 =v a across 3 =v a
56 Solutio of the RP i pills RP solvers ca be either exact (obtaied through a iterative procedure) or approximate (usually obtaied through a liearizatio of the Euler equatios). I approximate Riema solvers the rarefactio wave is usually approximated as a sigle state. Approximate structure of the Riema solutio (W=(,m, )T ): t S L S S R W L W L W R W R W L W R remember that oly the desity chages across the cotact! x
57 Solutio of the RP i pills The HLLC approximate Riema solver: P =max [0,0.5 P L P R 0.5 v R v L a ], =0.5 L R, a=0.5 a L a R S L = v L a L q L, S R = v R a R q R, q L, R = 1 P P [ L, R 2 1 P 1 1 ]1/ P 2 P P L, R L, R S = P R P L L v L S L v L R v R S R v R L S L v L R S R v R W L, R = R S L, R v L, R L, S L, R S 1 S [ L, R S v S L, R L, R ] P L, R L, R S L, R v L, R F L, R = F L, R S L, R W L, R W L, R
58 WAF methods last remarks 1. I spherical symmetry, as i the case of scalar problems, the Euler equatios ca be solved as: W i 1 =W i 3 t F i 1 / 2 F i 1 / r i 1 / 2 r i 1 / 2 ad the WAF fluxes aalyzed so far must be multiplied by r 2, i.e. F WAF i 1 / 2 F WAF 2 i 1 / 2 r i 1 /2 alteratively, oe ca solve the plaar Euler equatio with (geometric) source terms W t [ F W ] x =S W W = m, F W = m m v P u P, S W = 2 r v v 2 v P
59 WAF methods last remarks 2. I global solutios of the Euler equatio it is atural to defie scalar ad vector quatities (like m) at the grid ceter sice here is where the state W is defied > it is difficult to defie staggered grid as we did for serial solutios. However, this is ot importat because we shifted the defiitio of velocities i order to calculate fluxes at iterfaces! This is what we did also with WAF methods. t t 1,, P 1 /2 t m, v, t x i 1 /2 x i x i 1 /2 x
60 set up (grid based) Hydro code evolutio grid geeratio iitial values iput output storage solutio advacemet boudary values time stepper hydro solver source terms solutio update
61 Source terms t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė If you are solvig the equatios serially, you ca defie source terms where the correspodet quatities are defied (grid ceters for desity ad eergy, itercell boudaries for mometum). If you are solvig globally, the Euler equatio with source terms look like that: W t [ F W ] x =S W
62 Source terms W t [ F W ] x =S W Basic strategy: time splittig method! We first solve the pure advectio hyperbolic problem: W t [ F W ] x =0 W the we solve the system of Ordiary differetial equatios d W dt =S W usig the solutio of the pure advectio problem as iitial coditio. Ay of the methods leared i the first two discretizatio lectures could be theoretically adopted to solve this system
63 Source terms Example: scalar advectio equatio with source term u t au x = u the aalytical solutio is: u x, t =u 0 x at e t where u 0 represets the iitial coditios. The solutio of the pure advectio problem is: u t au x =0, u x, t =u 0 x at If we ow solve the ODE u'= u usig as iitial coditios the solutio of the pure advectio problem we obtai: u x, t = u x, t e t =u 0 x at e t
64 Source terms W t [ F W ] x =S W W t [ F W ] x =0 W i 1 =W i t x [ F i 1/ 2 F i 1/ 2 ] At this poit the system of ODEs d W dt =S W ca be solved, for istace with the Euler method W i 1 = W i 1 t S t, W i 1
65 W t [ F W ] x =0 Source terms W i 1 =W i t x [ F i 1/ 2 F i 1/ 2 ] W i 1 = W i 1 t S t, W i 1 t v =0 m m v P =0 t t v = P v Oce the pure advectio equatios have bee solved, the coolig ca be thus dealt with i this way: 1 i = 1 i t i i 1 m H c V i m H 2 T
66 Source terms W t [ F W ] x =0 W i 1 =W i t x [ F i 1/ 2 F i 1/ 2 ] W i 1 = W i 1 t S t, W i 1 This is sufficiet for first order (e.g. Goduov) methods. For (almost) secod order methods like the WAF methods it is more balaced to look for a secod order ODE solver. Two possibilities: 1. Heu method or similar schemes W i 1 = W i 1 t S t, W i 1 S [t 1, W i 1 t S t, W i 1 ] 2
67 Source terms 1. Heu method or similar schemes W 1 i = W 1 i t S t, W 1 i S [t 1, W 1 i t S t, W 1 i ] 2 2. A secod order accurate splittig scheme i which first the ODE for the source term is solved advacig the solutio by half a timestep (with a first order method) the the solutio is advected, the the solutio of the ODE is advaced by the remaiig half of the timestep. This ca be represeted by the sequece W i 1 =S t / 2 A t S t / 2 W i S: operator solutio of the ODE; A: operator solutio of the advectio
68 Source terms W i 1 =S t / 2 A t S t / 2 W i I order to update the eergy accordig to the coolig fuctio, we thus eed the followig 3 steps: 1. 1 /2 i = i t i 2 m H 2 i c V i 2. i 1 = i 1 / 2 t x [ f i 1/ 2 f i 1 / 2 ] 3. 1 i = 1 i t 1 i 2 m H 2 1 i 1 c V i
69 Source terms ad time steppig t v = m m v P = g ṁ t E t [ E P v ]= v q v g Ė Each source term itroduces a ew typical timescale (i additio to the Courat timescale), amely: t =, t m= ṁ m, t E t c = Ė E All these timescales must be take ito accout i the time steppig routie!
70 Source terms ad time steppig t =, t m= ṁ m, t E t c = Ė E All these timescales must be take ito accout i the time steppig routie! I order ot to be boud by these timesteps, a implicit method (for istace backward Euler) could be employed (remember: implicit methods are ucoditioally stable!). I the case of the coolig fuctio: 1 i = 1 i t i i 1 m H c V i
71 Thak you for the attetio!
Chapter 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 informationNumerical Astrophysics: hydrodynamics
Numerical Astrophysics: hydrodyamics Part 1: Numerical solutios to the Euler Equatios Outlie Numerical eperimets are a valuable tool to study astrophysical objects (where we ca rarely do direct eperimets).
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
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 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 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 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 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 information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
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 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 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 informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
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 informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
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 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 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 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 informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:
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 informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationMAT1026 Calculus II Basic Convergence Tests for Series
MAT026 Calculus II Basic Covergece Tests for Series Egi MERMUT 202.03.08 Dokuz Eylül Uiversity Faculty of Sciece Departmet of Mathematics İzmir/TURKEY Cotets Mootoe Covergece Theorem 2 2 Series of Real
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 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 informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
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 informationINFINITE SEQUENCES AND SERIES
INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES I geeral, it is difficult to fid the exact sum of a series. We were able to accomplish this for geometric series ad the series /[(+)]. This is
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 informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationEE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course
Sigal-EE Postal Correspodece Course 1 SAMPLE STUDY MATERIAL Electrical Egieerig EE / EEE Postal Correspodece Course GATE, IES & PSUs Sigal System Sigal-EE Postal Correspodece Course CONTENTS 1. SIGNAL
More informationLesson 10: Limits and Continuity
www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals
More informationQuestion 1: The magnetic case
September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS
SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS I a exothermic reactio the temperature will cotiue to rise as oe moves alog a plug flow reactor util all of the limitig reactat is exhausted. Schematically
More informationOnce we have a sequence of numbers, the next thing to do is to sum them up. Given a sequence (a n ) n=1
. Ifiite Series Oce we have a sequece of umbers, the ext thig to do is to sum them up. Give a sequece a be a sequece: ca we give a sesible meaig to the followig expressio? a = a a a a While summig ifiitely
More informationChapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:
Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets
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 informationMath 113 Exam 4 Practice
Math Exam 4 Practice Exam 4 will cover.-.. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for
More informationLinear Regression Demystified
Liear Regressio Demystified Liear regressio is a importat subject i statistics. I elemetary statistics courses, formulae related to liear regressio are ofte stated without derivatio. This ote iteds to
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
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 informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
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 informationChapter 7: Numerical Series
Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals
More informationAdvanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology
Advaced Aalysis Mi Ya Departmet of Mathematics Hog Kog Uiversity of Sciece ad Techology September 3, 009 Cotets Limit ad Cotiuity 7 Limit of Sequece 8 Defiitio 8 Property 3 3 Ifiity ad Ifiitesimal 8 4
More informationIntegrable Functions. { f n } is called a determining sequence for f. If f is integrable with respect to, then f d does exist as a finite real number
MATH 532 Itegrable Fuctios Dr. Neal, WKU We ow shall defie what it meas for a measurable fuctio to be itegrable, show that all itegral properties of simple fuctios still hold, ad the give some coditios
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationTrue Nature of Potential Energy of a Hydrogen Atom
True Nature of Potetial Eergy of a Hydroge Atom Koshu Suto Key words: Bohr Radius, Potetial Eergy, Rest Mass Eergy, Classical Electro Radius PACS codes: 365Sq, 365-w, 33+p Abstract I cosiderig the potetial
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 informationMATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and
MATH01 Real Aalysis (2008 Fall) Tutorial Note #7 Sequece ad Series of fuctio 1: Poitwise Covergece ad Uiform Covergece Part I: Poitwise Covergece Defiitio of poitwise covergece: A sequece of fuctios f
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationBoundary layer problem on conveyor belt. Gabriella Bognár University of Miskolc 3515 Miskolc-Egyetemváros, Hungary
Boudary layer problem o coveyor belt Gabriella Bogár Uiversity of Miskolc 355 Miskolc-Egyetemváros, Hugary e-mail: matvbg@ui-miskolc.hu Abstract: A techologically importat source of the boudary layer pheomeo
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationΩ ). Then the following inequality takes place:
Lecture 8 Lemma 5. Let f : R R be a cotiuously differetiable covex fuctio. Choose a costat δ > ad cosider the subset Ωδ = { R f δ } R. Let Ωδ ad assume that f < δ, i.e., is ot o the boudary of f = δ, i.e.,
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES 11.4 The Compariso Tests I this sectio, we will lear: How to fid the value of a series by comparig it with a kow series. COMPARISON TESTS
More informationCHAPTER I: Vector Spaces
CHAPTER I: Vector Spaces Sectio 1: Itroductio ad Examples This first chapter is largely a review of topics you probably saw i your liear algebra course. So why cover it? (1) Not everyoe remembers everythig
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More information10.6 ALTERNATING SERIES
0.6 Alteratig Series Cotemporary Calculus 0.6 ALTERNATING SERIES I the last two sectios we cosidered tests for the covergece of series whose terms were all positive. I this sectio we examie series whose
More informationOptimization Methods MIT 2.098/6.255/ Final exam
Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short
More informationENGI Series Page 6-01
ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,
More informationIntroduction to Extreme Value Theory Laurens de Haan, ISM Japan, Erasmus University Rotterdam, NL University of Lisbon, PT
Itroductio to Extreme Value Theory Laures de Haa, ISM Japa, 202 Itroductio to Extreme Value Theory Laures de Haa Erasmus Uiversity Rotterdam, NL Uiversity of Lisbo, PT Itroductio to Extreme Value Theory
More informationMechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter
Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationSeries III. Chapter Alternating Series
Chapter 9 Series III With the exceptio of the Null Sequece Test, all the tests for series covergece ad divergece that we have cosidered so far have dealt oly with series of oegative terms. Series with
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationAn Analysis of a Certain Linear First Order. Partial Differential Equation + f ( x, by Means of Topology
Iteratioal Mathematical Forum 2 2007 o. 66 3241-3267 A Aalysis of a Certai Liear First Order Partial Differetial Equatio + f ( x y) = 0 z x by Meas of Topology z y T. Oepomo Sciece Egieerig ad Mathematics
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
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 informationMa 530 Infinite Series I
Ma 50 Ifiite Series I Please ote that i additio to the material below this lecture icorporated material from the Visual Calculus web site. The material o sequeces is at Visual Sequeces. (To use this li
More information(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3
MATH 337 Sequeces Dr. Neal, WKU Let X be a metric space with distace fuctio d. We shall defie the geeral cocept of sequece ad limit i a metric space, the apply the results i particular to some special
More informationCALCULUS BASIC SUMMER REVIEW
CALCULUS BASIC SUMMER REVIEW NAME rise y y y Slope of a o vertical lie: m ru Poit Slope Equatio: y y m( ) The slope is m ad a poit o your lie is, ). ( y Slope-Itercept Equatio: y m b slope= m y-itercept=
More informationSolution. 1 Solutions of Homework 1. Sangchul Lee. October 27, Problem 1.1
Solutio Sagchul Lee October 7, 017 1 Solutios of Homework 1 Problem 1.1 Let Ω,F,P) be a probability space. Show that if {A : N} F such that A := lim A exists, the PA) = lim PA ). Proof. Usig the cotiuity
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More informationMATH4822E FOURIER ANALYSIS AND ITS APPLICATIONS
MATH48E FOURIER ANALYSIS AND ITS APPLICATIONS 7.. Cesàro summability. 7. Summability methods Arithmetic meas. The followig idea is due to the Italia geometer Eresto Cesàro (859-96). He shows that eve if
More informationSection 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations
Differece Equatios to Differetial Equatios Sectio. Calculus: Areas Ad Tagets The study of calculus begis with questios about chage. What happes to the velocity of a swigig pedulum as its positio chages?
More informationChapter 3. Strong convergence. 3.1 Definition of almost sure convergence
Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i
More informationTest One (Answer Key)
CS395/Ma395 (Sprig 2005) Test Oe Name: Page 1 Test Oe (Aswer Key) CS395/Ma395: Aalysis of Algorithms This is a closed book, closed otes, 70 miute examiatio. It is worth 100 poits. There are twelve (12)
More informationInverse Matrix. A meaning that matrix B is an inverse of matrix A.
Iverse Matrix Two square matrices A ad B of dimesios are called iverses to oe aother if the followig holds, AB BA I (11) The otio is dual but we ofte write 1 B A meaig that matrix B is a iverse of matrix
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 15
CS 70 Discrete Mathematics ad Probability Theory Summer 2014 James Cook Note 15 Some Importat Distributios I this ote we will itroduce three importat probability distributios that are widely used to model
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationStatistical Pattern Recognition
Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig
More informationChapter 8. Euler s Gamma function
Chapter 8 Euler s Gamma fuctio The Gamma fuctio plays a importat role i the fuctioal equatio for ζ(s) that we will derive i the ext chapter. I the preset chapter we have collected some properties of the
More informationOptimally Sparse SVMs
A. Proof of Lemma 3. We here prove a lower boud o the umber of support vectors to achieve geeralizatio bouds of the form which we cosider. Importatly, this result holds ot oly for liear classifiers, but
More informationThe Riemann Zeta Function
Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationNUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK
NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK For this piece of coursework studets must use the methods for umerical itegratio they meet i the Numerical Methods module
More informationIP Reference guide for integer programming formulations.
IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
More information