Numerical schemes for scalar conservation laws with stiff source terms
|
|
- Tiffany Roberts
- 5 years ago
- Views:
Transcription
1 XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Maemáica Aplicada Ciudad Real, 2-25 sepiembre 29 (pp. 8) Numerical schemes for scalar conservaion laws wih siff source erms R. Dona, I. Higueras 2, A. Marínez-Gavara 3 Dp. de Maemàica Aplicada, Univ. de València, Apdo. 46, Burjasso. dona@uv.es. 2 Dpo. de Ingeniería Maemáica e Informáica, Univ. Pública de Navarra, Apdo. 36, Pamplona. higueras@unavarra.es. 3 Dpo. de Ecuaciones Diferenciales y Análisis Numérico, Univ. de Sevilla, Apdo. 42, Sevilla. gavara@us.es. Palabras clave: siff source erm, conservaion laws Resumen We consider a simple model problem of a hyperbolic PDE wih a siff source erm, inroduced by Le Veque and Yee in [5], as he basic es case o explore he abiliy of Implici-Explici Runge-Kua schemes o produce and mainain non-oscillaory reacion frons. Since Srong Sabiliy conceps canno be applied in a direc manner o he model problem, we develop a Weak Sabiliy heory ha provides a convenien framework o sudy he non-oscillaory properies of numerically compued reacion frons.. Inroducion Many physical problems are governed by hyperbolic conservaion laws wih non vanishing siff source erms. These problems could describe he effec of relaxaion as in he kineic heory of gases, chemical reacions, elasiciy wih memory, waer waves, raffic flows, ec. In some problems he source erms depend only on he soluion, i.e. s(x, u) = s(u) and ye he soluion naurally develops srucures in which he source erms are nonzero, and possibly large, only over very small regions in space. This ofen happens if he source erms model chemical reacions beween differen species (reacing flow) in cases where he reacions happen on ime scales much faser han he fluid dynamic ime scales. Then, soluions can develop hin reacion zones where he chemical-kineics aciviy is concenraed. Such problems are said o have siff source erms, in analogy wih he classical case of siff ordinary differenial equaions (ODEs).
2 R. Dona, I. Higueras, A. Marínez-Gavara Numerical difficulies ofen appear when he fas reacions are in near-equilibrium during mos of he compuaion. Some ime scales, ypically hose driving he reacion erms, are several orders of magniude faser han he scale on which he soluion is evolving and on which one would like o compue. Wih many numerical mehods, including all explici mehods, aking a ime sep appropriae for he slower scale of ineres can resul in violen numerical insabiliy, caused by he fas scales. Sabiliy, meaning he absence of violen oscillaory behavior, can be achieved by using implici mehods. A variey of excellen implici mehods have been developed for solving siff sysems of ODEs, and many of he same echniques can be applied when he source erms are siff in order o obain sable resuls. A differen ype of numerical difficuly is, however, also encounered in numerical simulaions concerning hyperbolic PDEs wih siff source erms: he occurrence of frons propagaing a he wrong speeds. I is observed in [5] ha he essenial numerical difficulies o be encounered in he numerical approximaion of convecion-reacion problems wih siff reacion erms can be idenified and sudied mos easily by looking a he equaion u + u x = µu(u )(u ) < x < >, () 2 where he parameer µ > conrols he siffness of he problem. The soluion wih piecewise consan iniial daa {, if x < xd, u(x, ) = (2), if x > x d, is simply u(x, ) = u(x, ). Here we revisi LeVeque and Yee s model problem: a linear advecion equaion wih a parameer dependen source erm. In order o make more precise he relaion beween he occurrence of he pahological reacion frons, he specific form of he source erm and he various parameers involved in he numerical simulaions, we sar by considering very basic firs order echniques which allow us o analyze he source of he pahology. We obain an explici expression for he discrepancy beween he rue speed of propagaion of he reacion fron and he speed of he numerical fron produced by he scheme, and his allows us o quanify he observaion, made by Le Veque and Yee, ha i is precisely he inclusion of non-equilibrium values, a a discree shock profile, wha induces he pahology. Our analysis highlighs he imporance of ensuring ha he compued numerical profiles are monoone. 2. Mehod of lines discreizaions The applicaion of he mehod of lines o he model problem of he previous secion reduces he PDE o an iniial value problem for a sysem of ordinary differenial equaions (ODEs), U = L(U()) + S(U()), U() = (u(x, ), u(x 2, ),..., u(x N, )) T, (3) for he vecor U() = (U (), U 2 (),..., U N ()) T wih componens U i () u(x i, ). Due o he naure of he problem, differen operaors are assigned o he convecive derivaive and he source erm. 2
3 Numerical schemes for scalar conservaion laws wih a siff source erm The erm L(U) in (3) is he spaial discreizaion operaor of he convecive derivaive erm, u x. The erm S(U) represens he discree approximaion of he source erm, which will always be defined in his work as S(U) i = s(u i ). I is widely acceped ha siff source erms should be handled in an implici fashion, in order o avoid sabiliy problems relaed wih he fas scales. On he oher hand, here are a number of robus, and raher specialized, numerical flux funcions ha can be used if disconinuous, or nearly disconinuous, soluions need o be compued. These observaions lead, in a raher naural way, o consider IMplici-EXplici (IMEX) Runge-Kua (RK) schemes for he ime inegraion of MOL discreizaions (see [6]). The abiliy o rea he convecive par in an explici fashion, while sill mainaining an implici handling of he source erms gives a disinc advanage when designing a general purpose high order, high resoluion numerical scheme. As expeced, he occurrence of frons moving a wrong speeds canno be avoided by using any paricular numerical flux funcions in he discree convecive derivaive. The cause of he delay mechanism can be undersood, and quanified, by considering firs he simples numerical echniques for he model problem. In [3], we examine he limiaions ha are encounered when using he simples Eulerype firs order schemes for he ime discreizaion of he MOL sysem (3). For he sake of numerical illusraion, we shall consider he linear upwind discreizaion of he convecive derivaive in () L(U) j = U j U j. (4) x This echnique is paricularly simple because he only boundary condiions needed for he ime sepping process are he physical inflow boundary condiions, which, for he model problem, we ake as U n =, n. The Explici Euler mehod All explici mehods aking a ime sep appropriae for he slower scale of ineres can resul in violen numerical insabiliy caused by he faser scales. Typically, he compuaions end o become very inefficien because he imesep sizes dicaed by he sabiliy requiremens are much smaller han hose required by accuracy consideraions for he slowly varying soluion. The numerical insabiliy become larger when µ increases, o he poin of rendering he useless. In order o obain a free of unwaned oscillaions, he ime sep has o be reduced (on he same mesh). However, he reducion of only guaranees an effecive conrol over he oscillaions developed in he, bu he approximae soluion obained migh sill be qualiaively wrong. In figure, we display resuls for he numerical simulaion wih µ = on a mesh wih x =,. In figure -(a) we show a ypically oscillaory behavior, which is obained for CF L = / x =,4 (i.e. µ = 4, for larger values of µ more violen oscillaions are observed). In figure -(b) we have considered CF L =,2 (i.e. µ = 2). Here we see ha he soluion looks reasonable, however he disconinuous profile is delayed wih respec o ha of he rue soluion. This delay persiss when we reduce he ime sep even furher, as can be observed in figure -(c), where we show he numerical profile for CF L =, ( µ =,). 3
4 R. Dona, I. Higueras, A. Marínez-Gavara (a) CF L =,4 (b) CF L =,2 (c) CF L =, Figura : Numerical soluion for model problem a =,3 for µ = using he firs order explici scheme. x =,. The Implici Euler mehod In he numerical reamen of siff ODEs, i is cusomary o apply implici echniques in order o bypass he sric requiremens on he ime sep imposed by he siffness of he problem. This echnique was proposed by Ahmad and Berzins in [], as a firs sep in order o sudy he effec of neglecing various erms in he nonlinear solvers involved in higher order implici discreizaions of (). The vecor of unknowns U n+ = (U n+, U2 n+,..., U n+ N ) can be compued from he known soluion a ime n, U n = (U n,..., U N n ) by applying a Newon procedure on he sysem ( Uj n+ = Uj n + U j n+ Uj n+ ) + s(uj n+ ), j =,, N (5) x The Jacobian marix of he sysem is of he form JG(V ) = I L(V ) V S(V ) V. (6) The expression of JL := L depends on he advecive erms. Clearly, i migh be very V complicaed or even impossible o compue if a nonlinear high order scheme is used (see [] for a discussion of alernaives in his case). As expeced, he s are non-oscillaory. When he siffness of he model is increased even furher, he occurrence of frons moving a he wrong speeds is observed again. In figure 2, a fron moving a he wrong speed is displayed, due o lack of proper spaial resoluion. In boh cases, x =, is used and CF L =,9 for he lef plo while CF L =,4 for he righ plo. The numerical ess carried ou in his secion clearly display he wo main problems ha may arise when he balance law () is solved numerically: he presence of oscillaions and a numerical delay in he shock profile. Our numerical experimens wih firs order discreizaions indicae ha implici or semi-implici reamens of he siff source erms can conrol he numerical oscillaions in a sraighforward manner, bu will no solve he problem of he numerical delay in he shock profiles, which sill demands an adequae spaial resoluion. The source of he observed delay can be fully analyzed for hese firs order schemes, as we show in he following secion. 4
5 Numerical schemes for scalar conservaion laws wih a siff source erm (a) CFL=,9 (b) CFL=,4 Figura 2: Numerical soluion for model problem wih µ = a =,3, using he firs order fully implici scheme. x =,. CFL= x. 3. Wave speed analysis The analysis below, which is based on simple consideraions on he wave speed of he fron, shows ha he incorrec speed of propagaion of he numerical fron is, in fac, a by-produc of he spaial discreizaion, ha will always be presen as long as here are values of he unknown ha lie sricly beween and. The area under he disconinuous soluion for he model problem () wih iniial condiion (2) φ() := u(x, )dx (7) Hence φ () =, which is he speed of propagaion of he rue soluion. On he oher hand, if he soluion o he sysem of ODEs (3) is a monoone profile, i is naural o define he area under he discree profile, a ime, as φ x () := x N U i (). (8) For a MOL discreizaion, as esablished in (3), he ime variaion of his quaniy can also be easily compued: d d φ x() = x N i= i= d d U i() = µ x N i= U i (U i )(U i ). (9) 2 Clearly, d d φ x() represens he velociy of propagaion of he discree fron, hence he relaion above implies ha he discree profile can move a a speed ha can be quie differen from ha of he rue profile. In fac, he discrepancy is equal o he delay facor α(µ x, U) = µ x N i= U i (U i )(U i ). () 2 A he fully discree level, where he is given by a sequence (U n i ) i, 5
6 R. Dona, I. Higueras, A. Marínez-Gavara (a) µ x =, CFL=, (b) µ x =, CFL=,2 Figura 3: α(µ x, U n ) as a funcion of. Fully explici scheme (a) µ x =, CFL=,7 (b) µ x =, CFL=,2 (c) µ x =, CFL=,7 Figura 4: α(µ x, U n ) as a funcion of. Fully implici scheme. a each ime sep, he corresponding discree equivalen o (8) is N φ n x := x Ui n, () so ha, he discree wave speed defined by LeVeque and Yee in [5] becomes ws(n) = x Uj n U n j = φn+ x φn x, n >. (2) x j j The discree wave speed, ws(n) in (2) can easily be compued for he fully discree Explici Euler (EE) and Implici Euler (IE). A sraighforward compuaion easily leads o ws(n) EE = + α(µ x, U n ), ws(n) IE = + α(µ x, U n ). The funcion s(u) = u(u )(u,5) saisfies s(u) 5 2 for u [, ] and in a ypical simulaion here are only a small number of poins ha conribue o he sum, precisely he poins a he discree disconinuous profile, hence i is o be expeced ha α(µ x, U) = O(µ x 2 ), bu non-zero. Thus, he simples way o ensure a correc speed of propagaion for he fron is o ensure ha µ x is below a securiy hreshold. In figures 3 and 4, we show plos of α(µ x, U n ) wih respec o n for he firs order schemes (EE, IE). We noice ha he behavior of his funcion depends on he CFL number, and ha i migh be oscillaory, bu in all cases α(µ x, U n ) = O(5 2 µ x), i= 6
7 Numerical schemes for scalar conservaion laws wih a siff source erm as expeced and provided ha he soluion is a monoone profile wih only a finie (and small) number of non-equilibrium poins. Hence, for hese schemes we can only aemp o obain a qualiaively correc speed of propagaion for he numerical fron if we ensure ha µ x is below a securiy hreshold. If he is a monoone profile wih values in [, ], wih a finie, and small, number of ransiion values, we have ha for µ x, he delay becomes O( 2 ), which would render i undeecable upon visual examinaion. I should be noiced ha in order o be compuaionally meaningful, ws(n) in (2) needs o be compued on a monoone (i.e. non-oscillaory) discree profile. For homogeneous conservaion laws, he use of TVD consrains is a general procedure ha can be used o ensure a sable, and oscillaion free,. However, for balance laws, TVD consrains migh be oo srong, as i is indeed he case for he model problem. Hence, we seek o enforce weaker forms of sabiliy, and his issue will be sudied in he following secion. 4. Weak sabiliy In general, when a PDE is solved numerically, i is naural o require ha he numerical soluion saisfies as many qualiaive properies of he as possible. Sabiliy requiremens sem from he desire o have numerical schemes ha preserve, a he discree level, cerain properies of he analyic soluion of he problem o be solved. For homogeneous hyperbolic conservaion laws, he fac ha he Toal Variaion (TV) of he soluion does no increase in ime has lead o he developmen of he so-called TVD (Toal-Variaion Diminishing) schemes, a highly successful sabiliy framework ha emerged from imposing his condiion on he numerical scheme. In he case of he balance law, he heorem 5. in [3] (see also [2]) esablishes a monooniciy propery ha can be used o enforce a weaker ype of sabiliy in he numerical schemes. In paricular in he LeVeque and Yee s model problem () wih he iniial profile (2), he applicaion of he heorem gives ha u(x, ), (3) and hence, i makes sense o require similar inequaliies for he, U n i, i =,, N n, (4) perhaps under cerain sepsize resricions. The numerical experimens done [3] show ha numerical oscillaions are always associaed o violaions of propery (4) and hence we may consider i as a key propery o avoid he producion of numerical oscillaions. We, hus, inroduce a weaker form of sabiliy ha seeks o comply wih propery (3). We refer o his propery as weak sabiliy (WS), as opposed o he srong sabiliy (SS) conceps, such as in [4], which prevens growh in a given norm (or semi-norm), ha canno be applied o he model problem due o he paricular form of he source erm. A numerical scheme is hen ermed weakly sable if i saisfies (4) provided U i. Wih he aim o provide useful ools ha lead us o expec a monoone behavior in he s, we analyze in [3] he weak sabiliy properies of firs order schemes 7
8 R. Dona, I. Higueras, A. Marínez-Gavara of he ype considered in he previous secion. As in he case of Srong Sabiliy, we expec o obain Weak Sabiliy only under cerain resricions on x and. The weak sabiliy propery of he EE scheme for he model problem, can be wrien in erms of CFL condiion as CF L = x µ x. (5) This is a sufficien condiion, however i seems o be raher sharp, according o our numerical experimens. For example, for µ x =, we obain weak sabiliy for CF L 2/3,66, while for µ x = weak sabiliy is guaraneed only for CF L,66. In he fully implici firs order scheme, we are led o expec monoone profiles for s provided hese soluions can be compued. Noice ha if s(u), he implici reamen of he convecive derivaive leads o a linear sysem wih a non-singular bi-diagonal marix. However, here is no guaranee in general ha U n+ can be compued independenly of for s(u) in he model problem. In pracice, we have observed ha an increase in he CFL number seems o lead o convergence problems in he Newon-like ieraive process. Conclusion. We have analyzed he phenomenon of waves propagaing a nonphysical speeds for MOL schemes obained from firs order upwind spaial discreizaion of he convecive derivaive in a model problem proposed by LeVeque and Yee in [5]. The use of he simple upwind discreizaion allows us o carry ou a deailed sudy of he delay facor. I is seen ha he delay or advance of he discree wave profile in a MOL discreizaion is conrolled mainly by µ x. Even hough we only showed numerical resuls for he firs order upwind discreizaion of he convecive erm, he behavior obained in his case is ypical of wha would be obained wih more sophisicaed, sae of he ar conservaive discreizaions for he soluion of he model IVP. The advanage of he firs order upwind discreizaion is ha i allows for a raher sysemaic analysis of he numerical soluion, which allowed us o esablish cerain bounds ha, for all pracical purposes, ensure monooniciy of he numerical profiles. Acknowledgemen The auhors acknowledge suppor from MTM and MTM Referencias [] I. Ahmad and M. Berzins. MOL solvers for hyperbolic PDEs wih source erms. Mah. Compu. Simulaion, 56(2):5 25, 2. Mehod of lines (Ahens, GA, 999). [2] A. Chalabi. On convergence of numerical schemes for hyperbolic conservaion laws wih siff source erms. Mah. Compu., 66(28): , 997. [3] R. Dona, I. Higueras, and A. Marínez-Gavara. On sabiliy issues for IMEX schemes applied o hyperbolic equaions wih siff reacion erms. (in preparaion). [4] S. Golieb, W. C.Shu, and E. Tadmor. Srong sabiliy preserving high-order ime discreizaion mehods. SIAM Rev., 43()(TR-2-5):89 2, 2. [5] R. J. LeVeque and H. C. Yee. A sudy of numerical mehods for hyperbolic conservaion laws wih siff source erms. J. Compu. Phys., 86():87 2, 99. [6] L. Pareschi and G. Russo. Implici-Explici Runge-Kua schemes and applicaions o hyperbolic sysems wih relaxaion. J. Sci. Compu., 25(-2):29 55, 25. 8
Unsteady Flow Problems
School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationThe expectation value of the field operator.
The expecaion value of he field operaor. Dan Solomon Universiy of Illinois Chicago, IL dsolom@uic.edu June, 04 Absrac. Much of he mahemaical developmen of quanum field heory has been in suppor of deermining
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationu(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x
. 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationWall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.
MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationMA 366 Review - Test # 1
MA 366 Review - Tes # 1 Fall 5 () Resuls from Calculus: differeniaion formulas, implici differeniaion, Chain Rule; inegraion formulas, inegraion b pars, parial fracions, oher inegraion echniques. (1) Order
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationAn Iterative Method for Solving Two Special Cases of Nonlinear PDEs
Conemporary Engineering Sciences, Vol. 10, 2017, no. 11, 55-553 HIKARI Ld, www.m-hikari.com hps://doi.org/10.12988/ces.2017.7651 An Ieraive Mehod for Solving Two Special Cases of Nonlinear PDEs Carlos
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationLecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples
EE263 Auumn 27-8 Sephen Boyd Lecure 1 Overview course mechanics ouline & opics wha is a linear dynamical sysem? why sudy linear sysems? some examples 1 1 Course mechanics all class info, lecures, homeworks,
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationKinematics and kinematic functions
Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationChapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...
More informationEE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?
EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of
More informationarxiv:quant-ph/ v1 5 Jul 2004
Numerical Mehods for Sochasic Differenial Equaions Joshua Wilkie Deparmen of Chemisry, Simon Fraser Universiy, Burnaby, Briish Columbia V5A 1S6, Canada Sochasic differenial equaions (sdes) play an imporan
More informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More informationAdvanced Organic Chemistry
Lalic, G. Chem 53A Chemisry 53A Advanced Organic Chemisry Lecure noes 1 Kineics: A racical Approach Simple Kineics Scenarios Fiing Experimenal Daa Using Kineics o Deermine he Mechanism Doughery, D. A.,
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationCHBE320 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang
CHBE320 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Spring 208 Dep. of Chemical and Biological Engineering CHBE320 Process Dynamics and Conrol 4- Road Map of he Lecure
More informationVelocity-based boundary integral equation formulation in the time domain
Mesh Reducion Mehods 271 Velociy-based boundary inegral equaion formulaion in he ime domain G. D. Manolis & C. G. Panagioopoulos Deparmen of Civil Engineering, Arisole Universiy of Thessaloniki, Thessaloniki,
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More informationON THE BEAT PHENOMENON IN COUPLED SYSTEMS
8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationThis document was generated at 1:04 PM, 09/10/13 Copyright 2013 Richard T. Woodward. 4. End points and transversality conditions AGEC
his documen was generaed a 1:4 PM, 9/1/13 Copyrigh 213 Richard. Woodward 4. End poins and ransversaliy condiions AGEC 637-213 F z d Recall from Lecure 3 ha a ypical opimal conrol problem is o maimize (,,
More informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More informationHamilton Jacobi equations
Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Noes for EE7C Spring 018: Convex Opimizaion and Approximaion Insrucor: Moriz Hard Email: hard+ee7c@berkeley.edu Graduae Insrucor: Max Simchowiz Email: msimchow+ee7c@berkeley.edu Ocober 15, 018 3
More informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationIntroduction to Numerical Analysis. In this lesson you will be taken through a pair of techniques that will be used to solve the equations of.
Inroducion o Nuerical Analysis oion In his lesson you will be aen hrough a pair of echniques ha will be used o solve he equaions of and v dx d a F d for siuaions in which F is well nown, and he iniial
More informationExplaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015
Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become
More informationOrdinary differential equations. Phys 750 Lecture 7
Ordinary differenial equaions Phys 750 Lecure 7 Ordinary Differenial Equaions Mos physical laws are expressed as differenial equaions These come in hree flavours: iniial-value problems boundary-value problems
More informationCHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang
CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS Professor Dae Ryook Yang Fall 200 Dep. of Chemical and Biological Engineering Korea Universiy CHE302 Process Dynamics and Conrol Korea Universiy
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationCHEMICAL KINETICS: 1. Rate Order Rate law Rate constant Half-life Temperature Dependence
CHEMICL KINETICS: Rae Order Rae law Rae consan Half-life Temperaure Dependence Chemical Reacions Kineics Chemical ineics is he sudy of ime dependence of he change in he concenraion of reacans and producs.
More informationA Shooting Method for A Node Generation Algorithm
A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan
More informationThe L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay
Advances in Dynamical Sysems and Applicaions ISSN 973-5321, Volume 6, Number 2, pp. 177 184 (211) hp://campus.ms.edu/adsa The L p -Version of he Generalized Bohl Perron Principle for Vecor Equaions wih
More informationPade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol
Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne
More informationSliding Mode Controller for Unstable Systems
S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationA car following model for traffic flow simulation
Inernaional Journal of Applied Mahemaical Sciences ISSN 0973-076 Volume 9, Number (206), pp. -9 Research India Publicaions hp://www.ripublicaion.com A car following model for raffic flow simulaion Doudou
More informationChapter 8 The Complete Response of RL and RC Circuits
Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior
More information