Hyperbolic partial differential equation








In mathematics, a hyperbolic partial differential equation of order n is a partial differential equation (PDE) that, roughly speaking, has a well-posed initial value problem for the first n − 1 derivatives. More precisely, the Cauchy problem can be locally solved for arbitrary initial data along any non-characteristic hypersurface. Many of the equations of mechanics are hyperbolic, and so the study of hyperbolic equations is of substantial contemporary interest. The model hyperbolic equation is the wave equation. In one spatial dimension, this is


∂2u∂t2=c2∂2u∂x2{displaystyle {frac {partial ^{2}u}{partial t^{2}}}=c^{2}{frac {partial ^{2}u}{partial x^{2}}}}{displaystyle {frac {partial ^{2}u}{partial t^{2}}}=c^{2}{frac {partial ^{2}u}{partial x^{2}}}}

The equation has the property that, if u and its first time derivative are arbitrarily specified initial data on the line t = 0 (with sufficient smoothness properties), then there exists a solution for all time t.


The solutions of hyperbolic equations are "wave-like". If a disturbance is made in the initial data of a hyperbolic differential equation, then not every point of space feels the disturbance at once. Relative to a fixed time coordinate, disturbances have a finite propagation speed. They travel along the characteristics of the equation. This feature qualitatively distinguishes hyperbolic equations from elliptic partial differential equations and parabolic partial differential equations. A perturbation of the initial (or boundary) data of an elliptic or parabolic equation is felt at once by essentially all points in the domain.


Although the definition of hyperbolicity is fundamentally a qualitative one, there are precise criteria that depend on the particular kind of differential equation under consideration. There is a well-developed theory for linear differential operators, due to Lars Gårding, in the context of microlocal analysis. Nonlinear differential equations are hyperbolic if their linearizations are hyperbolic in the sense of Gårding. There is a somewhat different theory for first order systems of equations coming from systems of conservation laws.




Contents






  • 1 Definition


  • 2 Examples


  • 3 Hyperbolic system of partial differential equations


  • 4 Hyperbolic system and conservation laws


  • 5 See also


  • 6 Notes


  • 7 Bibliography


  • 8 External links





Definition


A partial differential equation is hyperbolic at a point P provided that the Cauchy problem is uniquely solvable in a neighborhood of P for any initial data given on a non-characteristic hypersurface passing through P.[1] Here the prescribed initial data consist of all (transverse) derivatives of the function on the surface up to one less than the order of the differential equation.



Examples


By a linear change of variables, any equation of the form


A∂2u∂x2+2B∂2u∂x∂y+C∂2u∂y2+(lower order derivative terms)=0{displaystyle A{frac {partial ^{2}u}{partial x^{2}}}+2B{frac {partial ^{2}u}{partial xpartial y}}+C{frac {partial ^{2}u}{partial y^{2}}}+{text{(lower order derivative terms)}}=0}{displaystyle A{frac {partial ^{2}u}{partial x^{2}}}+2B{frac {partial ^{2}u}{partial xpartial y}}+C{frac {partial ^{2}u}{partial y^{2}}}+{text{(lower order derivative terms)}}=0}

with


B2−AC>0{displaystyle B^{2}-AC>0}{displaystyle B^{2}-AC>0}

can be transformed to the wave equation, apart from lower order terms which are inessential for the qualitative understanding of the equation.[2] This definition is analogous to the definition of a planar hyperbola.


The one-dimensional wave equation:


∂2u∂t2−c2∂2u∂x2=0{displaystyle {frac {partial ^{2}u}{partial t^{2}}}-c^{2}{frac {partial ^{2}u}{partial x^{2}}}=0}frac{partial^2 u}{partial t^2} - c^2frac{partial^2 u}{partial x^2} = 0

is an example of a hyperbolic equation. The two-dimensional and three-dimensional wave equations also fall into the category of hyperbolic PDE. This type of second-order hyperbolic partial differential equation may be transformed to a hyperbolic system of first-order differential equations.[3]



Hyperbolic system of partial differential equations


The following is a system of s{displaystyle s}s first order partial differential equations for s{displaystyle s}s unknown functions u→=(u1,…,us){displaystyle {vec {u}}=(u_{1},ldots ,u_{s})} vec u = (u_1, ldots, u_s) , u→=u→(x→,t){displaystyle {vec {u}}={vec {u}}({vec {x}},t)} vec u =vec u (vec x,t), where x→∈Rd{displaystyle {vec {x}}in mathbb {R} ^{d}}vec x in mathbb{R}^d:


(∗)∂u→∂t+∑j=1d∂∂xjfj→(u→)=0,{displaystyle (*)quad {frac {partial {vec {u}}}{partial t}}+sum _{j=1}^{d}{frac {partial }{partial x_{j}}}{vec {f^{j}}}({vec {u}})=0,}(*) quad frac{partial vec u}{partial t}<br /> + sum_{j=1}^d frac{partial}{partial x_j}<br /> vec {f^j} (vec u) = 0,<br />

where fj→∈C1(Rs,Rs),j=1,…,d{displaystyle {vec {f^{j}}}in C^{1}(mathbb {R} ^{s},mathbb {R} ^{s}),j=1,ldots ,d}vec {f^j} in C^1(mathbb{R}^s, mathbb{R}^s), j = 1, ldots, d are once continuously differentiable functions, nonlinear in general.


Next, for each fj→{displaystyle {vec {f^{j}}}}vec {f^j} define the s×s{displaystyle stimes s}s times s Jacobian matrix


Aj:=(∂f1j∂u1⋯∂f1j∂us⋮⋱⋮∂fsj∂u1⋯∂fsj∂us), for j=1,…,d.{displaystyle A^{j}:={begin{pmatrix}{frac {partial f_{1}^{j}}{partial u_{1}}}&cdots &{frac {partial f_{1}^{j}}{partial u_{s}}}\vdots &ddots &vdots \{frac {partial f_{s}^{j}}{partial u_{1}}}&cdots &{frac {partial f_{s}^{j}}{partial u_{s}}}end{pmatrix}},{text{ for }}j=1,ldots ,d.}A^j:=<br /> begin{pmatrix} frac{partial f_1^j}{partial u_1} & cdots & frac{partial f_1^j}{partial u_s} \<br /> vdots & ddots & vdots \<br /> frac{partial f_s^j}{partial u_1} & cdots &<br /> frac{partial f_s^j}{partial u_s}<br /> end{pmatrix}<br /> ,text{ for }j = 1, ldots, d.

The system (∗){displaystyle (*)}(*) is hyperbolic if for all α1,…,αd∈R{displaystyle alpha _{1},ldots ,alpha _{d}in mathbb {R} }alpha_1, ldots, alpha_d in mathbb{R} the matrix A:=α1A1+⋯+αdAd{displaystyle A:=alpha _{1}A^{1}+cdots +alpha _{d}A^{d}}A := alpha_1 A^1 + cdots + alpha_d A^d
has only real eigenvalues and is diagonalizable.


If the matrix A{displaystyle A}A has s distinct real eigenvalues, it follows that it is diagonalizable. In this case the system (∗){displaystyle (*)}(*) is called strictly hyperbolic.


If the matrix A{displaystyle A}A is symmetric, it follows that it is diagonalizable and the eigenvalues are real. In this case the system (∗){displaystyle (*)}(*) is called symmetric hyperbolic.



Hyperbolic system and conservation laws


There is a connection between a hyperbolic system and a conservation law. Consider a hyperbolic system of one partial differential equation for one unknown function u=u(x→,t){displaystyle u=u({vec {x}},t)}u = u(vec x, t). Then the system (∗){displaystyle (*)}(*) has the form


(∗∗)∂u∂t+∑j=1d∂∂xjfj(u)=0.{displaystyle (**)quad {frac {partial u}{partial t}}+sum _{j=1}^{d}{frac {partial }{partial x_{j}}}{f^{j}}(u)=0.}(**) quad frac{partial u}{partial t}<br /> + sum_{j=1}^d frac{partial}{partial x_j}<br /> {f^j} (u) = 0.<br />

Here, u{displaystyle u}u can be interpreted as a quantity that moves around according to the flux given by f→=(f1,…,fd){displaystyle {vec {f}}=(f^{1},ldots ,f^{d})}vec f = (f^1, ldots, f^d). To see that the quantity u{displaystyle u}u is conserved, integrate (∗∗){displaystyle (**)}(**) over a domain Ω{displaystyle Omega }Omega


∫Ω∂u∂td⁡Ω+∫Ω∇⋅f→(u)d⁡Ω=0.{displaystyle int _{Omega }{frac {partial u}{partial t}}operatorname {d} Omega +int _{Omega }nabla cdot {vec {f}}(u)operatorname {d} Omega =0.}{displaystyle int _{Omega }{frac {partial u}{partial t}}operatorname {d} Omega +int _{Omega }nabla cdot {vec {f}}(u)operatorname {d} Omega =0.}

If u{displaystyle u}u and f→{displaystyle {vec {f}}}vec f are sufficiently smooth functions, we can use the divergence theorem and change the order of the integration and ∂/∂t{displaystyle partial /partial t}partial / partial t to get a conservation law for the quantity u{displaystyle u}u in the general form


dd⁡t∫Ωud⁡Ω+∫∂Ωf→(u)⋅n→d⁡Γ=0,{displaystyle {frac {operatorname {d} }{operatorname {d} t}}int _{Omega }uoperatorname {d} Omega +int _{partial Omega }{vec {f}}(u)cdot {vec {n}}operatorname {d} Gamma =0,}{displaystyle {frac {operatorname {d} }{operatorname {d} t}}int _{Omega }uoperatorname {d} Omega +int _{partial Omega }{vec {f}}(u)cdot {vec {n}}operatorname {d} Gamma =0,}

which means that the time rate of change of u{displaystyle u}u in the domain Ω{displaystyle Omega }Omega is equal to the net flux of u{displaystyle u}u through its boundary ∂Ω{displaystyle partial Omega }partialOmega. Since this is an equality, it can be concluded that u{displaystyle u}u is conserved within Ω{displaystyle Omega }Omega .



See also



  • Elliptic partial differential equation

  • Hypoelliptic operator

  • Parabolic partial differential equation

  • Relativistic heat conduction



Notes





  1. ^ Rozhdestvenskii


  2. ^ Evans 1998, p.400


  3. ^ Evans 1998, p.402




Bibliography




  • Evans, Lawrence C. (2010) [1998], Partial differential equations (PDF), Graduate Studies in Mathematics, 19 (2nd ed.), Providence, R.I.: American Mathematical Society, doi:10.1090/gsm/019, ISBN 978-0-8218-4974-3, MR 2597943.mw-parser-output cite.citation{font-style:inherit}.mw-parser-output q{quotes:"""""""'""'"}.mw-parser-output code.cs1-code{color:inherit;background:inherit;border:inherit;padding:inherit}.mw-parser-output .cs1-lock-free a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/6/65/Lock-green.svg/9px-Lock-green.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-lock-limited a,.mw-parser-output .cs1-lock-registration a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/d/d6/Lock-gray-alt-2.svg/9px-Lock-gray-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-lock-subscription a{background:url("//upload.wikimedia.org/wikipedia/commons/thumb/a/aa/Lock-red-alt-2.svg/9px-Lock-red-alt-2.svg.png")no-repeat;background-position:right .1em center}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration{color:#555}.mw-parser-output .cs1-subscription span,.mw-parser-output .cs1-registration span{border-bottom:1px dotted;cursor:help}.mw-parser-output .cs1-hidden-error{display:none;font-size:100%}.mw-parser-output .cs1-visible-error{font-size:100%}.mw-parser-output .cs1-subscription,.mw-parser-output .cs1-registration,.mw-parser-output .cs1-format{font-size:95%}.mw-parser-output .cs1-kern-left,.mw-parser-output .cs1-kern-wl-left{padding-left:0.2em}.mw-parser-output .cs1-kern-right,.mw-parser-output .cs1-kern-wl-right{padding-right:0.2em}

  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002.
    ISBN 1-58488-299-9


  • Rozhdestvenskii, B.L. (2001) [1994], "H/h048300", in Hazewinkel, Michiel, Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4



External links




  • Hazewinkel, Michiel, ed. (2001) [1994], "Hyperbolic partial differential equation, numerical methods", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4


  • Linear Hyperbolic Equations at EqWorld: The World of Mathematical Equations.


  • Nonlinear Hyperbolic Equations at EqWorld: The World of Mathematical Equations.




-

Popular posts from this blog

Xamarin.iOS Cant Deploy on Iphone

Image
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ height:90px;width:728px;box-sizing:border-box; } -1 I have a windows laptop and a mac mini my problem is that It wont deploy on iphone if I use Visual Studio Xamarin on windows to install my app, but it works fine with Visual Studio For Mac. Here's what I get after build succeded on Visual Studio Xamarin on windows : 1>------ Build started: Project: FinalCustomerApp.iOS, Configuration: Debug iPhone ------ 1> Connecting to Mac server 192.168.8.100... 1> FinalCustomerApp.iOS -> C:UsersJeremy PaulDesktopFinalCustomerAppFinalCustomerAppFinalCustomerApp.iOSbiniPhoneDebugFinalCustomerApp.iOS.exe 1> Detected signing identity: 1> Code Signing Key: &quo
Read more

Glorious Revolution

Image
"The Bloodless Revolution" redirects here. For a history of the vegetarian movement, see The Bloodless Revolution (book). This article is about the English revolution of 1688. For the revolution of 1868 in Spain, see Glorious Revolution (Spain). For other uses, see Glorious Revolution (disambiguation). Glorious Revolution The Prince of Orange lands at Torbay Date 1688–1689 Location British Isles Also known as Revolution of 1688 War of the English Succession Bloodless Revolution Participants English, Welsh and Scottish society, Dutch forces Outcome Replacement of James II by William III and Mary II Jacobite rising of 1689 Williamite War in Ireland War with France; England and Scotland join Grand Alliance Drafting of the Bill of Rights 1689 Part of a series on the History of England Timeline Prehistoric Britain Roman Britain Sub-Roman Britain Medieval period Economy in the Middle Ages Anglo-Saxon period
Read more

Dulmage-Mendelsohn matrix decomposition in Python

Image
0 Matlab has a function called dmperm that computes the so-called Dulmage–Mendelsohn decomposition of a n x n matrix. From wikipedia, the Dulmage–Mendelsohn is a partition of the vertices of a bipartite graph into subsets, with the property that two adjacent vertices belong to the same subset if and only if they are paired with each other in a perfect matching of the graph. Looking both on scipy and numpy, I could not find this function, nor some similar version. Is it possible to implement it using basic linear algebra operations? Any idea if this is implemented in some Python package? python matlab linear-algebra graph-theory matrix-decomposition share | improve this question
Read more