A large class of initial value (time-evolution) PDEs in one space dimension can be cast into the form of a flux-conservative equation, u (19.1.
|
|
- Pearl Stevens
- 5 years ago
- Views:
Transcription
1 834 Chapter 19. Partial Differetial Equatios egieerig; these methods allow cosiderable freedom i puttig computatioal elemets where you wat them, importat whe dealig with highly irregular geometries. Spectral methods [13-15] are preferred for very regular geometries ad smooth fuctios; they coverge more rapidly tha fiite-differece methods (cf. 19.4), but they do ot work well for problems with discotiuities. CITED REFERENCES AND FURTHER READING: Ames, W.F. 1977, Numerical Methods for Partial Differetial Equatios, 2d ed. (New York: Academic Press). [1] Richtmyer, R.D., ad Morto, K.W. 1967, Differece Methods for Iitial Value Problems, 2d ed. (New York: Wiley-Itersciece). [2] Roache, P.J. 1976, Computatioal Fluid Dyamics (Albuquerque: Hermosa). [3] Mitchell, A.R., ad Griffiths, D.F. 1980, The Fiite Differece Method i Partial Differetial Equatios (New York: Wiley) [icludes discussio of fiite elemet methods]. [4] Dorr, F.W. 1970, SIAM Review, vol. 12, pp [5] Meierik, J.A., ad va der Vorst, H.A. 1977, Mathematics of Computatio, vol. 31, pp [6] va der Vorst, H.A. 1981, Joural of Computatioal Physics, vol. 44, pp [review of sparse iterative methods]. [7] Kershaw, D.S. 1970, Joural of Computatioal Physics, vol. 26, pp [8] Stoe, H.J. 1968, SIAM Joural o Numerical Aalysis, vol. 5, pp [9] Jesshope, C.R. 1979, Computer Physics Commuicatios, vol. 17, pp [10] Strag, G., ad Fix, G. 1973, A Aalysis of the Fiite Elemet Method (Eglewood Cliffs, NJ: Pretice-Hall). [11] Burett, D.S. 1987, Fiite Elemet Aalysis: From Cocepts to Applicatios (Readig, MA: Addiso-Wesley). [12] Gottlieb, D. ad Orszag, S.A. 1977, Numerical Aalysis of Spectral Methods: Theory ad Applicatios (Philadelphia: S.I.A.M.). [13] Cauto, C., Hussaii, M.Y., Quarteroi, A., ad Zag, T.A. 1988, Spectral Methods i Fluid Dyamics (New York: Spriger-Verlag). [14] Boyd, J.P. 1989, Chebyshev ad Fourier Spectral Methods (New York: Spriger-Verlag). [15] 19.1 Flux-Coservative Iitial Value Problems A large class of iitial value (time-evolutio) PDEs i oe space dimesio ca be cast ito the form of a flux-coservative equatio, u t = F(u) (19.1.1) x where u ad F are vectors, ad where (i some cases) F may deped ot oly o u but also o spatial derivatives of u. The vector F is called the coserved flux. For example, the prototypical hyperbolic equatio, the oe-dimesioal wave equatio with costat velocity of propagatio v 2 u t 2 = v2 2 u x 2 (19.1.2)
2 19.1 Flux-Coservative Iitial Value Problems 835 ca be rewritte as a set of two first-order equatios where r t = v s x s t = v r x (19.1.3) r v u x s u (19.1.4) t I this case r ad s become the two compoets of u, ad the flux is give by the liear matrix relatio F(u) = ( ) 0 v u (19.1.5) v 0 (The physicist-reader may recogize equatios (19.1.3) as aalogous to Maxwell s equatios for oe-dimesioal propagatio of electromagetic waves.) We will cosider, i this sectio, a prototypical example of the geeral fluxcoservative equatio (19.1.1), amely the equatio for a scalar u, u t = v u (19.1.6) x with v a costat. As it happes, we already kow aalytically that the geeral solutio of this equatio is a wave propagatig i the positive x-directio, u = f(x vt) (19.1.7) where f is a arbitrary fuctio. However, the umerical strategies that we develop will be equally applicable to the more geeral equatios represeted by (19.1.1). I some cotexts, equatio (19.1.6) iscalled a advective equatio, because the quatity u is trasported by a fluid flow with a velocity v. How do we go about fiite differecig equatio (19.1.6) (or, aalogously, )? The straightforward approach is to choose equally spaced poits alog both the t- adx-axes. Thus deote x = x 0 +, t =t 0 +, =0,1,...,J =0,1,...,N (19.1.8) Let u deote u(t,x ). We have several choices for represetig the time derivative term. The obvious way is to set u t = u+1 u + O() (19.1.9), This is called forward Euler differecig (cf. equatio ). While forward Euler is oly first-order accurate i, it has the advatage that oe is able to calculate
3 836 Chapter 19. Partial Differetial Equatios t or FTCS x or Figure Represetatio of the Forward Time Cetered Space (FTCS) differecig scheme. I this ad subsequet figures, the ope circle is the ew poit at which the solutio is desired; filled circles are kow poits whose fuctio values are used i calculatig the ew poit; the solid lies coect poits that are used to calculate spatial derivatives; the dashed lies coect poits that are used to calculate time derivatives. The FTCS scheme is geerally ustable for hyperbolic problems ad caot usually be used. quatities at timestep +1i terms of oly quatities kow at timestep. Forthe space derivative, we ca use a secod-order represetatio still usig oly quatities kow at timestep : u x = u +1 u 1 + O( 2 ) ( ), 2 The resultig fiite-differece approximatio to equatio (19.1.6) is called the FTCS represetatio (Forward Time Cetered Space), u +1 u ( u +1 u ) 1 = v 2 ( ) which ca easily be rearraged to be a formula for u +1 i terms of the other quatities. The FTCS scheme is illustrated i Figure It s a fie example of a algorithm that is easy to derive, takes little storage, ad executes quickly. Too bad it does t work! (See below.) The FTCS represetatio is a explicit scheme. This meas that u +1 for each ca be calculated explicitly from the quatities that are already kow. Later we shall meet implicit schemes, which require us to solve implicit equatios couplig the u +1 for various. (Explicit ad implicit methods for ordiary differetial equatios were discussed i 16.6.) The FTCS algorithm is also a example of a sigle-level scheme, sice oly values at time level have to be stored to fid values at time level +1. vo Neuma Stability Aalysis Ufortuately, equatio ( ) is of very limited usefuless. It is a ustable method, which ca be used oly (if at all) to study waves for a short fractio of oe oscillatio period. To fid alterative methods with more geeral applicability, we must itroduce the vo Neuma stability aalysis. The vo Neuma aalysis is local: We imagie that the coefficiets of the differece equatios are so slowly varyig as to be cosidered costat i space ad time. I that case, the idepedet solutios, or eigemodes, of the differece equatios are all of the form u = ξ e ik ( )
4 19.1 Flux-Coservative Iitial Value Problems 837 t or Lax x or Figure Represetatio of the Lax differecig scheme, as i the previous figure. The stability criterio for this scheme is the Courat coditio. where k is a real spatial wave umber (which ca have ay value) ad ξ = ξ(k) is a complex umber that depeds o k. The key fact is that the time depedece of a sigle eigemode is othig more tha successive iteger powers of the complex umber ξ. Therefore, the differece equatios are ustable (have expoetially growig modes) if ξ(k) > 1 for some k. The umber ξ is called the amplificatio factor at a give wave umber k. To fid ξ(k), we simply substitute ( ) back ito ( ). Dividig by ξ, we get ξ(k)=1 i v si k ( ) whose modulus is > 1 for all k; so the FTCS scheme is ucoditioally ustable. If the velocity v were a fuctio of t ad x, the we would write v i equatio ( ). I the vo Neuma stability aalysis we would still treat v as a costat, theideabeigthatforvslowly varyig the aalysis is local. I fact, eve i the case of strictly costat v, the vo Neuma aalysis does ot rigorously treat the ed effects at =0ad = N. More geerally, if the equatio s right-had side were oliear i u, thea vo Neuma aalysis would liearize by writig u = u 0 + δu, expadig to liear order i δu. Assumig that the u 0 quatities already satisfy the differece equatio exactly, the aalysis would look for a ustable eigemode of δu. Despite its lack of rigor, the vo Neuma method geerally gives valid aswers ad is much easier to apply tha more careful methods. We accordigly adopt it exclusively. (See, for example, [1] for a discussio of other methods of stability aalysis.) Lax Method The istability i the FTCS method ca be cured by a simple chage due to Lax. Oe replaces the term u i the time derivative term by its average (Figure ): u 1 2 ( u +1 + u 1) ( ) This turs ( ) ito u +1 = 1 ( u u v ( 1) u 2 +1 u 1) ( )
5 838 Chapter 19. Partial Differetial Equatios stable ustable t or (a) x or Figure Courat coditio for stability of a differecig scheme. The solutio of a hyperbolic problem at a poit depeds o iformatio withi some domai of depedecy to the past, show here shaded. The differecig scheme ( ) has its ow domai of depedecy determied by the choice of poits o oe time slice (show as coected solid dots) whose values are used i determiig a ew poit (show coected by dashed lies). A differecig scheme is Courat stable if the differecig domai of depedecy is larger tha that of the PDEs, as i (a), ad ustable if the relatioship is the reverse, as i (b). For more complicated differecig schemes, the domai of depedecy might ot be determied simply by the outermost poits. Substitutig equatio ( ), we fid for the amplificatio factor ξ =cosk i v si k ( ) The stability coditio ξ 2 1 leads to the requiremet v 1 ( ) This is the famous Courat-Friedrichs-Lewy stability criterio, ofte called simply the Courat coditio. Ituitively, the stability coditio ca be uderstood as follows (Figure ): The quatity u +1 i equatio ( ) is computed from iformatio at poits 1 ad +1at time. I other words, x 1 ad x +1 are the boudaries of the spatial regio that is allowed to commuicate iformatio to u +1. Now recall that i the cotiuum wave equatio, iformatio actually propagates with a maximum velocity v. If the poit u +1 is outside of the shaded regio i Figure , the it requires iformatio from poits more distat tha the differecig scheme allows. Lack of that iformatio gives rise to a istability. Therefore, caot be made too large. The surprisig result, that the simple replacemet ( ) stabilizes the FTCS scheme, is our first ecouter with the fact that differecig PDEs is a art as much as a sciece. To see if we ca demystify the art somewhat, let us compare the FTCS ad Lax schemes by rewritig equatio ( ) so that it is i the form of equatio ( ) with a remaider term: u +1 u ( u +1 u ) 1 = v + 1 ( u +1 2u + ) u 1 ( ) 2 2 But this is exactly the FTCS represetatio of the equatio (b) u t = v u x + ()2 2 2 u ( )
6 19.1 Flux-Coservative Iitial Value Problems 839 where 2 = 2 / x 2 i oe dimesio. We have, i effect, added a diffusio term to the equatio, or, if you recall the form of the Navier-Stokes equatio for viscous fluid flow, a dissipative term. The Lax scheme is thus said to have umerical dissipatio, or umerical viscosity. We ca see this also i the amplificatio factor. Uless v is exactly equal to, ξ < 1 ad the amplitude of the wave decreases spuriously. Is t a spurious decrease as bad as a spurious icrease? No. The scales that we hope to studyaccurately are those that ecompass may grid poits, so that they have k 1. (The spatial wave umber k is defied by equatio ) For these scales, the amplificatio factor ca be see to be very close to oe, i both the stable ad ustable schemes. The stable ad ustable schemes are therefore about equally accurate. For the ustable scheme, however, short scales with k 1, which we are ot iterested i, will blow up ad swamp the iterestig part of the solutio. Much better to have a stable scheme i which these short wavelegths die away iocuously. Both the stable ad the ustable schemes are iaccurate for these short wavelegths, but the iaccuracy is of a tolerable character whe the scheme is stable. Whe the idepedet variable u is a vector, the the vo Neuma aalysis is slightly more complicated. For example, we ca cosider equatio (19.1.3), rewritte as The Lax method for this equatio is r +1 s +1 [ ] r = t s x [ ] vs vr = 1 2 (r +1 + v r 1 )+ 2 (s +1 s 1 ) = 1 2 (s +1 + s 1)+ v 2 (r +1 r 1) ( ) ( ) The vo Neuma stability aalysis ow proceeds by assumig that the eigemode is of the followig (vector) form, [ ] [ ] r = ξ e ik r 0 ( ) s Here the vector o the right-had side is a costat (both i space ad i time) eigevector, ad ξ is a complex umber, as before. Substitutig ( ) ito ( ), ad dividig by the power ξ, gives the homogeeous vector equatio (cos k) ξ iv si k i v r0 = 0 ( ) si k (cos k) ξ s 0 0 This admits a solutio oly if the determiat of the matrix o the left vaishes, a coditio easily show to yield the two roots ξ s 0 ξ =cosk±i v si k ( ) The stability coditio is that both roots satisfy ξ 1. This agai turs out to be simply the Courat coditio ( ).
7 840 Chapter 19. Partial Differetial Equatios Other Varieties of Error Thus far we have bee cocered with amplitude error, because of its itimate coectio with the stability or istability of a differecig scheme. Other varieties of error are relevat whe we shift our cocer to accuracy, rather tha stability. Fiite-differece schemes for hyperbolic equatios ca exhibit dispersio, or phase errors. For example, equatio ( ) ca be rewritte as ( ξ = e ik + i 1 v ) si k ( ) A arbitrary iitial wave packet is a superpositio of modes with differet k s. At each timestep the modes get multiplied by differet phase factors ( ), depedig o their value of k. If =/v, the the exact solutio for each mode of a wave packet f(x vt) is obtaied if each mode gets multiplied by exp( ik). For this value of, equatio ( ) shows that the fiite-differece solutio gives the exact aalytic result. However, if v/ is ot exactly 1, the phase relatios of the modes ca become hopelessly garbled ad the wave packet disperses. Note from ( ) that the dispersio becomes large as soo as the wavelegth becomes comparable to the grid spacig. A third type of error is oe associated with oliear hyperbolic equatios ad is therefore sometimes called oliear istability. For example, a piece of the Euler or Navier-Stokes equatios for fluid flow looks like v t v = v +... ( ) x The oliear term i v ca cause a trasfer of eergy i Fourier space from log wavelegths to short wavelegths. This results i a wave profile steepeig util a vertical profile or shock develops. Sice the vo Neuma aalysis suggests that the stability ca deped o k, a scheme that was stable for shallow profiles ca become ustable for steep profiles. This kid of difficulty arises i a differecig scheme where the cascade i Fourier space is halted at the shortest wavelegth represetable o the grid, that is, at k 1/. If eergy simply accumulates i these modes, it evetually swamps the eergy i the log wavelegth modes of iterest. Noliear istability ad shock formatio is thus somewhat cotrolled by umerical viscosity such as that discussed i coectio with equatio ( ) above. I some fluid problems,however, shock formatio is ot merely a aoyace, but a actual physical behavior of the fluid whose detailed study is a goal. The, umerical viscosity aloe may ot be adequate or sufficietly cotrollable. This is a complicated subect which we discuss further i the subsectio o fluid dyamics, below. For wave equatios, propagatio errors (amplitude or phase) are usually most worrisome. For advective equatios, o the other had, trasport errors are usually of greater cocer. I the Lax scheme, equatio ( ), a disturbace i the advected quatity u at mesh poit propagates to mesh poits +1ad 1 at the ext timestep. I reality, however, if the velocity v is positive the oly mesh poit +1 should be affected.
8 19.1 Flux-Coservative Iitial Value Problems 841 v t or x or upwid Figure Represetatio of upwid differecig schemes. The upper scheme is stable whe the advectio costat v is egative, as show; the lower scheme is stable whe the advectio costat v is positive, also as show. The Courat coditio must, of course, also be satisfied. The simplest way to model the trasport properties better is to use upwid differecig (see Figure ): u +1 u u u 1, v = v > 0 u +1 ( ) u, v < 0 Note that this scheme is oly first-order, ot secod-order, accurate i the calculatio of the spatial derivatives. How ca it be better? The aswer is oe that aoys the mathematicias: The goal of umerical simulatios is ot always accuracy i a strictly mathematical sese, but sometimes fidelity to the uderlyig physics i a sese that is looser ad more pragmatic. I such cotexts, some kids of error are much more tolerable tha others. Upwid differecig geerally adds fidelity to problems where the advected variables are liable to udergo sudde chages of state, e.g., as they pass through shocks or other discotiuities. You will have to be guided by the specific ature of your ow problem. For the differecig scheme ( ), the amplificatio factor (for costat v)is ξ =1 v (1 cos k) iv si k ( ) ( ) ξ 2 =1 2 v 1 v (1 cos k) ( ) v So the stability criterio ξ 2 1 is (agai) simply the Courat coditio ( ). There are various ways of improvig the accuracy of first-order upwid differecig. I the cotiuum equatio, material origially a distace v away
9 842 Chapter 19. Partial Differetial Equatios t or x or staggered leapfrog Figure Represetatio of the staggered leapfrog differecig scheme. Note that iformatio from two previous time slices is used i obtaiig the desired poit. This scheme is secod-order accurate i both space ad time. arrives at a give poit after a time iterval. I the first-order method, the material always arrives from away. If v (to isure accuracy), this ca cause a large error. Oe way of reducig this error is to iterpolate u betwee 1 ad before trasportig it. This gives effectively a secod-order method. Various schemes for secod-order upwid differecig are discussed ad compared i [2-3]. Secod-Order Accuracy i Time Whe usig a method that is first-order accurate i time but secod-order accurate i space, oe geerally has to take v sigificatly smaller tha to achieve desired accuracy, say, by at least a factor of 5. Thus the Courat coditio is ot actually the limitig factor with such schemes i practice. However, there are schemes that are secod-order accurate i both space ad time, ad these ca ofte be pushed right to their stability limit, with correspodigly smaller computatio times. For example, the staggered leapfrog method for the coservatio equatio (19.1.1) is defied as follows (Figure ): Usig the values of u at time t, compute the fluxes F. The compute ew values u+1 usig the time-cetered values of the fluxes: u +1 u 1 = (F +1 F 1 ) ( ) The ame comes from the fact that the time levels i the time derivative term leapfrog over the time levels i the space derivative term. The method requires that u 1 ad u be stored to compute u +1. For our simple model equatio (19.1.6), staggered leapfrog takes the form u +1 u 1 = v (u +1 u 1) ( ) The vo Neuma stability aalysis ow gives a quadratic equatio for ξ, rather tha a liear oe, because of the occurrece of three cosecutive powers of ξ whe the
10 19.1 Flux-Coservative Iitial Value Problems 843 form ( ) for a eigemode is substituted ito equatio ( ), ξ 2 1= 2iξ v si k ( ) whose solutio is ξ = i v si k ± 1 ( ) 2 v si k ( ) Thus the Courat coditio is agai required for stability. I fact, i equatio ( ), ξ 2 =1for ay v. This is the great advatage of the staggered leapfrog method: There is o amplitude dissipatio. Staggered leapfrog differecig of equatios like ( ) is most trasparet if the variables are cetered o appropriate half-mesh poits: r+1/2 v u x s +1/2 u t +1/2 +1/2 = v u +1 u = u+1 u ( ) This is purely a otatioal coveiece: we ca thik of the mesh o which r ad s are defied as beig twice as fie as the mesh o which the origial variable u is defied. The leapfrog differecig of equatio ( ) is r +1 +1/2 r +1/2 s +1/2 s 1/2 = s+1/2 +1 s +1/2 = v r +1/2 r 1/2 ( ) If you substitute equatio ( ) i equatio ( ), you will fid that oce agai the Courat coditio is required for stability, ad that there is o amplitude dissipatio whe it is satisfied. If we substitute equatio ( ) i equatio ( ), we fid that equatio ( ) is equivalet to u +1 2u + u 1 () 2 = v 2 u +1 2u + u 1 () 2 ( ) This is ust the usual secod-order differecig of the wave equatio (19.1.2). We see that it is a two-level scheme, requirig both u ad u 1 to obtai u +1. I equatio ( ) this shows up as both s 1/2 ad r beig eeded to advace the solutio. For equatios more complicated tha our simple model equatio, especially oliear equatios, the leapfrog method usually becomes ustable whe the gradiets get large. The istability is related to the fact that odd ad eve mesh poits are completely decoupled, like the black ad white squares of a chess board, as show
11 844 Chapter 19. Partial Differetial Equatios Figure Origi of mesh-drift istabilities i a staggered leapfrog scheme. If the mesh poits are imagied to lie i the squares of a chess board, the white squares couple to themselves, black to themselves, but there is o couplig betwee white ad black. The fix is to itroduce a small diffusive mesh-couplig piece. i Figure This mesh driftig istability is cured by couplig the two meshes through a umerical viscosity term, e.g., addig to the right side of ( ) a small coefficiet ( 1) times u +1 2u + u 1. For more o stabilizig differece schemes by addig umerical dissipatio, see, e.g., [4]. The Two-Step Lax-Wedroff scheme is a secod-order i time method that avoids large umerical dissipatio ad mesh driftig. Oe defies itermediate values u +1/2 at the half timesteps t +1/2 ad the half mesh poits x +1/2.These are calculated by the Lax scheme: u +1/2 +1/2 = 1 2 (u +1 + u ) 2 (F +1 F ) ( ) Usig these variables, oe calculates the fluxes F +1/2 +1/2. The the updated values are calculated by the properly cetered expressio u +1 u +1 = u ( ) F +1/2 +1/2 F +1/2 1/2 ( ) The provisioal values u +1/2 +1/2 are ow discarded. (See Figure ) Let us ivestigate the stability of this method for our model advective equatio, where F = vu. Substitute ( ) i ( ) to get u +1 [ 1 = u α 2 (u +1 + u ) 1 2 α(u +1 u ) 1 ] ( ) 2 (u + u 1 )+1 2 α(u u 1 )
12 19.1 Flux-Coservative Iitial Value Problems 845 two-step Lax Wedroff halfstep poits t or x or Figure Represetatio of the two-step Lax-Wedroff differecig scheme. Two halfstep poits ( ) are calculated by the Lax method. These, plus oe of the origial poits, produce the ew poit via staggered leapfrog. Halfstep poits are used oly temporarily ad do ot require storage allocatio o the grid. This scheme is secod-order accurate i both space ad time. where The so α v ( ) ξ =1 iα si k α 2 (1 cos k) ( ) ξ 2 =1 α 2 (1 α 2 )(1 cos k) 2 ( ) The stability criterio ξ 2 1 is therefore α 2 1, or v as usual. Icidetally, you should ot thik that the Courat coditio is the oly stability requiremet that ever turs up i PDEs. It keeps doig so i our model examples ust because those examples are so simple i form. The method of aalysis is, however, geeral. Except whe α =1, ξ 2 <1i ( ), so some amplitude dampig does occur. The effect is relatively small, however, for wavelegths large compared with the mesh size. If we expad ( ) for small k, wefid ξ 2 =1 α 2 (1 α 2 ) (k) ( ) The departure from uity occurs oly at fourth order i k. This should be cotrasted with equatio ( ) for the Lax method, which shows that ξ 2 =1 (1 α 2 )(k) ( ) for small k.
13 846 Chapter 19. Partial Differetial Equatios I summary, our recommedatio for iitial value problems that ca be cast i flux-coservative form, ad especially problems related to the wave equatio, is to use the staggered leapfrog method whe possible. We have persoally had better success with it tha with the Two-Step Lax-Wedroff method. For problems sesitive to trasport errors, upwid differecig or oe of its refiemets should be cosidered. Fluid Dyamics with Shocks As we alluded to earlier, the treatmet of fluid dyamics problems with shocks has become a very complicated ad very sophisticated subect. All we ca attempt to do here is to guide you to some startig poits i the literature. There are basically three importat geeral methods for hadlig shocks. The oldest ad simplest method, iveted by vo Neuma ad Richtmyer, is to add artificial viscosity to the equatios, modelig the way Nature uses real viscosity to smooth discotiuities. A good startig poit for tryig out this method is the differecig scheme i 12.11of [1]. This scheme is excellet for early all problems i oe spatial dimesio. The secod method combies a high-order differecig scheme that is accurate for smooth flows with a low order scheme that is very dissipative ad ca smooth the shocks. Typically, various upwid differecig schemes are combied usig weights chose to zero the low order scheme uless steep gradiets are preset, ad also chose to eforce various mootoicity costraits that prevet ophysical oscillatios from appearig i the umerical solutio. Refereces [2-3,5] are a good place to start with these methods. The third, ad potetially most powerful method, is Goduov s approach. Here oe gives up the simple liearizatio iheret i fiite differecig based o Taylor series ad icludes the oliearity of the equatios explicitly. There is a aalytic solutio for the evolutio of two uiform states of a fluid separated by a discotiuity, the Riema shock problem. Goduov s idea was to approximate the fluid by a large umber of cells of uiform states, ad piece them together usig the Riema solutio. There have bee may geeralizatios of Goduov s approach, of which the most powerful is probably the PPM method [6]. Readable reviews of all these methods, discussig the difficulties arisig whe oe-dimesioal methods are geeralized to multidimesios, are give i [7-9]. CITED REFERENCES AND FURTHER READING: Ames, W.F. 1977, Numerical Methods for Partial Differetial Equatios, 2d ed. (New York: Academic Press), Chapter 4. Richtmyer, R.D., ad Morto, K.W. 1967, Differece Methods for Iitial Value Problems, 2d ed. (New York: Wiley-Itersciece). [1] Cetrella, J., ad Wilso, J.R. 1984, Astrophysical Joural Supplemet, vol. 54, pp , Appedix B. [2] Hawley, J.F., Smarr, L.L., ad Wilso, J.R. 1984, Astrophysical Joural Supplemet, vol. 55, pp , 2c. [3] Kreiss, H.-O. 1978, Numerical Methods for Solvig Time-Depedet Problems for Partial Differetial Equatios (Motreal: Uiversity of Motreal Press), pp. 66ff. [4] Harte, A., Lax, P.D., ad Va Leer, B. 1983, SIAM Review, vol. 25, pp [5] Woodward, P., ad Colella, P. 1984, Joural of Computatioal Physics, vol. 54, pp [6]
14 19.2 Diffusive Iitial Value Problems 847 Roache, P.J. 1976, Computatioal Fluid Dyamics (Albuquerque: Hermosa). [7] Woodward, P., ad Colella, P. 1984, Joural of Computatioal Physics, vol. 54, pp [8] Rizzi, A., ad Egquist, B. 1987, Joural of Computatioal Physics, vol. 72, pp [9] 19.2 Diffusive Iitial Value Problems Recall the model parabolic equatio, the diffusio equatio i oe space dimesio, u t = ( D u ) (19.2.1) x x where D is the diffusio coefficiet. Actually, this equatio is a flux-coservative equatio of the form cosidered i the previous sectio, with F = D u (19.2.2) x the flux i the x-directio. We will assume D 0, otherwise equatio (19.2.1) has physically ustable solutios: A small disturbace evolves to become more ad more cocetrated istead of dispersig. (Do t make the mistake of tryig to fid a stable differecig scheme for a problem whose uderlyig PDEs are themselves ustable!) Eve though (19.2.1) is of the form already cosidered, it is useful to cosider it as a model i its ow right. The particular form of flux (19.2.2), ad its direct geeralizatios, occur quite frequetly i practice. Moreover, we have already see that umerical viscosity ad artificial viscosity ca itroduce diffusive pieces like the right-had side of (19.2.1) i may other situatios. Cosider first the case whe D is a costat. The the equatio u t = u D 2 x 2 (19.2.3) ca be differeced i the obvious way: u +1 u [ u +1 2u = D + ] u 1 () 2 (19.2.4) This is the FTCS scheme agai, except that it is a secod derivative that has bee differeced o the right-had side. But this makes a world of differece! The FTCS scheme was ustable for the hyperbolic equatio; however, a quick calculatio shows that the amplificatio factor for equatio (19.2.4) is ξ =1 4D () 2 si2 ( ) k 2 (19.2.5) The requiremet ξ 1leads to the stability criterio 2D () 2 1 (19.2.6)
Partial 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 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 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 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 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 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 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 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 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 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 informationHydro code II. Simone Recchi
Hydro code II Simoe Recchi 16.06.2011 Literature E.F. Toro: Riema Solvers ad Numerical Methods for Fluid Dyamics, Spriger R.J. LeVeque: Fiite Volume Methods for Hyperbolic Problems, Cambridge Uiversity
More 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 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 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 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 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 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 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 informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
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 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 informationRoberto s Notes on Series Chapter 2: Convergence tests Section 7. Alternating series
Roberto s Notes o Series Chapter 2: Covergece tests Sectio 7 Alteratig series What you eed to kow already: All basic covergece tests for evetually positive series. What you ca lear here: A test for series
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 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 informationNumerical Methods for PDEs
Numerical Methods for PDEs Hyperbolic PDEs: Coupled system/noliear coservatio laws/a oliear Lax-Wedroff scheme (Lecture 18, Week 6 Markus Schmuck Departmet of Mathematics ad Maxwell Istitute for Mathematical
More 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 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 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 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 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 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 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 informationREGRESSION (Physics 1210 Notes, Partial Modified Appendix A)
REGRESSION (Physics 0 Notes, Partial Modified Appedix A) HOW TO PERFORM A LINEAR REGRESSION Cosider the followig data poits ad their graph (Table I ad Figure ): X Y 0 3 5 3 7 4 9 5 Table : Example Data
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
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 informationSeunghee Ye Ma 8: Week 5 Oct 28
Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value
More informationFall 2013 MTH431/531 Real analysis Section Notes
Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters
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 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 information4.1 Sigma Notation and Riemann Sums
0 the itegral. Sigma Notatio ad Riema Sums Oe strategy for calculatig the area of a regio is to cut the regio ito simple shapes, calculate the area of each simple shape, ad the add these smaller areas
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 information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More informationmx bx kx F t. dt IR I LI V t, Q LQ RQ V t,
Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio
More 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 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 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 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 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 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 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 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 information( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!
.8,.9: Taylor ad Maclauri Series.8. Although we were able to fid power series represetatios for a limited group of fuctios i the previous sectio, it is ot immediately obvious whether ay give fuctio has
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 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 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 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 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 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 informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationLecture 6: Integration and the Mean Value Theorem. slope =
Math 8 Istructor: Padraic Bartlett Lecture 6: Itegratio ad the Mea Value Theorem Week 6 Caltech 202 The Mea Value Theorem The Mea Value Theorem abbreviated MVT is the followig result: Theorem. Suppose
More informationSECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES
SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES Read Sectio 1.5 (pages 5 9) Overview I Sectio 1.5 we lear to work with summatio otatio ad formulas. We will also itroduce a brief overview of sequeces,
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 informationThe Phi Power Series
The Phi Power Series I did this work i about 0 years while poderig the relatioship betwee the golde mea ad the Madelbrot set. I have fially decided to make it available from my blog at http://semresearch.wordpress.com/.
More informationNumber of fatalities X Sunday 4 Monday 6 Tuesday 2 Wednesday 0 Thursday 3 Friday 5 Saturday 8 Total 28. Day
LECTURE # 8 Mea Deviatio, Stadard Deviatio ad Variace & Coefficiet of variatio Mea Deviatio Stadard Deviatio ad Variace Coefficiet of variatio First, we will discuss it for the case of raw data, ad the
More informationBertrand s Postulate
Bertrad s Postulate Lola Thompso Ross Program July 3, 2009 Lola Thompso (Ross Program Bertrad s Postulate July 3, 2009 1 / 33 Bertrad s Postulate I ve said it oce ad I ll say it agai: There s always a
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 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 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 informationOutput Analysis and Run-Length Control
IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α%
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 informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationAn efficient time integration method for extra-large eddy simulations
A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler Departmet of Mathematics Master s Thesis A efficiet time itegratio method for extra-large eddy simulatios M.A. Scheibeler
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator  is defied by its actio
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Statistics
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER 1 018/019 DR. ANTHONY BROWN 8. Statistics 8.1. Measures of Cetre: Mea, Media ad Mode. If we have a series of umbers the
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 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 informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
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 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 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 information(3) If you replace row i of A by its sum with a multiple of another row, then the determinant is unchanged! Expand across the i th row:
Math 5-4 Tue Feb 4 Cotiue with sectio 36 Determiats The effective way to compute determiats for larger-sized matrices without lots of zeroes is to ot use the defiitio, but rather to use the followig facts,
More informationSequences of Definite Integrals, Factorials and Double Factorials
47 6 Joural of Iteger Sequeces, Vol. 8 (5), Article 5.4.6 Sequeces of Defiite Itegrals, Factorials ad Double Factorials Thierry Daa-Picard Departmet of Applied Mathematics Jerusalem College of Techology
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
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 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 information7 Sequences of real numbers
40 7 Sequeces of real umbers 7. Defiitios ad examples Defiitio 7... A sequece of real umbers is a real fuctio whose domai is the set N of atural umbers. Let s : N R be a sequece. The the values of s are
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 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 informationRiemann Sums y = f (x)
Riema Sums Recall that we have previously discussed the area problem I its simplest form we ca state it this way: The Area Problem Let f be a cotiuous, o-egative fuctio o the closed iterval [a, b] Fid
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
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 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 informationAlternating Series. 1 n 0 2 n n THEOREM 9.14 Alternating Series Test Let a n > 0. The alternating series. 1 n a n.
0_0905.qxd //0 :7 PM Page SECTION 9.5 Alteratig Series Sectio 9.5 Alteratig Series Use the Alteratig Series Test to determie whether a ifiite series coverges. Use the Alteratig Series Remaider to approximate
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 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 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 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 informationLecture 6: Integration and the Mean Value Theorem
Math 8 Istructor: Padraic Bartlett Lecture 6: Itegratio ad the Mea Value Theorem Week 6 Caltech - Fall, 2011 1 Radom Questios Questio 1.1. Show that ay positive ratioal umber ca be writte as the sum of
More information