Laplaces Ligning i Polære Koordinater Forklaret

Laplaces ligning er en fundamental delvis differentialligning, der optræder i mange grene af fysik og ingeniørvidenskab, herunder elektromagnetisme, varmeledning og fluidmekanik. Den beskriver systemer i en steady-state eller ligevægtstilstand. Mens den ofte præsenteres i kartesiske koordinater (x, y), bliver mange virkelige problemer, især dem med cirkulær eller cylindrisk symmetri, meget lettere at løse, når ligningen udtrykkes i polære koordinater (r, θ). Denne artikel vil guide dig gennem processen med at omskrive Laplaces ligning til dens polære form og demonstrere, hvordan den bruges til at løse randværdiproblemer i cirkulære domæner.

Hvad er Laplaces Ligning?

I to dimensioner og kartesiske koordinater er Laplaces ligning for en funktion u(x, y) givet ved:

∇²u = ∂²u/∂x² + ∂²u/∂y² = 0

Operatoren ∇² (nabla i anden) kaldes Laplace-operatoren. En funktion, der opfylder Laplaces ligning, kaldes en harmonisk funktion. Fysisk set betyder det, at værdien af funktionen i et punkt er gennemsnittet af værdierne i en lille omkreds omkring punktet. Dette er grunden til, at den beskriver ligevægtstilstande, hvor der ikke er nogen lokale kilder eller dræn.

Transformation til Polære Koordinater

For at transformere ligningen til polære koordinater skal vi relatere de kartesiske koordinater (x, y) til de polære (r, θ). Relationerne er:

x = r cos(θ)

y = r sin(θ)

Og omvendt:

r = √(x² + y²)

θ = arctan(y/x)

Ved hjælp af kædereglen for differentiation kan vi udtrykke de partielle afledede med hensyn til x og y ved hjælp af de partielle afledede med hensyn til r og θ. Efter en del algebraiske manipulationer, som involverer at finde ∂u/∂x, ∂u/∂y, ∂²u/∂x² og ∂²u/∂y², når man til den endelige og meget mere elegante form af Laplaces ligning i polære koordinater:

u_rr + (1/r)u_r + (1/r²)u_θθ = 0

Her repræsenterer u_rr den anden partielle afledede af u med hensyn til r, u_r er den første partielle afledede med hensyn til r, og u_θθ er den anden partielle afledede med hensyn til θ. Denne form er afgørende, når man arbejder med problemer på en cirkulær disk eller en annulus (et ringformet område).

Løsning af Randværdiproblemer på en Cirkulær Disk

Et klassisk problem er at finde temperaturen (eller det elektriske potentiale) u(r, θ) inde i en cirkulær disk med radius ρ, når temperaturen på kanten af disken er givet ved en funktion f(θ). Dette er kendt som et Dirichlet-problem for en disk. Vi skal løse:

u_rr + (1/r)u_r + (1/r²)u_θθ = 0, for 0 < r < ρ

med randbetingelsen u(ρ, θ) = f(θ).

Metoden: Separation af Variable

Vi antager, at en løsning kan skrives som et produkt af en funktion, der kun afhænger af r, og en funktion, der kun afhænger af θ: u(r, θ) = R(r)Θ(θ). Ved at indsætte dette i den polære Laplace-ligning og omarrangere leddene, kan vi adskille de variable:

(r²R'' + rR') / R = -Θ'' / Θ = λ

Her er λ en separationskonstant. Dette giver os to separate ordinære differentialligninger:

1. Θ'' + λΘ = 0
2. r²R'' + rR' - λR = 0

For den angulære del, Θ(θ), skal løsningen være periodisk, da θ og θ + 2π repræsenterer den samme position. Dette pålægger betingelserne Θ(-π) = Θ(π) og Θ'(-π) = Θ'(π). Disse periodiske randbetingelser fører til, at separationskonstanten λ skal være egenværdier af formen λ_n = n² for n = 0, 1, 2, ... De tilhørende løsninger (egenfunktioner) er konstanten 1 (for n=0) og funktionerne cos(nθ) og sin(nθ) (for n > 0).

Den Radiale Ligning

For hver egenværdi λ_n = n² løser vi den radiale ligning. For n = 0 bliver ligningen r²R'' + rR' = 0, hvilket har den generelle løsning R(r) = c₁ + c₂ ln(r). For at løsningen skal være fysisk meningsfuld og begrænset i centrum af disken (r=0), må vi sætte c₂ = 0, da ln(r) går mod uendelig. Så for n=0 er den radiale løsning en konstant.

For n > 0 er den radiale ligning r²R'' + rR' - n²R = 0. Dette er en Cauchy-Euler-ligning med den generelle løsning R(r) = c₁rⁿ + c₂r⁻ⁿ. Igen, for at løsningen skal være begrænset ved r=0, må vi sætte c₂ = 0. Den acceptable løsning er derfor R(r) = c₁rⁿ.

Den Generelle Løsning som en Fourier-række

Ved at kombinere de radiale og angulære løsninger og summere over alle mulige n (superpositionsprincippet), får vi den generelle løsning for u(r, θ) inde i disken:

u(r, θ) = a₀ + Σ [ (r/ρ)ⁿ (a_n cos(nθ) + b_n sin(nθ)) ] (sum fra n=1 til ∞)

Konstanterne a₀, a_n og b_n bestemmes ud fra randbetingelsen u(ρ, θ) = f(θ). De viser sig at være koefficienterne i Fourier-rækken for funktionen f(θ). Denne løsning er kendt som Poisson's integralformel for en disk.

Sammenligning af Koordinatsystemer

Valget af koordinatsystem er afgørende for at simplificere et matematisk problem. Her er en sammenligning for Laplaces ligning.

Generaliseringer: Sfæriske Koordinater og Videre

Konceptet kan udvides til tre dimensioner. I sfæriske koordinater (ρ, φ, θ) bliver Laplace-operatoren mere kompleks:

∇²u = (1/ρ²) ∂/∂ρ (ρ² ∂u/∂ρ) + (1/(ρ²sinφ)) ∂/∂φ (sinφ ∂u/∂φ) + (1/(ρ²sin²φ)) ∂²u/∂θ²

Denne form er essentiel for at løse problemer med sfærisk symmetri, såsom at finde det elektriske potentiale omkring en ladet kugle. Den generelle formel for Laplace-operatoren i vilkårlige ortogonale koordinatsystemer er et kraftfuldt værktøj i teoretisk fysik og er tæt forbundet med begrebet Laplace-Beltrami operatoren på Riemannske mangfoldigheder, som er en central del af den generelle relativitetsteori.

Ofte Stillede Spørgsmål