Lack of BV Bounds for Approximate Solutions to the p-system with Large Data

Size: px
Start display at page:

Download "Lack of BV Bounds for Approximate Solutions to the p-system with Large Data"

Transcription

1 Lack of BV Bounds for Approimae Soluions o he p-sysem wih Large Daa Albero Bressan (, Geng Chen (, and Qingian Zhang ( (* Deparmen of Mahemaics, Penn Sae Universiy Universiy Park, Pa. 680, U.S.A. (** School of Mahemaics Georgia Insiue of Technology Alana, Ga. 0, U.S.A.. s: bressan@mah.psu.edu, gchen7@mah.gaech.edu, zhang q@mah.psu.edu January, 04 Absrac We consider fron racking approimae soluions o he p-sysem of isenropic gas dynamics. A ineracion imes, he ougoing wave frons have he same srengh as in he eac soluion of he Riemann problem, bu some error is allowed in heir speed. For large BV iniial daa, we consruc eamples showing ha he oal variaion of hese approimae soluions can become arbirarily large, or even blow up in finie ime. This happens even if he densiy of he gas remains uniformly posiive. Inroducion For hyperbolic sysems of conservaion laws in one space dimension, a saisfacory eisenceuniqueness heory is currenly available for enropy weak soluions wih small oal variaion [, 5]. The well-posedness of he Cauchy problem holds also in he case of large daa, as long as he oal variaion remains bounded [, 4]. A major remaining open problem is wheher, for large BV iniial daa, he oal variaion remains uniformly bounded or can blow up in finie ime. Eamples of finie ime blow up have been consruced in [, ]. However, hese sysems do no come from physical models and do no admi a sricly conve enropy. For iniial daa having small oal variaion, regardless of he order in which wave frons cross each oher, he Glimm ineracion esimaes [9] show ha he oal srengh of waves remains small for all imes. There are few eamples of hyperbolic sysems where uniform BV esimaes hold also for soluions wih large daa [7, 0]. In he presen paper we sudy BV bounds for he p-sysem of isenropic gas dynamics in Lagrangean variables: { v u = 0, (. u + p(v = 0. Here u is he velociy, ρ is he densiy, v = is specific volume, while p = p(v is he pressure. Our main concern is wheher, for fron racking approimae soluions o he p- sysem, uniform BV esimaes can be esablished. More precisely, we sudy he following quesion:

2 (Q Consider a fron racking approimaion for (. wih large BV iniial daa. Assume ha, a each ineracion ime, he ougoing wave frons have he same srengh as in he eac soluion of he Riemann problem bu some error is allowed in heir speed. Can he oal variaion of such approimae soluion become arbirarily large? In his paper, some eamples will be consruced, showing ha he oal srengh of waves in a fron racking approimaion can indeed approach infiniy in finie ime. This confirms he non-eisence of a Lyapunov funcional which decreases a every wave-fron ineracion, as proved in [8]. For ineracions occurring near vacuum, i was already noiced in [6] ha uniform Glimmype esimaes were no longer valid. I hus comes as no surprise ha an approimae fron racking soluion can be consruced, where he oal srengh of waves (measured by he change in Riemann invarians blows up in finie ime. Remarkably, our las wo eamples show ha an arbirarily large amplificaion of he oal variaion is sill possible even if he gas densiy remains uniformly posiive. I should be clear ha he presen counereamples do no prove ha, for large BV soluions of he p-sysem, he oal variaion can blow up in finie ime. Indeed, we sill conjecure ha global BV bounds do hold. Our analysis only shows ha such BV bounds canno be proved by wave ineracion esimaes alone, and addiional properies of soluions mus be aken ino accoun. Apparenly, he decay of rarefacion waves due o genuine nonlineariy [, 4, 6, 0] should be used in a crucial way. In he las secion we revisi wo of he earlier eamples and show ha, if such decay is aken ino accoun, hese specific ineracion paerns do no produce blow up. Wave ineracions for he p-sysem We review here some sandard properies of characerisic curves and of shock curves. For deails we refer o [8, 9]. To simplify he compuaions, we assume ha in (. he pressure has he special form p(v = v = ρ, p (v = v 4. Smooh soluions of (. saisfy he quasilinear sysem ρ + ρ u = 0, u + ρ ρ = 0, (. wih characerisic speeds The variables λ = ± ρ, (. w = u, w = ρ + u, (. provide a coordinae sysem of Riemann invarians, in he u-ρ plane.

3 ρ Q Q η(ε σ ε P P ε P Q P η(ε Q Figure : Ineracion of a -shock and a -rarefacion u A shock wih lef sae (u, and righ sae (u +, ρ +, raveling wih speed λ, saisfies he Rankine-Hugonio equaions ( λ = u u +, ρ + λ(u + u = p(v + p(v = ρ +. The La admissibiliy condiion here yields u + < u. Hence u + u = ( ρ + ( ρ ρ +, (.4 Seing from (.4 i follows λ = ± ρ ρ +. (.5 ρ + s = u u +, θ = ρ +, ( θ( θ s = ρ. (.6 θ For he p-sysem, he ineracion of wave frons has been horoughly sudied [9, 8]. For reader s convenience, we derive here wo elemenary esimaes which will be used in he sequel. Consider a -shock wih lef and righ saes (see Fig. P = (u,, Q = (u +, ρ +, wih srengh σ = (ρ + u + ( u

4 measured by he change in he -Riemann invarian. Assume ha his shock crosses a small -wave (compression or rarefacion of srengh σ = ε. We seek an esimae on he size of he ougoing waves, up o leading order. Le ( P = (u ε, ε, Q = u + η(ε, ρ + η(ε be he lef and righ saes across he -shock afer he ineracion. Seing s(ε = s ε + η(ε, θ(ε = θ η(ε ε, ψ(θ = ( θ( θ θ (.7 from (.6 i follows Differeniaing w.r.. ε, a ε = 0 we obain Two cases are relevan o our analysis. s(ε = ( ε ψ(θ(ε. (.8 η = ψ(θ + ψ(θ ψ (θθ = s + s θ + θ 4 θ θ η. (.9 CASE : The -shock has small ampliude. In his case, since shock and rarefacion curves coincide up o second order, for s small we have he epansion Insering (.0 in (.9 we obain θ = + s + O(s. (.0 ( s + O(s 4 (. η = + CASE : The -shock has a fied srengh σ, while he densiy approaches zero. By definiion, he srengh is compued by θ σ = (θ + s = (θ + θ + θ. (. As 0, from (. i follows ha θ. Indeed, dropping lower order erms we find θ / σ, θ ( / σ. (. Insering (. in (.9 and reaining only leading order erms, we obain η σ + σ θ (θ η, 4

5 η σ + σ + σ ( σ 4/ / ( σ /, η as 0. (.4 In oher words, when an infiniesimal -wave crosses a -shock of fied srengh σ, is size is amplified by a facor η κ / which becomes arbirarily large as 0, i.e. as he densiy approaches vacuum. As a special case, if a -shock of srengh σ = crosses a small -wave (eiher compression or rarefacion of srengh ε 0, he srengh of he shock remains consan, while he srengh of he -wave is amplified by a facor η(ε ε > ρ /. (.5 The above esimae is valid as soon as he densiy of he lef sae (i.e., ahead of he -shock is sufficienly small. One more case will be of relevance. Consider a small -shock wih righ sae (u +, ρ + = (0,. By (.4 is lef sae (u, saisfies ( u = ρ. (.6 Taking = s, we obain ( ( u = s s + s4 ( + s + s + = s + s 6 + s 6 + o(s. (.7 Ne, consider a small -shock, again wih righ sae (ũ +, ρ + = (0,. Le (ũ, be he lef sae. Taking = + r, by (.6 we now have ( ( ũ = r + r + r4 ( r + r = r + r 6 r 6 + o(r. (.8 Imposing u = ũ yields r = s + s4 + o(s4. Referring o Fig. 5, consider he poins and le B = A = (0,, C = (u,, D = (u,, ( u +, + Then he slope of he segmen A B is = s o(s4 s + s 6 + s4 + o(s4 ( s + s 6 + s4 + o(s4, + s4 6 + o(s4 5 = s + o(s. (.9

6 Head-on ineracions 0 δ 0 Figure : The head-on ineracion of a -shock wih a rain of smooh -waves, in Eample. Eample. Consider an iniial daa ( ρ, ū consising of a -shock of srengh σ =, approaching a rain of smooh -waves, near vacuum (Fig.. In erms of Riemann invarians (w, w, assume ha he iniial daa is given by { w (0, = w ( = 0, w (0, = w ( = α ( + sin β for [0, δ], (. for suiable consans α, β > 0 and some δ > 0 suiably small. In addiion, assume ha his iniial daa has a -shock of size σ = a =, and is consan on he half lines where < 0 and >. For his iniial daa we consruc an approimae soluion such ha, a every ineracion, he srengh of ougoing waves is he same as in an eac soluion. However, insead of (. or (.5, we le he waves ravel wih consan speeds, say, c for he -shock and c for he -waves, for some consan c > 0. Then, as T = δ/c, all he -waves cross he shock and he oal variaion of he soluion approaches infiniy. Indeed, by (.5, choosing δ > 0 sufficienly small he following holds. A any ime 0 < < T we have w, (, + c ρ( / w,( for all [0, 0 ] such ha + c > c. (. ( α / / α α ( + sin β β α β cos β, If we now choose 0 < α/ < β < α <, hen he iniial daa has bounded variaion, because 0 w, ( d = 0 α α ( + sin β β α β cos β d <. On he oher hand, as T he oal variaion blows up, because 0 lim /c c w, (, + c d = ( 0 α (α/ ( + sin β β (α/ β cos β d =

7 In he previous eample he iniial daa conains vacuum. Moreover, in he erminal profile he large variaion is achieved a very low gas densiy. The ne eample shows ha his arbirarily large amplificaion of he oal variaion can be achieved also wih an iniial daum having uniformly posiive densiy. As before, we require ha a each ineracion he srengh of ougoing waves is he same as in an eac soluion, bu we allow a small error in he wave speed. Namely, in our approimae soluion all -waves ravel wih speed c while -waves ravel wih speed c. ρ = 0 ρ = ρ = ρ = a b Figure : Evoluion of he densiy profile in Eample. Top lef: a ime = 0 he iniial densiy is uniformly posiive. Top righ: afer wo rarefacion waves of he opposie families cross, in he middle secion he same configuraion in Eample is recreaed. Boom lef: when he -shock crosses he rain of -waves a low densiy, he oal variaion grows wihou bounds. Boom righ: afer crossing a -compression, an arbirarily large amoun of oal variaion occurs wihin he inerval [a, b], a uniformly posiive densiy. Eample. As shown in Fig., consider an iniial daa similar o (., bu wih he inserion of wo addiional rarefacion waves, and a compression. A ime, afer crossing rarefacion waves of he opposie family, a ime = he rain of -waves and he -shock recreae he iniial daa in Eample. A ime, he oal variaion of he second Riemann invarian w (, becomes infinie. A ime, afer crossing a -compression, his infinie oal variaion occurs a uniformly posiive densiy. In he - plane, he soluion is described in Fig. 4. Remark. Eample shows ha here eiss an iniial profile (ū, ρ BV wih ρ( ρ 0 > 0, and an approimae soluion wih frons moving a consan speed ±c such ha he following holds. If a each ineracion he srengh of ougoing waves is he same as in he eac soluion of he Riemann problem, hen he oal variaion blows up in finie ime. More 7

8 precisely, a some ime { } ρ(, for 0, To.Var. ρ(, ; [0, [ =. (. In paricular, is is no possible o pu a coninuous weigh on wave srenghs, possibly approaching infiniy as ρ 0, in order o conrol he oal variaion. Indeed, for any choice of he weighs, coninuous on he region where ρ > 0, he iniial weighed oal variaion will be finie, while he final weighed oal variaion will be infinie. To.Var.= + 8 ρ > ρ ~ 0 compression small waves raref. raref. shock Figure 4: The wave paern described in Eample. The bo shows he region where he ineracions in Eample ake place. 4 An eample wih uniformly posiive densiy In he previous eamples, he blow up of he oal variaion was achieved because waves crossing a shock of uni srengh were amplified by an arbirarily large facor as he gas densiy approached vacuum. The following eample shows ha he oal variaion can become arbirarily large even if he gas densiy remains uniformly bounded away from zero, a all imes. Eample. STEP. We begin by consrucing a symmeric ineracion paern conaining four wave frons, as shown in Fig. 5. We choose he srenghs of he wo large shocks S, S and of he wo inermediae waves in such a way ha, afer a whole round of ineracions, hese srenghs are he same as a he iniial ime. Working in he u-ρ plane, his is done as follows. (i Sar by consrucing wo symmeric shocks: he -shock A C and he -shock A C, approaching each oher. (ii Deermine he wo ougoing shocks DA and DA, resuling from he crossing of he above wo shocks. 8

9 D S S 6 Ur A B B C A U l 5 4 Figure 5: A periodic ineracion paern. (iii Consruc a square having wo opposie verices a C and D. Call B, B he remaining wo verices. (iv Choose U l so ha he wo poins B and A are on he same -shock curve wih lef sae U l. Symmerically, choose U r so ha he wo poins B and A are on he same -shock curve wih righ sae U r. The eisence of saes U l, U r saisfying he condiions (iv is now proved (see Fig. 6. Lemma. In he u-ρ plane, consider wo poins B = (u, ρ and A = (u, ρ. Assume ha (i u < u, and ρ > ρ. (ii Calling A = (u, ρ he poin on he -shock curve wih righ sae B having he same ρ-componen as A, one has u < u. Then here eiss a sae U l = (u l, ρ l, wih 0 < ρ l < ρ, such ha boh B and A lie on he -shock curve wih lef sae sae U l. Proof. We shall use (.4 wih (u, = (u l, ρ l while (u +, ρ + = (u, ρ or (u +, ρ + = (u, ρ. To prove he lemma we need o find (u l, ρ l such ha ( u l u = ( ρ l ρ ρ l ρ (, u l u = ( ρ l ρ ρ l ρ. (4. This will be achieved if we can find ρ l such ha u u = G(ρ l, (4. where G is he funcion defined as. ( G(ρ = ( ρ ρ ρ ρ ( ρ ρ ( ρ ρ. The assumpion (ii implies G(ρ = ( ρ ( ρ ρ ρ < u u. 9

10 Moreover, a direc compuaion yields ρ l G(ρ l < 0, lim G(ρ ρl l = +. 0 Therefore here eiss a unique value of ρ l such ha G(ρ l = u u. ρ ρ B 0 u A U A u U = (u, ρ l l l u Figure 6: By moving he poin U along he -shock curve wih righ sae B, we evenually reach a lef sae U l such ha he -shock curve hrough U l conains A as well. Having consruced he above wave curves, consider he following ineracion paern, shown in Fig. 5, righ: A ime he iniial daum consising of four shocks, connecing he saes U l, A, C, A, U r. A ime he profile sill consiss of four shocks (he wo middle ones have crossed each oher, connecing he saes U l, A, D, A, U r. A ime he profile consiss of wo large shocks and wo rarefacions (he wo middle shocks have joined he big ones, generaing wo rarefacions, connecing he saes U l, B, D, B, U r. A ime 4 he profile sill consiss of wo large shocks and wo rarefacions (he wo rarefacions have crossed each oher, connecing he saes U l, B, C, B, U r. A ime 5 he iniial daum consising of four shocks (he wo rarefacions have impinged on he big shocks, generaing wo inermediae shocks of he opposie families, connecing he saes U l, A, C, A, U r, eacly he same as a ime. The paern hus repeas iself. Remark. I is clear ha for he above eample i is essenial o have large oal variaion. Indeed, if we choose he middle shock A C small, hen he line A B will be almos horizonal and he poin U r mus be far away, wih densiy close o zero. On he oher hand, in his soluion obained by fron racking he densiy rivially remains uniformly bounded away from zero. Since he ineracion paern is periodic, we conclude ha Even under he assumpion ha he densiy remains uniformly posiive, here is no way o consruc a Lyapunov funcional conrolling he oal variaion, which is sricly decreasing a every ineracion. 0

11 In he above approimae soluion he wave srenghs follow a periodic paern. To achieve an arbirarily large amplificaion of he oal variaion, a furher consrucion is needed. B P P A Figure 7: A periodic paern ha amplifies a small wave fron. STEP : As shown in Fig. 7, righ, on op of he previous paern we add a very small wave fron. To fi he ideas, consider a -rarefacion of srengh ε > 0, locaed a A. Wihin a ime period, his fron will Cross he inermediae -shock. Inerac wih he large -shock a P producing a -compression. Cross he inermediae -shock. Cross he inermediae -rarefacion. Inerac wih he large -shock a P producing a -rarefacion. Cross he inermediae -rarefacion. We analyze he case where he wo middle shocks are small and he densiy of heir lef sae is. We claim ha, when he addiional fron reaches B, is size will be increased by a facor κ >. Indeed, when a small wave of srengh ε crosses a shock of he opposie family of srengh s a densiy =, by (. he srengh of he ougoing fron is ( ε + = + s + o(s ε. (4. When he fron crosses a rarefacion of he opposie family, is srengh does no change. Finally, when he small wave impinges on a large shock a P or a P, we need o esimae he relaive size of he refleced wave fron. Toward his goal, le ρ = ψ(u be he equaion of he shock curve wih righ sae U r, passing hrough boh A and B, as consruced in

12 Lemma. Calling ε, ε + he srenghs of he fron before and afer ineracion, o leading order we have ( ε + = ( ψ (uε = s 6 + o(s ε. (4.4 Indeed, by (.9 we have ψ (u = s / + o(s. Calling ε A and ε B respecively he srenghs of he small wave-fron a A and a B, we hus have ( ( ( ε B = + s + o(s s 6 + o(s ε A = + s + o(s ε A. (4.5 B σ ε ε P ε ε σ P A Figure 8: Righ: a periodic paern ha amplifies infiniesimal waves. Cener: if a fron of srengh σ crosses a compression (no a shock! of size ε and hen a rarefacion of size ε, is final srengh is no changed. The ougoing frons have srenghs ε, ε. Here ε > ε if he fron σ is a shock, oherwise ε = ε. Lef: o consruc a ineracion paern ha yields an arbirarily large oal variaion, we replace he single infiniesimally small fron in Fig. 7 by counably many pairs rarefacion + compression, of opposie size. The oal srengh of hese waves is finie, each fron having srengh ε. STEP. Consider a periodic paern ha amplifies an infiniesimal wave, as in Sep. By coninuiy, here eiss λ > and ε > 0 such ha any wave-fron (rarefacion or compression of size < ε, raveling from A o B along he pah in Fig. 7, is amplified by a facor λ. We now consruc an iniial se of wave frons where he infiniesimal fron is replaced by counably many pairs rarefacion + compression, whose sizes eacly cancel each oher (Fig. 8, lef. A ime = 0, he oal srengh of all hese small frons can be aken o be =. The key observaion is ha each of hese pairs leaves no fooprin on he underlying soluion consruced in Sep. Indeed, if a fron of size σ crosses a rarefacion and a compression of eacly opposie sizes, afer he wo crossings he size of he fron is sill σ (Fig. 8, cener. As a resul, he paern of four large frons reains is periodiciy. By consrucion, afer each period each pair of opposie small wavefrons is enlarged by a facor λ. When a pair grows o size > ε, we can perform a parial cancellaion so ha is size remains [ε/, ε]. Since he oal number of small wave-frons is infinie, afer several periods a larger and larger number of pairs (compression + rarefacion reaches size > ε/. Hence, as, he oal variaion of his approimae soluion grows wihou bounds.

13 5 Blow up in finie ime The previous consrucion shows ha, if he oal variaion is iniially sufficienly large, hen here eiss an ineracion paern ha renders he oal variaion arbirarily large as. This can be achieved even wih a uniform lower bound on he densiy. The ne quesion is wheher one can arrange he order of wave-fron ineracions so ha he oal variaion blows up in finie ime. Noice ha his is no he case in he previous eample. Indeed, if -frons ravel wih speed ẋ [ C, 0[ and -frons ravel wih speed ẋ ]0, C], i akes a uniformly posiive amoun of ime for each inermediae fron o bounce back and forh beween wo large shocks. Hence he arbirarily large amplificaion of he oal variaion is only achieved in he limi as +. T 4 S S S S Figure 9: A paern yielding finie ime blow up of he oal variaion. A each sep he oal amoun of small waves bouncing back and forh beween he wo large shocks keeps increasing. In his secion we briefly indicae how he previous consrucion can be modified, providing finie ime blow up of he oal variaion. The main idea is illusraed in Fig. 9. We consider a sequence of imes 0 = 0 < < < < T. During each ime inerval J i = [ i, i ], a counable number of pairs of small waves (compression + rarefacion is amplified by a very large facor. Before ime i, he large -shock S i is compleely canceled by impinging -rarefacions, and a new -shock S (i+ of he same srengh is recreaed a a locaion closer o he large -shock S. Figures 0 and show how his can be achieved, saring wih very many pairs of small waves (compression + rarefacion. By leing each compression fron collapse o a shock, and hen canceling his shock wih a rarefacion fron of he same family, we obain a rain of pairs of small waves (compression + rarefacion in he opposie family (Fig. 0, lef. By varying he locaions of hese ineracions, insead of many pairs of small waves we can achieve a large compression followed by a large rarefacion fron (Fig. 0, righ. The basic sep is illusraed in Fig.. A large number of small compression+rarefacion pairs produces a large -rarefacion, which sars depleing he -shock along he line AB, and a

14 Figure 0: Saring wih a large number of pairs of small waves, one can generae a large number of similar pairs in he oher characerisic family (lef, or one single large pair of frons (righ. τ D B τ U l S ρ U U l τ C S U r U r U 4 U S S U u A Figure : Lef: beween A and B he shock S is compleely canceled by impinging -rarefacions. Beween C and D, he new -shock S is formed by impinging compression waves. Noice ha he -compression frons emerging from S are used o compleely cancel he -rarefacions ha would oherwise be produced by he ineracions of -compressions wih S. Righ: if he shock S is very large, i canno be canceled by one single large rarefacion fron. Therefore, we need o produce several small rarefacions a subsequen imes, so ha he densiy ρ remains uniformly posiive. 4

15 large compression, which builds up a new -shock along he line CD. Since we canno allow he densiy o become negaive, i may no be possible o cancel he large shock S wih one single large -rarefacion. For his reason, his cancellaion may be accomplished in several sages. For eample, he firs se of -rarefacions reduce he size of he -shock S from (U l, U r o (U, U r. A a ime τ, he profile u(τ, hus conains he -shock (U l, U, he -compression (U, U, and he new -shock (U, U r. A a laer ime τ, he profile u(τ, conains he (shrinking -shock (U l, U 4, he -compression (U 4, U, and he (growing -shock (U, U r. A a laer ime τ, he original -shock has been compleely depleed by -rarefacions. A -shock connecing eacly he same wo saes (U l, U r is formed a a differen locaion, as desired. By canceling he large -shock and reconsrucing i a a differen locaion, shifed o he lef, we can reproduce he paern in Fig. 9. Since a each sep he oal srengh of he small inermediae waves can be amplified by an arbirarily large facor, as T he oal variaion of our approimae soluion blows up o +. 6 Concluding remarks The eamples presened in his paper show ha, if he srengh of wave-frons is compued eacly bu some error is allowed heir speeds, hen he oal variaion of approimae soluions can blow up. I is ineresing o revisi some of he previous eamples, aking ino accoun he decay of rarefacion waves due o genuine nonlineariy. Looking a eac soluions, i becomes clear ha hese paricular ineracion paerns do no yield an arbirarily large amplificaion of he oal variaion. 6. Head-on ineracions, near vacuum. Consider an eac soluion of he sysem (., wih iniial daa as in Eample. We show ha here is no way o choose α, β in (. so ha he following requiremens are simulaneously saisfied: (i The -shock crosses all -waves in finie ime. (ii The -waves do no break before crossing he shock. (iii The sum of srenghs of he -waves is iniially finie, and becomes infinie as hey all cross he -shock. (iv The -waves do no break immediaely afer crossing he shock. For a -shock of uni srengh, assume ha he lef sae (ahead of he shock has densiy 0. Then by (.4 he righ sae (behind he shock has densiy ρ + saisfying ρ +. 5

16 Hence he righ sae and he speed of he shock are given respecively by ρ + ( /, ẋ ρ + = / ρ /. (6. If he iniial profile is given by (., so ha w (0, = 0, w (0, = α ( + sin β [0, δ], (6. hen he requiremen (i will be saisfied provided ha α <. (6. Ne, o make sure ha he -waves do no break before crossing he -shock, we look a he evoluion of w, along characerisics. From ( w w, + w, = 0 i follows ( w w, + w, = w 4 w,. By a comparison argumen, we conclude ha he gradien w, will no blow up before ime T > 0 provided ha he iniial daa saisfy w, (0, w (0, < T. Recalling (6., he condiion (ii is hus saisfied if α + (α β > 0. (6.4 As shown in he discussion of Eample, condiion (iii is saisfied provided ha 0 < α < β < α <. (6.5 To check wheher (iv can be saisfied, le T be he ime when he -shock reaches he origin, crossing all -waves. Denoe by (, y he posiion of a -characerisic saring a (0, y = y, Calling τ(y he ime where his -characerisic crosses he -shock, we find T τ(y y α/. We consider he evoluion of w, (, (, y along his -characerisic. For < τ(y we have w, (, (, y y α β. When his characerisic crosses he -shock a ime τ = τ(y, by (.5, his gradien is amplified by a facor w, (τ+ w, (τ y α/. Moreover, w (τ+ y α/. To make sure ha his gradien remains bounded during he ime inerval [τ(y, T ], we need Therefore we should have y α/ w, (τ w (τ+ (T = O(. y α/ y α β y α/ y α/ = y β = O(. (6.6 This condiion is incompaible wih he requiremen (6.5 ha β > 0. 6

17 T γ 0 Figure : Lef: Before crossing he large -shock, he small -waves do no break, because when he gas densiy ρ > 0 is very small, he sysem is almos linearly degenerae. Cener: afer crossing he -shock, he -compression waves break and a large amoun of cancellaion beween -rarefacions and -shocks occurs. A he erminal ime T he oal srengh of waves is sill finie, due o hese cancellaions. Righ: in he - plane his paern produces an infinie oal variaion only along he (dashed ime-like curve γ. 6. Waves bouncing back and forh beween wo large shocks From our earlier analysis, his should be he paern ha achieves he greaes amplificaion of wave srenghs (Fig.. If he size of he shocks S, S is sufficienly large, he srengh of a refleced -fron σ is almos he same as he srengh of he impinging -fron σ. Aferwards, as his -fron crosses oher -shocks, is srengh increases by a large facor. Repeaing his paern, i may appear ha an arbirarily large amplificaion of wave srenghs can be achieved. The following analysis shows ha his is no he case, if we ake ino accoun he decay of rarefacion waves due o genuine nonlineariy. For some consan c 0, he wo large shocks will have speeds ẋ ( c 0 < 0 < c 0 ẋ (. (6.7 Consider a -rarefacion wave (Fig. emerging from he large -shock a some ime τ and impinging on he opposie -shock a ime τ. The upper and lower esimaes on he velociy yield an esimae of he form τ κτ. (6.8 By wave decay esimaes, he densiy of such -rarefacion a ime τ is C τ τ Cκ κ Therefore, he oal amoun of -rarefacions ha impinge on he large -shock wihin a ime inerval [T 0, T ] is T Cκ κ τ dτ. T 0 An enirely similar esimae holds for he -rarefacions impinging on he large -shock. Ne, fi any large ime T. As in Figure, consider he maimal backward -characerisic hrough he poin (T, (T. This will cross he large -shock a an earlier ime T 0. The oal τ. 7

18 S S T τ T 0 τ ( ( σ σ Figure : The oal amoun of compression waves (including shocks a ime T is bounded in erms of he oal amoun of rarefacion waves ha impinge on he large -shock during he inerval [T 0, T ]. amoun of -shocks (ogeher wih -compressions a a given ime T can be esimaed as: Toal amoun of -shocks a ime T conained in he inerval [ (T, (T ] = O( amoun of -rarefacions impinging on he -shock for [T 0, T ] = O( T T 0 Cκ κ τ dτ = O( ln T T 0 = O(. (6.9 Here O( denoes a quaniy ha remans uniformly bounded (provided ha some upper and lower bounds on he densiy ρ are given. An enirely similar esimae of course holds for -shocks (ogeher wih -compressions. In addiion, he rarefacion waves can be esimaed as oal amoun of -rarefacions a ime T conained in he inerval [ (T, (T ] = O( amoun of -shocks impinging on he -shock for [T 0, T ] oal amoun of -shocks a ime T 0 conained in he inerval [ (T 0, (T 0 ] cons. (6.0 Here he las inequaliy follows from (6.9, wih T replaced by T 0. This yields a uniform a priori bound on he oal srengh of waves produced by his paricular wave-ineracion paern. 8

19 References [] P. Baii and H. K. Jenssen, Blowup in L for a class of genuinely nonlinear hyperbolic sysems of conservaion laws. Discree Conin. Dynam. Sysems 7 (00, [] A. Bressan, Hyperbolic Sysems of Conservaion Laws. The One Dimensional Cauchy Problem. Oford Universiy Press, 000. [] A. Bressan and R. M. Colombo, Unique soluions of conservaion laws wih large daa, Indiana Univ. Mah. J. 44 (995, [4] A. Bressan and R. M. Colombo, Decay of posiive waves in nonlinear sysems of conservaion laws, Ann. Scuola Normale Superiore Pisa IV - 6 (998, 60. [5] A. Bressan, T. P. Liu and T. Yang, L sabiliy esimaes for n n conservaion laws, Arch. Raional Mech. Anal. 49 (999,. [6] A. Bressan and T. Yang, A sharp decay esimae for posiive nonlinear waves, SIAM Jour. Mah. Anal. 6 (004, [7] T. Chang and L. Hsiao, The Riemann problem and ineracion of waves in gas dynamics, Longman Scienific & Technical, Harlow, 989. [8] G. Chen and H. K. Jenssen, No TVD fields for -d isenropic gas flow, Comm. Parial Differenial Equaions, 8 (0, [9] J. Glimm, Soluions in he large for nonlinear hyperbolic sysems of equaions, Comm. Pure Appl. Mah. 8 (965, [0] J. Glimm and P. La, Decay of soluions of sysems of nonlinear hyperbolic conservaion laws, Amer. Mah. Soc. Memoir 0 (970. [] D. Hoff, Invarian regions for sysems of conservaion laws. Trans. Amer. Mah. Soc. 89 (985, [] H. Holden and N. H. Risebro, Fron Tracking for Hyperbolic Conservaion Laws. Springer- Verlag, New York, 00. [] H. K. Jenssen, Blowup for sysems of conservaion laws, SIAM J. Mah. Anal. (000, [4] M. Lewicka, Well-posedness for hyperbolic sysems of conservaion laws wih large BV daa. Arch. Raional Mech. Anal. 7 (004, [5] L. W. Lin, On he vacuum sae for he equaions of isenropic gas dynamics. J. Mah. Analysis Appl. (987, [6] T. P. Liu and J. Smoller, On he vacuum sae for he isenropic gas dynamics equaions, Advances Pure Appl. Mah. (980, [7] T. Nishida, Global soluion for an iniial boundary value problem of a quasilinear hyperbolic sysem. Proc. Japan Acad. 44 (968,

20 [8] O. Oleinik, Disconinuous soluions of nonlinear differenial equaions, Amer. Mah. Soc. Transl. 6, [9] J. Smoller, Shock waves and reacion-diffusion equaions, Second ediion. Springer-Verlag, New York, 994. [0] B. Temple, Sysems of conservaion laws wih invarian submanifolds, Trans. Amer. Mah. Soc. 80 (98,

NO BV BOUNDS FOR APPROXIMATE SOLUTIONS TO P-SYSTEM WITH GENERAL PRESSURE LAW

NO BV BOUNDS FOR APPROXIMATE SOLUTIONS TO P-SYSTEM WITH GENERAL PRESSURE LAW NO BV BOUNDS FOR APPROXIMATE SOLUTIONS TO P-SYSTEM WITH GENERAL PRESSURE LAW ALBERTO BRESSAN, GENG CHEN, QINGTIAN ZHANG, AND SHENGGUO ZHU Absrac. For he p-sysem wih large BV iniial daa, an assumpion inroduced

More information

2. Nonlinear Conservation Law Equations

2. 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 information

Optima and Equilibria for Traffic Flow on a Network

Optima and Equilibria for Traffic Flow on a Network Opima and Equilibria for Traffic Flow on a Nework Albero Bressan Deparmen of Mahemaics, Penn Sae Universiy bressan@mah.psu.edu Albero Bressan (Penn Sae) Opima and equilibria for raffic flow 1 / 1 A Traffic

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Solutions from Chapter 9.1 and 9.2

Solutions from Chapter 9.1 and 9.2 Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is

More information

Vehicle Arrival Models : Headway

Vehicle 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 information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.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 information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 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 information

An Introduction to Malliavin calculus and its applications

An Introduction to Malliavin calculus and its applications An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214

More information

Theory of! Partial Differential Equations!

Theory 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 information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. 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 information

arxiv:math/ v1 [math.nt] 3 Nov 2005

arxiv:math/ v1 [math.nt] 3 Nov 2005 arxiv:mah/0511092v1 [mah.nt] 3 Nov 2005 A NOTE ON S AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION D. A. GOLDSTON AND S. M. GONEK Absrac. Le πs denoe he argumen of he Riemann zea-funcion a he poin 1 + i. Assuming

More information

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity

L p -L q -Time decay estimate for solution of the Cauchy problem for hyperbolic partial differential equations of linear thermoelasticity ANNALES POLONICI MATHEMATICI LIV.2 99) L p -L q -Time decay esimae for soluion of he Cauchy problem for hyperbolic parial differenial equaions of linear hermoelasiciy by Jerzy Gawinecki Warszawa) Absrac.

More information

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Hamilton- 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 information

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

The 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 information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

1. VELOCITY AND ACCELERATION

1. VELOCITY AND ACCELERATION 1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under

More information

Class Meeting # 10: Introduction to the Wave Equation

Class 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 information

The motions of the celt on a horizontal plane with viscous friction

The motions of the celt on a horizontal plane with viscous friction The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and

More information

Differential Equations

Differential 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 information

Theory of! Partial Differential Equations-I!

Theory 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 information

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan Ground Rules PC11 Fundamenals of Physics I Lecures 3 and 4 Moion in One Dimension A/Prof Tay Seng Chuan 1 Swich off your handphone and pager Swich off your lapop compuer and keep i No alking while lecure

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

Matlab and Python programming: how to get started

Matlab 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 information

IB Physics Kinematics Worksheet

IB Physics Kinematics Worksheet IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?

More information

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Hamilton- 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 information

Hamilton Jacobi equations

Hamilton 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 information

Chapter 2. First Order Scalar Equations

Chapter 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 information

4. Advanced Stability Theory

4. Advanced Stability Theory Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium

More information

15. Vector Valued Functions

15. Vector Valued Functions 1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,

More information

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - Prey Model Trajectories and the nonlinear conservation law Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories

More information

5.1 - Logarithms and Their Properties

5.1 - Logarithms and Their Properties Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We

More information

Lecture 20: Riccati Equations and Least Squares Feedback Control

Lecture 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 information

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS

CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien

More information

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4. PHY1 Elecriciy Topic 7 (Lecures 1 & 11) Elecric Circuis n his opic, we will cover: 1) Elecromoive Force (EMF) ) Series and parallel resisor combinaions 3) Kirchhoff s rules for circuis 4) Time dependence

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

Some Basic Information about M-S-D Systems

Some 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 information

Traveling Waves. Chapter Introduction

Traveling Waves. Chapter Introduction Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from

More information

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems. di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x

More information

4.6 One Dimensional Kinematics and Integration

4.6 One Dimensional Kinematics and Integration 4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of

More information

EXERCISES FOR SECTION 1.5

EXERCISES 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 information

Chapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws

Chapter 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 information

A NOTE ON S(t) AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION

A NOTE ON S(t) AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION Bull. London Mah. Soc. 39 2007 482 486 C 2007 London Mahemaical Sociey doi:10.1112/blms/bdm032 A NOTE ON S AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION D. A. GOLDSTON and S. M. GONEK Absrac Le πs denoe he

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)

(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 information

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

!!#$%&#'()!#&'(*%)+,&',-)./0)1-*23) "#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5

More information

t 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...

t 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 information

A Shooting Method for A Node Generation Algorithm

A 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 information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

d 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 information

Convergence of the Neumann series in higher norms

Convergence of the Neumann series in higher norms Convergence of he Neumann series in higher norms Charles L. Epsein Deparmen of Mahemaics, Universiy of Pennsylvania Version 1.0 Augus 1, 003 Absrac Naural condiions on an operaor A are given so ha he Neumann

More information

Cash Flow Valuation Mode Lin Discrete Time

Cash Flow Valuation Mode Lin Discrete Time IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion

More information

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15.

SMT 2014 Calculus Test Solutions February 15, 2014 = 3 5 = 15. SMT Calculus Tes Soluions February 5,. Le f() = and le g() =. Compue f ()g (). Answer: 5 Soluion: We noe ha f () = and g () = 6. Then f ()g () =. Plugging in = we ge f ()g () = 6 = 3 5 = 5.. There is a

More information

Two Coupled Oscillators / Normal Modes

Two 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 information

8. Basic RL and RC Circuits

8. Basic RL and RC Circuits 8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics

More information

arxiv: v1 [math.fa] 9 Dec 2018

arxiv: v1 [math.fa] 9 Dec 2018 AN INVERSE FUNCTION THEOREM CONVERSE arxiv:1812.03561v1 [mah.fa] 9 Dec 2018 JIMMIE LAWSON Absrac. We esablish he following converse of he well-known inverse funcion heorem. Le g : U V and f : V U be inverse

More information

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation

An Invariance for (2+1)-Extension of Burgers Equation and Formulae to Obtain Solutions of KP Equation Commun Theor Phys Beijing, China 43 2005 pp 591 596 c Inernaional Academic Publishers Vol 43, No 4, April 15, 2005 An Invariance for 2+1-Eension of Burgers Equaion Formulae o Obain Soluions of KP Equaion

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0. Time-Domain Sysem Analysis Coninuous Time. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 1. J. Robers - All Righs Reserved. Edied by Dr. Rober Akl 2 Le a sysem be described by a 2 y ( ) + a 1

More information

CONTRIBUTION TO IMPULSIVE EQUATIONS

CONTRIBUTION TO IMPULSIVE EQUATIONS European Scienific Journal Sepember 214 /SPECIAL/ ediion Vol.3 ISSN: 1857 7881 (Prin) e - ISSN 1857-7431 CONTRIBUTION TO IMPULSIVE EQUATIONS Berrabah Faima Zohra, MA Universiy of sidi bel abbes/ Algeria

More information

1. Kinematics I: Position and Velocity

1. Kinematics I: Position and Velocity 1. Kinemaics I: Posiion and Velociy Inroducion The purpose of his eperimen is o undersand and describe moion. We describe he moion of an objec by specifying is posiion, velociy, and acceleraion. In his

More information

Oscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson

Oscillation 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 information

Linear Response Theory: The connection between QFT and experiments

Linear 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 information

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

Vanishing 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 information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

IMPACT OF AN OBLIQUE BREAKING WAVE ON A WALL

IMPACT OF AN OBLIQUE BREAKING WAVE ON A WALL Source: Physics of Fluids Vol 6 No pp 6-64 4 DOI: 6/64445 IMPACT OF AN OIQUE REAKING WAVE ON A WA Jian-Jun SHU School of Mechanical & Aerospace Engineering Nanyang Technological Universiy 5 Nanyang Avenue

More information

( ) = Q 0. ( ) R = R dq. ( t) = I t

( ) = Q 0. ( ) R = R dq. ( t) = I t ircuis onceps The addiion of a simple capacior o a circui of resisors allows wo relaed phenomena o occur The observaion ha he ime-dependence of a complex waveform is alered by he circui is referred o as

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-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 information

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990),

SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F. Trench. SIAM J. Matrix Anal. Appl. 11 (1990), SPECTRAL EVOLUTION OF A ONE PARAMETER EXTENSION OF A REAL SYMMETRIC TOEPLITZ MATRIX* William F Trench SIAM J Marix Anal Appl 11 (1990), 601-611 Absrac Le T n = ( i j ) n i,j=1 (n 3) be a real symmeric

More information

Ordinary Differential Equations

Ordinary 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 information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

Effects of Coordinate Curvature on Integration

Effects of Coordinate Curvature on Integration Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold

More information

RC, RL and RLC circuits

RC, RL and RLC circuits Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.

More information

1 Differential Equation Investigations using Customizable

1 Differential Equation Investigations using Customizable Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for

More information

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

dt = 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 information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Math 334 Fall 2011 Homework 11 Solutions

Math 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 information

t dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP

t dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP Coninuous Linear Programming. Separaed Coninuous Linear Programming Bellman (1953) max c () u() d H () u () + Gsusds (,) () a () u (), < < CLP (Danzig, yndall, Grinold, Perold, Ansreicher 6's-8's) Anderson

More information

The L p -Version of the Generalized Bohl Perron Principle for Vector Equations with Infinite Delay

The 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 information

4.1 - Logarithms and Their Properties

4.1 - Logarithms and Their Properties Chaper 4 Logarihmic Funcions 4.1 - Logarihms and Their Properies Wha is a Logarihm? We define he common logarihm funcion, simply he log funcion, wrien log 10 x log x, as follows: If x is a posiive number,

More information

Numerical Dispersion

Numerical 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

Solutions to Assignment 1

Solutions 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 information

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion

More information

Expert Advice for Amateurs

Expert Advice for Amateurs Exper Advice for Amaeurs Ernes K. Lai Online Appendix - Exisence of Equilibria The analysis in his secion is performed under more general payoff funcions. Wihou aking an explici form, he payoffs of he

More information

Electrical and current self-induction

Electrical and current self-induction Elecrical and curren self-inducion F. F. Mende hp://fmnauka.narod.ru/works.hml mende_fedor@mail.ru Absrac The aricle considers he self-inducance of reacive elemens. Elecrical self-inducion To he laws of

More information

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE

MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS

More information

MEI STRUCTURED MATHEMATICS 4758

MEI STRUCTURED MATHEMATICS 4758 OXFORD CAMBRIDGE AND RSA EXAMINATIONS Advanced Subsidiary General Cerificae of Educaion Advanced General Cerificae of Educaion MEI STRUCTURED MATHEMATICS 4758 Differenial Equaions Thursday 5 JUNE 006 Afernoon

More information

arxiv: v1 [math.fa] 3 Jan 2019

arxiv: v1 [math.fa] 3 Jan 2019 DAMPED AND DIVERGENCE EXACT SOLUTIONS FOR THE DUFFING EQUATION USING LEAF FUNCTIONS AND HYPERBOLIC LEAF FUNCTIONS A PREPRINT arxiv:9.66v [mah.fa] Jan 9 Kazunori Shinohara Deparmen of Mechanical Sysems

More information

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

u(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

u(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 information

Lab 10: RC, RL, and RLC Circuits

Lab 10: RC, RL, and RLC Circuits Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in

More information

Math 333 Problem Set #2 Solution 14 February 2003

Math 333 Problem Set #2 Solution 14 February 2003 Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability 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 information