Matematyka i fizyka od podstaw
PLEN

Matematyka › Część III. Funkcje wielu zmiennych › 3.4 szkic · PDF do druku

Zanim zaczniesz 3.1 Reguła łańcuchowa, ekstrema i mnożniki Lagrange'a w wielu zmiennych3.3 Całki podwójne i potrójne

Nowe pojęcia w tej lekcji (4)

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 mm na wysokość hh wyżej. Ile pracy wykonujesz przeciw grawitacji? Oba szlaki dają mghmgh. 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]F = [y;\ x] albo F=[−y; x]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]\vec F(x,\ y,\ z) = [P;\ Q;\ R], gdzie PP, QQ, RR to funkcje położenia. Przykłady z fizyki: siła grawitacji −GMm r⃗/r3-GMm\,\vec r/r^3, pole elektryczne, prędkość wiatru albo wody w rzece, siła sprężyny −kr⃗-k\vec r. 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][x;\ y] to strzałki skierowane od środka, coraz dłuższe. Pole [−y; x][-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 AA do BB. Wybierz pole, a obok zobaczysz pracę ∫F⃗⋅dr⃗\int\vec F\cdot\mathrm d\vec r po każdej z trzech dróg. Dla pola [y; x][y;\ x] wszystkie są równe. Dla wiru [−y; x][-y;\ x] różne. Dla wiru [−y; x]x2+y2\tfrac{[-y;\ x]}{x^2 + y^2} 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⃗\vec F\cdot\Delta\vec 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⃗\vec F po krzywej CC. Liczymy ją, opisując krzywą parametrycznie, r⃗(t)\vec r(t) dla a≤t≤ba \le t \le b, tak jak tor w lekcji 1.3 Pochodne. Wtedy dr⃗=r⃗ ′(t) dt\mathrm d\vec r = \vec r\,'(t)\,\mathrm dt i

Całka krzywoliniowa (praca)

W=∫CF⃗⋅dr⃗=∫abF⃗(r⃗(t))⋅r⃗ ′(t) dt=∫CP dx+Q dy+R dz.W = \int_C \vec F\cdot\mathrm d\vec r = \int_a^b \vec F\big(\vec r(t)\big)\cdot\vec r\,'(t)\,\mathrm dt = \int_C P\,\mathrm dx + Q\,\mathrm dy + R\,\mathrm dz .

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 F⃗1=[y; x]\vec F_1 = [y;\ x] i F⃗2=[−y; x]\vec F_2 = [-y;\ x] od (0; 0)(0;\ 0) do (1; 1)(1;\ 1) po trzech drogach: prostej y=xy = x, paraboli y=x2y = x^2 i łamanej przez (1; 0)(1;\ 0).

  1. Prosta: r⃗=[t; t]\vec r = [t;\ t], r⃗ ′=[1; 1]\vec r\,' = [1;\ 1].   F⃗1⋅r⃗ ′=t+t\;\vec F_1\cdot\vec r\,' = t + t, ∫012t dt=1\int_0^1 2t\,\mathrm dt = 1.   F⃗2⋅r⃗ ′=−t+t=0\;\vec F_2\cdot\vec r\,' = -t + t = 0.
  2. Parabola: r⃗=[t; t2]\vec r = [t;\ t^2], r⃗ ′=[1; 2t]\vec r\,' = [1;\ 2t].   F⃗1\;\vec F_1: ∫01(t2+2t2) dt=1\int_0^1(t^2 + 2t^2)\,\mathrm dt = 1.   F⃗2\;\vec F_2: ∫01(−t2+2t2) dt=13\int_0^1(-t^2 + 2t^2)\,\mathrm dt = \tfrac13.
  3. Łamana: najpierw [t; 0][t;\ 0] (dy=0\mathrm dy = 0), potem [1; t][1;\ t] (dx=0\mathrm dx = 0).   F⃗1\;\vec F_1: ∫010 dt+∫011 dt=1\int_0^1 0\,\mathrm dt + \int_0^1 1\,\mathrm dt = 1.   F⃗2\;\vec F_2: 0+∫011 dt=10 + \int_0^1 1\,\mathrm dt = 1.

Wynik: F⃗1\vec F_1 daje zawsze 11; F⃗2\vec F_2 daje 00, 13\tfrac13, 11, zależnie od drogi.

Sprawdzenie: F⃗1=∇(xy)\vec F_1 = \nabla(xy), a xyxy zmienia się od 00 do 11: praca to różnica wartości, niezależnie od drogi (część 3) ✓. Dla F⃗2\vec F_2 na prostej y=xy = x pole jest prostopadłe do drogi ([−t; t]⋅[1; 1]=0[-t;\ t]\cdot[1;\ 1] = 0), więc praca zero ✓.

Bywa też potrzebna całka po długości łuku z funkcji skalarnej, ∫Cf ds\int_C f\,\mathrm ds z ds=∣r⃗ ′(t)∣ dt\mathrm ds = |\vec r\,'(t)|\,\mathrm dt: na przykład masa drutu o gęstości liniowej ff. Ta całka nie zależy od kierunku przebiegu (zadanie 2.5).

Zadania · część 2

6 zadań · ok. 40 min

  1. Oblicz pracę pola [−y; x][-y;\ x] po okręgu jednostkowym przebieganym raz przeciwnie do ruchu wskazówek zegara, a potem zgodnie z nim.

    Odpowiedź

    r⃗=[cos⁡t; sin⁡t]\vec r = [\cos t;\ \sin t]: F⃗⋅r⃗ ′=sin⁡2t+cos⁡2t=1\vec F\cdot\vec r\,' = \sin^2t + \cos^2t = 1, praca 2π2\pi. W drugą stronę −2π-2\pi.

  2. Oblicz ∫Cx dy−y dx\int_C x\,\mathrm dy - y\,\mathrm dx po odcinku od (1; 0)(1;\ 0) do (0; 1)(0;\ 1).

    Odpowiedź

    r⃗=[1−t; t]\vec r = [1 - t;\ t]: (1−t)⋅1−t⋅(−1)=1(1 - t)\cdot1 - t\cdot(-1) = 1, całka 11.

  3. Pokaż, że zmiana parametryzacji t=2st = 2s (dwa razy szybszy przebieg) nie zmienia całki z przykładu 2.1 po paraboli dla F⃗2\vec F_2.

    Odpowiedź

    r⃗=[2s; 4s2]\vec r = [2s;\ 4s^2], 0≤s≤120 \le s \le \tfrac12, r⃗ ′=[2; 8s]\vec r\,' = [2;\ 8s]: ∫01/2(−4s2⋅2+2s⋅8s) ds=∫01/28s2 ds=13\int_0^{1/2}(-4s^2\cdot2 + 2s\cdot8s)\,\mathrm ds = \int_0^{1/2}8s^2\,\mathrm ds = \tfrac13 ✓. Szybszy przebieg daje większe r⃗ ′\vec r\,', ale krótszy przedział ss.

  4. Wynosisz satelitę o masie 11 t z powierzchni Ziemi (R=6371R = 6371 km) na wysokość 400400 km. Jaką pracę wykonujesz przeciw grawitacji F⃗=−GMmr2r^\vec F = -\tfrac{GMm}{r^2}\hat r (z GM=3,986⋅1014 m3/s2GM = 3{,}986\cdot10^{14}\ \mathrm{m^3/s^2})? Porównaj z mghmgh.

    Odpowiedź

    W=GMm(1R−1R+h)≈3,70⋅109W = GMm\left(\tfrac1R - \tfrac1{R + h}\right) \approx 3{,}70\cdot10^9 J.   mgh≈3,92⋅109\;mgh \approx 3{,}92\cdot10^9 J, o 6% za dużo, bo gg maleje z wysokością.

  5. Drut ma kształt półokręgu x2+y2=1x^2 + y^2 = 1 (w metrach), y≥0y \ge 0, i gęstość liniową λ=1+y\lambda = 1 + y (kg/m). Oblicz jego masę.

    Odpowiedź

    ds=dt\mathrm ds = \mathrm dt dla [cos⁡t; sin⁡t][\cos t;\ \sin t]: ∫0π(1+sin⁡t) dt=π+2≈5,14\int_0^\pi(1 + \sin t)\,\mathrm dt = \pi + 2 \approx 5{,}14 kg. Zmiana kierunku przebiegu nic nie zmienia, bo ds>0\mathrm ds > 0.

  6. Pchasz skrzynię z siłą pokonującą tarcie 1010 N (stały opór skierowany przeciw ruchowi). Ile pracy wykonasz na drodze prostej 55 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ź

    5050 J i 10⋅π⋅2,5≈78,510\cdot\pi\cdot2{,}5 \approx 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)\vec F(x,\ y). Feynman podkreśla, że tarcie tylko pozornie łamie zasadę zachowania energii: energia przechodzi w ciepło.FLP I-14-4

3 Pola potencjalne

Jeśli pole jest gradientem, F⃗=∇φ\vec F = \nabla\varphi, to wzdłuż drogi F⃗⋅dr⃗=∇φ⋅dr⃗=dφ\vec F\cdot\mathrm d\vec r = \nabla\varphi\cdot\mathrm d\vec r = \mathrm d\varphi (lekcja 3.1 Reguła łańcuchowa, ekstrema i mnożniki Lagrange'a w wielu zmiennych), a całka z różniczki to przyrost funkcji:

Podstawowe twierdzenie dla całek krzywoliniowych

∫C∇φ⋅dr⃗=φ(B)−φ(A).\int_C \nabla\varphi\cdot\mathrm d\vec r = \varphi(B) - \varphi(A) .

Praca pola gradientowego zależy tylko od końców drogi. Pole F⃗=∇φ\vec F = \nabla\varphi nazywamy potencjalnym (fizycy mówią: zachowawczym), a φ\varphi jego potencjałem. W fizyce piszemy F⃗=−∇U\vec F = -\nabla U z energią potencjalną U=−φU = -\varphi.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 11 do 22 to U(1)−U(2)U(1) - U(2), a F⃗=−∇U\vec F = -\nabla U, bo dla małego przesunięcia FxΔx=−ΔUF_x\Delta x = -\Delta 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:

  1. F⃗=∇φ\vec F = \nabla\varphi dla pewnej funkcji φ\varphi;
  2. praca ∫CF⃗⋅dr⃗\int_C\vec F\cdot\mathrm d\vec r zależy tylko od końców drogi;
  3. praca po każdej drodze zamkniętej (cyrkulacja ∮CF⃗⋅dr⃗\oint_C\vec F\cdot\mathrm d\vec r) jest zerowa.

Równoważność 2 i 3: drogę zamkniętą rozcinamy na dwie drogi z AA do BB, jedną przebytą w przeciwną stronę. Z 2 wynika 1: definiujemy φ(P)\varphi(P) jako pracę od ustalonego punktu do PP, jak Feynman definiuje UU.

3.1 Test, który nie wymaga sprawdzania dróg

Jeśli F⃗=[P; Q]=∇φ\vec F = [P;\ Q] = \nabla\varphi, to P=φxP = \varphi_x, Q=φyQ = \varphi_y, a z twierdzenia Schwarza o pochodnych mieszanych (lekcja 1.4 Rachunek różniczkowy w praktyce) ∂P/∂y=φxy=φyx=∂Q/∂x\partial P/\partial y = \varphi_{xy} = \varphi_{yx} = \partial Q/\partial x. To warunek konieczny. Na obszarze bez „dziur” wystarcza:

Kryterium pola potencjalnego

Na płaszczyźnie: jeśli PP, QQ 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  ⟺  ∂Q∂x=∂P∂y.\vec F = [P;\ Q]\ \text{jest potencjalne} \iff \frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y} .

W przestrzeni warunek to trzy równości ∂yR=∂zQ\partial_yR = \partial_zQ,   ∂zP=∂xR\;\partial_zP = \partial_xR,   ∂xQ=∂yP\;\partial_xQ = \partial_yP, czyli ∇×F⃗=0⃗\nabla\times\vec F = \vec 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\mathbb R^3 bez początku układu jest jednospójne, dlatego test działa dla pola Coulomba.

Wróćmy do prologu. Dla [y; x][y;\ x]: ∂xQ=1=∂yP\partial_xQ = 1 = \partial_yP, pole potencjalne z φ=xy\varphi = xy. Dla [−y; x][-y;\ x]: ∂xQ=1\partial_xQ = 1, ∂yP=−1\partial_yP = -1, różnica 22: pole nie jest potencjalne, a praca zależy od drogi.

Przykład 3.1 Szukanie potencjału

Sprawdź, że F⃗=[2xy+z; x2; x]\vec F = [2xy + z;\ x^2;\ x] jest potencjalne, znajdź potencjał i pracę od (0; 0; 0)(0;\ 0;\ 0) do (1; 2; 3)(1;\ 2;\ 3).

  1. Test: ∂yR=0=∂zQ\partial_yR = 0 = \partial_zQ;   ∂zP=1=∂xR\;\partial_zP = 1 = \partial_xR;   ∂xQ=2x=∂yP\;\partial_xQ = 2x = \partial_yP ✓. Przestrzeń nie ma dziur, więc pole jest potencjalne.
  2. φx=2xy+z\varphi_x = 2xy + z, więc φ=x2y+xz+g(y, z)\varphi = x^2y + xz + g(y,\ z). Z φy=x2+gy=x2\varphi_y = x^2 + g_y = x^2 mamy gy=0g_y = 0. Z φz=x+gz=x\varphi_z = x + g_z = x: gz=0g_z = 0, więc gg jest stałą.
  3. φ=x2y+xz\varphi = x^2y + xz; praca φ(1; 2; 3)−φ(0; 0; 0)=2+3=5\varphi(1;\ 2;\ 3) - \varphi(0;\ 0;\ 0) = 2 + 3 = 5.

Wynik: φ=x2y+xz\varphi = x^2y + xz (z dokładnością do stałej), praca 55.

Sprawdzenie: ∇(x2y+xz)=[2xy+z; x2; x]\nabla(x^2y + xz) = [2xy + z;\ x^2;\ x] ✓. Po odcinku r⃗=[t; 2t; 3t]\vec r = [t;\ 2t;\ 3t]: F⃗⋅r⃗ ′=(4t2+3t)+2t2+3t=6t2+6t\vec F\cdot\vec r\,' = (4t^2 + 3t) + 2t^2 + 3t = 6t^2 + 6t, ∫01=2+3=5\int_0^1 = 2 + 3 = 5 ✓.

3.2 Pułapka: wir wokół przewodu

Pole F⃗=[−y; x]x2+y2\vec F = \dfrac{[-y;\ x]}{x^2 + y^2} (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\partial_xQ = \partial_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π2\pi, a nie zero. Nie ma sprzeczności: pole nie jest określone w (0; 0)(0;\ 0), więc obszar ma dziurę i kryterium nie obowiązuje. Lokalnie to jest gradient kąta, F⃗=∇θ\vec F = \nabla\theta, ale kąt po obejściu początku rośnie o 2π2\pi 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

  1. Które pola są potencjalne: [3x2y; x3+2y][3x^2y;\ x^3 + 2y],   [ycos⁡x; sin⁡x+1]\;[y\cos x;\ \sin x + 1],   [xy; xy]\;[xy;\ xy],   [y2; 2xy+ey]\;[y^2;\ 2xy + e^y]? Znajdź potencjały.

    Odpowiedź

    Pierwsze: φ=x3y+y2\varphi = x^3y + y^2. Drugie: φ=ysin⁡x+y\varphi = y\sin x + y. Trzecie nie (∂xQ−∂yP=y−x\partial_xQ - \partial_yP = y - x). Czwarte: φ=xy2+ey\varphi = xy^2 + e^y.

  2. Oblicz pracę pola [3x2y; x3+2y][3x^2y;\ x^3 + 2y] od (0; 0)(0;\ 0) do (2; 1)(2;\ 1) po dowolnej drodze.

    Odpowiedź

    φ(2; 1)−φ(0; 0)=8+1=9\varphi(2;\ 1) - \varphi(0;\ 0) = 8 + 1 = 9.

  3. Pokaż, że siła kulombowska (lub grawitacyjna) F⃗=k r⃗/r3\vec F = k\,\vec r/r^3 jest potencjalna, i znajdź jej potencjał.

    Wskazówka

    Policz ∇(1/r)\nabla(1/r) we współrzędnych sferycznych (lekcja 3.2 Współrzędne biegunowe, walcowe i sferyczne).

    Odpowiedź

    ∇1r=−1r2r^=−r⃗r3\nabla\tfrac1r = -\tfrac{1}{r^2}\hat r = -\tfrac{\vec r}{r^3}, więc F⃗=∇(−kr)\vec F = \nabla\left(-\tfrac kr\right). Energia potencjalna U=−φ=krU = -\varphi = \tfrac kr (dla przyciągania k<0k < 0, np. k=−GMmk = -GMm: U=−GMmrU = -\tfrac{GMm}{r}).

  4. Sprawdź, że dla F⃗=[−y; x]x2+y2\vec F = \tfrac{[-y;\ x]}{x^2 + y^2} zachodzi ∂xQ=∂yP\partial_xQ = \partial_yP poza zerem, i oblicz cyrkulację po okręgu o promieniu 33 wokół zera.

    Odpowiedź

    ∂xxx2+y2=y2−x2(x2+y2)2=∂y−yx2+y2\partial_x\tfrac{x}{x^2 + y^2} = \tfrac{y^2 - x^2}{(x^2 + y^2)^2} = \partial_y\tfrac{-y}{x^2 + y^2} ✓. Na okręgu F⃗=1rθ^\vec F = \tfrac1r\hat\theta, a dr⃗=r dθ θ^\mathrm d\vec r = r\,\mathrm d\theta\,\hat\theta, więc cyrkulacja to ∫02πdθ=2π\int_0^{2\pi}\mathrm d\theta = 2\pi, niezależnie od promienia.

  5. Energia potencjalna to U=x2+xy+y2U = x^2 + xy + y^2. Oblicz siłę i pracę tej siły od (1; 0)(1;\ 0) do (0; 1)(0;\ 1).

    Odpowiedź

    F⃗=−∇U=[−2x−y; −x−2y]\vec F = -\nabla U = [-2x - y;\ -x - 2y]. Praca =U(A)−U(B)=1−1=0= U(A) - U(B) = 1 - 1 = 0.

  6. Pokaż, że pole stałe F⃗=[Fx; Fy; Fz]\vec F = [F_x;\ F_y;\ F_z] jest potencjalne i wywnioskuj stąd, że praca grawitacji w pobliżu Ziemi to −mgΔz-mg\Delta z po każdej drodze. Dlaczego zakosy na szlaku nic nie zmieniają w pracy, a jednak pomagają turyście?

    Odpowiedź

    φ=Fxx+Fyy+Fzz\varphi = F_xx + F_yy + F_zz. Dla F⃗=[0; 0; −mg]\vec F = [0;\ 0;\ -mg]: praca −mgΔz-mg\Delta 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 P dx+Q dyP\,\mathrm dx + Q\,\mathrm dy nazywamy różniczką zupełną, jeśli jest różniczką pewnej funkcji: P dx+Q dy=dφP\,\mathrm dx + Q\,\mathrm dy = \mathrm d\varphi. To dokładnie warunek, że pole [P; Q][P;\ Q] jest potencjalne, więc test jest ten sam: ∂yP=∂xQ\partial_yP = \partial_xQ (na obszarze bez dziur). Całka z różniczki zupełnej po drodze to przyrost φ\varphi, niezależny od drogi. Gdy test nie zachodzi, piszemy często δ\delta zamiast d\mathrm d, żeby pamiętać, że wynik zależy od drogi.

4.1 Termodynamika: co jest funkcją stanu

Stan gazu doskonałego opisują np. temperatura TT i objętość VV. Pierwsza zasada termodynamiki (Dział 9: Termodynamika i entropia) mówi, że ciepło dostarczone przy małej, quasi-statycznej (odwracalnej) zmianie to δQ=CV dT+p dV\delta Q = C_V\,\mathrm dT + p\,\mathrm dV, gdzie p=nRT/Vp = nRT/V i CVC_V to pojemność cieplna przy stałej objętości (dla gazu doskonałego zależy co najwyżej od TT; dla jednoatomowego stała, 32nR\tfrac32nR). Zastosujmy test do δQ=CV dT+nRTV dV\delta Q = C_V\,\mathrm dT + \tfrac{nRT}{V}\,\mathrm dV:

∂CV∂V=0≠∂∂TnRTV=nRV.\frac{\partial C_V}{\partial V} = 0 \neq \frac{\partial}{\partial T}\frac{nRT}{V} = \frac{nR}{V} .

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 TT:

δQT=CVT dT+nRV dV,∂∂VCVT=0=∂∂TnRV.\frac{\delta Q}{T} = \frac{C_V}{T}\,\mathrm dT + \frac{nR}{V}\,\mathrm dV,\qquad \frac{\partial}{\partial V}\frac{C_V}{T} = 0 = \frac{\partial}{\partial T}\frac{nR}{V} .

To już różniczka zupełna funkcji S=CVln⁡T+nRln⁡V+constS = C_V\ln T + nR\ln V + \text{const}: entropii. Czynnik 1/T1/T, który zamienia δQ\delta 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/T1/T.

Przykład 4.1 Dwie drogi między tymi samymi stanami

Jeden mol jednoatomowego gazu doskonałego (CV=32RC_V = \tfrac32R, R=8,314 J/(mol K)R = 8{,}314\ \mathrm{J/(mol\,K)}) przechodzi ze stanu AA (300300 K, 1010 L) do stanu BB (400400 K, 2020 L). Droga 1: najpierw ogrzanie przy stałej objętości do 400400 K, potem rozprężanie izotermiczne. Droga 2: najpierw rozprężanie izotermiczne w 300300 K, potem ogrzanie przy stałej objętości. Oblicz QQ, WW, ΔU\Delta U i ΔS\Delta S.

  1. ΔU=CVΔT=32⋅8,314⋅100≈1247\Delta U = C_V\Delta T = \tfrac32\cdot8{,}314\cdot100 \approx 1247 J na obu drogach (energia wewnętrzna to funkcja stanu).
  2. Praca gazu przy rozprężaniu izotermicznym W=nRTln⁡VBVAW = nRT\ln\tfrac{V_B}{V_A}: droga 1 w 400400 K: W1≈2305W_1 \approx 2305 J; droga 2 w 300300 K: W2≈1729W_2 \approx 1729 J. Przy stałej objętości praca jest zerowa.
  3. Q=ΔU+WQ = \Delta U + W: Q1≈3552Q_1 \approx 3552 J, Q2≈2976Q_2 \approx 2976 J.
  4. ΔS=CVln⁡TBTA+nRln⁡VBVA=12,47⋅0,2877+8,314⋅0,6931≈9,35 J/K\Delta S = C_V\ln\tfrac{T_B}{T_A} + nR\ln\tfrac{V_B}{V_A} = 12{,}47\cdot0{,}2877 + 8{,}314\cdot0{,}6931 \approx 9{,}35\ \mathrm{J/K} na obu drogach.

Wynik: QQ i WW zależą od drogi (różnica 576576 J), a ΔU\Delta U i ΔS\Delta S nie.

Sprawdzenie: Q1−Q2=W1−W2=R⋅100⋅ln⁡2≈576Q_1 - Q_2 = W_1 - W_2 = R\cdot100\cdot\ln2 \approx 576 J ✓. Dla trzeciej drogi, prostej na wykresie p(V)p(V), numeryczna całka ∫δQ/T\int\delta Q/T też daje 9,359{,}35 J/K (symulacja 2).

Symulacja 2 Trzy drogi między stanami AA i BB na wykresie p(V)p(V). Pole pod krzywą to praca, różna dla każdej drogi. Odczyt pokazuje, że QQ i WW zależą od drogi, a ΔU\Delta U i ΔS=∫δQ/T\Delta S = \int\delta Q/T nie. Zmieniaj temperaturę i objętość stanu BB.

Zadania · część 4

5 zadań · ok. 35 min

  1. Czy 2xy dx+x2 dy2xy\,\mathrm dx + x^2\,\mathrm dy jest różniczką zupełną? A y dx−x dyy\,\mathrm dx - x\,\mathrm dy?

    Odpowiedź

    Tak: d(x2y)\mathrm d(x^2y). Nie: ∂y(y)=1≠∂x(−x)=−1\partial_y(y) = 1 \neq \partial_x(-x) = -1.

  2. Pokaż, że 1x2\tfrac{1}{x^2} jest czynnikiem całkującym dla y dx−x dyy\,\mathrm dx - x\,\mathrm dy (dla x≠0x \neq 0) i znajdź funkcję, której różniczką jest wynik.

    Odpowiedź

    yx2 dx−1x dy\tfrac{y}{x^2}\,\mathrm dx - \tfrac1x\,\mathrm dy: ∂yyx2=1x2=∂x(−1x)\partial_y\tfrac{y}{x^2} = \tfrac1{x^2} = \partial_x\left(-\tfrac1x\right) ✓. To d(−yx)\mathrm d\left(-\tfrac yx\right).

  3. Oblicz QQ i WW dla drogi z przykładu 4.1 biegnącej po prostej na wykresie p(V)p(V), z pA=249,4p_A = 249{,}4 kPa do pB=166,3p_B = 166{,}3 kPa.

    Odpowiedź

    WW to pole trapezu: 12(pA+pB)(VB−VA)≈12⋅415,7 kPa⋅0,01 m3=2078,5\tfrac12(p_A + p_B)(V_B - V_A) \approx \tfrac12\cdot415{,}7\ \mathrm{kPa}\cdot0{,}01\ \mathrm{m^3} = 2078{,}5 J.   Q=ΔU+W≈3326\;Q = \Delta U + W \approx 3326 J. Inne niż na drogach 1 i 2.

  4. Gaz przechodzi cykl zamknięty. Ile wynoszą ∮dU\oint\mathrm dU i ∮dS\oint\mathrm dS? Czy ∮δQ\oint\delta Q musi być zerowe? Co oznacza niezerowe ∮δQ\oint\delta Q?

    Odpowiedź

    ∮dU=∮dS=0\oint\mathrm dU = \oint\mathrm dS = 0 (funkcje stanu, cyrkulacja zerowa). ∮δQ\oint\delta Q nie musi być zerowe; z pierwszej zasady ∮δQ=∮δW\oint\delta Q = \oint\delta W, czyli pole cyklu na wykresie p(V)p(V): tak pracuje silnik cieplny.

  5. Pokaż, że dp=nRV dT−nRTV2 dV\mathrm dp = \tfrac{nR}{V}\,\mathrm dT - \tfrac{nRT}{V^2}\,\mathrm dV dla gazu doskonałego jest różniczką zupełną, i wyjaśnij, dlaczego musiało tak wyjść.

    Odpowiedź

    ∂VnRV=−nRV2=∂T(−nRTV2)\partial_V\tfrac{nR}{V} = -\tfrac{nR}{V^2} = \partial_T\left(-\tfrac{nRT}{V^2}\right) ✓. Musiało, bo p=nRT/Vp = 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][-y;\ x] daje na trzech drogach 00, 13\tfrac13 i 11 (przykład 2.1).
„Jeśli ∂xQ=∂yP\partial_xQ = \partial_yP, to pole jest potencjalne.”
Tylko na obszarze bez dziur. Wir [−y; x]x2+y2\tfrac{[-y;\ x]}{x^2 + y^2} spełnia test poza zerem, a ma cyrkulację 2π2\pi 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\delta 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⃗=∇φ\vec F = \nabla\varphi, fizycy F⃗=−∇U\vec F = -\nabla U. Siła pcha w stronę spadku energii potencjalnej.

6 Sprawdź się

1 Jak obliczyć ∫CF⃗⋅dr⃗\int_C\vec F\cdot\mathrm d\vec r dla krzywej r⃗(t)\vec r(t)?

∫abF⃗(r⃗(t))⋅r⃗ ′(t) dt\int_a^b\vec F(\vec r(t))\cdot\vec r\,'(t)\,\mathrm dt.

2 Ile wynosi praca pola ∇φ\nabla\varphi od AA do BB?

φ(B)−φ(A)\varphi(B) - \varphi(A), po każdej drodze.

3 Jaki test sprawdza, czy [P; Q][P;\ Q] jest potencjalne?

∂xQ=∂yP\partial_xQ = \partial_yP na obszarze jednospójnym.

4 Ile wynosi cyrkulacja pola potencjalnego?

Zero po każdej drodze zamkniętej.

5 Dlaczego wir [−y; x]x2+y2\tfrac{[-y;\ x]}{x^2 + y^2} nie jest potencjalny?

Bo ma cyrkulację 2π2\pi wokół zera; test nie obowiązuje, bo obszar ma dziurę w zerze.

6 Czy δQ\delta Q jest różniczką zupełną?

Nie. Dla procesów odwracalnych różniczką zupełną jest δQ/T=dS\delta Q/T = \mathrm dS.

7 Jaki jest związek siły z energią potencjalną?

F⃗=−∇U\vec F = -\nabla U.

7 Karta wzorów i co dalej

Całki krzywoliniowe

∫CF⃗⋅dr⃗=∫abF⃗(r⃗(t))⋅r⃗ ′(t) dt=∫CP dx+Q dy+R dz\int_C\vec F\cdot\mathrm d\vec r = \int_a^b\vec F(\vec r(t))\cdot\vec r\,'(t)\,\mathrm dt = \int_C P\,\mathrm dx + Q\,\mathrm dy + R\,\mathrm dz

zmiana kierunku: znak minus; parametryzacja: bez znaczenia

∫Cf ds\int_C f\,\mathrm ds,   ds=∣r⃗ ′∣ dt\;\mathrm ds = |\vec r\,'|\,\mathrm dt (masa drutu, długość)

cyrkulacja: ∮CF⃗⋅dr⃗\oint_C\vec F\cdot\mathrm d\vec r

Pola potencjalne i różniczki

∫C∇φ⋅dr⃗=φ(B)−φ(A)\int_C\nabla\varphi\cdot\mathrm d\vec r = \varphi(B) - \varphi(A); fizyka: F⃗=−∇U\vec F = -\nabla U

potencjalne   ⟺  \iff niezależność od drogi   ⟺  \iff zerowe cyrkulacje

test 2D: ∂xQ=∂yP\partial_xQ = \partial_yP; 3D: ∇×F⃗=0⃗\nabla\times\vec F = \vec 0; obszar jednospójny!

potencjał: całkuj φx=P\varphi_x = P po xx, dobierz g(y, z)g(y,\ z)

P dx+Q dy=dφP\,\mathrm dx + Q\,\mathrm dy = \mathrm d\varphi (zupełna)   ⟺  \iff ten sam test

δQ\delta Q nie zupełna;   δQ/T=dS\;\delta Q/T = \mathrm dS (proces odwracalny; czynnik całkujący 1/T1/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\vec F = -\nabla 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.

Cyrkulację po małej pętli zamienimy w lekcji 3.5 Strumień, dywergencja i rotacja w rotację, a strumień przez małą zamkniętą powierzchnię w dywergencję. Twierdzenia Greena i Stokesa (3.6 Twierdzenia Greena, Stokesa i Gaussa) powiedzą, że cyrkulacja po brzegu to suma rotacji wewnątrz, co wyjaśni, dlaczego dziura w obszarze tak wszystko zmienia. Różniczki zupełne wrócą w termodynamice (Dział 9: Termodynamika i entropia), a pola potencjalne w elektrostatyce (Dział 14: Potencjał, przewodniki i kondensatory).

Korzystają z tej lekcji 3.5 Strumień, dywergencja i rotacjaFizyka 9 Termodynamika i entropia