Accuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics
|
|
- Barrie Davis
- 5 years ago
- Views:
Transcription
1 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, there are sigiicat quatitative diereces. The derivatio o the umerical approximatios or the derivatives showed that the error depeds o the size o h ad. First we test or dieret. Evolutio or U=; D=.5; k= N= =.5 m=; time= Number o time steps= T/ Exact Numerical Repeat with a smaller time-step =.5.5 m=; time=.5 Repeat with a smaller time-step =.5.5 m=4; time=.5 N=.5 N= Exact - Exact - Numerical Numerical
2 How accurate solutio ca we obtai? Take =.5 ad N= Very ie spatial resolutio ad a small time step U=; D=.5; k= N= = m=; time=.5 Exact - Numerical Quatiyig the Error Order o Examie the spatial accuracy by takig a very small time step, =.5 ad vary the umber o grid poits, N, used to resolve the spatial directio. The grid size is h = L/N where L = or our case Exact Numerical N E = h ( exact ) = time=.5 N=; E =.633 N =; E =.43 N =4; E =.96 N =6; E =.4 N =8; E =. N =; E =.5 N =; E =. N =6; E = 9.6e-4. Eect o spatial resolutio dt=.5 N= to N=6 - time=.5 I the error is o secod order: E = Ch = C h E - -3 Takig the log: l E = l C h = lc l h -4 / h O a log-log plot, the E versus (/h) curve should thereore have a slope -
3 . Eect o spatial resolutio dt=.5 N= to N=6 E time=.5-4 / h Why is does the error deviate rom the lie or the highest values o N? Trucatio Error Roudo Total Number o Steps, N = /h Summary Fiite dierece approximatios by Taylor expasio Approximatig a partial dieretial equatio Showed, by umerical experimets, that accuracy icreases as the resolutio is icreased Showed that the error behaves i a way that should be predictable As log as accuracy is reasoable, itegratio at larger time steps is more eiciet ad desirable. Ca we icrease the time step ideiitely? Let s repeat the -D advectio-diusio equatio with larges time step. Use =., istead o =.5 Evolutio or U=; D=.5; k= N= =. Exact Numerical
4 Istead o decayig as it should, the amplitude o the umerical solutio keeps icreasig. Ideed, i we cotiued the calculatios, we would evetually produce umbers larger tha the computer ca hadle. This results i a overlow or NaN (Not a Number). Ordiary Dieretial Equatio Take: d dt = The exact solutio is ( t) = e t Forward Euler ODE Example + = with iitial coditio () = + oly i + = ( ) + = ( ) =.5 =.5 =.5 ODE Example + oly i However Backward Euler + = + ODE Example + Obviously, oscillates uless + = ( ) = ( ) = = ( ) or all + = ( + ) Aalysis o the Advectio-Diusio Equatio: Vo Neuma Method
5 Geerally, stability aalysis o the ull oliear system o equatios is too ivolved to be practical, ad we study a model problem that i some way mimics the ull equatios. The liear advectio-diusio equatio is oe such model equatio, ad we will apply vo Neuma's method to check the stability o a simple iite dierece approximatio to that equatio. Cosider the -D advectio-diusio equatio: I iite-dierece orm: + + U = D t x x + U + + = D h h + Look at the evolutio o a small perturbatio = The evolutio o the perturbatio is govered by: + + U + h = D + + h Write the error as a wave (expad as a Fourier series): = ( x ) = Droppig the subscript ikx k e k = ikx = e Recall: e ikx = coskx + isikx The error at ode is: = e ikx + The error at + ad - ca be writte as = e ikx + = e ik (x + h) = e ikx e ikh = e ikx = e ik (x h) = e ikx e ikh Substitutig = e ikx ito yields e ikx e ikx = e ikx e ikh + U + h + U h (eikh e ikx e ikh e ikx ) = D h (eikh e ikx e ikx + e ikh e ikx ) = e ikx e ikh = D + + h + The equatio or the error is: + + U h (eikh e ikh ) = D Solvig or the ratio o the errors: U ikh ikh D = ( e e ) + ( e h h h (eikh + e ikh ) ikh + e ikh )
6 + Dividig by the error amplitude at : Usig: e ikh + e ikh = coskh; ampliicatio actor U ikh ikh D ikh = ( e e ) + ( e h h = U D isikh + (coskh ) h h = 4 D si k h h i U h sikh e ikh e ikh = isikh; + e ikh ) si θ = cosθ + = 4 D h requires that The ratio o the error amplitude at + ad is: si k h i U h sikh + Sice the ampliicatio actor is a complex umber, ad k, the wave umber o the error, ca be aythig, the determiatio o the stability limit is slightly ivolved. We will look at two special cases: (a) U = ad (b) D = (a) Cosider irst the case whe U =, so the problem reduces to a pure diusio + = 4 D h si k h Sice si () the ampliicatio actor is always less tha, ad we id that it is bigger tha - i 4 D h D h (b) Cosider ow the other limit where D = ad we have a pure advectio problem. + = i U h sikh Sice the ampliicatio actor has the orm +i() the absolute value o this complex umber is always larger tha uity ad the method is ucoditioally ustable or this case. + i U h sikh For the geeral case we must ivestigate the stability coditio i more detail. We will ot do so here, but simply quote the results: For a two-dimesioal problem, assume a error o the orm A stability aalysis gives: i, = i(kxi +ly ) e D h ad U D D h 4 ( U + V ) ad 4 D Notice that high velocity ad low viscosity lead to istability accordig to the secod restrictio. For a three-dimesioal problem we get: D h 6 ad ( U + V + W ) 8 D
7 Now you kow! Covergece the solutio to the iite-dierece equatio approaches the true solutio to the PDE havig the same iitial ad boudary coditios as the mesh is reied. Lax s Equivalece Theorem Give a properly posed iitial value problem ad a iitedierece approximatio to it that satisies the cosistecy coditio, stability is the ecessary ad suiciet coditio or covergece. The Modiied Equatio Cosistecy Usig the iite dierece approximatio, we are eectively solvig a equatio that is slightly dieret tha the origial partial dieretial equatios. Does the iite dierece equatio approach the partial dieretial equatio i the limit o zero ad h? Cosistecy Cosider the -D advectio-diusio equatio ad its iite-dierece approximatio + t + U x = D x +U + h = D + + h The discrepacy betwee the two equatios ca be oud by derivig the modiied equatio. Substitutig + Cosistecy + = + = h + h (t) t + (t) t + (x) + 3 (x) h x x = (x) x ito the iite dierece equatio + + U + = D + h + 4 (x) x 4 h + + h t + U x D x Cosistecy Results i the Modiied Equatio: Origial Equatio = (t) t U 3 (x) h x 3 Error terms Shorthad: t + U x D x = O(,h ) 6 + D 4 (x) x 4 h + I this case, the error goes to zero as h ad, so the approximatio is said to be CONSISTENT
8 Cosistecy Although most iite dierece approximatios are cosistet, iocet-lookig modiicatios ca sometimes lead to approximatios that are ot! The Frakel-Duort is a example o a o-cosistet scheme. You will examie it i the homework Cosistecy HW: Examie the Frakel-Duort method Solve the diusio equatio t = D x Usig the Leaprog time itegratio method ad stadard iite-dierece approximatio or the spatial derivative gives: + = D h + + [ ] Leaprog method Now modiy it slightly: Cosistecy + = D h + + [ ] Replace by: = + ( + ) This gives: + = + D h + + ( + ). Which is easily solved or at the ew time step Cosistecy The MODIFIED EQUATION is obtaied by substitutig the expressio or the iite dierece approximatios, icludig the error terms, ito the iite dierece equatio. For a CONSISTENT iite dierece approximatio the error terms go to zero as h ad. The modiied equatio ca ote be used to ier the ature o the error o the iite dierece scheme. More about that later. I the HW, you will examie the error! Multidimesioal Equatios Two-Dimesioal Advectio- Diusio Equatio We will use the model equatio: +U t x + V y = D x + y to demostrate how to solve a partial equatio (iitial value problem) umerically. The extesio to two-dimesios is relatively straight orward, oce the oe-dimesioal problem is ully uderstood.
9 Multidimesioal Equatios For a two-dimesioal low discretize the variables o a two-dimesioal grid + - (x, y) i - i i+ i, + = (x, y + h) i, = (x,y) i +, = (x + h, y) + i, i, Multidimesioal Equatios t + U x + V y = D x The discrete equatio is: = U i+, h D i, i+, + y V i, + h i, + i, i, + i, + h + + i, h i, + Solve or i, or + i, = i, + i, = i, Multidimesioal Equatios + U i+, h D i+, U h : O,h i, + i, V i, + h h + i, + h i, + i, ( i+, i, ) V ( i, + i, ) ( ) + D h A stability aalysis gives: + 4 i, i+, + i, + i, + + ( i, 4 i, ) D h 4 ad ( U + V ) 4 D Example Multidimesioal Equatios t + U x + V y = D x + y + i, = i, U h Multidimesioal Equatios ( i+, i, ) + D h i +, =NY + i, + i, + + ( i, 4 i, ) x = y = =. =. Uiorm low through the domai U=- V= i, stored at each grid poit or give o the boudary =. = = i= i= i=nx
10 Boudary coditios Multidimesioal Equatios Where is give, we simply speciy its value Where the ormal derivative is speciied, we approximate the value at the boudary by oesided diereces At the i= boudary, or example, y = y i, ad by usig i, = h we id that: i, = i, Multidimesioal Equatios % two-dimesioal usteady diusio by the FTCS scheme % =3;m=3;step=;D=.5;legth=.;h=legth/(-); dt=.*.5*h*h/d;=zeros(,m);o=zeros(,m);time=.; u=-.; v=-.; (:,)=.; or l=:step,l,time hold o;mesh(); axis([ m.5]);pause; o=; or i=:-, or =:m- (i,)=o(i,)-(.5*dt*u/h)*(o(i+,)-o(i-,))-... (.5*dt*v/h)*(o(i,+)-o(i,-))+... (D*dt/h^)*(o(i+,)+ o(i,+)+o(i-,)+o(i, -)-4*o(i,)); ed,ed or i=:, (i,)=(i,);ed;or =:m,(,)=(,);(m,)=(m-,);ed; time=time+dt; ed; The usteady evolutio o the solutio Multidimesioal Equatios Multidimesioal Boudary Value Problems (Steady-State) Boudary Value Problems Cosider the Poisso Equatio: + = S x y This equatio has a solutio i or is speciied o the boudary Use stadard iite diereces to discretize: i +, + i, i, + + i, i, + i, = S h h i, For uiorm grids: i+, + i, i, Solve or i, : Boudary Value Problems h ca be writte as + i, + i, + i, h i+, + i, + i, + i, + 4 i, = S h = S i, i, = ( 4 i +, + i, + i, + i, + h S i, )
11 Boudary Value Problems Boudary Value Problems Solve or i, ad use the right had side to compute a ew value. Deote the old values by α ad the ew oes with α+ i, ( ) α + = 4 α i +, + α i, + α i, + α i, + h S i, This iteratio process Jacobi iteratio is very robust but may iteratios are required to reach a accurate solutio. The iteratio must be carried out util the solutio is suicietly accurate. To measure the error, deie the residual: R i, = i+, + i, + i, + i, + 4 i, h S i, At steady-state the residual should be zero. The poitwise residual or the average absolute residual ca be used, depedig o the problem. Ote, simpler criteria, such as the chage rom oe iteratio to the ext is used Boudary Value Problems Although the Jacobi iteratio is a very robust iteratio techique, it coverges VERY slowly. We thereore seek a way to ACCELERATE the covergece to steady-state, makig use o the act that it is oly the steady-state that is o iterest. Here we itroduce the Gauss-Seidler method ad the Successive Over-Relaxatio (SOR) method. Boudary Value Problems The Jacobi iteratio ca be improved somewhat by usig ew values as soo as they become available. + - i- i i+ α + i, = 4 ( α i +, + i, α + α + i, + or =:m or i=: iterate ed ed α + + i, h S i, ) From a programmig poit o view, Gauss-Seidler iteratio is eve simpler tha Jacobi iteratio sice oly oe vector with values is eeded. Boudary Value Problems The Gauss-Seidler iteratio ca be accelerated eve urther by various acceleratio techiques. The simplest oe is the Successive Over-Relaxatio (SOR) iteratio α+ i, = β 4 ( α i+, + α+ α i, + i, + + α+ i, h S i, ) α + ( β) i, Example The SOR iteratio is very simple to program, ust as the Gauss-Seidler iteratio. The user must select the coeiciet. It must be bouded by <β<. β=.5 is usually a good startig value.
12 Boudary Value Problems Boudary Value Problems x + y = = = % two-dimesioal steady-state problem by SOR =4;m=4;iteratios=5;legth=.;h=legth/(-); T=zeros(,m);bb=.7; T(:-,)=.; or l=:iteratios, or i=:-, or =:m- T(i,)=bb*.5*(T(i+,)+... T(i,+)+T(i-,)+T(i,-))+(.-bb)*T(i,); ed,ed % id residual res=; or i=:-, or =:m- res=res+abs(t(i+,)+... T(i,+)+T(i-,)+T(i,-)-4*T(i,))/h^; ed,ed l,res/((m-)*(-)) % Prit iteratio ad residual i (res/((m-)*(-)) <.), break,ed ed; cotour(t); Boudary Value Problems The program is easily modiied or the Jacobi ad the Gauss-Seidler iteratio: Average absolute error:. Number o iteratios Jacobi: 989 Gauss-Seidler: 986 SOR (β =.5): 3 SOR (β =.7): 6 SOR (β =.9): 9 SOR (β =.95): The coverged solutio: 3 Boudary Value Problems 3 4 x + y = Summary smaller time step ad ier resolutio should get us the exact solutio itroduced the vo Neuma method. Fairly mechaical process, we will provide more isight by the iite volume poit o view The modiied equatios helps us see how the approximate ad the exact equatio dier ad i the ormer is cosistet with the latter Multidimesioal advectio-diusio equatio. Essetially the same as the oe-dimesioal problem Iterative methods or boudary value problems. Elemetary approaches to steady state problems
More 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 informationThis chapter describes different methods to discretize the diffusion equation. f z 2 = 0. y ) x f
Chapter 8 Diusio Equatio This chapter describes dieret methods to discretize the diusio equatio 2 t α x 2 + 2 y 2 + 2 z 2 = 0 which represets a combied boudary ad iitial value problem, i.e., requires to
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 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 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 informationSolutions Numerical Simulation - Homework, 4/15/ The DuFort-Frankel method is given by. f n. x 2. (a) Truncation error: Taylor expansion
Solutios Numerical Simulatio - Homework, 4/15/211 36. The DuFort-Frakel method is give by +1 1 = αl = α ( 2 t 2 1 1 +1 + +1 (a Trucatio error: Taylor epasio i +1 t = m [ m ] m= m! t m. (1 i Sice we cosider
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 informationNumerical Methods! for! Elliptic Equations-II! Multigrid! Methods! f Analytic Solution of! g = bk. da k dt. = α. ( t)
http://www.d.edu/~gtryggva/cfd-course/! umerical Methods! or! Elliptic Equatios-II! Grétar Tryggvaso! Sprig 13! Examples o elliptic equatios! Direct Methods or 1D problems! Elemetary Iterative Methods!
More informationNumerical Methods! for! Elliptic Equations-II! Multigrid! Methods! f Analytic Solution of! g = bk. da k dt. = α. ( t)
http://www.d.edu/~gtryggva/cfd-course/! Numerical Methods! or! Elliptic Equatios-II! Grétar Tryggvaso! Sprig! Examples o elliptic equatios! Direct Methods or D problems! Elemetary Iterative Methods! Note
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 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 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 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 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 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 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 informationLecture 11. Solution of Nonlinear Equations - III
Eiciecy o a ethod Lecture Solutio o Noliear Equatios - III The eiciecy ide o a iterative ethod is deied by / E r r: rate o covergece o the ethod : total uber o uctios ad derivative evaluatios at each step
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 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 informationTaylor Polynomials and Approximations - Classwork
Taylor Polyomials ad Approimatios - Classwork Suppose you were asked to id si 37 o. You have o calculator other tha oe that ca do simple additio, subtractio, multiplicatio, or divisio. Fareched\ Not really.
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 informationNonlinear regression
oliear regressio How to aalyse data? How to aalyse data? Plot! How to aalyse data? Plot! Huma brai is oe the most powerfull computatioall tools Works differetly tha a computer What if data have o liear
More informationOn Strongly Consistent Finite Dierence Approximations
D.Michels et al. (KAUST,JINR,SSU) Strogly Cosistet Approximatios 17 April 2018 1 / 25 O Strogly Cosistet Fiite Dierece Approximatios Domiik Michels 1, Vladimir Gerdt 2, Dmitry Lyakhov 1, ad Yuri Blikov
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 informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationNumerical Methods for Elliptic Equations-I. Computational Fluid Dynamics. φ f = S. x + 2 f. x a f x + b f x + cf = s. f y b f.
http://www.d.edu/~gtryggva/cfd-course/ http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture 13 March 6, 17 umerical Methods or Elliptic Equatios-I Grétar Tryggvaso Examples o elliptic
More information5. Fast NLMS-OCF Algorithm
5. Fast LMS-OCF Algorithm The LMS-OCF algorithm preseted i Chapter, which relies o Gram-Schmidt orthogoalizatio, has a compleity O ( M ). The square-law depedece o computatioal requiremets o the umber
More informationx 2 x x x x x + x x +2 x
Math 5440: Notes o particle radom walk Aaro Fogelso September 6, 005 Derivatio of the diusio equatio: Imagie that there is a distributio of particles spread alog the x-axis ad that the particles udergo
More informationSome Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations
Copyright, Darbose Iteratioal Joural o Applied Mathematics ad Computatio Volume (), pp -6, 9 http//: ijamc.darbose.com Some Variats o Newto's Method with Fith-Order ad Fourth-Order Covergece or Solvig
More informationPhys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12
Physics Departmet, Yarmouk Uiversity, Irbid Jorda Phys. Mathematical Physics Dr. Nidal M. Ershaidat Doc. Fourier Series Deiitio A Fourier series is a expasio o a periodic uctio (x) i terms o a iiite sum
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 informationWhere do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?
Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)
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 informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationFIR Filter Design: Part I
EEL3: Discrete-Time Sigals ad Systems FIR Filter Desig: Part I. Itroductio FIR Filter Desig: Part I I this set o otes, we cotiue our exploratio o the requecy respose o FIR ilters. First, we cosider some
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 informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
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 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 informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
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 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 informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationChapter 2 The Solution of Numerical Algebraic and Transcendental Equations
Chapter The Solutio of Numerical Algebraic ad Trascedetal Equatios Itroductio I this chapter we shall discuss some umerical methods for solvig algebraic ad trascedetal equatios. The equatio f( is said
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 informationENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 4 Solutions [Numerical Methods]
ENGI 3 Advaced Calculus or Egieerig Facult o Egieerig ad Applied Sciece Problem Set Solutios [Numerical Methods]. Use Simpso s rule with our itervals to estimate I si d a, b, h a si si.889 si 3 si.889
More informationThe natural exponential function
The atural expoetial fuctio Attila Máté Brookly College of the City Uiversity of New York December, 205 Cotets The atural expoetial fuctio for real x. Beroulli s iequality.....................................2
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 informationME NUMERICAL METHODS Fall 2007
ME - 310 NUMERICAL METHODS Fall 2007 Group 02 Istructor: Prof. Dr. Eres Söylemez (Rm C205, email:eres@metu.edu.tr ) Class Hours ad Room: Moday 13:40-15:30 Rm: B101 Wedesday 12:40-13:30 Rm: B103 Course
More informationSequences. Notation. Convergence of a Sequence
Sequeces A sequece is essetially just a list. Defiitio (Sequece of Real Numbers). A sequece of real umbers is a fuctio Z (, ) R for some real umber. Do t let the descriptio of the domai cofuse you; it
More information10-701/ Machine Learning Mid-term Exam Solution
0-70/5-78 Machie Learig Mid-term Exam Solutio Your Name: Your Adrew ID: True or False (Give oe setece explaatio) (20%). (F) For a cotiuous radom variable x ad its probability distributio fuctio p(x), it
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 informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
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 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 informationIterative Techniques for Solving Ax b -(3.8). Assume that the system has a unique solution. Let x be the solution. Then x A 1 b.
Iterative Techiques for Solvig Ax b -(8) Cosider solvig liear systems of them form: Ax b where A a ij, x x i, b b i Assume that the system has a uique solutio Let x be the solutio The x A b Jacobi ad Gauss-Seidel
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
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 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 informationPractice Problems: Taylor and Maclaurin Series
Practice Problems: Taylor ad Maclauri Series Aswers. a) Start by takig derivatives util a patter develops that lets you to write a geeral formula for the -th derivative. Do t simplify as you go, because
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 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 informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationExplicit Group Methods in the Solution of the 2-D Convection-Diffusion Equations
Proceedigs of the World Cogress o Egieerig 00 Vol III WCE 00 Jue 0 - July 00 Lodo U.K. Explicit Group Methods i the Solutio of the -D Covectio-Diffusio Equatios a Kah Bee orhashidah Hj. M. Ali ad Choi-Hog
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 informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationNotes for Lecture 11
U.C. Berkeley CS78: Computatioal Complexity Hadout N Professor Luca Trevisa 3/4/008 Notes for Lecture Eigevalues, Expasio, ad Radom Walks As usual by ow, let G = (V, E) be a udirected d-regular graph with
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
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 informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationTR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT
TR/46 OCTOBER 974 THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION by A. TALBOT .. Itroductio. A problem i approximatio theory o which I have recetly worked [] required for its solutio a proof that the
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
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 informationj=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)
Problem. Compute the itegrals C r d for Z, where C r = ad r >. Recall that C r has the couter-clockwise orietatio. Solutio: We will use the idue Theorem to solve this oe. We could istead use other (perhaps
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 informationPhysics 116A Solutions to Homework Set #1 Winter Boas, problem Use equation 1.8 to find a fraction describing
Physics 6A Solutios to Homework Set # Witer 0. Boas, problem. 8 Use equatio.8 to fid a fractio describig 0.694444444... Start with the formula S = a, ad otice that we ca remove ay umber of r fiite decimals
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
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 informationCaputo s Implicit Solution of Time-Fractional Diffusion Equation Using Half-Sweep AOR Iteration
Global Joural o Pure ad Applied Mathematics. ISSN 0973-768 Volume Number 4 (06) pp. 3469-3479 Research Idia Publicatios http://www.ripublicatio.com/gpam.htm Caputo s Implicit Solutio o Time-Fractioal Diusio
More informationENGI 9420 Lecture Notes 3 - Numerical Methods Page 3.01
ENGI 940 Lecture Notes 3 - Numerical Methods Page 3.0 3. Numerical Methods The majority of equatios of iterest i actual practice do ot admit ay aalytic solutio. Eve equatios as simple as = e ad I = e d
More informationNumerical Method for Blasius Equation on an infinite Interval
Numerical Method for Blasius Equatio o a ifiite Iterval Alexader I. Zadori Omsk departmet of Sobolev Mathematics Istitute of Siberia Brach of Russia Academy of Scieces, Russia zadori@iitam.omsk.et.ru 1
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 informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationECON 3150/4150, Spring term Lecture 3
Itroductio Fidig the best fit by regressio Residuals ad R-sq Regressio ad causality Summary ad ext step ECON 3150/4150, Sprig term 2014. Lecture 3 Ragar Nymoe Uiversity of Oslo 21 Jauary 2014 1 / 30 Itroductio
More informationTopic 9 - Taylor and MacLaurin Series
Topic 9 - Taylor ad MacLauri Series A. Taylors Theorem. The use o power series is very commo i uctioal aalysis i act may useul ad commoly used uctios ca be writte as a power series ad this remarkable result
More informationFourier Series and the Wave Equation
Fourier Series ad the Wave Equatio We start with the oe-dimesioal wave equatio u u =, x u(, t) = u(, t) =, ux (,) = f( x), u ( x,) = This represets a vibratig strig, where u is the displacemet of the strig
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 information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More information-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION
NEW NEWTON-TYPE METHOD WITH k -ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds West Yorkshire, LS7 JS, ENGLAND ABSTRACT The objective
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 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 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 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 informationIN many scientific and engineering applications, one often
INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS, VOL 3, NO, FEBRUARY 07 5 Secod Degree Refiemet Jacobi Iteratio Method for Solvig System of Liear Equatio Tesfaye Kebede Abstract Several
More information9.3 Taylor s Theorem: Error Analysis for Series. Tacoma Narrows Bridge: November 7, 1940
9. Taylor s Theorem: Error Aalysis or Series Tacoma Narrows Bridge: November 7, 940 Last time i BC So the Taylor Series or l x cetered at x is give by ) l x ( ) ) + ) ) + ) ) 4 Use the irst two terms o
More information