Laplace denklemi

Matematikte Laplace denklemi, özellikleri ilk defa Pierre-Simon Laplace tarafından çalışılmış bir kısmi diferansiyel denklemdir. Laplace denkleminin çözümleri, elektromanyetizma, astronomi ve akışkanlar dinamiği gibi birçok bilim alanında önemlidir çünkü çözümler bilhassa elektrik ve yerçekim potansiyeli ile akışkan potansiyelinin davranışını açıklar. Laplace denkleminin çözümlerinin genel teorisi aynı zamanda potansiyel teorisi olarak da bilinmektedir.

Tanım

Üç boyutta, problem x, y ve z gibi üç gerçel değişkene sahip ve iki kere türevlenebilir.Laplace denkleminin çözümlerine aynı zamanda harmonik fonksiyonlar da denmektedir. Laplace ve Poisson denklemleri eliptik kısmi diferansiyel denklemlerin en basit örnekleridir. Kısmi diferansiyel operatörü olan ve herhangi bir boyutta tanımlanabilen \scriptstyle\nabla^2'ye veya \scriptstyle\Delta'ya Laplasyen işlemcisi veya kısaca Laplasyen denmektedir.

Kartezyen koordinatlar da

 \Delta f = \frac{\partial^2 f}{\partial x^2 } + \frac{\partial^2 f}{\partial y^2 } + \frac{\partial^2 f}{\partial z^2 } = 0.

Silindirik koordinatlar da

\Delta f=\frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial f}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 f}{\partial \phi^2} + \frac{\partial^2 f}{\partial z^2} =0

Küresel koordinatlar da,

 \Delta f = \frac{1}{\rho^2}\frac{\partial}{\partial \rho} \left(\rho^2 \frac{\partial f}{\partial \rho}\right) + \frac{1}{\rho^2 \sin\theta} \frac{\partial}{\partial \theta} \left(\sin\theta \frac{\partial f}{\partial \theta}\right) + \frac{1}{\rho^2 \sin^2\theta} \frac{\partial^2 f}{\partial \varphi^2} =0.

Eğrisel koordinatlar da,

 \Delta f =\frac{\partial}{\partial \xi^j}\left(\frac{\partial f}{\partial \xi^k}g^{ki}\right) + \frac{\partial f}{\partial \xi^j} g^{jm}\Gamma^n_{mn} =0,

veya

 \Delta f = \frac{1}{\sqrt{|g|}} \frac{\partial}{\partial \xi^i}\!\left(\sqrt{|g|}g^{ij} \frac{\partial f}{\partial \xi^j}\right) =0, \qquad (g=\mathrm{det}\{g_{ij}\}).

Bu sıklıkla

\nabla^2 f = 0

olarak yazılır veya,daha genel kavramlar içinde özel olarak,

\Delta f = 0,

burada ∆ = 2 Laplace operatörü veya "Laplasyen"dir

\Delta f = \nabla^2 f =\nabla \cdot \nabla f =\operatorname{div}\operatorname{grad} f,

burada ∇   diverjans işlemcidir ("div" ile sembolize edilir) bu skalerler için vektörler gönderme, ve ∇ gradyan işlemcidir ("grad" ile ifade edilir) bu vektörler için skaler göndermedir. (bu nedenle, Laplasyen Δf ≝ div grad f, bir skaler büyüklük için skaler f fonksiyonu göndermedir ; özellikle bu bir skaler için (fonksiyon) f'in vektör grad (kısmi türevler ) göndermesidir.)

Eğer sağ-el taraf verilen özel bir fonksiyon ise, h(x, y, z), yani, Eğer böyle bir denklem

\Delta f = h

ise "Poisson denklemi" denir.

Laplace denklemi ayrıca Helmholtz denkleminin özel bir durumudur.

Not:Delta sembolü, Δ, ayrıca "bir değişiklik içinde" bazı çokluklar gösterim için yaygın olarak kullanılıyor,yani ∆Q ≝ Q + δ veya ∆Q ≝ Q + εQ bazı çok küçük skalerler δ veya ε için. Laplasyen temsil kullanımı bu kullanım ile karıştırılmamalıdır. Laplace ve Poisson denklemleri eliptik kısmi diferansiyel denklemlerin en basit örnekleridir. Kısmi diferansiyel operatörü olan ve herhangi bir boyutta tanımlanabilen \scriptstyle\nabla^2'ye veya \scriptstyle\Delta'ya Laplace operatörü veya kısaca Laplasyen denmektedir.

Sınır koşulları

Bir halka üzerinde Laplace Denklemi (iç yarıçap r=2 ve dış yarıçap R=4) ile Dirichlet sınır Koşulları: u(r=2)=0 ve u(R=4)=4sin(5*θ)

Laplace denklemi için Dirichlet problemi bir D bölgesi üzerinde tanımlı ve verilmiş başka bir fonksiyona D 'nin sınırı üzerinde eşit olan bir \varphi fonksiyonu bulmaktan ibarettir. Laplace operatörü ısı denkleminde yer aldığı için, problemin bir diğer yorumu da şöyledir: Bölgenin sınırındaki sıcaklık sabit tutulur ve bölgenin iç tarafındaki sıcaklık artık değişmeyecek şekilde beklenilir. İç bölgedeki sıcaklık dağılımı artık ilişkin Dirichlet probleminin çözümü tarafından verilecektir.

Laplace denklemi için Neumann sınır koşulları D'nin sınırında \varphi fonsiyonunu belirtmez ancak bu fonksiyonun normal türevini belirtir. Fiziksel olarak bu durum, yalnız D'nin sınırında etkisi bilinen bir vektör alanı için olan bir potansiyelin inşasına (oluşturulmasına) denk gelmektedir.

Laplace denkleminin çözümlerine harmonik fonksiyonlar denilmektedir ve bu fonksiyonların hepsi denklemin sağlandığı bölge içinde analitiktir. Eğer iki fonksiyon Laplace denkleminin (veya herhangi doğrusal homojen diferansiyel denklemin) çözümüyse, toplamları (veya herhangi doğrusal kombinasyonları) da ayrıca bir çözümdür. Süperpozisyon ilkesi de denilen bu özellik özellikle karmaşık problemlerin basit çözümlerin toplanılması yoluyla yapılan çözümlerinde çok yararlıdır.

İki boyutta Laplace denklemi

İki değişkenli Laplace denklemi

\varphi_{xx} + \varphi_{yy} = 0\,

formuna sahiptir.

Analitik fonksiyonlar

Karmaşık analitik bir fonksiyonun gerçel ve sanal kısmının her ikisi de Laplace denklemini sağlar. Eğer z=x+iy ise ve

f(z) = u(x,y) + iv(x,y)\,

ise, o zaman f(z) 'nin analitik olması için gerekli koşul aşağıdaki Cauchy-Riemann denklemlerinin sağlanmasıdır:

u_x = v_y, \quad v_x = -u_y.\,

Takip eden ifade ise

u_{yy} = (-v_x)_y = -(v_y)_x = -(u_x)_x.\,

olacaktır. Bu yüzden u Laplace denklemini sağlar. Benzer bir hesaplama yine v 'nin de Laplace denklemini sağladığını gösterir.

Aksine diğer taraftan bir harmonik fonksiyon verilirse, bu fonksiyon analitik bir f(z) fonksiyonunun gerçel kısmı olur (en azından yerel olarak). Eğer

f(z) = \varphi(x,y) + i \psi(x,y)\,

olarak alınırsa ve

\psi_x = -\varphi_y, \quad \psi_y = \varphi_x\,

şartı konulursa, o zaman Cauchy-Riemann denklemleri sağlanacaktır.

Bu ilişki ψ'yi belirlemese de artışlarını belirler:

d \psi = -\varphi_y\, dx + \varphi_x\, dy.\,

φ için Laplace denklemi ψ'nin integrallenebilme koşulunun sağlandığını gösterir:

\psi_{xy} = \psi_{yx}\,

ve bu yüzden ψ bir çizgi integrali yoluyla tanımlanabilir. İntegrallenebilme koşulu ve Stoke teoremi iki noktayı birleştiren çizgi integralinin değerinin takip edilen yoldan bağımsız olduğunu gösterir. Laplace denkleminin sonucunda çıkan çözüm çiftine eşlenik harmonik fonksiyonlar adı verilir. Bu inşa sadece yerel olarak veya takip edilen yolun bir tekilliği çevrelememesi koşuluyla geçerlidir. Örneğin, r ve θ kutupsal koordinatlar olursa ve

\varphi = \log r \,

ise, o zaman karşılık gelen analitik fonksiyon

f(z) = \log z = \log r + i\theta \,

fonksiyonudur. Bununla birlikte, θ açısı orijini çevrelemeyen bir bölge içinde tek (bir) değerlidir.

Laplace denklemi ve analitik fonksiyonlar arasındaki yakın ilişki Laplace denkleminin çözümünün her mertebeden türevi olduğunu gösterir ve bu çözüm en azından bir tekilliği çevrelemeyen bir çember içinde kuvvet serilerine genişletilebilir. Bu durum, daha az düzenliliğe sahip ısı denklemi çözümleriyle tezat bir haldedir.

Kuvvet serileri ve Fourier serileri arasında sıkı bir ilişki vardır. Bir f fonksiyonu R yarıçaplı bir çember içinde kuvvet serisine genişletilirse, bu gerçel ve sanal kısımları

c_n = a_n + i b_n.\,

şeklinde olan uygun katsayıların olduğu

f(z) = \sum_{n=0}^\infty c_n z^n,\,

ifadesi anlamına gelir. Bu yüzden,

f(z) = \sum_{n=0}^\infty \left[ a_n r^n \cos n \theta - b_n r^n \sin n \theta\right] + i \sum_{n=1}^\infty \left[ a_n r^n \sin n\theta + b_n r^n \cos n \theta\right]\,

olur ki bu da f 'nin Fourier seridir.

Akışkan akımı

u ve v nicelikleri durağan sıkıştırılamaz, dönmez bir akımın iki boyutta yatay ve dikey bileşenleri olsun. Akımın sıkıştırılamaz olmasının koşulu,

u_x + v_y=0\,

olmasıdır ve akımın dönmez olmasının şartı da

\nabla \times \mathbf{V}=v_x - u_y =0 \,

olmasıdır. Bir ψ fonksiyonunun diferansiyeli

d \psi = v\, dx - u\, dy\,

olarak tanımlanırsa, o zaman sıkıştırılamama şartı bu diferansiyel için integrallenebilme koşulu olur: Sonuçtaki fonksiyona akış fonksiyonu adı verilir çünkü bu fonksiyon akım çizgileri boyunca sabittir. ψ'nin birinci türevi

\psi_x = v, \quad \psi_y=-u \,

ile verilir ve sıkıştırılamama şartı ψ 'nin Laplace denklemini sağladığını gösterir. ψ 'ye eşlenik olan harmonik φ fonksiyonuna hız potansiyeli denilir. Cauchy-Riemann denklemleri

\varphi_x=-u, \quad \varphi_y=-v \,

ifadesini verir.

Bu yüzden her analitik fonksiyon düzlemde durağan sıkıştırılamaz, dönmez bir akışkan akıma karşılık gelir. Gerçel kısım hız potansiyeli olurken sanal kısım akış fonksiyonu olur.

Elektrostatik

Maxwell denklemleri'ne göre, iki uzay boyutunda yer alan ve zamandan bağımsız olan bir elektrik alanı (u,v),

\nabla \times (u,v) = v_x -u_y =0\,

ifadesini ve ρ'nun yük yoğunluğu olduğu

\nabla \cdot (u,v) = \rho\,

ifadesini sağlar.

Birinci Maxwell denklemi

d \varphi = -u\, dx -v\, dy\,

diferansiyeli için integrallenebilme koşuludur. Böylece elektrik potansiyeli olan φ

\varphi_x = -u, \quad \varphi_y = -v\,

ifadesini sağlayacak şekilde inşa edilebilir.

İkinci Maxwell denklemi o zaman Poisson denklemi olarak ifade edilen

\varphi_{xx} + \varphi_{yy} = -\rho\,

denklemini verir.

İki boyutta kullanılana benzer olarak, Laplace denklemi elektrostatik ve akışkan akımının üç boyutlu problemlerinde de kullanılabilir.

Üç boyutta Laplace denklemi

Temel çözüm

Laplace denkleminin temel çözümü, Dirac delta fonksiyonu \delta'nın  (x',\, y', \, z') noktasında toplanmış bir birim kaynağı gösterdiği

 \Delta u = u_{xx} + u_{yy} + u_{zz} = -\delta(x-x',y-y',z-z') \,

denklemini sağlar. Hiçbir fonksiyon bu özelliğe sahip değildir ancak yine de bu, integralleri uzay üzerinde birlik olan ve desteği (fonksiyonun sıfır olmadığı bölge) bir noktaya küçülen bir fonksiyonlar limiti olarak düşünülebilir. Temel çözümün tanımı bu yüzden, u 'nun Laplasyeninin kaynak noktasını çevreleyen herhengi bir hacim üzerinde integrali alındığında, o zaman

 \iiint_V  \operatorname{div} \nabla u \, dV =-1 \,

olduğunu gösterir.

Laplace denklemi koordinatların rotasyonuyla değişmez kalır ve bu yüzden bir temel çözümün, sadece (kaynak noktasından uzaklığı gösteren) r 'ye bağımlı olan çözümler arasından elde edilebileceği beklenir. Hacim kaynak noktası etrafında a yarıçaplı bir top olarak düşünülürse, o zaman Gauss diverjans teoremi

 -1= \iiint_V \operatorname{div} \nabla u \, dV = \iint_S u_r dS = 4\pi a^2 u_r(a)\,

ifadesini verir. O zaman takip eden ifade ise, kaynak noktası etrafında r yarıçaplı bir küre üzerindeki

  u_r(r) = -\frac{1}{4\pi r^2}\,

ifadesidir ve bu yüzden

 u = \frac{1}{4\pi r}\,

olur. Benzer bir hesap ise iki boyutta

 u = \frac{-\log r}{2\pi} \,

olduğunu gösterir.

Green fonksiyonu

Bir Green fonksiyonu da bir V hacminin S sınırındaki uygun şartı sağlayan temel bir çözümdür. Örneğin, \scriptstyle G(x,y,z;x',y',z')\, ,

 \quad V'de \quad \nabla \cdot \nabla G = -\delta(x-x',y-y',z-z') , \,
 \quad (x,y,z) \quad S'de \ ise\, \quad G = 0  \,

ifadelerini sağlayabilir.

Eğer u, V üzerinde Poisson denkleminin herhangi bir çözümüyse

 \nabla \cdot \nabla u = -f \,

olur ve u, S üzerinde g sınır değerlerini alır. O zaman aşağıdaki eşitlikleri veren (diverjans teoreminin bir sonucu olan) Green özdeşliğine başvurulabilir:

 \iiint_V \left[ G \, \nabla \cdot \nabla u - u \, \nabla \cdot \nabla G \right]\, dV  = \iiint_V \nabla \cdot \left[ G \nabla u - u \nabla G \right]\, dV = \iint_S \left[ G u_n -u G_n \right] \, dS. \,

un ve Gn gösterimleri S üzerindeki normal türevleri ifade etmektedir. u ve G 'nin sağladığı şartlar bağlamında, bu sonuç

 u(x',y',z') =  \iiint_V G f \, dV + \iint_S G_n g \, dS \,

haline gelir.

Bu yüzden, Green fonksiyonu f ve g 'nin \scriptstyle (x',y',z')\, noktalarındaki etkisini açıklar. a yarıçaplı kürenin içi düşünüldüğünde ise, Green fonksiyonu yansıtma yoluyla elde edilebilir (Sommerfeld, 1949): Kürenin merkezinden ρ kadar uzaklıkta olan P kaynak noktası,

 \rho' = \frac{a^2}{\rho} \,

uzaklıkta bulunan bir N noktasına yarıçapsal doğru boyunca yansıtılır.

Unutulmaması gereken nokta P küre içindeyse, N 'nin küre dışında olması gerektiğidir. O zaman Green fonksiyonu R 'nin P kaynak noktasına uzaklığı ve T 'nin yansıtılmış N noktasına olan uzaklığı gösterdiği

 \frac{1}{4 \pi R} - \frac{a}{4 \pi \rho T} \,

ifadesi tarafından verilir. Green fonksiyonu için olan bu ifadenin bir sonucu ise Poisson integral formülüdür. ρ, θ, ve φ, P kaynak noktası için küresel koordinatlar olsun. Burada θ dikey eksenle olan açıyı göstermektedir. (Amerikan matematik gösterimine uymaz ancak standard Avrupa ve fiziksel uygulamalarına uyum gösteren bir gösterimdir) O zaman, küre içindeki Laplace denkleminin çözümü

 \cos \Theta = \cos \varphi \cos \varphi' + \sin\varphi \sin\varphi'\cos(\theta -\theta') \,

olarak alınırsa

 u(P) = \frac{1}{4\pi} a^3\left( 1 - \frac{\rho^2}{a^2} \right) \iint \frac{g(\theta',\varphi') \sin \varphi' \, d\theta' \, d\varphi'}{(a^2 + \rho^2 - 2 a \rho \cos \Theta)^{3/2} }  \,

tarafından verilir.

Bu formülün basit bir sonucu ise şudur: u harmonikse, o zaman u 'nun kürenin merkezindeki değerleri, u 'nun küre üzerindeki değerlerinin ortalama değerleridir. Bu ortalama değer özelliği ise ivedilikle sabit olmayan bir fonksiyonun maksimum değerini kürenin içinde alamayacağı sonucunu verir.

Ayrıca bakınız

Kaynakça

  • L.C. Evans, Partial Differential Equations, American Mathematical Society, Providence, 1998. ISBN 0-8218-0772-2Kısmi Diferansiyel Denklemler.
  • I. G. Petrovsky, Partial Differential Equations, W. B. Saunders Co., Philadelphia, 1967.--Kısmi Diferansiyel Denklemler.
  • A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9Mühendisler ve Biliminsanları için Doğrusal Kısmi Diferansiyel Denklemler El Kitabı.
  • A. Sommerfeld, Partial Differential Equations in Physics, Academic Press, New York, 1949.--Fizikte Kısmi Diferansiyel Denklemler.

Dış bağlantılar

This article is issued from Vikipedi - version of the 5/21/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.