arxiv: v1 [math-ph] 8 Oct 2015

Size: px
Start display at page:

Download "arxiv: v1 [math-ph] 8 Oct 2015"

Transcription

1 Jordan form, parabolicity and other featres of change of type transition for hydrodynamic type systems. arxi:5.97 [math-ph] 8 Oct 5 B.G. Konopelchenko Dipartimento di Matematica e Fisica Ennio de Giorgi, Uniersità del Salento INFN, Sezione di Lecce, 7 Lecce, Italy konopel@le.infn.it G. Ortenzi Dipartimento di Matematica Pra ed Applicazioni, Uniersità di Milano Bicocca, 5 Milano, Italy gioanni.ortenzi@nimib.it Jly, 8 Abstract Changes of type transitions for the two-component hydrodynamic type systems are discssed. It is shown that these systems generically assme the Jordan form (with Jordan block) on the transition line with hodograph eqations becoming parabolic. Conditions which allow or forbid the transition from hyperbolic domain to elliptic one are discssed. Hamiltonian systems and their special sbclasses and eqations, like dispersionless nonlinear Schrödinger, dispersionless Bossinesq, one-dimensional isentropic gas dynamics eqations and nonlinear wae eqations are stdied. Nmerical reslts concerning the crossing of transition line for the dispersionless Bossinesq eqation are presented too. Introdction Differential eqations and systems of mixed type always hae attracted a great interest de to the presence of both hyperbolic and elliptic regimes, possibility of transition to each other and nmeros interpretations of sch transitions in arios fields of mathematics, physics and applied science (see e.g. []-[7]). Stdy of the properties of the systems of qasi-linear eqations near the transition line (referred also as sonic line, parabolic line or hyperbolic-elliptic bondary) is of particlar interest de to the connection with the problem of nonlinear stability of systems of mixed type. The reslts obtained in the recent papers [8, 9, ] contribted significantly to nderstanding and clarifying the sitation. On the other hand some assmptions made, for instance, in [] seems to be rather restrictie. In particlar, the calclations made in [] are based on the hypothesis that the matrix V for the hydrodynamic type systems t = V( ) x () is degenerated yet diagonalizable on sonic line []. There are, howeer, nmber of systems for which it is not the case. The simplest example is proided by the well-known one-layer Benney system (or dipersionless nonlinear Schrödinger eqation (dnls)) t = x + x, t = x + x. ()

2 Characteristic speeds for () are λ ± = ± and the transition line is gien by eqation (x,t) =. On the transition line the matrix V takes the form ( ) V = () which is obiosly non-diagonalizable. Another example is proided by the dispersionless Bossinesq (db) eqation or the system In this case the matrix V on the transition line is tt = ( ) xx () V = t = x, t = x. ( (5) ), (6) i.e. the Jordan block with zero eigenale. The appearance of Jordan blocks in the these examples is a clear manifestation of the generic strctre for the hydrodynamics type systems on the transition line. These two systems represent two different classes of hydrodynamic systems of mixed type. For dnls eqation () hyperbolic domain ( > ) is separated from the elliptic one ( < ). For the db eqation () the transition line can be crossed and soltions from the hyperbolic domain ( > ) can pass to the elliptic domain ( < ). In the present paper these phenomenons are stdied for the two-component systems () of mixed type. It is shown that generically a two-component system () at the transition line is of Jordan form. Hodograph eqations are manifestly parabolic on the transition line. This parabolic regime separates the hyperbolic domain describing wae propagation and elliptic domain containing qasi-conformal mapping. Conditions nder which soltions of the system (8) may belong to both hyperbolic and elliptic domain or aoid the crossing of the transition line are discssed. Hamiltonian systems are considered in detail as illstratie examples. It is shown that the presence of the Jordan block on the transition line is a typical behaior of Hamiltonian systems. It is also that in the generic case the characteristics in (, ) plane (simple waes) hae niersal behaior ( ) / (7) near the point (, ) of contact with the transition line. The db eqation is characteristic representatie of sch a behaior. Particlar classes of Hamiltonian systems, inclding gasdynamics eqations and nonlinear wae eqations are considered. Nmerical reslts for the db eqation showing the particlarities of crossing of the transition line are presented too. The paper is organized as follows. Some basic well-known reslts for the system (8), inclding hodograph eqations are gien in Section. Behaior of the system (8) on the transition line, its Jordan form, and parabolic character are considered in Section. In section it is shown the behaior of the system near the transition line from the elliptic side. Necessary and sfficient conditions which allows or forbid the crossing of the transition line are discssed in Section 5. These reslts applied to general Hamiltonian systems are presented in Section 6. Special classes of Hamiltonian systems and, in particlar, eqations of motion for isentropic gas eqations and nonlinear wae eqations are considered in Sections 7 and and 8. Some nmerical reslts for the db eqation near to the transition line are presented in Section 9. General formlae We will consider two component qasi-linear system of mixed type of first order ( ) ( )( ) t A B x = C D t x (8)

3 where A,B,C,D are certain real fnctions of and and sbscript denotes deriaties. For conenience we will recall here some basic known facts (see e.g. [,, ]). Generically the matrix ( ) A B V = (9) C D has two distinct eigenales gien by λ ± = A+D± (A D) +BC. () If Ω (A D) + BC > the system is hyperbolic, while at Ω < it is elliptic. In this paper we will assme that Ω(,) is a smooth fnction of and. So, the hyperbolic and elliptic domains are separated by the transition line gie by the eqation Classical hodograph eqations for the system (8) is ( ) ( x D C = B A x Ω((x,t),(x,t)) =. () )( t As a conseqence, the ariables t and x obey the second order eqations and t ). () Ct + A D t Bt (B +D )t +(A +C )t = () ( ) Ax +Cx AD BC ( ) Bx +Dx =. () AD BC If the system (8) has a conseration eqation Q t = P x then Q obeys the eqation CQ +(A D)Q BQ (B A )Q +( D +C )Q =. (5) In the hyperbolic domain there are two real Riemann inariants r + and r sch that the system (8) is eqialent to r ± t = λ ±r ± x (6) with two distinct characteristic speeds λ + and λ. In the elliptic domain λ + and λ are complex-conjgate to each-other and one has the single complex eqation r + t = λ +r + x (7) with r being the complex conjgate to r +. Riemann inariants obey to the system Only two among these eqations are independent say or (A λ ± )r ± +Cr ± =, Br ± +(D λ ±)r ± =. (8) Br ± +(D λ ±)r ± = (9) r ± = A D± Ω r ± B. () Eqations () are diagonal form of the hodograph eqations () rewritten as ( t x ) = ( B A BC AD D )( t x ). ()

4 Indeed the eigenales µ ± of the matrix present in () are µ ± = A D± Ω B. () The characteristics for eqation () are defined by the eqation ( ) d = µ ±. () d ± Riemann inariants are constants along these characteristics ± in the hodograph space. Eqations() are those which define simple waes for the system(8). The simple waes are the hodograph conterpart of the sal characteristics ( ) dx dt ± = λ ± in the space (x,t). We note also that the components y and y of the eigenector y of the matrix V in (9) on the transition line obey the eqation (A D)y +By =. () Hydrodynamic type systems and, in particlar, the system (8) exhibit one more important phenomenon, the so-called gradient catastrophe, i.e. nbondedness of deriaties of and at finite x and t while (x,t) and (x, t) remain bonded (see e.g. []). Interference of the gradient catastrophe and crossing the transition line is a rather complicated problem. To simplify the analysis we will assme in the rest of the paper that the soltions of the system (8) aoid gradient catastrophe, at least, before the crossing of the transition line (if so). Transition line and Jordan form Stdy of the behaior of systems of mixed type on the transition line and nearby is fndamental for nderstanding their properties. The system (8) can be of mixed type only when BC <. In the case BC (inclding symmetric matrices) it is hyperbolic except the degeneration at the set of points defined by the eqation A = D and B = (or C = ) in generic case or at the line in the case A = D and B = C =. It was already shown in the Introdction that on the transition lines the dnls and db eqations assme the special form with Jordan blocks. These reslts can be obtained as the limit, performed accrately, of the eqations (6) for Riemann inariants when a soltion (, ) approaches the transition line. In the dnls eqation () case λ ± = ±, r ± = ± (5) and the transition line is gien by the eqation (x,t) =. Eqations for the Riemann inariants in terms of and are (± ) t = (± )(± ) x (6) or t ± t = x ± x ± x + x. (7) In the limit the two leading order terms / and gie the system t = x + x, t = x. (8) This system has matrix V gien by (), that is the Jordan form with λ =. In the db eqation (5) case λ ± = ± /, r ± = ± /, (9) and the transition line is gien by the eqation (x,t) =. So eqations (6) in terms of and are ( ± / ) t = ± / ( ± / ) x () or t ± / t = ± / x + x. ()

5 When x and t approaches the transition line (x,t) = in the leading orders and / one gets t =, t = x, () i.e. the Jordan form with λ =. In the general case the matrix V for the system (8) apparently is not of the Jordan block form on the transition line Ω (). It can be parameterized at C as ( ) λ+ BC B V V Ω= = C λ () BC where λ = λ + = λ = (A+D)/. Sch a matrix has the form ( ) V = λ +N () where N is the general nilpotent matrix. It is straightforward to check that there exists a two parameter family of inertible matrices P sch that ( ) PV P λ =. (5) λ The family P at C is gien by P = a b + BC C C BC b, (6) where a = a(,) and b = b(,) are two arbitrary fnctions. In the case of C = and B the matrix P becomes ( ) b P = a, (7) B where a = a(,) and b = b(,) are still two arbitrary fnctions. Ths is all non diagonal cases the matrix V is eqialent to a Jordan block, and the system (8) is eqialent to ( ) ( ) ( ) t λ x P = P. (8) λ t In or constrction the systems in the Jordanforms or the system (8) arise on the transition line only. In the paper [] the system (8) with matrix V gien by () on the whole plane (x,t) and its mlti-component analogs with Jordan blocks has been deried ia the conflence process for the Laricella-type fnctions associated with Grassmannians Gr(, 5) and Gr(, n). Let s consider the system (8) with the matrix V = V gien by (). It is parabolic on the plane (x,t). In this case there are ariables, sch the system (8) takes the form ( ) t = t ( λ λ )( x x x ). (9) which will be referred as the Jordan form. The relation between the Jordan ariables, and the original ones, can be fond soling the eqations ( ) ( ) ( ) ( ) t t = P, x x = P. () t t These eqations can be soled only if one finds the sitable integrating factors a and b which mst be chosen sch that ( ) BC d = abd+a + C b d, d = ad a BCd, () 5 x x

6 which fix a and b thanks to the compatibility conditions ( (ab) = a+ a ) BC b, C a = (a BC). () Therefore the Jordan ariables in parabolic systems play a role similar to the Riemann inariants for the standard diagonalizable case. Elliptic domain At the hyperbolic domain soltion of the system (8) describe wae motions. Properties of soltions of the system (8) in the elliptic domain are qite different. Their treatment as the fnction defining qasi-conformal mappings is one of the possible interpretations []. Indeed eqation (6) can be rewritten as [5] r z = +iλ iλ r z () where λ = λ + = (A+D+ Ω)/, r = r +, z = x+it and the oerline stands for complex conjgate. Sch eqations are known as Beltrami nonlinear eqations [6, 7]. A soltion of () defines a qasi-conformal mapping r : (z,z) (r,r) if the complex dilation µ = +iλ iλ obeys the condition µ = +iλ iλ <. () Using the explicit form of λ, in the region Ω < we hae µ = +i(a+d) Ω i(a+d)+ Ω. (5) It is easy to see that condition () is always satisfied when Ω <. So any soltion of the system (8) in the elliptic domain defines qasi-conformal mapping. At the transition line Ω = µ = +i(a+d) i(a+d) = (6) and the qasi conformal mappings degenerates. For instance it maps the nit circle in the plane (z,z) in the degenerate ellipsis in the plane (r,r) with ratio of major and minor axes going to infinity (see e.g. [6]). In the hodograph space, eqations () are eqialent to the linear Beltrami eqation r w = +i λ(w,w) i λ(w,w) r w (7) where w = + i and λ = A D+i Ω B. Similar to the calclations presented before, one shows that the condition () is always satisfied in the elliptic domain and so any soltion of the eqation (7) defines a qasi-conformal mapping (w,w) (r,r). On the transition line Ω = again +i λ = and qasi-conformal i λ mappings become singlar. So, both in the elliptic and hyperbolic domains soltions of the mixed system (8) exhibit particlar behaior when they approach the transition line Ω =. Namely, approaching the transition line from the hyperbolic side waes become nstable (see e.g. [8]) conerting into mess, goerned by parabolic eqations and transforming into qasi-conformal mappings dynamics beyond the transition line. Approaching the transition line from the elliptic side, the qasi-conformal mappings degenerate into singlar ones with µ = which maps the two dimensional domains in C into qasi one-dimensional ones. Beyond the transition line these qasi one-dimensional objects are transformed into moing waes. 6

7 5 Transition line and its crossing Since Riemann inariants are constants along characteristics (real or complex) the problem whether or not a transition line can be crossed is redced to the stdy of respectie properties of characteristics and transition lines (see e.g. [8, 9, ]). Comparison of the formlae for characteristics and transition line in the original ariables (t, x) and hodograph ariables (, ) (see e.g. formlae () and () clearly indicates that the latter ones is more appropriate for or prpose. The se of simple waes in [9] proides s another spport of sch obseration. In the hodograph space (, ) the characteristics and transition lines are gien by formla () and () respectiely. Let s begins with hyperbolic domain and let s assme that the deriaties inoled are bonded. Ths we hae two families of plane characteristic lines (ChL) in the hodograph space (B ) and a single transition line (TL) ChL : ( ) d = A D ± Ω, (8) d ± B TL : Ω(,) =. (9) The two simplest cases are: ) the two families (8) do not hae common points with (9) and ) they coincide at least on some interal. In the latter case, on the transition line one has the eqations which shold be eqialent to each other. Since on the transition line the necessary condition for this is gien by (if Ω ) d d = D A, Ω =, (5) B dω = Ω d+ω d =, (5) D A B + Ω Ω =. (5) Obiosly in both cases the transition from the hyperbolic domain to the elliptic one is impossible. Another simple case corresponds to the transersal intersections of ChL and TL. To derie the corresponding condition it is sfficient to consider these lines at points of intersection. Two characteristic toch each other and at the point on TL their tangents are (assming that both cres are smooth) Tangent to the TL at the same point is gien by (at Ω ) Characteristic and transition line cross transersally (with angle ) if i.e. d = D A d ChL B. (5) d = Ω TL. (5) d TL Ω d d, (55) d ChL d TL ( D A B + Ω ). (56) Ω TL Ths, if condition (56) is satisfied, the transition from the hyperbolic domain to the elliptic one is not forbidden. There are eight other possibilities. First for aregien by the figre and its reflections each oftwo cres with respect the straight line of common tangent at the point (, ). In these for cases characteristic lines 7

8 TL elliptic domain hyperbolic domain contact point ChL Figre : A configration of characteristic (solid) and transition line (dashed) in the (, )-space where the hyperbolic-elliptic transition in the contact point is forbidden. TL hyperbolic domain elliptic domain contact point ChL Figre : A configration of characteristic (solid) and transition line (dashed) in the (, )-space where the hyperbolic-elliptic transition in the contact point reqires a deeper analysis and it cold not be forbidden. In the elliptic region the characteristic line (here depicted as dotted) becomes complex. toches the transition line at the point (, ) and then trns back to the hyperbolic domain. So the transition is forbidden. At the point (, ) tangents of both side coincides and so ( ) d ( ) d ( d D A (, ) d = (, ) B + Ω ) =. (57) ChL TL Ω (, ) The fact of non-crossing is inariant nder the transformation of coordinates. Choosing the coordinates(, ) near to the point (, ) in sch a way that the axes = coincides with the common tangent, it is not difficlt to show in all for cases that the difference ( ) ( ) d d T (58) d ChL d TL changes sign passing the point (, ). Second for cases are gien by figre () and three others possible are obtained by reflection each of two cres with respect the straight line of common tangent at the point (, ). In these cases characteristic toches the transition line at the point (, ) and then pass into elliptic domain with characteristic speeds becoming complex. Ths, the transition is not forbidden. Ths in the eight cases considered aboe the behaior of T near the point of toch (, ) distingishes the cases of crossing and non-crossing. So we conclde that the transition from the hyperbolic to the elliptic domain is not possible if either T Ω= = on some interal of the transition line of T Ω(, )= = at 8

9 some point on TL and T changes sign passing from one side of the toch point (, ) to another. In particlar, the comparison of the condition (57) rewritten as ((D A)Ω +BΩ ) (, ) = (59) and the relation () shows that the eigenector of the matrix V corresponding to the doble eigenale λ is tangent to the transition line at the point (, ) in agreement with necessary condition of non-crossing (nonlinear stability) proposed in [8]. On the other hand if T Ω= or T Ω(, )= = and T does not change sign at the toching point (, ) the conditions of non-crossing are not satisfied and the transition from the hyperbolic domain to the elliptic one is not forbidden. In the analysis presented aboe it was assmed that all deriaties inclding d d are bonded. The cases of possible nbondedness reqire special consideration. To clarify the point let s consider the system ( ) ( )( ) t x = (6) f() t where f() = α and α = n+ or α = /(n+) with n =,,,,.... In this case the transition line is gien by = and the characteristics in the hodograph space are defined by the eqation So, the characteristics are gien by the lines x d d = α/. (6) ± α ( α)/ = = const, (6) with arbitrary andby the straightline =. The transition line = is then (degenerate) characteristics. The behaior of other characteristics is qite different for α > and < α <. In the case α = n+ and n, () has a singlarity at = and it it may toch the transition line only at the infinity. So the transition from hyperbolic to elliptic domain is not possible. For α = n+, the characteristics (6) toch the transition line at finite point (,) where d =. Since d, d ( ) α (α ) d = α ( ) = α α (6) the characteristics hae completely different behaior in the cases α = and α <. For α = one has d d (,),α= = const. (6) So characteristics approach smoothly the transition line and, hence, the transition is not possible. Note that at α = the system (6) represent the eqation (log) tt = xx. (65) In contrast, for α < (n =,,,...) the normal to the characteristic, i.e. elocity d d with which characteristic approaches (in normal direction) the transition line = grows to infinity (at the point (,)). Sch a behaior allows s to sggest that the characteristics may jmp across the line = and transition from the hyperbolic domain to the elliptic one ( < ) cold be possible. The nmerical reslts of the system (6) with f = / (α = /) presented in the paper [8] spport this obseration. Another indication that the system (6) with f = α, α < has has rather special properties is proided by eqations () and (), i.e. α t t +a α t = (66) and similar eqation for x. At the transition line = they are singlar eqations of parabolic type. Analysisofthesystemswhichhaepropertiessimilartothoseofsystem(6)with f = α, α < reqiresa separate stdy which will be performed elsewhere. Possibility of transition from elliptic to hyperbolic domain, corresponding conditions and associated qasi-conformal mapping are of interest too. These problems will be considered in a separate pblication. 9

10 6 Hamiltonian systems In this and next sections we will consider some concrete classes of eqations (8). We begin with the Hamiltonian systems which for two component case can always be locally pt in the form [9, ] ( ) ( )( ) t x h = (67) x t with Hamiltonian H = h(,)dx. So, the system (8) is ( ) ( )( ) t h h = x. (68) h h x In this case Ω = h h. Eqation () is t h h t h t h t +h t =, (69) while characteristics in the hodograph space (simple waes) are defined by the eqation ( ) d h =. (7) d ± h Transition line is gien by h h =, (7) assming that h h may change sign. In order to deal with generic non-diagonalizable case we defined the transition line as h =, h h=. (7) ( ) λ On the transition line the matrix V is eqialent to the Jordan block with λ = h λ h=. The transformation matrix P is gien by (7) with B = h h=. We note that at points where h = and h = on the transition line the matrix V is degenerated to a constant diagonal matrix. When it holds (7) one also has Ths in generic case with h h= we hae T = h h h=. (7) T (7) and transition from the hyperbolic domain h h > to the elliptic one is not forbidden. Both dnls and db eqations are Hamiltonian ones (see e.g. [, 9, ] ) with respectiely h dnls = +, (75) and h db = + 6. (76) Now let s stdy the behaior of characteristics in (,) plane (simple waes) near the point (, ) of contact with the transition line h = and h for general Hamiltonian systems (68). Expanding the right hand side of (7) near the point (, ) and assming that h and the deriaties enoled are finite, one obtains ( ) d d ± = ± a( )+b( ), with a = h h b = h h, (77)

11 where f f(, ). For infinitesimal δ = and δ = eqation (77) takes the form ( ) δ = aδ+bδ. (78) δ Hence δ = b b(δ) ± (δ) +a(δ). (79) So at δ, δ one has at the leading order h ± h ( ) /. (8) This formla gies s the niersal behaior of (, ) characteristics near the transition line for general Hamiltonian system (68) in the generic case h. The simplest and characteristic example of sch a behaior is proided by the db eqation for which h =. It shold be also noted that the fact, that for the general stationary plane motion of compressible gas (described by Chaplygin eqation) the behaior of characteristics near sonic line is gien by the formla (8), has been known for a long time (see e.g. [], 8). In particlar cases the behaior of characteristics near to the transition line is qite different. If at the transition point (, ) also h =, then, instead of eqation (78) one has ( ) δ = bδ +c(δ) +dδδ +f(δ) (8) δ where c,d and f are certain constants depending on h and its deriaties ealated at (, ) gien by b = h h, c = h h, d = h h h h (h ), f = h h h h (h ). (8) Soling this eqation and considering the limit of infinitesimal δ and δ one gets h h ( ). (8) In this case which changes sign at the point and, hence, transition is forbidden. In a similar way one can show that in the case when all deriaties k h k h= the behaior of near the transition line is of the type T ( ) (8) for k =,,5,...,n (85) ( ) n+. (86) So, for odd n, the transition is forbidden allowed while for odd n its is not. Finally if the transition line is gien by the eqation h = with h h= then one has the reslts presented aboe with exchange. For the dnls eqation the transition line = is a characteristic and d d =. Hence, the transition = is not possible. We remark that this case is in some sense degenerate becase the Hamiltonian is qadratic in and all the partial deriaties of k h/ n =, k are zero. For the db eqation, in contrast, the transition line is =, characteristic cross the transition line orthogonally d =, i.e. = T = h h =, (87) and conseqently the transition from the hyperbolic domain to elliptic one is not forbidden. d

12 7 Special classes of Hamiltonian systems: Gas dynamics eqations Expressions (75) and (76) sggest s to consider two special classes of systems with Hamiltonian densities h = F ()+, h = F ()+F () (88) where F,F,F are fnctions of a single ariable. For the Hamiltonian density of the form h one has the system ( ) ( )( ) t F x =. (89) t and Ω = F. Simple waes are defined by eqation F d = ±d. (9) The relations h = and h = imply T =. There are two qite different sitations. First corresponds to the case when he transition line is gien by = and F > near =. For the system of mixed type one has F = the characteristic in (, )-plane are gien by the eqation and so x n(n ) n, n =,,,..., (9) d d = ± n(n ) n (9) d d ( n). (9) Ths if n the acceleration d d dierges as the characteristic approaches the transition line = (which is also a characteristic) and, hence, the transition is not forbidden. If eerywhere the properties of the system (89) are qite different. Introdcing the fnction P() defined by the eqation F = P (), (9) one rewrites the system (89) as t = x + P x(ρ) ρ ρ t =(ρ) x., (95) It is the general isentropic one-dimensional gas-dynamic eqation with being elocity, being density ρ, P(ρ) is the pressre and t t (see e.g. [, ]). The TL-line is defined by Ω = ρf = P (ρ) = (96) and the characteristics in the space (,) are defined by the eqation dρ d = ± ρ P (ρ). (97) For ordinary media dp = c, (98) dρ S

13 i.e. sqared sond speed c in the medim. Ths, for normal cases the system is hyperbolic eerywhere. So the system (95) is of mixed type for particlar macroscopic systems for which the deriatie dp can anish S at some ale of density ρ (zero sond speed point) and change sign passing throgh this ale (see e.g. []-[5]). Sch a sitation is realized, for instance, for the fnctions P which for small ρ ρ are of the form a) : P (ρ ρ ) n+, b) : P (ρ ρ ) n+ n+, (99) where n =,±,±,.... Near the transition line ρ = ρ one has d ρ d (ρ ρ ) (+γ) () where γ = n+, or γ = n+ n+. In both cases d ρ d as ρ ρ and, hence, the transition is not forbidden. 8 Nonlinear wae type eqations dρ The system with Hamiltonian density h (88) is of the form ( ) ( t F = () t F () )( x x ). () This system is, in fact, the system of two conseration laws and it is eqialent to the single eqation t = (F ()) x, t = (F ()) x () (A(w t )) t = (B(w x )) x () where w x =, w t = F (), B(y) = F (y) and A(y) = (F ) (y). In the particlar case F () = / eqation () takes the form w tt = (F (w x)) x () or tt = (F ()) xx, (5) i.e. the standard form of the nonlinear wae eqations (see e.g. [] with F = P). For the db eqation For the system () one has and eqations for (, ) characteristics are gien by F ()d = ± F db = 6 + c. (6) Ω = F ()F () (7) F First we consider the case when the transition line is gien by ()d. (8) F ( ) =, F (), (9) which incldes the nonlinear wae case (5) where F () =. For the system of mixed type near the point the fnction F () shold be of the form F () const( ) n+, n =,,,.... ()

14 Figre : Simple wae strctre () for db eqation (). For (,) characteristics near the transition line = one has d d ( ) n. () So d d as and, hence,the characteristic approaches the transition line orthogonally. Ths the change of type is not forbidden. It is clearly so for nonlinear wae eqations with F () = n+, n =,,,.... Note that nonlinear wae eqations with F () = n, n =,,,... are hyperbolic eerywhere. The same is alid for the dispersionless Toda eqation tt = exp() xx for which F () = exp() (see e.g. []). If, instead of (9), the transition line is gien by one has the same reslts with the exchange. F (), F ( ) =, () 9 Nmerical example of transitions for the db eqation Here we present some nmerical reslts for the db eqation as the characteristic representatie of the generic class of Hamiltonian systems. Let s consider the class of periodic soltions with fixed bondary ales and with initial conditions ((x,) c) + =, (x,) = c+sin(x) () where c > is assmed in order to start in the hyperbolic sector. The simple waes for db eqation are (see figre ) ± / = k, k R. () For eery ale of c there are for simple waes tangent to the circle. The two lower ones satisfy the system ( ) / (k ± c) = c c, (k c ) / = c c. (5) where c is the ale of at contact point. For eery ale of c > the preios systems admits two soltions corresponding to two different simple waes symmetric with respect to the axis. Only aboe a minimm

15 Figre : In the example of initial conditions class () we shown that, in dependence of the parameter ales the hyperbolic-elliptic transition is forbidden (solid line) or allowed (dotted line). ale of c, called c crit, the cople of tangent simple waes intersects each other. Therefore only aboe c crit the hyperbolic elliptic transition is forbidden becase the initial condition (see e.g. [9]) is separated from the transition line by the two tangent simple waes. The critical ale c crit can be estimated as follows. At c crit the intersection of the tangent simple waes, becase of the problem symmetry, is at the origin =, =, i.e. with k = in (5). Soling therefore (5) with k = we obtain the ale of this critical constant (in case ()) which is c crit.77. In figre the hyperbolic elliptic transition possibility as the fnction of the parameter c of the family of initial conditions () is shown. The solid circle are the initial data at c = which is greaterthan the critical ale. In this case the initial data are bonded from below by two simple waes (two solid open cres) which are tangent to the initial data and intersect each other aboe the transition line =. This behaior preents the transition. The dashed circle (initial conditions with c = c crit ) are the critical conditions in the family which forbid the transition. Actally the two tangent simple waes (dashed open cres) intersect exactly at the transition line. Finally the dotted circle are initial data at c =. which is below the critical ale: the tangent simple waes (two dotted cres) hae no intersection and the initial data cold reach the transition line. In figres 5 (left) and 6 it is shown the eoltion of in dependence of (x,t) and respectiely in the non transition case of c = > c crit. In figres 5 (right) and 7 is shown the eoltion of in dependence of (x,t) and respectiely in the transition case of c =. < c crit. At the second to last step the cre is tangent to the transition line. Howeer this line is not a characteristic in the (,) space and the cre can cross the transition line as can be seen in the last plot. The nmerical eoltions in the figres 5, 6 and 7 are obtained sing Mathematica. Acknowledgments G.O. thanks R. Camassa for sefl discssions and informations. Partial spport by MIUR Cofin - project JJKPA is acknowledged. This work was carried ot nder the aspices of the GNFM Section of INdAM. 5

16 Hx,tL Hx,tL. x x t..5. t Figre 5: On the right plot we show the eoltion of with c = in the initial conditions () while on the left plot with c =.. The motion remains in the hyperbolic sector if c = > ccrit, while in the case c =. < ccrit we see the transition in the elliptic sector (pper-right region of the plot) Figre 6: Circle hyperbolic eoltion with c = > ccrit, t from zero to t = with eqispaced steps. The cre eoltion does not admits the transition below the intersection point. 6

17 Figre 7: Circle eoltion with c =. < c crit, t from zero to t = with eqispaced steps. The cre can cross the transition line =. References [] R. Corant, and D. Hilbert, Methods of mathematical physics. Vol. II: Partial differential eqations, Interscience Pblishers (John Wiley & Sons), New York-London (96) [] L. Bers, Mathematical aspects of sbsonic and transonic gas dynamics. Sreys in Applied Mathematics, Vol. John Wiley & Sons, Inc., New York; Chapman & Hall, Ltd., London (958) [] A. Bitsadze, V. Eqations of the mixed type. A Pergamon Press Book The Macmillan Co., New York (96) [] M. Schneider, Introdction to partial differential eqations of mixed type. Lectre Notes, No.. Institte for Mathematical Sciences, Uniersity of Delaware, Newark, Del., (977) [5] J. M. Rassias, Lectre notes on mixed type partial differential eqations. World Scientific Pblishing Co., Inc., Teaneck, NJ, (99) [6] G.-Q G. Chen, On Nonlinear Partial Differential Eqations of Mixed Type, Talk gien at UK-Japan Winter School on Nonlinear Analysis London 7- janary ( berndt/conferences/uk-japan/chen.pdf) [7] T. H. Otway Elliptic-Hyperbolic Partial Differential Eqations A Mini-Corse in Geometric and Qasilinear Methods, Springer (5) [8] P. Milewski, E. Tabak, C. Trner, R. Rosales and F. Menzaqe Nonlinear stability of two-layer flows. Commn. Math. Sci., no., 7- () [9] L. Chmakoa, F. E. Menzaqe, P. A. Milewski, R. R. Rosales, E. G. Tabak, and C. V. Trner Stability properties and nonlinear mappings of two and three-layer stratified flows, SAPM :-7 (9) 7

18 [] L. Chmakoa and E. Tabak Simple Waes Do Not Aoid Eigenale Crossings Comm. Pre Appl. Math. Vol. LXIII, 9- () [] G. B. Whitham, Linear and nonlinear waes. Pre and Applied Mathematics. Wiley-Interscience-John Wiley & Sons (97) [] L. D. Landa and E. M. Lifshitz Flid mechanics. Translated from the Rssian by J. B. Sykes and W. H. Reid. Corse of Theoretical Physics, Vol. 6 Pergamon Press, London-Paris-Frankfrt; Addison-Wesley Pblishing Co., Inc., Reading, Mass. 959 [] Y. Kodama, B. G. Konopelchenko, Conflence of hypergeometric fnctions and integrable hydrodynamic type systems arxi:5.5 5 [] M. Larente, A general problem of the theory of qasi-conformal representation of plane regions, (Rssian) Mat. Sbornik N.S. :85- (97). [5] B. G. Konopelchenko, G. Ortenzi, Qasi-classical approximation in ortex filament dynamics. Integrable systems, gradient catastrophe, and fltter Std. Appl. Math. (), no., [6] L. V Ahlfors Lectres on qasiconformal mappings, in Van Nostrand Mathematical Stdies, No. D. Van Nostrand Co., Inc., Toronto, Ontario, New York, London, 966. [7] B. Bojarski, Qasiconformal mappings and general strctral properties of systems of non linear eqations elliptic in the sense of Larente, Symp. Math. 8:85-99 (976). [8] A. D. D. Craik, Wae interactions and flid flows. Cambridge Monographs on Mechanics and Applied Mathematics. Cambridge Uniersity Press, Cambridge, (988) [9] B. A. Dbroin and S. P. Noiko Poisson brackets of hydrodynamic type. (Rssian) Dokl. Akad. Nak SSSR 79, no., 9-97 (98) (English Translation: Soiet Math. Dokl., no., (98) ) [] B. A. Dbroin On niersality of critical behaior in Hamiltonian PDEs. Geometry, topology, and mathematical physics, 59-9, Amer. Math. Soc. Transl. Ser.,, Amer. Math. Soc., Proidence, RI, (8) [] S. M. Carroll, M. Hoffman, and M. Trodden Can the dark energy eqation-of-state parameter w be less than? Phys. Re. D 68, 59, () [] L. P. Chimento Extended tachyon field, Chaplygin gas, and solable k-essence cosmologies. (English smmary) Phys. Re. D () 69, no., 57 () [] P. Creminelli, G. D Amico, J. Norẽa, L. Senatore and F. Vernizzi, Spherical collapse in qintessence models with zero speed of sond Jornal of Cosmology and Astroparticle Physics Volme, Isse, Article nmber 7 () [] L.R. Gómez, A.M. Trner, M. Van Hecke and V. Vitelli, Shocks near jamming Phys. Re. Lett. Volme 8, Isse 5, Article nmber 58 () [5] O.Longo and H.Qeedo, A nified dark energy model from a anishing speed of sond with emergent cosmological constant, Int. J. Mod. Phys. D, 5 () 8

Reduction of over-determined systems of differential equations

Reduction of over-determined systems of differential equations Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical

More information

Direct linearization method for nonlinear PDE s and the related kernel RBFs

Direct linearization method for nonlinear PDE s and the related kernel RBFs Direct linearization method for nonlinear PDE s and the related kernel BFs W. Chen Department of Informatics, Uniersity of Oslo, P.O.Box 1080, Blindern, 0316 Oslo, Norway Email: wenc@ifi.io.no Abstract

More information

1 The space of linear transformations from R n to R m :

1 The space of linear transformations from R n to R m : Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations

More information

The Brauer Manin obstruction

The Brauer Manin obstruction The Braer Manin obstrction Martin Bright 17 April 2008 1 Definitions Let X be a smooth, geometrically irredcible ariety oer a field k. Recall that the defining property of an Azmaya algebra A is that,

More information

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems

A scalar nonlocal bifurcation of solitary waves for coupled nonlinear Schrödinger systems INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 5 (22) 265 292 NONLINEARITY PII: S95-775(2)349-4 A scalar nonlocal bifrcation of solitary waes for copled nonlinear Schrödinger systems Alan R Champneys and

More information

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently

Relativity II. The laws of physics are identical in all inertial frames of reference. equivalently Relatiity II I. Henri Poincare's Relatiity Principle In the late 1800's, Henri Poincare proposed that the principle of Galilean relatiity be expanded to inclde all physical phenomena and not jst mechanics.

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

We automate the bivariate change-of-variables technique for bivariate continuous random variables with

We automate the bivariate change-of-variables technique for bivariate continuous random variables with INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,

More information

Formation and Transition of Delta Shock in the Limits of Riemann Solutions to the Perturbed Chromatography Equations

Formation and Transition of Delta Shock in the Limits of Riemann Solutions to the Perturbed Chromatography Equations Jornal of Mathematical Research with Applications Jan., 29, Vol.39, No., pp.6 74 DOI:.377/j.issn:295-265.29..7 Http://jmre.dlt.ed.cn Formation and Transition of Delta Shock in the Limits of Riemann Soltions

More information

ON THE PERFORMANCE OF LOW

ON THE PERFORMANCE OF LOW Monografías Matemáticas García de Galdeano, 77 86 (6) ON THE PERFORMANCE OF LOW STORAGE ADDITIVE RUNGE-KUTTA METHODS Inmaclada Higeras and Teo Roldán Abstract. Gien a differential system that inoles terms

More information

Pulses on a Struck String

Pulses on a Struck String 8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a

More information

Numerical Model for Studying Cloud Formation Processes in the Tropics

Numerical Model for Studying Cloud Formation Processes in the Tropics Astralian Jornal of Basic and Applied Sciences, 5(2): 189-193, 211 ISSN 1991-8178 Nmerical Model for Stdying Clod Formation Processes in the Tropics Chantawan Noisri, Dsadee Skawat Department of Mathematics

More information

The Linear Quadratic Regulator

The Linear Quadratic Regulator 10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.

More information

EE2 Mathematics : Functions of Multiple Variables

EE2 Mathematics : Functions of Multiple Variables EE2 Mathematics : Fnctions of Mltiple Variables http://www2.imperial.ac.k/ nsjones These notes are not identical word-for-word with m lectres which will be gien on the blackboard. Some of these notes ma

More information

arxiv: v1 [math.co] 25 Sep 2016

arxiv: v1 [math.co] 25 Sep 2016 arxi:1609.077891 [math.co] 25 Sep 2016 Total domination polynomial of graphs from primary sbgraphs Saeid Alikhani and Nasrin Jafari September 27, 2016 Department of Mathematics, Yazd Uniersity, 89195-741,

More information

Modelling, Simulation and Control of Quadruple Tank Process

Modelling, Simulation and Control of Quadruple Tank Process Modelling, Simlation and Control of Qadrple Tan Process Seran Özan, Tolgay Kara and Mehmet rıcı,, Electrical and electronics Engineering Department, Gaziantep Uniersity, Gaziantep, Trey bstract Simple

More information

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qalitatie featres of a discrete dynamical system of homogeneos difference

More information

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer 6. Similarit Soltions Chapter 6 Momentm Transfer in an Eternal Laminar Bondar Laer Consider a laminar incompressible bondar laer with constant properties. Assme the flow is stead and two-dimensional aligned

More information

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary

More information

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad

Linear Strain Triangle and other types of 2D elements. By S. Ziaei Rad Linear Strain Triangle and other tpes o D elements B S. Ziaei Rad Linear Strain Triangle (LST or T6 This element is also called qadratic trianglar element. Qadratic Trianglar Element Linear Strain Triangle

More information

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE

THE DIFFERENTIAL GEOMETRY OF REGULAR CURVES ON A REGULAR TIME-LIKE SURFACE Dynamic Systems and Applications 24 2015 349-360 TH DIFFRNTIAL OMTRY OF RULAR CURVS ON A RULAR TIM-LIK SURFAC MIN OZYILMAZ AND YUSUF YAYLI Department of Mathematics ge Uniersity Bornoa Izmir 35100 Trkey

More information

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps

Minimizing Intra-Edge Crossings in Wiring Diagrams and Public Transportation Maps Minimizing Intra-Edge Crossings in Wiring Diagrams and Pblic Transportation Maps Marc Benkert 1, Martin Nöllenbrg 1, Takeaki Uno 2, and Alexander Wolff 1 1 Department of Compter Science, Karlsrhe Uniersity,

More information

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2

Vectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2 MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,

More information

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.

Krauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields. Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore

More information

Existence of HCMU metrics on higher genus Riemann surfaces

Existence of HCMU metrics on higher genus Riemann surfaces Existence of HCMU metrics on higher gens Riemann srfaces October 4, 0 bstract We will generalize the reslt in [CCW05] and roe the existence of HCMU metrics on higher gens K-srfaces, i.e. Riemann srfaces

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1

u P(t) = P(x,y) r v t=0 4/4/2006 Motion ( F.Robilliard) 1 y g j P(t) P(,y) r t0 i 4/4/006 Motion ( F.Robilliard) 1 Motion: We stdy in detail three cases of motion: 1. Motion in one dimension with constant acceleration niform linear motion.. Motion in two dimensions

More information

GRAY CODES FAULTING MATCHINGS

GRAY CODES FAULTING MATCHINGS Uniersity of Ljbljana Institte of Mathematics, Physics and Mechanics Department of Mathematics Jadranska 19, 1000 Ljbljana, Sloenia Preprint series, Vol. 45 (2007), 1036 GRAY CODES FAULTING MATCHINGS Darko

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

The Open Civil Engineering Journal

The Open Civil Engineering Journal Send Orders for Reprints to reprints@benthamscience.ae 564 The Open Ciil Engineering Jornal, 16, 1, 564-57 The Open Ciil Engineering Jornal Content list aailable at: www.benthamopen.com/tociej/ DOI: 1.174/187414951611564

More information

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany

Subcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2

More information

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007

Graphs and Networks Lecture 5. PageRank. Lecturer: Daniel A. Spielman September 20, 2007 Graphs and Networks Lectre 5 PageRank Lectrer: Daniel A. Spielman September 20, 2007 5.1 Intro to PageRank PageRank, the algorithm reportedly sed by Google, assigns a nmerical rank to eery web page. More

More information

Estimations of the Influence of the Non-Linearity of the Aerodynamic Coefficients on the Skewness of the Loading. Vincent Denoël *, 1), Hervé Degée 1)

Estimations of the Influence of the Non-Linearity of the Aerodynamic Coefficients on the Skewness of the Loading. Vincent Denoël *, 1), Hervé Degée 1) Estimations of the Inflence of the Non-Linearity of the Aerodynamic oefficients on the Skewness of the Loading Vincent enoël *, 1), Heré egée 1) 1) epartment of Material mechanics and Strctres, Uniersity

More information

MATH2715: Statistical Methods

MATH2715: Statistical Methods MATH275: Statistical Methods Exercises III (based on lectres 5-6, work week 4, hand in lectre Mon 23 Oct) ALL qestions cont towards the continos assessment for this modle. Q. If X has a niform distribtion

More information

Study on the Mathematic Model of Product Modular System Orienting the Modular Design

Study on the Mathematic Model of Product Modular System Orienting the Modular Design Natre and Science, 2(, 2004, Zhong, et al, Stdy on the Mathematic Model Stdy on the Mathematic Model of Prodct Modlar Orienting the Modlar Design Shisheng Zhong 1, Jiang Li 1, Jin Li 2, Lin Lin 1 (1. College

More information

Sources of Non Stationarity in the Semivariogram

Sources of Non Stationarity in the Semivariogram Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary

More information

Spring, 2008 CIS 610. Advanced Geometric Methods in Computer Science Jean Gallier Homework 1, Corrected Version

Spring, 2008 CIS 610. Advanced Geometric Methods in Computer Science Jean Gallier Homework 1, Corrected Version Spring, 008 CIS 610 Adanced Geometric Methods in Compter Science Jean Gallier Homework 1, Corrected Version Febrary 18, 008; De March 5, 008 A problems are for practice only, and shold not be trned in.

More information

Modelling by Differential Equations from Properties of Phenomenon to its Investigation

Modelling by Differential Equations from Properties of Phenomenon to its Investigation Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University

More information

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica

Visualisations of Gussian and Mean Curvatures by Using Mathematica and webmathematica Visalisations of Gssian and Mean Cratres by Using Mathematica and webmathematica Vladimir Benić, B. sc., (benic@grad.hr), Sonja Gorjanc, Ph. D., (sgorjanc@grad.hr) Faclty of Ciil Engineering, Kačićea 6,

More information

Restricted cycle factors and arc-decompositions of digraphs. J. Bang-Jensen and C. J. Casselgren

Restricted cycle factors and arc-decompositions of digraphs. J. Bang-Jensen and C. J. Casselgren Restricted cycle factors and arc-decompositions of digraphs J. Bang-Jensen and C. J. Casselgren REPORT No. 0, 0/04, spring ISSN 0-467X ISRN IML-R- -0-/4- -SE+spring Restricted cycle factors and arc-decompositions

More information

3.3 Operations With Vectors, Linear Combinations

3.3 Operations With Vectors, Linear Combinations Operations With Vectors, Linear Combinations Performance Criteria: (d) Mltiply ectors by scalars and add ectors, algebraically Find linear combinations of ectors algebraically (e) Illstrate the parallelogram

More information

Change of Variables. (f T) JT. f = U

Change of Variables. (f T) JT. f = U Change of Variables 4-5-8 The change of ariables formla for mltiple integrals is like -sbstittion for single-ariable integrals. I ll gie the general change of ariables formla first, and consider specific

More information

Flood flow at the confluence of compound river channels

Flood flow at the confluence of compound river channels Rier Basin Management VIII 37 Flood flow at the conflence of compond rier channels T. Ishikawa 1, R. Akoh 1 & N. Arai 2 1 Department of Enironmental Science and Technology, Tokyo Institte of Technology,

More information

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS

FREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS 7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,

More information

Shock wave structure for Generalized Burnett Equations

Shock wave structure for Generalized Burnett Equations Shock wave strctre for Generalized Brnett Eqations A.V. Bobylev, M. Bisi, M.P. Cassinari, G. Spiga Dept. of Mathematics, Karlstad University, SE-65 88 Karlstad, Sweden, aleander.bobylev@ka.se Dip. di Matematica,

More information

Analysis of Enthalpy Approximation for Compressed Liquid Water

Analysis of Enthalpy Approximation for Compressed Liquid Water Analysis of Entalpy Approximation for Compressed Liqid Water Milioje M. Kostic e-mail: kostic@ni.ed Nortern Illinois Uniersity, DeKalb, IL 60115-2854 It is cstom to approximate solid and liqid termodynamic

More information

Differential Geometry. Peter Petersen

Differential Geometry. Peter Petersen Differential Geometry Peter Petersen CHAPTER Preliminaries.. Vectors-Matrices Gien a basis e, f for a two dimensional ector space we expand ectors sing matrix mltiplication e e + f f e f apple e f and

More information

Turbulence and boundary layers

Turbulence and boundary layers Trblence and bondary layers Weather and trblence Big whorls hae little whorls which feed on the elocity; and little whorls hae lesser whorls and so on to iscosity Lewis Fry Richardson Momentm eqations

More information

Concept of Stress at a Point

Concept of Stress at a Point Washkeic College of Engineering Section : STRONG FORMULATION Concept of Stress at a Point Consider a point ithin an arbitraril loaded deformable bod Define Normal Stress Shear Stress lim A Fn A lim A FS

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2

Lecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE

More information

DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY

DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS ( < 20) AND CONSEQUENCES ON CRITICALITY SAFETY N. Leclaire, S. Evo, I.R.S.N., France Introdction In criticality stdies, the blk density

More information

Baroclinic Buoyancy-Inertia Joint Stability Parameter

Baroclinic Buoyancy-Inertia Joint Stability Parameter Jornal of Oceanography, Vol. 6, pp. 35 to 46, 5 Baroclinic Boyancy-Inertia Joint Stability Parameter HIDEO KAWAI* 3-8 Shibagahara, Kse, Joyo, Kyoto Pref. 6-, Japan (Receied 9 September 3; in reised form

More information

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS

NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS NUCLEATION AND SPINODAL DECOMPOSITION IN TERNARY-COMPONENT ALLOYS COLLEEN ACKERMANN AND WILL HARDESTY Abstract. The Cahn-Morral System has often been sed to model the dynamics of phase separation in mlti-component

More information

ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES

ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES ON TRANSIENT DYNAMICS, OFF-EQUILIBRIUM BEHAVIOUR AND IDENTIFICATION IN BLENDED MULTIPLE MODEL STRUCTURES Roderick Mrray-Smith Dept. of Compting Science, Glasgow Uniersity, Glasgow, Scotland. rod@dcs.gla.ac.k

More information

Group and Renormgroup Symmetry of Quasi-Chaplygin Media

Group and Renormgroup Symmetry of Quasi-Chaplygin Media Nonlinear Mathematical Physics 1996, V.3, N 3 4, 351 356. Group and Renormgroup Symmetry of Quasi-Chaplygin Media Vladimir F. KOVALEV Institute for Mathematical Modelling, Miusskaya 4-A, Moscow, 125047,

More information

Second-Order Wave Equation

Second-Order Wave Equation Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order

More information

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students

BLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they

More information

A Model-Free Adaptive Control of Pulsed GTAW

A Model-Free Adaptive Control of Pulsed GTAW A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,

More information

A Single Species in One Spatial Dimension

A Single Species in One Spatial Dimension Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,

More information

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.

Classify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled. Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

Gravitational Instability of a Nonrotating Galaxy *

Gravitational Instability of a Nonrotating Galaxy * SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy

More information

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use.

Change of Variables. f(x, y) da = (1) If the transformation T hasn t already been given, come up with the transformation to use. MATH 2Q Spring 26 Daid Nichols Change of Variables Change of ariables in mltiple integrals is complicated, bt it can be broken down into steps as follows. The starting point is a doble integral in & y.

More information

Five Basic Concepts of Axiomatic Rewriting Theory

Five Basic Concepts of Axiomatic Rewriting Theory Fie Basic Concepts of Axiomatic Rewriting Theory Pal-André Melliès Institt de Recherche en Informatiqe Fondamentale (IRIF) CNRS, Uniersité Paris Diderot Abstract In this inited talk, I will reiew fie basic

More information

PHASE STEERING AND FOCUSING BEHAVIOR OF ULTRASOUND IN CEMENTITIOUS MATERIALS

PHASE STEERING AND FOCUSING BEHAVIOR OF ULTRASOUND IN CEMENTITIOUS MATERIALS PHAS STRING AND FOCUSING BHAVIOR OF ULTRASOUND IN CMNTITIOUS MATRIALS Shi-Chang Wooh and Lawrence Azar Department of Civil and nvironmental ngineering Massachsetts Institte of Technology Cambridge, MA

More information

Remarks on strongly convex stochastic processes

Remarks on strongly convex stochastic processes Aeqat. Math. 86 (01), 91 98 c The Athor(s) 01. This article is pblished with open access at Springerlink.com 0001-9054/1/010091-8 pblished online November 7, 01 DOI 10.1007/s00010-01-016-9 Aeqationes Mathematicae

More information

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector Chapter 4: Flids Kinematics 4. Velocit and Description Methods Primar dependent ariable is flid elocit ector V V ( r ); where r is the position ector If V is known then pressre and forces can be determined

More information

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk

Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Porous Rotating Disk Effects of MHD Laminar Flow Between a Fixed Impermeable Disk and a Poros otating Disk Hemant Poonia * * Asstt. Prof., Deptt. of Math, Stat & Physics, CCSHAU, Hisar-54.. C. Chadhary etd. Professor, Deptt.

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method

Evaluation of the Fiberglass-Reinforced Plastics Interfacial Behavior by using Ultrasonic Wave Propagation Method 17th World Conference on Nondestrctive Testing, 5-8 Oct 008, Shanghai, China Evalation of the Fiberglass-Reinforced Plastics Interfacial Behavior by sing Ultrasonic Wave Propagation Method Jnjie CHANG

More information

Low-emittance tuning of storage rings using normal mode beam position monitor calibration

Low-emittance tuning of storage rings using normal mode beam position monitor calibration PHYSIAL REVIEW SPEIAL TOPIS - AELERATORS AND BEAMS 4, 784 () Low-emittance tning of storage rings sing normal mode beam position monitor calibration A. Wolski* Uniersity of Lierpool, Lierpool, United Kingdom

More information

Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen

Reaction-Diusion Systems with. 1-Homogeneous Non-linearity. Matthias Buger. Mathematisches Institut der Justus-Liebig-Universitat Gieen Reaction-Dision Systems ith 1-Homogeneos Non-linearity Matthias Bger Mathematisches Institt der Jsts-Liebig-Uniersitat Gieen Arndtstrae 2, D-35392 Gieen, Germany Abstract We describe the dynamics of a

More information

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL

FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL RESPONSE OF A CANTILEVER MODEL BBAA VI International Colloqim on: Blff Bodies Aerodynamics & Applications Milano, Italy, Jly, 2-24 28 FLUCTUATING WIND VELOCITY CHARACTERISTICS OF THE WAKE OF A CONICAL HILL THAT CAUSE LARGE HORIZONTAL

More information

Affine Invariant Total Variation Models

Affine Invariant Total Variation Models Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,

More information

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63

Kragujevac J. Sci. 34 (2012) UDC 532.5: :537.63 5 Kragjevac J. Sci. 34 () 5-. UDC 53.5: 536.4:537.63 UNSTEADY MHD FLOW AND HEAT TRANSFER BETWEEN PARALLEL POROUS PLATES WITH EXPONENTIAL DECAYING PRESSURE GRADIENT Hazem A. Attia and Mostafa A. M. Abdeen

More information

Homotopy Perturbation Method for Solving Linear Boundary Value Problems

Homotopy Perturbation Method for Solving Linear Boundary Value Problems International Jornal of Crrent Engineering and Technolog E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rights Reserved Available at http://inpressco.com/categor/ijcet Research Article Homotop

More information

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance

Efficiency Increase and Input Power Decrease of Converted Prototype Pump Performance International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power

More information

A NUMERICAL STUDY OF A PATHOLOGICAL EXAMPLE OF p-system

A NUMERICAL STUDY OF A PATHOLOGICAL EXAMPLE OF p-system SIAM J. NUMER. ANAL. Vol. 36, No. 5, pp. 1588 163 c 1999 Society for Indstrial and Applied Mathematics A NUMERICAL STUDY OF A PATHOLOGICAL EXAMPLE OF p-system PHILIPPE HOCH AND MICHEL RASCLE Abstract.

More information

The Heat Equation and the Li-Yau Harnack Inequality

The Heat Equation and the Li-Yau Harnack Inequality The Heat Eqation and the Li-Ya Harnack Ineqality Blake Hartley VIGRE Research Paper Abstract In this paper, we develop the necessary mathematics for nderstanding the Li-Ya Harnack ineqality. We begin with

More information

MAT389 Fall 2016, Problem Set 6

MAT389 Fall 2016, Problem Set 6 MAT389 Fall 016, Problem Set 6 Trigonometric and hperbolic fnctions 6.1 Show that e iz = cos z + i sin z for eer comple nmber z. Hint: start from the right-hand side and work or wa towards the left-hand

More information

Two-media boundary layer on a flat plate

Two-media boundary layer on a flat plate Two-media bondary layer on a flat plate Nikolay Ilyich Klyev, Asgat Gatyatovich Gimadiev, Yriy Alekseevich Krykov Samara State University, Samara,, Rssia Samara State Aerospace University named after academician

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model

Computational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi

More information

Linear System Theory (Fall 2011): Homework 1. Solutions

Linear System Theory (Fall 2011): Homework 1. Solutions Linear System Theory (Fall 20): Homework Soltions De Sep. 29, 20 Exercise (C.T. Chen: Ex.3-8). Consider a linear system with inpt and otpt y. Three experiments are performed on this system sing the inpts

More information

Self-induced stochastic resonance in excitable systems

Self-induced stochastic resonance in excitable systems Self-indced stochastic resonance in excitable systems Cyrill B. Mrato Department of Mathematical Sciences, New Jersey Institte of Technology, Newark, NJ 7 Eric Vanden-Eijnden Corant Institte of Mathematical

More information

4.2 First-Order Logic

4.2 First-Order Logic 64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore

More information

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane

Conditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching

More information

IN this paper we consider simple, finite, connected and

IN this paper we consider simple, finite, connected and INTERNATIONAL JOURNAL OF MATHEMATICS AND SCIENTIFIC COMPUTING (ISSN: -5), VOL., NO., -Eqitable Labeling for Some Star and Bistar Related Graphs S.K. Vaidya and N.H. Shah Abstract In this paper we proe

More information

Minimum-Latency Beaconing Schedule in Multihop Wireless Networks

Minimum-Latency Beaconing Schedule in Multihop Wireless Networks This fll text paper was peer reiewed at the direction of IEEE Commnications Society sbject matter experts for pblication in the IEEE INFOCOM 009 proceedings Minimm-Latency Beaconing Schedle in Mltihop

More information

Math 4A03: Practice problems on Multivariable Calculus

Math 4A03: Practice problems on Multivariable Calculus Mat 4A0: Practice problems on Mltiariable Calcls Problem Consider te mapping f, ) : R R defined by fx, y) e y + x, e x y) x, y) R a) Is it possible to express x, y) as a differentiable fnction of, ) near

More information

Chaotic and Hyperchaotic Complex Jerk Equations

Chaotic and Hyperchaotic Complex Jerk Equations International Jornal of Modern Nonlinear Theory and Application 0 6-3 http://dxdoiorg/0436/ijmnta000 Pblished Online March 0 (http://wwwscirporg/jornal/ijmnta) Chaotic and Hyperchaotic Complex Jerk Eqations

More information

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE

TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE TRANSONIC EVAPORATION WAVES IN A SPHERICALLY SYMMETRIC NOZZLE XIAOBIAO LIN AND MARTIN WECHSELBERGER Abstract. This paper stdies the liqid to vapor phase transition in a cone shaped nozzle. Using the geometric

More information

entropy ISSN by MDPI

entropy ISSN by MDPI Entropy, 007, 9, 113-117 Letter to the Editor entropy ISSN 1099-4300 007 by MDPI www.mdpi.org/entropy Entropy of Relativistic Mono-Atomic Gas and Temperatre Relativistic Transformation in Thermodynamics

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

Palindromes and local periodicity

Palindromes and local periodicity Palindromes and local periodicity A. Blondin Massé, S. Brlek, A. Garon, S. Labbé Laboratoire de Combinatoire et d Informatiqe Mathématiqe, Uniersité d Qébec à Montréal, C. P. 8888 Sccrsale Centre-Ville,

More information

Exercise 4. An optional time which is not a stopping time

Exercise 4. An optional time which is not a stopping time M5MF6, EXERCICE SET 1 We shall here consider a gien filtered probability space Ω, F, P, spporting a standard rownian motion W t t, with natral filtration F t t. Exercise 1 Proe Proposition 1.1.3, Theorem

More information

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms

Spanning Trees with Many Leaves in Graphs without Diamonds and Blossoms Spanning Trees ith Many Leaes in Graphs ithot Diamonds and Blossoms Pal Bonsma Florian Zickfeld Technische Uniersität Berlin, Fachbereich Mathematik Str. des 7. Jni 36, 0623 Berlin, Germany {bonsma,zickfeld}@math.t-berlin.de

More information

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications

Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that

More information