Fiziğin Muhteşem DünyasıEtkileşimli fizik simülasyonları
Ana sayfa › Elektromanyetik Teori I › Laplace Denklemi · Değişkenlere Ayırma

Laplace Denklemi: Sayısal Gevşeme ve Değişkenlere Ayırma

Yükten arınmış bölgede elektrostatik potansiyel $\nabla^2 V = 0$ denklemini sağlar; çözüm sınırdaki değerlerle tek olarak belirlenir. Kenar uzunluğu $h$ olan bir ızgarada sonlu farklarla yazılınca her iç noktadaki potansiyel, dört komşusunun ortalaması olur:

$$V_{i,j} = \tfrac14\left(V_{i+1,j} + V_{i-1,j} + V_{i,j+1} + V_{i,j-1}\right) \;\big(+\tfrac14 h^2\rho_{i,j}/\varepsilon_0\ \text{Poisson için}\big).$$

Gevşeme yöntemi bu ortalamayı tekrar tekrar uygular. Aşırı gevşemede (SOR) düzeltme $\omega$ katsayısıyla büyütülür: $V \leftarrow V + \omega\,(V_{\text{ort}} - V)$; $1<\omega<2$ yakınsamayı çok hızlandırır, en iyi değer $N$ noktalı kare ızgara için yaklaşık $\omega_{\text{opt}} = 2/(1+\sin(\pi/N))$'dir.

Dikdörtgen kutuda ($0\le x\le a$, $0\le y\le b$) üst kenarı $V_0$, diğer kenarları 0 olan problemin değişkenlere ayırma çözümü:

$$V(x,y) = \frac{4V_0}{\pi}\sum_{n=1,3,5,\dots}\frac{1}{n}\,\sin\frac{n\pi x}{a}\,\frac{\sinh(n\pi y/a)}{\sinh(n\pi b/a)} .$$
Ortalama değer özelliği ve tekillik

Laplace denkleminin çözümleri için bir çember üzerindeki ortalama, çemberin merkezindeki değere eşittir (çember içinde yük yoksa). Bu yüzden yükten arınmış bölgede potansiyelin yerel maksimumu ya da minimumu olamaz; uçdeğerler sınırdadır. Köşelerde sınır değeri süreksiz olduğundan seri orada yavaş yakınsar ve sınır çizgisi üzerinde Gibbs aşımı görülür: aşım sıçramanın ≈ %9'udur; tek uzatmada sıçrama $2V_0$ olduğundan seri köşe yanında ≈ $1.18\,V_0$ değerine çıkar.

Parametreler

Hazır sınır koşulu

Haritaya tıklayınca o noktadaki V ve E hesaplanır.

V > 0V < 0eşpotansiyelseri (kesit)

Hesaplama

Bunu dene

1. Kare kutuda (b/a = 1) "Üst kenar V₀" seçin ve gevşetin. Kutunun merkezindeki potansiyel kaç volt? Neden tam olarak V₀/4?

Açıklama
Dört kenarın her biri sırayla $V_0$ olan dört problemi üst üste koyarsanız tüm kenarlar $V_0$ olur ve çözüm her yerde $V_0$'dır. Simetri gereği merkezde her problemin katkısı eşittir, yani $V(\text{merkez}) = V_0/4$; $V_0 = 10$ V için $2.500$ V. Seri ve SOR bu değeri dört basamak doğrulukla verir.

2. ω = 1.00 (Gauss–Seidel) ile ω ≈ ωopt arasında yakınsama hızını karşılaştırın: artığın 10⁻⁶'ya inmesi için kaç yineleme gerekiyor?

Açıklama
$N = 81$ ızgarada $\omega = 1$ ile birkaç bin yineleme gerekir; $\omega_{\text{opt}} = 2/(1+\sin(\pi/80)) \approx 1.92$ ile yaklaşık 150–250 yinelemede yakınsar. SOR'da yineleme sayısı $N$ ile, Gauss–Seidel'de $N^2$ ile ölçeklenir. $\omega \ge 2$ ise yöntem ıraksar.

3. Kesit yüksekliğini y = b (üst kenar) yapın ve seri terim sayısını artırın. Kenar üzerinde seri $V_0$ değerini aşıyor mu? Aşım terim sayısıyla kayboluyor mu?

Açıklama
Kenar üzerinde seri bir kare dalganın Fourier serisidir; köşelerin yanında Gibbs aşımı oluşur: sıçramanın (tek uzatmada $-V_0$'dan $+V_0$'a, yani $2V_0$) ≈ %9'u kadar, yani seri ≈ $1.18\,V_0$'a çıkar ($V_0 = 10$ V için ≈ 11.8 V). Terim sayısı arttıkça aşım köşeye yaklaşır ama yüksekliği kaybolmaz. Kutunun içinde ise $\sinh$ çarpanı yüksek terimleri hızla bastırdığından seri çabuk yakınsar ve SOR ile farkı köşelerden uzakta %0.1'in altına iner.

Elektromanyetik Teori I dersindeki diğer simülasyonlar