całka krzywoliniowa, praca, cyrkulacjapole potencjalne i potencjałpole skalarne i wektoroweróżniczka zupełna i niezupełna
Matematyka 3.4
Całki krzywoliniowe, pola potencjalne i różniczki zupełne
Wnosisz plecak na szczyt jednym szlakiem albo drugim. Praca przeciw grawitacji jest ta sama, a praca przeciw tarciu nie. Czy patrząc na sam wzór siły, da się rozpoznać, która to sytuacja, bez sprawdzania wszystkich dróg? Da się, i ta lekcja pokazuje jak. Przy okazji dostaniemy energię potencjalną jako „funkcję pierwotną” siły, prosty test na pole potencjalne z ważną pułapką (wir wokół przewodu) i odpowiedź na pytanie, dlaczego w termodynamice energia wewnętrzna jest funkcją stanu, a ciepło nie.
Czas
ok. 3 h teorii + 2,5 h zadań
Zadania
17 zadań w treści, z podpowiedziami i odpowiedziami
Źródła
Feynmana wykłady z fizyki, t. I, 14-1, 14-3…14-5; J. Stewart, Multivariable Calculus, wyd. 7, §16.1–16.3
0 Prolog: dwa szlaki na szczyt
Ze schroniska na szczyt prowadzą dwa szlaki: krótki i stromy oraz długi zakosami. Wnosisz plecak o masie m na wysokość h wyżej. Ile pracy wykonujesz przeciw grawitacji? Oba szlaki dają mgh. Teraz pchasz skrzynię po podłodze z tarciem z jednego rogu pokoju do drugiego. Droga prosta i droga wzdłuż ścian dają różną pracę, bo tarcie działa przez całą drogę.
Feynman nazywa siłę zachowawczą, jeśli praca między dwoma punktami jest taka sama po każdej drodze. Dla takiej siły można zdefiniować energię potencjalną, a energia mechaniczna jest zachowana.FLP I-14-3 Ale „po każdej drodze” to nieskończenie wiele sprawdzeń.
Pytanie na start. Siła w płaszczyźnie ma składowe F=[y;x] albo F=[−y;x]. Która z nich jest zachowawcza? Czy da się to rozstrzygnąć bez liczenia pracy po wielu drogach? Odpowiedź w części 3, a liczby w przykładzie 2.1.
1 Pola wektorowe
Pole wektorowe przypisuje każdemu punktowi przestrzeni wektor: F(x,y,z)=[P;Q;R], gdzie P, Q, R to funkcje położenia. Przykłady z fizyki: siła grawitacji −GMmr/r3, pole elektryczne, prędkość wiatru albo wody w rzece, siła sprężyny −kr. Gradient z lekcji 1.4 Rachunek różniczkowy w praktyce to też pole wektorowe, zbudowane z funkcji skalarnej.
Pole rysujemy jako strzałki w wybranych punktach albo jako linie pola, styczne w każdym punkcie do wektora pola. Dwa przykłady warto mieć przed oczami. Pole [x;y] to strzałki skierowane od środka, coraz dłuższe. Pole [−y;x] to wir: strzałki okrążają środek przeciwnie do ruchu wskazówek zegara, jak prędkość obracającej się płyty (lekcja 3.2 Współrzędne biegunowe, walcowe i sferyczne, zadanie 2.4).
Symulacja 1 Pole wektorowe i praca po trzech drogach z A do B. Wybierz pole, a obok zobaczysz pracę ∫F⋅dr po każdej z trzech dróg. Dla pola [y;x] wszystkie są równe. Dla wiru [−y;x] różne. Dla wiru x2+y2[−y;x] dwie drogi dają to samo, a trzecia, która okrąża początek układu z drugiej strony, coś innego.
2 Całka krzywoliniowa i praca
Feynman definiuje pracę jako sumę F⋅Δs po małych odcinkach drogi, w granicy coraz krótszych: tylko składowa siły wzdłuż przesunięcia wykonuje pracę.FLP I-14-1 To jest całka krzywoliniowa (skierowana) pola F po krzywej C. Liczymy ją, opisując krzywą parametrycznie, r(t) dla a≤t≤b, tak jak tor w lekcji 1.3 Pochodne. Wtedy dr=r′(t)dt i
Całka krzywoliniowa (praca)
W=∫CF⋅dr=∫abF(r(t))⋅r′(t)dt=∫CPdx+Qdy+Rdz.
Wynik nie zależy od wyboru parametryzacji (szybkości przebiegania krzywej), ale zmienia znak przy zmianie kierunku przebiegu.por. Stewart §16.2
Przykład 2.1 Trzy drogi, dwa pola
Oblicz pracę pól F1=[y;x] i F2=[−y;x] od (0;0) do (1;1) po trzech drogach: prostej y=x, paraboli y=x2 i łamanej przez (1;0).
Łamana: najpierw [t;0] (dy=0), potem [1;t] (dx=0). F1: ∫010dt+∫011dt=1. F2: 0+∫011dt=1.
Wynik:F1 daje zawsze 1; F2 daje 0, 31, 1, zależnie od drogi.
Sprawdzenie:F1=∇(xy), a xy zmienia się od 0 do 1: praca to różnica wartości, niezależnie od drogi (część 3) ✓. Dla F2 na prostej y=x pole jest prostopadłe do drogi ([−t;t]⋅[1;1]=0), więc praca zero ✓.
Bywa też potrzebna całka po długości łuku z funkcji skalarnej, ∫Cfds z ds=∣r′(t)∣dt: na przykład masa drutu o gęstości liniowej f. Ta całka nie zależy od kierunku przebiegu (zadanie 2.5).
Zadania · część 2
6 zadań · ok. 40 min
Oblicz pracę pola [−y;x] po okręgu jednostkowym przebieganym raz przeciwnie do ruchu wskazówek zegara, a potem zgodnie z nim.
Odpowiedź
r=[cost;sint]: F⋅r′=sin2t+cos2t=1, praca 2π. W drugą stronę −2π.
Oblicz ∫Cxdy−ydx po odcinku od (1;0) do (0;1).
Odpowiedź
r=[1−t;t]: (1−t)⋅1−t⋅(−1)=1, całka 1.
Pokaż, że zmiana parametryzacji t=2s (dwa razy szybszy przebieg) nie zmienia całki z przykładu 2.1 po paraboli dla F2.
Odpowiedź
r=[2s;4s2], 0≤s≤21, r′=[2;8s]: ∫01/2(−4s2⋅2+2s⋅8s)ds=∫01/28s2ds=31 ✓. Szybszy przebieg daje większe r′, ale krótszy przedział s.
Wynosisz satelitę o masie 1 t z powierzchni Ziemi (R=6371 km) na wysokość 400 km. Jaką pracę wykonujesz przeciw grawitacji F=−r2GMmr^ (z GM=3,986⋅1014m3/s2)? Porównaj z mgh.
Odpowiedź
W=GMm(R1−R+h1)≈3,70⋅109 J. mgh≈3,92⋅109 J, o 6% za dużo, bo g maleje z wysokością.
Drut ma kształt półokręgu x2+y2=1 (w metrach), y≥0, i gęstość liniową λ=1+y (kg/m). Oblicz jego masę.
Odpowiedź
ds=dt dla [cost;sint]: ∫0π(1+sint)dt=π+2≈5,14 kg. Zmiana kierunku przebiegu nic nie zmienia, bo ds>0.
Pchasz skrzynię z siłą pokonującą tarcie 10 N (stały opór skierowany przeciw ruchowi). Ile pracy wykonasz na drodze prostej 5 m, a ile po półokręgu o tej samej cięciwie? Dlaczego siła tarcia nie jest polem wektorowym w sensie tej lekcji?
Odpowiedź
50 J i 10⋅π⋅2,5≈78,5 J. Siła tarcia zależy od kierunku ruchu, a nie tylko od położenia, więc nie da się jej zapisać jako F(x,y). Feynman podkreśla, że tarcie tylko pozornie łamie zasadę zachowania energii: energia przechodzi w ciepło.FLP I-14-4
Praca pola gradientowego zależy tylko od końców drogi. Pole F=∇φ nazywamy potencjalnym (fizycy mówią: zachowawczym), a φ jego potencjałem. W fizyce piszemy F=−∇U z energią potencjalną U=−φ.por. Stewart §16.3
To odpowiednik wzoru Newtona–Leibniza z lekcji 1.5 Całki: potencjał jest „funkcją pierwotną” pola. Feynman dochodzi do tego samego od strony fizyki: praca siły zachowawczej od 1 do 2 to U(1)−U(2), a F=−∇U, bo dla małego przesunięcia FxΔx=−ΔU.FLP I-14-3, 14-5
Dla pola ciągłego na zbiorze otwartym i spójnym (z każdego punktu da się dojść do każdego) trzy warunki są równoważne:
F=∇φ dla pewnej funkcji φ;
praca ∫CF⋅dr zależy tylko od końców drogi;
praca po każdej drodze zamkniętej (cyrkulacja∮CF⋅dr) jest zerowa.
Równoważność 2 i 3: drogę zamkniętą rozcinamy na dwie drogi z A do B, jedną przebytą w przeciwną stronę. Z 2 wynika 1: definiujemy φ(P) jako pracę od ustalonego punktu do P, jak Feynman definiuje U.
3.1 Test, który nie wymaga sprawdzania dróg
Jeśli F=[P;Q]=∇φ, to P=φx, Q=φy, a z twierdzenia Schwarza o pochodnych mieszanych (lekcja 1.4 Rachunek różniczkowy w praktyce) ∂P/∂y=φxy=φyx=∂Q/∂x. To warunek konieczny. Na obszarze bez „dziur” wystarcza:
Kryterium pola potencjalnego
Na płaszczyźnie: jeśli P, Q mają ciągłe pochodne na obszarze jednospójnym (bez dziur: każdą krzywą zamkniętą da się ściągnąć do punktu, nie wychodząc z obszaru), to
F=[P;Q]jest potencjalne⟺∂x∂Q=∂y∂P.
W przestrzeni warunek to trzy równości ∂yR=∂zQ, ∂zP=∂xR, ∂xQ=∂yP, czyli ∇×F=0 (rotacja, lekcja 3.5 Strumień, dywergencja i rotacja), znów na obszarze jednospójnym. Uwaga: w przestrzeni przeszkadza wycięta prosta (np. przewód), a nie wycięty punkt. R3 bez początku układu jest jednospójne, dlatego test działa dla pola Coulomba.
Wróćmy do prologu. Dla [y;x]: ∂xQ=1=∂yP, pole potencjalne z φ=xy. Dla [−y;x]: ∂xQ=1, ∂yP=−1, różnica 2: pole nie jest potencjalne, a praca zależy od drogi.
Przykład 3.1 Szukanie potencjału
Sprawdź, że F=[2xy+z;x2;x] jest potencjalne, znajdź potencjał i pracę od (0;0;0) do (1;2;3).
Test: ∂yR=0=∂zQ; ∂zP=1=∂xR; ∂xQ=2x=∂yP ✓. Przestrzeń nie ma dziur, więc pole jest potencjalne.
φx=2xy+z, więc φ=x2y+xz+g(y,z). Z φy=x2+gy=x2 mamy gy=0. Z φz=x+gz=x: gz=0, więc g jest stałą.
φ=x2y+xz; praca φ(1;2;3)−φ(0;0;0)=2+3=5.
Wynik:φ=x2y+xz (z dokładnością do stałej), praca 5.
Sprawdzenie:∇(x2y+xz)=[2xy+z;x2;x] ✓. Po odcinku r=[t;2t;3t]: F⋅r′=(4t2+3t)+2t2+3t=6t2+6t, ∫01=2+3=5 ✓.
3.2 Pułapka: wir wokół przewodu
Pole F=x2+y2[−y;x] (tak, z dokładnością do stałej, wygląda pole magnetyczne prostego przewodu z prądem, Dział 15: Prąd i pole magnetyczne) spełnia ∂xQ=∂yP wszędzie poza początkiem układu (sprawdź w zadaniu 3.4). A jednak jego cyrkulacja po okręgu wokół początku wynosi 2π, a nie zero. Nie ma sprzeczności: pole nie jest określone w (0;0), więc obszar ma dziurę i kryterium nie obowiązuje. Lokalnie to jest gradient kąta, F=∇θ, ale kąt po obejściu początku rośnie o 2π i nie jest jednoznaczną funkcją. Cyrkulacja po krzywej, która nie okrąża początku, jest zerowa. W fizyce dokładnie to mówi prawo Ampère'a: cyrkulacja pola magnetycznego zależy tylko od prądu objętego drogą.
Zadania · część 3
6 zadań · ok. 45 min
Które pola są potencjalne: [3x2y;x3+2y], [ycosx;sinx+1], [xy;xy], [y2;2xy+ey]? Znajdź potencjały.
Odpowiedź
Pierwsze: φ=x3y+y2. Drugie: φ=ysinx+y. Trzecie nie (∂xQ−∂yP=y−x). Czwarte: φ=xy2+ey.
Oblicz pracę pola [3x2y;x3+2y] od (0;0) do (2;1) po dowolnej drodze.
Odpowiedź
φ(2;1)−φ(0;0)=8+1=9.
Pokaż, że siła kulombowska (lub grawitacyjna) F=kr/r3 jest potencjalna, i znajdź jej potencjał.
∇r1=−r21r^=−r3r, więc F=∇(−rk). Energia potencjalna U=−φ=rk (dla przyciągania k<0, np. k=−GMm: U=−rGMm).
Sprawdź, że dla F=x2+y2[−y;x] zachodzi ∂xQ=∂yP poza zerem, i oblicz cyrkulację po okręgu o promieniu 3 wokół zera.
Odpowiedź
∂xx2+y2x=(x2+y2)2y2−x2=∂yx2+y2−y ✓. Na okręgu F=r1θ^, a dr=rdθθ^, więc cyrkulacja to ∫02πdθ=2π, niezależnie od promienia.
Energia potencjalna to U=x2+xy+y2. Oblicz siłę i pracę tej siły od (1;0) do (0;1).
Odpowiedź
F=−∇U=[−2x−y;−x−2y]. Praca =U(A)−U(B)=1−1=0.
Pokaż, że pole stałe F=[Fx;Fy;Fz] jest potencjalne i wywnioskuj stąd, że praca grawitacji w pobliżu Ziemi to −mgΔz po każdej drodze. Dlaczego zakosy na szlaku nic nie zmieniają w pracy, a jednak pomagają turyście?
Odpowiedź
φ=Fxx+Fyy+Fzz. Dla F=[0;0;−mg]: praca −mgΔz. Zakosy rozkładają tę samą pracę na dłuższą drogę, więc wymagana siła (i moc przy tej samej prędkości) jest mniejsza. Praca mięśni to inna sprawa: organizm zużywa energię także na samo napinanie mięśni, co Feynman wyraźnie odróżnia od pracy fizycznej.FLP I-14-1
4 Różniczki zupełne i funkcje stanu
Wyrażenie Pdx+Qdy nazywamy różniczką zupełną, jeśli jest różniczką pewnej funkcji: Pdx+Qdy=dφ. To dokładnie warunek, że pole [P;Q] jest potencjalne, więc test jest ten sam: ∂yP=∂xQ (na obszarze bez dziur). Całka z różniczki zupełnej po drodze to przyrost φ, niezależny od drogi. Gdy test nie zachodzi, piszemy często δ zamiast d, żeby pamiętać, że wynik zależy od drogi.
4.1 Termodynamika: co jest funkcją stanu
Stan gazu doskonałego opisują np. temperatura T i objętość V. Pierwsza zasada termodynamiki (Dział 9: Termodynamika i entropia) mówi, że ciepło dostarczone przy małej, quasi-statycznej (odwracalnej) zmianie to δQ=CVdT+pdV, gdzie p=nRT/V i CV to pojemność cieplna przy stałej objętości (dla gazu doskonałego zależy co najwyżej od T; dla jednoatomowego stała, 23nR). Zastosujmy test do δQ=CVdT+VnRTdV:
∂V∂CV=0=∂T∂VnRT=VnR.
Ciepło nie jest różniczką zupełną: nie istnieje „ilość ciepła zawarta w gazie”, a ciepło pobrane między dwoma stanami zależy od drogi. Ale po podzieleniu przez T:
TδQ=TCVdT+VnRdV,∂V∂TCV=0=∂T∂VnR.
To już różniczka zupełna funkcji S=CVlnT+nRlnV+const: entropii. Czynnik 1/T, który zamienia δQ w różniczkę zupełną, to czynnik całkujący. Dla dwóch zmiennych jakiś czynnik całkujący zawsze istnieje; treścią drugiej zasady jest to, że dla każdego układu działa ten sam, uniwersalny czynnik 1/T.
Przykład 4.1 Dwie drogi między tymi samymi stanami
Jeden mol jednoatomowego gazu doskonałego (CV=23R, R=8,314J/(molK)) przechodzi ze stanu A (300 K, 10 L) do stanu B (400 K, 20 L). Droga 1: najpierw ogrzanie przy stałej objętości do 400 K, potem rozprężanie izotermiczne. Droga 2: najpierw rozprężanie izotermiczne w 300 K, potem ogrzanie przy stałej objętości. Oblicz Q, W, ΔU i ΔS.
ΔU=CVΔT=23⋅8,314⋅100≈1247 J na obu drogach (energia wewnętrzna to funkcja stanu).
Praca gazu przy rozprężaniu izotermicznym W=nRTlnVAVB: droga 1 w 400 K: W1≈2305 J; droga 2 w 300 K: W2≈1729 J. Przy stałej objętości praca jest zerowa.
Q=ΔU+W: Q1≈3552 J, Q2≈2976 J.
ΔS=CVlnTATB+nRlnVAVB=12,47⋅0,2877+8,314⋅0,6931≈9,35J/K na obu drogach.
Wynik:Q i W zależą od drogi (różnica 576 J), a ΔU i ΔS nie.
Sprawdzenie:Q1−Q2=W1−W2=R⋅100⋅ln2≈576 J ✓. Dla trzeciej drogi, prostej na wykresie p(V), numeryczna całka ∫δQ/T też daje 9,35 J/K (symulacja 2).
Symulacja 2 Trzy drogi między stanami A i B na wykresie p(V). Pole pod krzywą to praca, różna dla każdej drogi. Odczyt pokazuje, że Q i W zależą od drogi, a ΔU i ΔS=∫δQ/T nie. Zmieniaj temperaturę i objętość stanu B.
Zadania · część 4
5 zadań · ok. 35 min
Czy 2xydx+x2dy jest różniczką zupełną? A ydx−xdy?
Odpowiedź
Tak: d(x2y). Nie: ∂y(y)=1=∂x(−x)=−1.
Pokaż, że x21 jest czynnikiem całkującym dla ydx−xdy (dla x=0) i znajdź funkcję, której różniczką jest wynik.
Odpowiedź
x2ydx−x1dy: ∂yx2y=x21=∂x(−x1) ✓. To d(−xy).
Oblicz Q i W dla drogi z przykładu 4.1 biegnącej po prostej na wykresie p(V), z pA=249,4 kPa do pB=166,3 kPa.
Odpowiedź
W to pole trapezu: 21(pA+pB)(VB−VA)≈21⋅415,7kPa⋅0,01m3=2078,5 J. Q=ΔU+W≈3326 J. Inne niż na drogach 1 i 2.
Gaz przechodzi cykl zamknięty. Ile wynoszą ∮dU i ∮dS? Czy ∮δQ musi być zerowe? Co oznacza niezerowe ∮δQ?
Odpowiedź
∮dU=∮dS=0 (funkcje stanu, cyrkulacja zerowa). ∮δQ nie musi być zerowe; z pierwszej zasady ∮δQ=∮δW, czyli pole cyklu na wykresie p(V): tak pracuje silnik cieplny.
Pokaż, że dp=VnRdT−V2nRTdV dla gazu doskonałego jest różniczką zupełną, i wyjaśnij, dlaczego musiało tak wyjść.
Odpowiedź
∂VVnR=−V2nR=∂T(−V2nRT) ✓. Musiało, bo p=nRT/V jest funkcją stanu, a to jej różniczka zupełna (lekcja 1.4 Rachunek różniczkowy w praktyce).
5 Błędne intuicje
„Praca zależy tylko od początku i końca drogi.”
Tylko dla pól potencjalnych. Pole [−y;x] daje na trzech drogach 0, 31 i 1 (przykład 2.1).
„Jeśli ∂xQ=∂yP, to pole jest potencjalne.”
Tylko na obszarze bez dziur. Wir x2+y2[−y;x] spełnia test poza zerem, a ma cyrkulację 2π wokół zera.
„Całka krzywoliniowa zależy od tego, jak szybko przebiegamy krzywą.”
Nie zależy od parametryzacji (zadanie 2.3). Zależy tylko od kierunku przebiegu: odwrócenie kierunku zmienia znak.
„Ciało ma określoną ilość ciepła.”
Ciepło nie jest funkcją stanu: δQ nie jest różniczką zupełną, a ciepło pobrane między dwoma stanami zależy od drogi (przykład 4.1). Funkcjami stanu są energia wewnętrzna i entropia.
„Potencjał jest jeden.”
Z dokładnością do stałej. Fizyczne są tylko różnice potencjałów, dlatego wybór zera energii potencjalnej jest umowny.
„W fizyce i w matematyce potencjał to to samo.”
Znak bywa inny: matematycy piszą F=∇φ, fizycy F=−∇U. Siła pcha w stronę spadku energii potencjalnej.
6 Sprawdź się
1 Jak obliczyć ∫CF⋅dr dla krzywej r(t)?
∫abF(r(t))⋅r′(t)dt.
2 Ile wynosi praca pola ∇φ od A do B?
φ(B)−φ(A), po każdej drodze.
3 Jaki test sprawdza, czy [P;Q] jest potencjalne?
∂xQ=∂yP na obszarze jednospójnym.
4 Ile wynosi cyrkulacja pola potencjalnego?
Zero po każdej drodze zamkniętej.
5 Dlaczego wir x2+y2[−y;x] nie jest potencjalny?
Bo ma cyrkulację 2π wokół zera; test nie obowiązuje, bo obszar ma dziurę w zerze.
6 Czy δQ jest różniczką zupełną?
Nie. Dla procesów odwracalnych różniczką zupełną jest δQ/T=dS.
7 Jaki jest związek siły z energią potencjalną?
F=−∇U.
7 Karta wzorów i co dalej
Całki krzywoliniowe
∫CF⋅dr=∫abF(r(t))⋅r′(t)dt=∫CPdx+Qdy+Rdz
zmiana kierunku: znak minus; parametryzacja: bez znaczenia
∫Cfds, ds=∣r′∣dt (masa drutu, długość)
cyrkulacja: ∮CF⋅dr
Pola potencjalne i różniczki
∫C∇φ⋅dr=φ(B)−φ(A); fizyka: F=−∇U
potencjalne ⟺ niezależność od drogi ⟺ zerowe cyrkulacje
test 2D: ∂xQ=∂yP; 3D: ∇×F=0; obszar jednospójny!
potencjał: całkuj φx=P po x, dobierz g(y,z)
Pdx+Qdy=dφ (zupełna) ⟺ ten sam test
δQ nie zupełna; δQ/T=dS (proces odwracalny; czynnik całkujący 1/T)
Do przeczytania u Feynmana
Obowiązkowo: 14-1, 14-3…14-5 (praca jako całka, siły zachowawcze i niezachowawcze, potencjał i F=−∇U).
Do przeczytania w podręcznikach matematyki
J. Stewart, Multivariable Calculus, wyd. 7: §16.1 (pola wektorowe), §16.2 (całki krzywoliniowe), §16.3 (podstawowe twierdzenie dla całek krzywoliniowych, pola zachowawcze).
K. A. Stroud, Engineering Mathematics: całki krzywoliniowe i różniczki zupełne.