Metody matematyczne fizyki/Operatory różniczkowe

Metody matematyczne fizyki
Metody matematyczne fizyki
Operatory różniczkowe

Licencja
Autor: Mirosław Makowiecki
Absolwent UMCS Fizyki Komputerowej Uniwersytetu Marii Curie-Skłodowskiej w Lublinie
Email: miroslaw(kropka)makowiecki(małpa)gmail(kropka)pl
Dotyczy: książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami.
Użytkownika książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami nie zwalnia z odpowiedzialności prawnoautorskiej nieprzeczytanie warunków licencjonowania.
Umowa prawna: Creative Commons: uznanie autorstwa, na tych samych warunkach, z możliwością obowiązywania dodatkowych ograniczeń.
Autor tej książki dołożył wszelką staranność, aby informacje zawarte w książce były poprawne i najwyższej jakości, jednakże nie udzielana jest żadna gwarancja, czy też rękojma. Autor nie jest odpowiedzialny za wykorzystanie informacji zawarte w książce, nawet jeśli wywołaby jakąś szkodę, straty w zyskach, zastoju w prowadzeniu firmy, przedsiębiorstwa lub spółki bądź utraty informacji, niezależnie czy autor (a nawet Wikibooks) został powiadomiony o możliwości wystąpienie szkód. Informacje zawarte w książce mogą być wykorzystane tylko na własną odpowiedzialność.


Operator Nabla edytuj

Operator nabla oznaczamy symbolem ∇, jest wektorem różniczkowania, każda jego składowa jest pochodną cząstkową za każdym razem po innej współrzędnej trójwymiarowego układu współrzędnej, tzn. po kolei x,y,z, jego definicja jest:

(7.1)

Baza ortogonalna we współrzędnych uogólnionych edytuj

Niech znamy funkcję transformująca współrzędne uogólnione trójwymiarowego układu współrzędnych transformujących się według do trójwymiarowego układu kartezjańskiego, wtedy możemy napisać definicję na wektory jednostkowe zależne od "i" w krzywoliniowym układzie współrzędnych jako:

, gdzie:
(7.2)

Niech operator ∇ działa na współrzędną uogólnioną qr, wtedy wynik tej operacji rozkładamy w bazie opisywaną wzorem (7.2), którego wersory z założenia są prostopadłe do siebie:

(7.3)

Ale wektory tworzą bazę wektorów bazowych unormowanych i ortogonalnych do jedności, to pomnóżmy ostatni wzór obustronnie przez wektor , to:

(7.4)

Jeśli zastąpimy wektor przez jego definicję (7.2) we wzorze (7.4), wtedy możemy napisać lewą stronę tego równania:

(7.5)

Wzór (7.3) na podstawie przeprowadzonych obliczeń (7.5) piszemy przy definicji wersora (7.2):

(7.6)

Operator Nabla we współrzędnych uogólnionych edytuj

Napiszmy operator ∇, który jest we współrzędnych trójwymiarowego układu współrzędnych kartezjańskich, który przekształcimy przy pomocy wzoru (7.6), by póżniej przepisać jego definicję we współrzędnych uogólnionych:

(7.7)

Na podstawie obliczeń przeprowadzonych w punkcie (7.7), operator ∇ napisany jest we współrzędnych uogólnionych przy pomocy wersorów uogólnionego układu współrzędnych:

(7.8)

Operator ∇ we współrzędnych cylindrycznych i kulistych edytuj

Operator Nabla we współrzędnych cylindrycznych edytuj

Niech mamy wektor wodzący danego punktu w trójwymiarowym układzie współrzędnych w zależności od współrzędnych cylindrycznych, tzn. od (ρ,θ,z), którego piszemy:

(7.9)

Wektory bazy cylindrycznego układu współrzędnych obliczamy, korzystając przy tym z definicji wektora wodzącego w układzie kartezjańskich (7.9) i przy wykorzystaniu definicji ortogonalnych wersorów (7.2), którego są definicjami wektorów bazowych uogólnionego układu współrzędnych:

(7.10)
(7.11)
(7.12)
(7.13)
(7.14)
(7.15)
(7.16)
(7.17)
(7.18)

Jeśli wykorzystamy wzory (7.12), (7.15) i (7.18), którego definicję operatora ∇ wyjdziemy z (7.8), którego ten nasz operator był początkowo zapisany we współrzędnych kartezjańskiego układu współrzędnych opisywanego wzorem (7.1), to we współrzędnych cylindrycznych piszemy go:

(7.19)

Operator Nabla we współrzędnych kulistych edytuj

We współrzędnych kulistych wektor wodzący , napisany we współrzędnych kartezjańskich układu współrzędnych w zależności od współrzędnych kulistych, jest zdefiniowany:

(7.20)

Korzystając z wektora wodzącego w układzie kartezjańskich (7.20) i przy wykorzystaniu wzoru (7.2), wtedy możemy wyznaczyć wersory uogólnionego układu współrzędnych:

(7.21)
(7.22)
(7.23)
(7.24)
(7.25)
(7.26)
(7.27)
(7.28)
(7.29)

Jeśli wykorzystamy wzory (7.23), (7.26) i (7.29), którego definicję operatora ∇ wyjdziemy z (7.8), którego ten nasz operator był początkowo zapisanej we współrzędnych kartezjańskiego układu współrzędnych opisywanej wzorem (7.1), a we współrzędnych kulistych piszemy go:

(7.30)

Operator Δ edytuj

Kwadrat operatora ∇, czyli operator Δ piszemy przy definicji operatora wspomnianego operatora ∇ (7.1), którego to operator tutaj rozważany wychodzi z definicji operatora ∇ i ortogonalności wersorów w kartezjańskim układzie współrzędnych:

(7.31)

Wedle obliczeń przeprowadzonych w punkcie (7.31) operator Δ piszemy we współrzędnych kartezjańskich wedle sposobu zapisanego poniżej, jako sumą drugich pochodnych względem każdej ze zmiennych, tzn. x, y ,z:

(7.32)

Operator Δ we współrzędnych uogólnionych edytuj

Operator Δ we współrzędnych uogólnionych zapisujemy jako kwadrat operatora ∇ we współrzędnych uogólnionych (7.8), wtedy możemy napisać tą tożsamość operatorową dla tego opisywanego tutaj operatora:

(7.33)

Zatem na podstawie obliczeń przeprowadzonych w punkcie (7.33) dla kwadratu operatora ∇, czyli operatora Δ, którego definicja we współrzędnych uogólnionych jest:

(7.34)

Operator Δ w układzie cylindrycznym i kulistym edytuj

Operator Δ we współrzędnych cylindrycznych edytuj

Wyznaczmy operator Δ we współrzędnych cylindrycznych (7.9), wykorzystując wzór na operator Δ zapisanego w sposób ogólny dla współrzędnych w dowolnym układzie krzywoliniowym (7.34):




(7.35)

Ostatecznie wzór na operator Δ zapisanego we współrzędnych cylindrycznych jako kwadrat operatora ∇ zapisujemy jako:

(7.36)

Operator Δ we współrzędnych kulistych edytuj

Wyznaczmy operator Δ opisany we kulistym układzie współrzędnych (7.20) wykorzystując wzór na operator Δ we współrzędnych uogólnionych (7.34):







(7.37)

Aby przestawić końcowy wzór wynikających z obliczeń (7.37) musimy przekształcić pewne tożsamości, zatem dokonajmy obliczeń na wzorach ogólnych:

(7.38)

A teraz wyznaczmy drugą tożsamość operatorową wykorzystując twierdzenie o pochodnej iloczynu i z faktu istnienia jedynki operatorowej, zatem to drugie przekształcenie operatorowe możemy napisać sposobem:

(7.39)

Operator Δ napisanej początkowo w punkcie (7.37) przy przekształceniach przeprowadzonych według (7.38) i (7.39) na pewnych wyrażeniach operatorowych, to wzór na operator Δ we współrzędnych kulistych piszemy wedle sposobu:

(7.40)

Powyższe równanie na Δ (7.40) można inaczej przedstawić przy pomocy operatora Λ, gdzie ten ostatni operator zależy tylko od współrzędnych kątowych w kulistym układzie współrzędnym:

(7.41)
  • gdzie definicja operatora Λ:
(7.42)

Rotacja we współrzędnych kartezjańskich, cylindrycznych i kulistych edytuj

Rotacja we współrzędnych kartezjańskiej edytuj

Rotacje operatora nabla przy działaniu na jej funkcję wektorową piszemy na podstawie definicji iloczynu wektorowego zapisanego w punkcie (1.16):

(7.43)

Rotacja we współrzędnych cylindrycznych edytuj

Aby wyznaczyć we współrzędnych cylindrycznych wektor , w którym to w układzie współrzędnych obowiązują wektory ortogonalne:

(7.44)

Bardzo ważnym dla nas spotkaniem jest policzenie współrzędnych kartezjańskich w zależności od współrzędnych cylindrycznego układu współrzędnych, zatem korzystając ze wzorów na wektory bazowe cylindrycznym układzie współrzędnym, tzn. wektory (7.12), (7.15) i (7.18), wtedy możemy wyznaczyć trzy równania tworzące układ równań:

(7.45)

Operatory jednostkowe w cylindrycznym układzie współrzędnych w zależności od wersorów w kartezjańskim układzie współrzędnych są:

(7.46)

Wyznaczmy wektory jednostkowe w kartezjańskiego układzie współrzędnych, wiedząc, że wyznacznik główny układu równań (7.46) jest równy jeden, wtedy możemy wyznaczyć te wektory, bo wspomniany układ jest układem Cramera:

(7.47)

Mając otrzymane wzory napisane w punkcie (7.47), które są przepisem na wektory jednostkowe w kartezjańskim układzie współrzędnych w zależności od wektorów ortogonalnej bazy w układzie cylindrycznym, a także korzystając z trzech tożsamości (7.45), wtedy możemy napisać rotację wektora używając jego formalnej definicji (7.43):







(7.48)

Rotacja wektora , która jest zdefiniowana przy pomocy argumentów cylindrycznego układu współrzędnych jest napisana w troszeczkę uproszczonej postaci, ale równoważnej do rotacji obliczonego w punkcie (7.48), to w końcu dochodzimy do wniosku:

(7.49)

Rotacja we współrzędnych kulistej edytuj

Aby wyznaczyć całkę we współrzędnych cylindrycznych, to wektor , w którym to układzie współrzędnych obowiązują wektory ortogonalizowane, piszemy wedle:

(7.50)

Bardzo ważnym dla spotkaniem jest policzenie współrzędnych kartezjańskich w zależności od współrzędnych kulistego układu współrzędnych, zatem korzystając ze wzorów na wektory bazowe w naszym układzie współrzędnym, tzn. wektory (7.23), (7.26) i (7.29), zatem możemy wyznaczyć trzy równania tworzące układ równań:

(7.51)

Wyznaczmy operatory jednostkowe w kartezjańskim układzie współrzędnych w zależności od wektorów bazy jednostkowej cylindrycznego układu współrzędnych:

(7.52)

Wyznaczmy wyznacznik główny układu równań określonych w punkcie (7.52), którego definicja jest dla powyższego układu równań:


(7.53)

Na podstawie obliczeń przeprowadzonych w punkcie (7.53) wyznacznik główny układu równań (7.52) jest równy minus jeden, czyli jego wartość bezwzględna jest równa jeden, co nam ułatwi później bardzo obliczenia. Następnym krokiem dla naszego rozważanego układu równań (7.52) jest wyznaczenie wyznacznika:



(7.54)

Następnym krokiem jest wyznaczenie wyznacznika Wy dla układu równań (7.52):



(7.55)

Następnym krokiem jest wyznaczenie wyznacznika Wz dla układu równań (7.52):



(7.56)

Wyznaczmy czemu są równe wektory jednostkowe w kartezjańskim układzie współrzędnych w zależności od wektorów w kulistym układzie odniesienia wykorzystując fakt, że macierz główna układu W, czyli układu (7.53) jest równa minus jeden, i wyznaczniki Wi, które są już wyznaczone w punktach (7.54), (7.55) i (7.56), zatem te wersory obowiązujące w układzie kartezjańskim w zależność od wersorów w układzie kulistymi są napisane:

(7.57)

Otrzymane wzory napisane w punkcie (7.57), które są przepisem na wektory jednostkowe w kartezjańskim układzie współrzędnych w zależności od wektorów ortogonalnej bazy w układzie kulistym, a także korzystając z trzech tożsamości (7.51), zatem możemy napisać rotację wektora używając jego formalnej definicji w układzie współrzędnych kartezjańskich (7.43):




(7.58)

Podamy teraz końcowy wzór na rotację wektora , wyprowadzonych przy pomocy wzorów (7.30), a także (7.51) i (7.57), to dochodzimy do końcowego wzoru wynikający z (7.58):

(7.59)

Pochodne iloczynów edytuj

Tożsamość pierwsza edytuj

Tutaj będziemy wyznaczać niektóre tożsamości, korzystając z właności operatora ∇, zatem przejdźmy do naszych dowodów:

(7.60)

Oto dowód pierwszej tożsamości przedstawiająca się:



Tożsamość druga edytuj

Tożsamość druga przedstawiająca się wzorem:

(7.61)

Oto dowód drugiej tożsamości przedstawiająca się:

Tożsamość trzecia edytuj

Tożsamość trzecia przedstawiająca się wzorem:

(7.62)

Oto dowód trzeciej tożsamości przedstawiająca się:

Tożsamości czwarta edytuj

Tożsamość czwarta przedstawiająca się wzorem:

(7.63)

Oto dowód czwartej tożsamości przedstawiająca się:

Tożsamość piąta edytuj

Tożsamość piąta przedstawiająca się wzorem:

(7.64)

Oto dowód piątej tożsamości przedstawiająca się:


Twierdzenia całkowe edytuj

Tutaj podamy pewne twierdzenia całkowe z dowodami, tzn. twierdzenie gradientów, prawo Gaussa i Stokesa, którego to dowody przestawimy poniżej.

Twierdzenie dla gradientów edytuj

Twierdzenie dla gradientów przedstawia całkowanie funkcji ∇f od punktu do punktu i otrzymujemy podobny wzór jak przy całkowaniu na lini prostej, zatem to twierdzenie przedstawiamy:

(7.65)

Aby udowodnić wzór (7.65) należy policzyć wyrażenie przy założeniu, że jest wektorem jednostkowym skierowanych w kierunku różniczkowania w przestrzeni trójwymiarowej:

(7.66)

Na podstawie końcowego wniosku wynikających z obliczeń (7.66), które całkujemy obustronnie całką oznaczoną od punktu do punktu , zatem stąd wniosek, że twierdzenie dla gradientów (7.65) jest spełnione.

Twierdzenie Gaussa (tzn. twierdzenie dywergencji) edytuj

Załóżmy, że mamy pewną objętość zamkniętą ograniczoną powierzchnią zamkniętą, gdzie jest infinitezymalnym elementem tejże powierzchni, zatem to nasze wspomniane wyżej twierdzenie przyjmuje postać:

(7.67)

Aby udowodnić wzór (7.67) musimy podać wzór na dywergencją wektora , którą nazywamy w analizie matematycznej wielkość określona:

(7.68)

Wyobraźmy sobie teraz pewien sześcian o objętości dążącej do zera, dla którego wyznaczymy dywergencję wektora , zatem niech w środku tego sześcianu istnieje pewien wektor , w takim przypadku:



(7.69)

Na podstawie obliczeń przeprowadzonych w punkcie (7.68) dywergencja wektora jest to działanie operatora ∇ na wspomnianą wcześniej funkcję wektorową, zatem na podstawie definicji dywergencji pewnej funkcji wektorowej (7.68) prawo Gaussa (7.67) jest twierdzeniem udowodnionym.

Twierdzenie Stokesa (tj. dla rotacji) edytuj

Załózmy, że mamy powierzchnię ograniczoną pewnym zamkniętym konturem, wtedy prawo Stokesa przedstawiamy:

(7.70)

Rzut na wektor jednostkowy rotacji pewnej funkcji wektorowej, tzn. ten rzut jest dodatni gdy pokrywa się ten rzut z kierunkiem wektora jednostkowego , w którym, to definicji wprowadzimy powierzchnię S nieskończenie małą, którego ogranicza pewny kontur obiegany odwrotnie względem ruchu wskazówek zegara, którego zwrot naszego wektora jednostkowego jest do góry tej powierzchni (nad zegarem), nazywamy wielkością określoną:

(7.71)

We wzorze (7.71) rotacją wektora , nazywamy rotacją wzdłuż wektora normalnego , którego wartość jest określona wspomnianym wyżej wzorem, jeśli będziemy określać współrzędną iksową rotacji, to należy obrać prostokąt prostopadły do osi iksowej, ale ten wektor normalny powinien mieć ten sam zwrot co ta nasza oś, wtedy iloczyn współrzędnej iksowej rotacji i współrzędnej malutkiego elementy powierzchni wektora prostopadłego do tego wycinka powierzchni względem osi iksowej określamy przez:


(7.72)

Aby wyznaczyć współrzędną igrekową i zetową rotacji wektora możemy policzyć zmieniając w sposób cykliczny zmienne z zestawu zmiennych x,y,z w powyższych obliczeniach, to rotacja wspomnianego wektora:

(7.73)

Mając już zdefiniowaną rotację naszego wektora, która jest równa iloczynowi wektorowemu operatora ∇ i naszego wektora. Znając wzór (7.71) definicją rotacji możemy zapisać dla bardzo małego prostokącika w przybliżeniu:

(7.74)

Całą powierzchnię ograniczoną konturem zamkniętym dzielimy na malutkie prostokąciki, by te prostokąciki były równoległe wobec siebie na tej zamkniętej powierzchni, zatem cyrkulacja dla krzywej ograniczającej wielką powierzchnię jest równa sumie cyrkulacji po małych powierzchniach wcześniej określonych, można tak powiedzieć, bo cyrkulacja po małym wycinku danej powierzchni redukuje się częściowo z całką, które różnią się znakiem na tym samym wspólnym boku, czyli należącego do dwóch malutkich powierzchni. Na podstawie powyższych rozważań piszemy wzór (7.74) dla S bardzo małego, którego to biorąc całkę na obu stronach naszej przybliżonej tożsamości po małych prostokącikach, co stąd wynika wzór (7.70). Co kończy dowód.