Partielle differensialligninger

(Omdirigert fra «Partiell differensialligning»)

Partielle differensialligninger er en type differensialligning som inneholder en ukjent funksjon (eller funksjoner) av flere uavhengige variabler og funksjonens partiellderiverte med hensyn på disse variablene. Partielle differensialligninger anvendes til å formulere og bidra til løsningen av problemer som inneholder flere variabler; spredning av lyd eller varme, elektrostatikk, elektrodynamikk, fluidmekanikk, elastisitet med flere. Et interessant poeng er at tilsynelatende veldig forskjellige fysiske fenomener kan baseres på den samme matematiske formuleringen og dermed være styrt av akkurat samme dynamikk.

Introduksjon

rediger

En enkel differensialligning er

 

Denne relasjonen impliserer at verdiene u(x,y) er uavhengige av x. Dette betyr at den generelle løsningen av ligningen er

 

hvor f er en vilkårlig funksjon av y. En analog ordinær differensialligning er

 

som har løsningen

 

hvor C er en vilkårlig konstant verdi (uavhengig av x). Disse to eksemplene viser at generelle løsninger av ordinære differensialligninger inneholder vilkårlige konstanter, mens løsninger av partielle differensialligninger inneholder vilkårlige funksjoner. En løsning av en partiell differensialligning er generelt ikke unik. Typisk må det spesifiseres randbetingelser på det topologiske området hvor løsningen er definert. Ta eksemplet ovenfor: her kan funksjonen   bestemmes fullstendig hvis   spesifiseres på linjen  .

Eksistens og entydighet

rediger

Problemet rundt eksistens og entydighet av løsningen av ordinære differensialligninger gis et tilfredsstillende svar av Picard-Lindelöfs teorem; dette er langt fra sannheten når man studerer partielle differensialligninger. Cauchy-Kovalevskaya-teoremet sier at Cauchy-problemet for en partiell differensialligning som er analytisk i den ukjente funksjonen og dens deriverte, har en entydig analytisk løsning. Selv om dette tilsynelatende fastslår eksistens og entydighet av løsningene, så fins det eksempler på lineære partielle differensialligninger hvis koeffisienter har deriverte av alle ordener, men som ikke har noen løsning overhodet. Se Lewy (1957).

Selv om løsningen av en partiell differensialligning eksisterer og er unik, kan den ha uønskede egenskaper. Et eksempel på patologisk oppførsel er rekken av Cauchyproblemer (avhengig av n) til Laplace-ligningen

 

med initialbetingelser

 
 

hvor n er et heltall. Den deriverte av u med hensyn på y går mot 0 uniformt i x mens n øker, men løsningen er

 

Denne løsningen går mot uendelig hvis nx ikke er et heltallsmultiplum av π for enhver verdi av y forskjellig fra null. Cauchyproblemet for Laplaceligningen kalles dårlig formulert (eng: not well posed), siden løsningen ikke avhenger kontinuerlig av dataene til problemet. Slike dårlig formulerte problemer har vanligvis ikke en tilfredsstillende anvendelse til fysiske problemstillinger.

Notasjon

rediger

For PDLer er det vanlig å benytte senket skrift for å angi de partiellderiverte. Det vil si:

 
 

Spesielt i matematisk fysikk foretrekkes ofte nabla-operator som i kartesiske koordinater skrives

  for romlig deriverte.

For tidsderiverte brukes gjerne en prikk  , men senket skrift brukes også. Man skriver da bølgeligningen (se nedenfor) slik

  (matematisk notasjon)
  (fysisk notasjon)

Eksempler

rediger

Varmeligningen i én romlig dimensjon

rediger

Ligningen for spredning av varme i én dimensjon for et homogent legeme er gitt som

 

hvor u(t,x) er temperatur, og α er en positiv konstant som beskriver raten til diffusjonen. Cauchy-problemet til denne ligningen består i å spesifisere  , hvor f(x) er en vilkårlig funksjon. Generelle løsninger av varmeligningen kan finnes ved hjelp av metoden separasjon av variabler. Ved å benytte Fouriertransformasjon finner man at løsningen av varmeligningen har formen

 

hvor F er en vilkårlig funksjon. For å oppfylle initialbetingelsene er F gitt ved Fouriertransformasjonen av f, dvs.

 

Hvis f er representert ved en veldig liten, men intens varmekilde, så kan det foregående integralet tilnærmes med deltafordelingen multiplisert med lengden av kilden. I tilfellet hvor kildens styrke er normalisert til 1, er resultatet

 

og løsningen på varmeligningen blir da

 

Dette er et gaussisk integral som kan løses og gir

 

Resultatet er det samme som normalfordelingen for x med forventning 0 og varians 2αt. Varmeligningen og tilsvarende diffusjonsligninger er nyttige verktøy i studiet av tilfeldige fenomener.

Bølgeligningen i én romlig dimensjon

rediger

I bølgeligningen skal den ukjente funksjonen u(t, x) tilfredsstille følgende ligning

 

Her angir u forflytningen av en strukket fjær fra likevekt, trykkforskjellen i et rør, eller størrelsen på et elektromagnetisk felt i et rør, og konstanten c er et tall som forteller noe om hastigheten til bølgen. Cauchy-problemet for denne ligningen består i å beskrive den opprinnelige forflytningen og hastigheten av en fjær eller et annet medium:

 
 

Her er f og g vilkårlige gitte funksjoner. Løsningen på problemet er gitt ved d’Alemberts formel

 

Formelen impliserer at løsningen i punktet (t,x) kun avhenger av dataene på segmentet av startlinjen som kuttes ut av de karakteristiske kurvene

 

som trekkes bakover fra det punktet. Disse kurvene tilsvarer signaler som propagerer med en hastighet c forover og bakover.

Omvendt, innflytelsen av dataene i et vilkårlig punkt på startlinjen propagerer med den endelige hastigheten c: det er ingen effekt på utsiden av et triangel gjennom punktet hvis sider er karakteristiske kurver. Denne oppførselen er svært forskjellig fra løsningen av varmeligningen hvor effekten av en punktkilde forekommer øyeblikkelig i ethvert punkt i rommet. Løsningen gitt over gjelder også hvis t er negativ og den eksplisitte formelen viser at løsningen har en glatt avhengighet i forhold til dataene: Cauchy-problemet for bølgeligningen er velformulert både forover og bakover for bølgeligningen.

Sfæriske bølger

rediger

Dette er bølger hvor amplituden kun avhenger av den radielle distansen r fra en sentral punktkilde. For slike bølger ser den tredimensjonale bølgeligningen slik ut

 

Dette er ekvivalent til

 

hvilket betyr at størrelsen ru oppfyller den 1-dimensjonale bølgeligningen. En generell løsning for sfæriske bølger har derfor formen

 

hvor F og G er helt vilkårlige funksjoner. Stråling fra en antenne tilsvarer tilfellet med G identisk lik 0. Bølgefronten som sendes fra antennen har altså ingen forvrengning i tid: den eneste forvrengningsfaktoren er 1/r. Dette særpreget av uforstyrret propagering til bølger er ikke å finne igjen i to dimensjoner.

Laplaceligningen i to dimensjoner

rediger

Laplace-ligningen for en ukjent funksjon av to variabler har denne formen:

 

Løsninger av denne ligningen kalles harmoniske funksjoner

Sammenheng med holomorfe funksjoner

rediger

Løsninger av Laplace-ligningen i to dimensjoner er nært knyttet til analytiske funksjoner med en kompleks variabel (kjent som holomorfe funksjoner). De reelle og imaginære delene av en analytisk funksjon er konjugert harmoniske funksjoner. De oppfyller begge Laplaceligningen og gradientene deres er ortogonale. Hvis f=u+iv, så sier Cauchy–Riemanns ligninger at

 

og det følger at

 

Omvendt, for en vilkårlig harmonisk funksjon i to dimensjoner så er den realdelen til en analytisk funksjon, i det minste lokalt. Se Laplace-ligningen for detaljer.

Et typisk randverdiproblem

rediger

Man ønsker å finne en løsning som tilfredsstiller vilkårlige verdier på randen av et område. Det søkes for eksempel etter en harmonisk funksjon som har verdiene u(θ) på enhetssirkelen. Poisson fant følgende løsning

 

Petrovsky (1967, p. 248) viser hvordan denne formelen oppnås ved å summere Fourierrekken til φ. Hvis r < 1, kan den deriverte av φ beregnes ved å derivere under integraltegnet. Det kan også verifiseres at φ er analytisk selv om u er kontinuerlig, men ikke nødvendigvis deriverbar. Dette er typisk for løsninger av elliptiske PDLer: løsningen er vesentlig glattere enn randdataene. Dette skiller seg fra løsningene av bølgeligningen, og mer generelle hyperbolske PDLer som typisk ikke har flere deriverte enn data.

Euler-Tricomi ligningen

rediger

Euler-Tricomi ligningen brukes til å studere strømninger i nærheten av lydhastigheten.

 

Adveksjonsligningen

rediger

Denne beskriver transport av en bevart skalar   i et hastighetsfelt  . Det vil si:

 

Hvis hastighetsfeltet har egenskapen,  , kan ligningen forenkles til

 

Den endimensjonale adveksjonsligningen for stabil strømning   (  konstant) blir vanligvis referert til som grisebingeproblemet. Hvis   ikke er konstant og lik  , blir ligningen referert til som Burgers’ ligning.

Ginzburg-Landau-ligningen

rediger

Ligningen brukes i modelleringen av superledning.

 

Her er   og   konstanter og   er den immaginære enheten.

Dyms ligning

rediger

Oppkalt etter Harry Dym og forekommer i studiet av solitoner.

 

Initial- og randverdiproblemer

rediger

Mange problemer i matematisk fysikk er formulert som initial- og randverdiproblemer.

Vibrerende streng

rediger

Hvis en svingende streng er strukket mellom to punkter x=0 og x=L, og u angir amplituden til forflytningen av strengen, så tilfredsstiller u den endimensjonale bølgeligningen i regionen 0<x<L og t er ubegrenset. Siden strengen er låst fast i endene, er også følgende randverdier oppfylt

 

i likhet med startverdiene

 

Metoden for separasjon av variabler for bølgeligningen

 

leder til løsningene

 

hvor

 

og konstanten k må bestemmes. Randverdiene gir da at X er et multiplum av sin kx, og k må være på formen

 

Her er n et heltall. Hert ledd i summen korresponderer til en bestemt vibrasjonsmodus av strengen. Modusen n=1 kalles den fundamentale, og frekvensen til de andre modiene er alle multiplum av denne. De danner overtonene til strengen, og de er grunnlaget for all musikalsk akustikk. Initialbetingelsene kan nå tilfredsstilles ved å representere f og g som uendelige summer av disse modiene. Blåseinstrumenter korresponderer til vibrasjoner av luftsøyler med en åpen og en lukket ende. Dette gir randbetingelsene

 

Metoden med separasjon av variabler leder i dette tilfellet til en løsning som består av en rekke med odde overtoner.

Det generelle problemet av denne typen løses med Sturm-Liouville-teori.

Vibrerende membran

rediger

Hvis en membran strekkes over en kurve C som utgjør randen av et område D i planet, så styres vibrasjonene til denne av bølgeligningen

 

gitt at t>0 og (x,y) ligger inni området D. Randbetingelsene er   hvis   ligger på kurven  . Metoden med separasjon av variabler leder til

 

som må tilfredsstille

 
 

Den siste ligningen kalles Helmholtz’ ligning. Konstanten k må bestemmes for å kunne tillate en ikke-triviell v som kan oppfylle randbetingelsene på C. Disse verdiene av k2 kalles egenverdiene til Laplaceoperatoren i D, og de assosierte løsningene er egenfunksjonene. Sturm-Liouville-teorien kan utvides til dette elliptiske egenverdiproblemet (Jost, 2002).

Det er ingen generelle metoder som kan anvendes for å løse ikkelineære PDLer. Likevel, eksistens og entydighetsresultater er ofte mulige, likeså bevis på viktige kvalitative og kvantitative egenskaper av løsningene. Enkelte ligninger har også beregnbare løsninger f.eks. den ikke-lineære Schrødingerligningen.

På tross av dette kan noen teknikker brukes på flere typer ligninger. h-prinsippet er den kraftigste metoden for å løse underbestemte ligninger. Riquier-Janet-teori er en effektiv metode for å finne informasjon om mange analytisk overbestemte systemer.

Karakteristikkmetoden, også kjent som similærtransformasjonsmetoden, kan brukes på noen svært spesielle tilfeller for å løse partielle differensialligninger.

I noen tilfeller kan en PDL bli løst via perturbasjonsanalyse hvor løsningen antas å være en korreksjon til en ligning med en kjent løsning. Alternativt kan en bruke teknikker fra numerisk analyse, fra enkle endelig-differanse-skjema til mer modne multigrid- og elementmetoder. Mange interessante problemer innen forskning og ingeniørvitenskap løses på denne måten med bruk av datamaskiner, gjerne superdatamaskiner.

Andre eksempler

rediger

Schrødingerligningen er en PDL i hjertet av ikke-relativistisk kvantemekanikk. I forbindelse med WKB-tilnærmingen gjenkjennes den som Hamilton-Jacobi-ligningen.

Med unntak av Dyms ligning og Ginzburg-Landau-ligningen er ligningene ovenfor lineære i den forstand at de kan skrives på formen Au = f for en gitt lineær operator A og en gitt funksjon f. Andre viktige ikke-lineære ligninger er Navier-Stokes-ligningene som beskriver væskebevegelse og Einsteins feltligninger knyttet til generell relativitet.

Klassifisering

rediger

Noen lineære andreordens partielle differensialligninger kan klassifiseres som parabolske, hyperbolske, ultrahyperbolske eller elliptiske. Euler-Tricomi-ligningen faller innenfor ulike klassifiseringer, avhengig av hvilket område man studerer/løser ligningen i. Klassifiseringen hjelper på å bestemme passende initial- og randbetingelser, og glattheten av løsningene.

Første ordens ligninger

rediger

Ligninger av annen orden

rediger

Under antakelsen  , har en generell andre-ordens PDL følgende form

 

hvor koeffisientene A, B, C osv. kan være avhengige av x og y. Denne formen er analog til ligningen for et kjeglesnitt:

 

På samme måte som kjeglesnitt klassifiseres som parabolske, hyperbolske eller elliptiske avhengig av diskriminanten  , kan man også klassifisere en andre-ordens PDL i et gitt punkt.

  1.   : løsninger av elliptiske differensialligninger er så glatte som koeffisientene tillater i det indre av området hvor ligningen og løsningene er definert. For eksempel er løsningene av Laplaces ligning, analytisk i området hvor de er definert, men løsninger kan anta randverdier som ikke er glatte. Væskebevegelse under lydens hastighet kan tilnærmes med elliptiske PDLer, og Euler-Triconi ligningen er elliptisk hvor x<0.
  1.   : Ligninger som er parabolske i ethvert punkt kan transformeres til en form som er analog til varmeligningen ved hjelp av en transformasjon av uavhengige variabler. Løsningene glattes ut under propageringen av den transformerte tidsvariabelen. Euler-Triconi-ligningen er parabolsk for x=0.
  1.   : Hyperbolske ligninger tar vare på enhver diskontinuitet til funksjonene eller de deriverte til initialdataene. Bølgeligningen er et eksempel. Supersoniske væskebevegelser kan tilnærmes med hyperbolske PDLer, og Euler-Tricomi-ligningen er hyperbolsk for x>0.

Gitt n uavhengige variable x1, x2, ..., xn, da har en generell lineær andreordens PDL formen:

 

Klassifikasjonen avhenger av fortegnene til egenverdiene til koeffisientmatrisen:

  1. Elliptisk: Egenverdiene har alle samme fortegn, enten positivt eller negativt.
  2. Parabolsk : Egenverdiene er alle positive eller alle negative, med unntak av én som er lik null.
  3. Hyperbolsk: Det er kun én egenverdi som er negativ og resten er positive, eller omvendt.
  4. Ultrahyperbolsk: Flere positive og negative egenverdier. Spesielt er ingen av egenverdiene null. Det er begrenset med teori på området (Courant og Hilbert, 1962).

Systemer av førsteordens ligninger og karakteristiske overflater

rediger

Klassifikasjonen av partielle differensialligninger kan utvides til systemer av førsteordens ligninger hvor den ukjente u blir en vektor med m elementer, og koeffisientmatrisene   er mxm matriser for  . Den partielle diffligningen blir da seende slik ut

 

hvor koeffesientmatrisene Aν og vektoren B kan avhenge av x og u. Hvis en hyperoverflate S er gitt på den implisitte formen

 

hvor φ har en gradient forskjellig fra null, da er S en karakteristisk overflate for operatoren L i et gitt punkt hvis den karakteristiske formen forsvinner:

 

Den geometriske tolkningen av denne betingelsen er som følger: hvis data for u er gitt på overflaten S, så kan det være mulig å bestemme den normalderiverte av uS fra differensialligningen. Hvis dataene på S og differensialligningen bestemmer den normalderiverte av uS, så er S ikke-karakteristisk. Hvis dataene på S og differensialligningen IKKE bestemmer den normalderiverte av uS, så er overflaten karakteristisk, og differensialligningen begrenser dataene på S: vi sier at den er intern til S.

  1. Et førsteordens system Lu=0 er elliptisk hvis ingen overflate er karakteristisk for L: verdiene av uS og differensialligningen bestemmer alltid den normalderiverte av uS.
  2. Et førsteordens system er hyperbolsk i et punkt hvis det fins en romlignende overflate S med normal ξ i punktet. Dette betyr at gitt en vilkårlig ikke-triviell vektor η ortogonal til ξ, og en skalar multiplikator λ, har ligningen
 

m reelle røtter λ1, λ2, ..., λm. Systemet er strengt hyperbolsk hvis disse røttene alltid er forskjellige. Den geometriske tolkningen av denne betingelsen er: karakteristiske formen Q(ζ)=0 definerer en kjegle med homogene koordinater ζ. I det hyperbolske tilfellet har kjeglen m flater og aksen ζ = λ ξ går inni disse. (Den skjærer ingen av dem.) Men en forflytning η fra origo, medfører at den skjærer alle sammen. I det elliptiske tilfellet har ikke kjeglen noen reelle flater.

Blandede ligninger

rediger

Hvis en PDL har koeffisienter som ikke er konstanter, vil den sannsynligvis ikke kunne karakteriseres som en av typene ovenfor, men heller være av en blandet type. Et enkelt eksempel er den viktige Euler-Tricomi-ligningen

 

som kalles elliptisk-hyperbolsk. Dette fordi den er elliptisk for x < 0, hyperbolsk for x > 0 og degenerert parabolsk på linjen x = 0;

Metoder for å løse PDLer

rediger

Separasjon av variabler

rediger

Metoden med separasjon av variabler gir spesielle løsninger av en lineær PDL på veldig enkle domener f.eks. rektangler som tilfredsstiller initial- eller randbetingelsene.

Bytte av variabler

rediger

En PDL kan ofte reduseres til en enklere form med en kjent løsning med et passende bytte av variabler. For eksempel er Black–Scholes-ligningen

 

mulig å redusere til varmeledningsligningen

 

med følgende variabelbytte (for komplette detaljer se løsning av Black Scholes ligningen):

 
 
 
 

Karakteristikkmetoden

rediger

Superposisjonsprinsippet

rediger

Siden enhver superposisjon av løsninger av en lineær PDE igjen er en løsning, kan spesielle løsninger kombineres til å gi en mer generell løsning.

Fourierrekker

rediger

Hvis domenet er endelig eller periodisk, er det passende med en uendelig sum av løsninger slik som Fourierrekker. På uendelige domener må man generelt sett ty til Fourierintegral. Løsningen for en punktkilde for varmeligningen gitt ovenfor er et eksempel på bruk av et Fourierintegral.

Litteratur

rediger
  • R. Courant and D. Hilbert, Methods of Mathematical Physics, vol II. Wiley-Interscience, New York, 1962.
  • L.C. Evans, Partial Differential Equations, American Mathematical Society, Providence, 1998. ISBN 0-8218-0772-2
  • F. John, Partial Differential Equations, Springer-Verlag, 1982.
  • J. Jost, Partial Differential Equations, Springer-Verlag, New York, 2002.
  • Hans Lewy (1957) An example of a smooth linear partial differential equation without solution. Annals of Mathematics, 2nd Series, 66(1),155-158.
  • I.G. Petrovskii, Partial Differential Equations, W. B. Saunders Co., Philadelphia, 1967.
  • 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
  • A. D. Polyanin and V. F. Zaitsev, Handbook of Nonlinear Partial Differential Equations, Chapman & Hall/CRC Press, Boca Raton, 2004. ISBN 1-58488-355-3
  • A. D. Polyanin, V. F. Zaitsev, and A. Moussiaux, Handbook of First Order Partial Differential Equations, Taylor & Francis, London, 2002. ISBN 0-415-27267-X
  • D. Zwillinger, Handbook of Differential Equations (3rd edition), Academic Press, Boston, 1997.
  • Y. Pinchover and J. Rubinstein, An Introduction to Partial Differential Equations, Cambridge University Press, Cambridge, 2005. ISBN 978-0-521-84886-2

Eksterne lenker

rediger
  Wikibøker: Partial Differential Equations – bøker