Partiel differentialligning - Partial differential equation

En visualisering af en løsning til den todimensionale varme ligning med temperatur repræsenteret af den lodrette retning og farve.

I matematik er en partiel differentialligning ( PDE ) en ligning, der pålægger relationer mellem de forskellige partielle derivater af en multivariabel funktion .

Funktionen betragtes ofte som en "ukendt", der skal løses for, på samme måde som x betragtes som et ukendt tal, der skal løses for i en algebraisk ligning som x 2 - 3 x + 2 = 0 . Det er imidlertid normalt umuligt at nedskrive eksplicitte formler til løsninger af partielle differentialligninger. Der findes tilsvarende en enorm mængde moderne matematisk og videnskabelig forskning om metoder til numerisk at tilnærme løsninger til visse partielle differentialligninger ved hjælp af computere. Partielle differentialligninger indtager også en stor sektor af ren matematisk forskning , hvor de sædvanlige spørgsmål stort set handler om identifikation af generelle kvalitative træk ved løsninger af forskellige partielle differentialligninger. Blandt de mange åbne spørgsmål er eksistensen og problemfriheden af løsninger til Navier -Stokes ligninger , der blev kaldt et af Millennium Prize -problemerne i 2000.

Delvise differentialligninger er allestedsnærværende på matematisk orienterede videnskabelige områder, såsom fysik og teknik . For eksempel er de grundlæggende i den moderne videnskabelige forståelse af lyd, varme, diffusion , elektrostatik , elektrodynamik , termodynamik , væskedynamik , elasticitet , generel relativitet og kvantemekanik ( Schrodinger -ligning , Pauli -ligning osv.). De stammer også fra mange rent matematiske overvejelser, såsom differentialgeometri og beregning af variationer ; blandt andre bemærkelsesværdige applikationer er de det grundlæggende værktøj i beviset på Poincaré -formodningen fra geometrisk topologi .

Dels på grund af denne mangfoldighed af kilder er der et bredt spektrum af forskellige typer af partielle differentialligninger, og der er udviklet metoder til at håndtere mange af de individuelle ligninger, der opstår. Som sådan anerkendes det normalt, at der ikke er nogen "generel teori" om partielle differentialligninger, idet specialistviden er noget opdelt mellem flere i det væsentlige adskilte underfelter.

Almindelige differentialligninger danner en underklasse af partielle differentialligninger, der svarer til funktioner i en enkelt variabel. Stokastiske partielle differentialligninger og ikke -lokale ligninger er fra 2020 særligt bredt studerede udvidelser af "PDE" -begrebet. Mere klassiske emner, hvor der stadig er meget aktiv forskning, omfatter elliptiske og parabolske partielle differentialligninger, fluidmekanik , Boltzmann -ligninger og dispersive partielle differentialligninger.

Introduktion

Man siger, at en funktion u ( x , y , z ) af tre variable er " harmonisk " eller "en opløsning af den Laplace-ligningen ", hvis det opfylder betingelsen

Sådanne funktioner blev bredt undersøgt i det nittende århundrede på grund af deres relevans for klassisk mekanik . Hvis eksplicit givet en funktion, er det normalt et spørgsmål om ligetil beregning for at kontrollere, om den er harmonisk eller ej. For eksempel
og
er begge harmoniske mens
er ikke. Det kan være overraskende, at de to givne eksempler på harmoniske funktioner har en så slående forskellig form fra hinanden. Dette er en afspejling af det faktum, at de ikke på nogen umiddelbar måde begge er særlige tilfælde af en "generel løsningsformel" i Laplace -ligningen. Dette er i slående kontrast til tilfældet med almindelige differentialligninger (ODE'er), der stort set ligner Laplace -ligningen, med det formål at mange indledende lærebøger er at finde algoritmer, der fører til generelle løsningsformler. For Laplace -ligningen, som for et stort antal partielle differentialligninger, eksisterer sådanne løsningsformler ikke.

Arten af ​​denne fejl kan ses mere konkret i tilfælde af følgende PDE: for en funktion v ( x , y ) af to variabler, overvej ligningen

Det kan kontrolleres direkte, at enhver funktion v af formen v ( x , y ) = f ( x ) + g ( y ) for enhver enkeltvariabel funktion f og g overhovedet vil opfylde denne betingelse. Dette er langt ud over de valgmuligheder, der er tilgængelige i ODE -løsningsformler, som typisk tillader frit valg af nogle tal. I studiet af PDE har man generelt det frie valg af funktioner.

Dette valgs art varierer fra PDE til PDE. For at forstå det for en given ligning er eksistens- og entydighedssætninger normalt vigtige organisatoriske principper. I mange indledende lærebøger kan eksistensens rolle og unikke sætninger for ODE være noget uigennemsigtig; eksistenshalvdelen er normalt unødvendig, da man direkte kan kontrollere enhver foreslået løsningsformel, mens den unikke halvdel ofte kun er til stede i baggrunden for at sikre, at en foreslået løsningsformel er så generel som muligt. I modsætning hertil er eksistens- og entydighedssætninger for PDE ofte det eneste middel til at navigere gennem overflod af forskellige løsninger. Af denne grund er de også grundlæggende, når der udføres en rent numerisk simulering, da man skal have en forståelse af, hvilke data der skal ordineres af brugeren, og hvad der skal overlades til computeren for at beregne.

For at diskutere en sådan eksistens og unikke sætninger er det nødvendigt at være præcis om domænet for den "ukendte funktion". Ellers er det umuligt at formulere resultaterne på en meningsfuld måde som f.eks. "En funktion af to variabler". Det vil sige, at domænet for den ukendte funktion skal betragtes som en del af selve PDE -strukturen.

Det følgende giver to klassiske eksempler på en sådan eksistens og unikke sætninger. Selvom de to PDE, der er tale om, er så ens, er der en slående forskel i adfærd: For den første PDE har man gratis recept på en enkelt funktion, mens for den anden PDE har man gratis recept på to funktioner.

  • Lad B betegne enhedsradius-disken omkring oprindelsen i flyet. For enhver kontinuerlig funktion U på enhedscirklen er der præcis en funktion uB, sådan at
    og hvis begrænsning til enhedscirklen er givet ved U .
  • For alle funktioner f og g på den reelle linje R er der præcis en funktion uR × (−1, 1), således at
    og med u ( x , 0) = f ( x ) og u/y( x , 0) = g ( x ) for alle værdier af x .

Endnu flere fænomener er mulige. For eksempel illustrerer følgende PDE , der opstår naturligt inden for differentialgeometri , et eksempel, hvor der er en enkel og fuldstændig eksplicit løsningsformel, men med frit valg af kun tre tal og ikke engang en funktion.

  • Hvis u er en funktion på R 2 med
    så er der tal a , b og c med u ( x , y ) = ax + med + c .

I modsætning til de tidligere eksempler er denne PDE ikke -lineær på grund af kvadratrødderne og firkanterne. En lineær PDE er en sådan, at hvis den er homogen, er summen af ​​to løsninger også en løsning, og alle konstante multipler af enhver løsning er også en løsning.

Godt stillet

Godt stillet refererer til en fælles skematisk pakke med information om en PDE. For at sige, at en PDE er godt stillet, skal man have:

  • en eksistens- og entydighedssætning, der hævder, at man ved foreskrivelse af nogle frit valgte funktioner kan udpege en specifik løsning af PDE
  • ved løbende at ændre de frie valg ændrer man løbende den tilsvarende løsning

Dette er af nødvendigheden af ​​at være gældende for flere forskellige PDE noget uklart. Kravet om "kontinuitet" er især tvetydigt, da der normalt er mange ækvivalente midler, hvormed det kan defineres nøje. Det er imidlertid noget usædvanligt at studere en PDE uden at angive en måde, hvorpå den er godt poseret.

Eksistens af lokale løsninger

I en lidt svag form siger Cauchy - Kowalevski -sætningen i det væsentlige, at hvis udtrykkene i en delvis differentialligning alle består af analytiske funktioner , så findes der nødvendigvis på visse områder løsninger af PDE, som også er analytiske funktioner. Selvom dette er et grundlæggende resultat, er det i mange situationer ikke nyttigt, da man ikke let kan kontrollere domænet for de producerede løsninger. Desuden er der kendte eksempler på lineære partielle differentialligninger, hvis koefficienter har derivater af alle ordrer (som dog ikke er analytiske), men som slet ikke har nogen løsninger: dette overraskende eksempel blev opdaget af Hans Lewy i 1957. Så Cauchy-Kowalevski-sætningen er nødvendigvis begrænset i sit omfang til analytiske funktioner. Denne kontekst udelukker mange fænomener af både fysisk og matematisk interesse.

Klassifikation

Notation

Når man skriver PDE'er, er det almindeligt at betegne partielle derivater ved hjælp af abonnementer. For eksempel:

I den generelle situation, at u er en funktion af n variabler, betegner u i det første partielle derivat i forhold til i -input, u ij angiver det andet partielle derivat i forhold til i -th og j -th input, og så på.

Det græske bogstav Δ betegner Laplace -operatøren ; hvis u er en funktion af n variabler, så

I fysiklitteraturen betegnes Laplace -operatøren ofte med 2 ; i matematiklitteraturen kan 2 u også betegne den hessiske matrix af u .

Ligninger af første orden

Lineære og ikke -lineære ligninger

En PDE kaldes lineær, hvis den er lineær i det ukendte og dets derivater. For eksempel for en funktion u af x og y er en andenordens lineær PDE af formen

hvor et jeg og f er funktioner af kun de uafhængige variable. (Ofte vil de blandede-partielle derivater u xy og u yx blive ligestillet, men dette er ikke påkrævet for diskussionen om linearitet.) Hvis a i er konstanter (uafhængigt af x og y ) kaldes PDE lineær med konstante koefficienter . Hvis f er nul overalt, så er den lineære PDE homogen , ellers er den inhomogen . (Dette er adskilt fra asymptotisk homogenisering , som studerer virkningerne af højfrekvente svingninger i koefficienterne på løsninger til PDE'er.)

Nærmest lineære PDE'er er semilinear PDE'er, hvor derivater i højeste række kun vises som lineære termer, med koefficienter, der kun er funktioner for de uafhængige variabler. De lavere ordens derivater og den ukendte funktion kan forekomme vilkårligt på anden måde. For eksempel er en generel andenordens semilinear PDE i to variabler

I en quasilinær PDE vises derivater af højeste orden ligeledes kun som lineære termer, men med koefficienter muligvis funktioner af de ukendte og lavere ordensderivater :

Mange af de grundlæggende PDE'er i fysik er quasilinære, såsom Einstein -ligningerne for generel relativitet og Navier -Stokes -ligningerne, der beskriver væskens bevægelse.

En PDE uden nogen linearitetsegenskaber kaldes fuldstændig ikke-lineær og besidder ikke-lineariteter på et eller flere af de højeste ordenderivater. Et eksempel er Monge -Ampère -ligningen , der opstår i differentialgeometri .

Lineære ligninger af anden orden

Elliptiske , parabolske og hyperboliske partielle differentialligninger af rækkefølge to er blevet bredt undersøgt siden begyndelsen af ​​det tyvende århundrede. Der er imidlertid mange andre vigtige typer PDE, herunder Korteweg – de Vries -ligningen . Der er også hybrider såsom Euler -Tricomi -ligningen , som varierer fra elliptisk til hyperbolsk for forskellige områder af domænet. Der er også vigtige udvidelser af disse grundlæggende typer til højere ordens PDE, men sådan viden er mere specialiseret.

Den elliptiske/parabolske/hyperbolske klassificering giver en vejledning til passende indledende og grænsebetingelser og til løsningernes glathed. Forudsat at u xy = u yx har den generelle lineære andenordens PDE i to uafhængige variabler formen

hvor koefficienterne A , B , C ... kan afhænge af x og y . Hvis A 2 + B 2 + C 2 > 0 over et område af xy -flyet, er PDE anden orden i dette område. Denne form er analog med ligningen for et keglesnit:

Mere præcist, ved at erstatte x med X og på samme måde for andre variabler (formelt gøres dette ved en Fouriertransform ), konverterer en konstant-koefficient PDE til et polynom af samme grad med vilkårene i højeste grad (et homogent polynom , her en kvadratisk form ) er mest signifikant for klassificeringen.

Ligesom en klassificerer keglesnit og kvadratiske former i parabolsk, hyperbolsk og elliptisk baseret på den diskriminant B 2 - 4 AC , kan det samme ske for en anden ordens PDE ved et givet punkt. Imidlertid diskriminant er i en PDE afgivet B 2 - AC grund konventionen af xy sigt er 2 B snarere end B ; formelt er diskriminanten (af den tilhørende kvadratiske form) (2 B ) 2 - 4 AC = 4 ( B 2 - AC ) , med faktoren 4 faldet for enkelhedens skyld.

  1. B 2 - AC <0 ( elliptisk partiel differentialligning ): Opløsninger af elliptiske PDE'er er lige så glatte som koefficienterne tillader inden for det indre af det område, hvor ligningen og løsningerne er defineret. For eksempel er løsninger til Laplaces ligning analytiske inden for det domæne, hvor de er defineret, men løsninger kan antage grænseværdier, der ikke er glatte. Bevægelsen af ​​en væske ved subsoniske hastigheder kan tilnærmes med elliptiske PDE'er, og Euler -Tricomi -ligningen er elliptisk, hvor x <0 .
  2. B 2 - AC = 0 ( parabolisk partiel differentialligning ): Ligninger, der er parabolske på hvert punkt, kan omdannes til en form, der er analog med varmelegningen ved en ændring af uafhængige variabler. Løsninger udjævnes, når den transformerede tidsvariabel stiger. Euler -Tricomi -ligningen har parabolsk type på linjen, hvor x = 0 .
  3. B 2 - AC > 0 ( hyperbolsk partiel differentialligning ): hyperboliske ligninger bevarer eventuelle diskontinuiteter af funktioner eller derivater i de indledende data. Et eksempel er bølgelegningen . Bevægelsen af ​​en væske ved supersonisk hastighed kan tilnærmes med hyperboliske PDE'er, og Euler -Tricomi -ligningen er hyperbolisk, hvor x > 0 .

Hvis der er n uafhængige variabler x 1 , x 2 ,…, x n , har en generel lineær partiel differentialligning af anden orden formen

Klassificeringen afhænger af signaturen af egenværdierne i koefficientmatrixen a i , j .

  1. Elliptisk: egenværdierne er alle positive eller alle negative.
  2. Parabolsk: egenværdierne er alle positive eller alle negative, undtagen en der er nul.
  3. Hyperbolisk: der er kun en negativ egenværdi, og resten er positiv, eller der er kun en positiv egenværdi, og resten er negativ.
  4. Ultrahyperbolisk: der er mere end én positiv egenværdi og mere end én negativ egenværdi, og der er ingen nul -egenværdier. Der er kun en begrænset teori for ultrahyperboliske ligninger (Courant og Hilbert, 1962).

Systemer af første ordens ligninger og karakteristiske overflader

Klassifikationen af ​​partielle differentialligninger kan udvides til systemer med førsteordensligninger, hvor det ukendte u nu er en vektor med m- komponenter, og koefficientmatricerne A ν er m ved m matricer for ν = 1, 2,…, n . Den delvise differentialligning har form

hvor koefficienten matricer A ν og vektoren B kan afhænge af x og u . Hvis en hypersurface S er givet i den implicitte form
hvor φ har en ikke-nul gradient, så er S en karakteristisk overflade for operatoren L på et givet punkt, hvis den karakteristiske form forsvinder:

Den geometriske fortolkning af denne betingelse er som følger: Hvis data for u er foreskrevet på overfladen S , kan det være muligt at bestemme det normale derivat af uS ud fra differentialligningen. Hvis dataene om S og differentialligningen bestemmer det normale derivat af uS , er S ikke-karakteristisk. Hvis dataene på S og differentialligningen ikke fastsætte den normale derivat af uS , så er overfladen karakteristisk , og differentialligningen begrænser dataene på S : differentialligningen er internt til S .

  1. En første-ordens system Lu = 0 er elliptisk hvis ingen overflade er karakteristisk for L : værdierne af uS og differentialligningen altid bestemme den normale derivat af uS .
  2. Et førsteordens system er hyperbolsk på et tidspunkt, hvis der er en rumlignende overflade S med normal ξ på det tidspunkt. Dette betyder , at ligningen Q ( λξ + η ) = 0 givet m reelle rødder λ 1 , λ 2 ,…, λ m givet enhver ikke-triviel vektor η ortogonal til ξ og en skalarmultiplikator λ . Systemet er strengt hyperbolsk, hvis disse rødder altid er forskellige. Den geometriske fortolkning af denne tilstand er som følger: Den karakteristiske form Q ( ζ ) = 0 definerer en kegle (den normale kegle) med homogene koordinater ζ. I det hyperboliske tilfælde har denne kegle m -ark, og aksen ζ = λξ løber inde i disse ark: den skærer ikke nogen af ​​dem. Men når den forskydes fra oprindelsen med η, skærer denne akse hvert ark. I det elliptiske tilfælde har den normale kegle ingen rigtige plader.

Analytiske løsninger

Adskillelse af variabler

Lineære PDE'er kan reduceres til systemer med almindelige differentialligninger ved hjælp af den vigtige teknik til adskillelse af variabler. Denne teknik hviler på en egenskab ved løsninger til differentialligninger: hvis man kan finde en løsning, der løser ligningen og opfylder randbetingelserne, så er det løsningen (dette gælder også for ODE'er). Vi antager som en ansatz, at en løsnings afhængighed af parametrene rum og tid kan skrives som et produkt af termer, der hver afhænger af en enkelt parameter, og derefter se om dette kan gøres for at løse problemet.

I metoden til adskillelse af variabler reducerer man en PDE til en PDE i færre variabler, hvilket er en almindelig differentialligning, hvis den er i en variabel - disse er til gengæld lettere at løse.

Dette er muligt for simple PDE'er, som kaldes separable partielle differentialligninger , og domænet er generelt et rektangel (et produkt af intervaller). Adskilbare PDE'er svarer til diagonale matricer - tænker på "værdien for fast x " som en koordinat, kan hver koordinat forstås separat.

Dette generaliserer til karakteristikmetoden og bruges også i integrerede transformationer .

Metode til kendetegn

I særlige tilfælde kan man finde karakteristiske kurver, hvorpå ligningen reduceres til en ODE -ændrende koordinater i domænet for at rette disse kurver, tillader adskillelse af variabler, og kaldes metoden for karakteristika .

Mere generelt kan man finde karakteristiske overflader.

Integreret transformation

En integreret transformation kan transformere PDE'en til en enklere, især en adskillelig PDE. Dette svarer til diagonalisering af en operatør.

Et vigtigt eksempel på dette er Fourier -analyse , som diagonaliserer varmeligningen ved hjælp af sinbase af sinusformede bølger.

Hvis domænet er begrænset eller periodisk, er en uendelig sum af løsninger, såsom en Fourier -serie , passende, men et integral af løsninger, såsom en Fourier -integral, er generelt påkrævet for uendelige domæner. Løsningen for en punktkilde til varmeudligningen ovenfor er et eksempel på brugen af ​​et Fourier -integral.

Ændring af variabler

Ofte kan en PDE reduceres til en enklere form med en kendt løsning ved en passende ændring af variabler . For eksempel Black -Scholes -ligningen

kan reduceres til varmeligningen
ved ændring af variabler

Grundlæggende løsning

Inhomogene ligninger kan ofte løses (for konstant koefficient PDE'er skal altid løses) ved at finde den grundlæggende løsning (løsningen for en punktkilde) og derefter tage konvolutionen med randbetingelserne for at få løsningen.

Dette er analogt i signalbehandling til at forstå et filter ved dets impulsrespons .

Superposition princip

Superpositionsprincippet gælder for ethvert lineært system, herunder lineære systemer af PDE'er. En almindelig visualisering af dette koncept er interaktionen mellem to bølger i fase, der kombineres for at resultere i en større amplitude, for eksempel sin x + sin x = 2 sin x . Det samme princip kan observeres i PDE'er, hvor løsningerne kan være reelle eller komplekse og additive. Hvis u 1 og u 2 er løsninger af lineær PDE i et funktionsrum R , så er u = c 1 u 1 + c 2 u 2 med alle konstanter c 1 og c 2 også en løsning af denne PDE i det samme funktionsrum.

Metoder til ikke-lineære ligninger

Der er ingen generelt anvendelige metoder til at løse ikke -lineære PDE'er. Alligevel er eksistens og unikke resultater (såsom Cauchy - Kowalevski -sætningen ) ofte mulige, ligesom beviser for vigtige kvalitative og kvantitative egenskaber ved løsninger (at få disse resultater er en vigtig del af analysen ). Beregningsløsning til de ikke-lineære PDE'er, split-step metoden , findes for specifikke ligninger som ikke-lineær Schrödinger ligning .

Ikke desto mindre kan nogle teknikker bruges til flere typer ligninger. Den h -princippet er den mest effektive metode til at løse underbestemte ligninger. Den Riquier-Janet teorien er en effektiv metode til at opnå information om mange analytiske overbestemte systemer.

Den fremgangsmåde egenskaber kan anvendes i nogle meget specielle tilfælde at løse lineære partielle differentialligninger.

I nogle tilfælde kan en PDE løses via forstyrrelsesanalyse, hvor løsningen anses for at være en korrektion til en ligning med en kendt løsning. Alternativer er numeriske analyseteknikker fra simple begrænsede forskelsordninger til de mere modne multigrid- og endelige elementmetoder . Mange interessante problemer inden for videnskab og teknik løses på denne måde ved hjælp af computere , undertiden højtydende supercomputere .

Lie group -metode

Fra 1870 satte Sophus Lies arbejde teorien om differentialligninger på et mere tilfredsstillende grundlag. Han viste, at integrationsteorierne fra de ældre matematikere ved indførelsen af ​​det, der nu kaldes Lie -grupper , kan henvises til en fælles kilde; og at almindelige differentialligninger, der indrømmer de samme uendelige transformationer, frembyder sammenlignelige vanskeligheder med integration. Han lagde også vægt på emnet transformationer af kontakt .

En generel tilgang til løsning af PDE'er bruger symmetriegenskaben ved differentialligninger, de kontinuerlige infinitesimale transformationer af løsninger til løsninger ( Lie -teori ). Kontinuerlig gruppeteori , Lie -algebraer og differentialgeometri bruges til at forstå strukturen af ​​lineære og ikke -lineære partielle differentialligninger til generering af integrerbare ligninger, for at finde dens Lax -par , rekursionsoperatorer, Bäcklund -transformationer og endelig finde præcise analytiske løsninger til PDE.

Symmetri metoder er blevet anerkendt til at studere differentialligninger, der opstår i matematik, fysik, teknik og mange andre discipliner.

Semianalytiske metoder

Den Adomian nedbrydning metode , den Lyapunov kunstige lille parameter metode, og hans homotopi forstyrrelse metode er alle specialtilfælde af den mere generelle Homotopiteori analysemetode . Disse er serieudvidelsesmetoder, og bortset fra Lyapunov -metoden er de uafhængige af små fysiske parametre i forhold til den velkendte forstyrrelsesteori , hvilket giver disse metoder større fleksibilitet og løsningsgeneralitet.

Numeriske løsninger

De tre mest anvendte numeriske metoder til at løse PDE'er er metoden med finite element (FEM), finite volume metoder (FVM) og finite difference metoder (FDM) samt andre former for metoder kaldet Meshfree -metoder , som blev lavet til at løse problemer, hvor de førnævnte metoder er begrænsede. FEM har en fremtrædende position blandt disse metoder og især den usædvanligt effektive version af højere ordre hp-FEM . Andre hybridversioner af FEM- og Meshfree-metoder inkluderer den generaliserede endelige elementmetode (GFEM), metoden med forlænget endelig element (XFEM), metoden med spektralfinite element (SFEM), meshfree finite element metode , diskontinuerlig Galerkin finite element metode (DGFEM), Element- Free Galerkin Method (EFGM), Interpolating Element-Free Galerkin Method (IEFGM) osv.

Endelig element metode

Den endelige element metode (FEM) (dens praktiske anvendelse ofte kendt som finite element analyse (FEA)) er en numerisk teknik til at finde omtrentlige løsninger af partielle differentialligninger (PDE) samt integrale ligninger. Løsningsmetoden er enten baseret på at eliminere differentialligningen fuldstændigt (steady state -problemer) eller gøre PDE til et tilnærmelsesværdigt system af almindelige differentialligninger, som derefter numerisk integreres ved hjælp af standardteknikker såsom Eulers metode, Runge – Kutta osv.

Endelig forskel metode

Finite-difference metoder er numeriske metoder til at tilnærme løsningerne til differentialligninger ved hjælp af endelige differensligninger til omtrentlige derivater.

Endelig volumenmetode

I lighed med metoden med begrænset forskel eller endelig elementmetode beregnes værdier på diskrete steder på en masket geometri. "Endelig volumen" refererer til den lille volumen, der omgiver hvert knudepunkt på et net. I metoden med begrænset volumen konverteres overfladeintegraler i en partiel differentialligning, der indeholder et divergensudtryk, til volumenintegraler ved hjælp af divergenssætningen . Disse udtryk vurderes derefter som flux på overfladerne af hvert begrænset volumen. Fordi fluxen, der kommer ind i et givet volumen, er identisk med den, der forlader det tilstødende volumen, bevarer disse metoder masse ved design.

Energimetoden

Energimetoden er en matematisk procedure, der kan bruges til at verificere, at de oprindelige grænseværdi-problemer er velstillede. I det følgende eksempel bruges energimetoden til at bestemme, hvor og hvilke randbetingelser der skal pålægges, således at den resulterende IBVP er veloplagt. Overvej den endimensionelle hyperbolske PDE givet af

hvor er en konstant og er en ukendt funktion med starttilstand . Multiplikation med og integration over domænet giver

Brug det

hvor integration af dele er blevet brugt til det andet forhold, får vi

Her betegner standarden L2-norm. For at være godt stillet kræver vi, at energien i løsningen er ikke-stigende, dvs. den , som opnås ved at specificere på if og at if . Dette svarer til kun at pålægge randbetingelser ved indstrømningen. Bemærk, at velstilling giver mulighed for vækst med hensyn til data (indledende og grænse), og det er derfor tilstrækkeligt at vise, at det holder, når alle data er sat til nul.

Se også

Nogle almindelige PDE'er

Typer af randbetingelser

Forskellige emner

Noter

Referencer

Yderligere læsning

eksterne links