O005 · C120 · Bab 12

Penyegaran: Kalkulus Vektor

Joceline Lega · Edisi Bahasa Indonesia

Informasi di bawah ini merupakan ikhtisar elementer mengenai sifat-sifat dasar gradien suatu fungsi, serta divergensi dan rotasi (curl) suatu medan vektor. Topik-topik ini biasanya dibahas dalam mata kuliah kalkulus semester ketiga.

Gradien

Gradien suatu fungsi tiga variabel yang dapat didiferensialkan f(x,y,z)f(x,y,z) didefinisikan oleh

gradf=f=fxı+fyȷ+fzk\text{grad}f=\vec \nabla f = \displaystyle \frac{\partial f}{\partial x} \vec \imath + \frac{\partial f}{\partial y} \vec \jmath + \frac{\partial f}{\partial z} \vec k

dan memiliki sifat-sifat berikut.

Integral garis dan medan gradien

Integral garis

Integral garis suatu medan vektor kontinu F\vec F di sepanjang lintasan berorientasi 𝒞\mathcal C dituliskan sebagai

𝒞F(r)dr,\displaystyle \int_{\mathcal C} \vec F(\vec r) \cdot d\vec r,

dengan r\vec r sebagai vektor posisi dan drd\vec r merupakan vektor tangen terhadap lintasan 𝒞\mathcal C. Jika kita mengetahui suatu parametrisasi 𝒞\mathcal C, yaitu jika 𝒞\mathcal C ditelusuri oleh r(t)\vec r(t) saat tt bervariasi antara t0t_0 dan t1t_1, integral garis di atas menjadi

𝒞Fdr=t0t1F(r(t))drdtdt,\displaystyle \int_{\mathcal C} \vec F \cdot d\vec r = \int_{t_0}^{t_1} \vec F(\vec r(t)) \cdot \frac{d \vec r}{d t}\ dt,

sehingga dapat dihitung.

Contoh 1

Perhatikan medan vektor F(x,y)=exı+eyȷ\vec F(x,y) = e^x \vec \imath + e^y \vec \jmath. Misalkan 𝒞\mathcal C merupakan bagian elips x2+4y2=4x^2+4 y^2 = 4 yang menghubungkan (0,1)(0,1) ke (2,0)(2,0) searah jarum jam. Suatu parametrisasi 𝒞\mathcal C adalah x=2sin(t)x = 2 \sin(t), y=cos(t)y=\cos(t), dengan tt bervariasi antara 0 dan π/2\pi/2. Oleh karena itu,

𝒞Fdr=0π/2[e2sin(t)ı+ecos(t)ȷ][2cos(t)ısin(t)ȷ]dt=0π/2[2cos(t)e2sin(t)sin(t)ecos(t)]dt=[e2sin(t)+ecos(t)]0π/2=e(e1).\begin{align*} \int_{\mathcal C} \vec F \cdot d\vec r & = \int_0^{\pi/2} \left[e^{2 \sin(t)} \vec \imath + e^{\cos(t)} \vec \jmath \right] \cdot \left[2 \cos(t) \vec \imath - \sin(t) \vec \jmath \right] dt \\ & = \int_0^{\pi/2} \left[2 \cos(t) e^{2 \sin(t)} - \sin(t) e^{\cos(t)} \right]\ dt \\ & = \left[ e^{2 \sin(t)} + e^{\cos(t)} \right]_0^{\pi/2} = e (e -1).\end{align*}

Contoh 2

Usaha yang dilakukan oleh gaya F\vec F di sepanjang lintasan 𝒞\mathcal C diberikan oleh W=𝒞FdrW = \int_{\mathcal C} \vec F \cdot d\vec r. Jika r\vec r bergerak menurut hukum Newton, F=md2r/dt2\vec F = m \ d^2 \vec r / d t^2, maka

W=t0t1md2rdt2drdtdt=[12mdrdt2]t0t1=[12mv2]t0t1,\displaystyle W = \int_{t_0}^{t_1} m\ \frac{d^2 \vec r}{d t^2} \cdot \frac{d \vec r}{d t}\ dt = \left[\frac{1}{2} m \left\|\frac{d \vec r}{d t}\right\|^2 \right]_{t_0}^{t_1} = \left[ \frac{1}{2} m v^2 \right]_{t_0}^{t_1},

di mana vv menyatakan kelajuan; dengan demikian, usaha yang dilakukan oleh F\vec F ketika suatu massa titik bergerak di sepanjang lintasan 𝒞\mathcal C sama dengan selisih energi kinetik massa titik tersebut di antara kedua titik ujung 𝒞\mathcal C.

Medan gradien

Suatu medan vektor F\vec F disebut medan gradien jika terdapat fungsi f(x,y,z)f(x,y,z) sedemikian sehingga F=f\vec F = \vec \nabla f. Medan vektor F\vec F tersebut bersifat konservatif, tidak bergantung pada lintasan, dan bebas sirkulasi.

Menurut teorema dasar kalkulus untuk integral garis, suatu medan gradien tidak bergantung pada lintasan, yaitu

𝒞fdr=f(Q)f(P),\displaystyle \int_{\mathcal C} \vec \nabla f \cdot d \vec r = f(Q) - f(P),

dengan 𝒞\mathcal C sebagai lintasan yang mulus sepotong demi sepotong dan menghubungkan PP ke QQ, serta f\vec \nabla f kontinu pada 𝒞\mathcal C. Sebaliknya, setiap medan vektor yang tidak bergantung pada lintasan merupakan medan gradien.

Contoh

Perhatikan gaya F=md2rdt2\vec F = m \frac{d^2 \vec r}{d t^2} yang dibahas di atas. Jika F\vec F memiliki fungsi potensial, yaitu jika terdapat fungsi V(r)V(\vec r) sedemikian sehingga F=V\vec F = - \vec \nabla V, maka

W=𝒞Fdr=𝒞Vdr=V(r(t0))V(r(t1))=[12mv2]t0t1.\begin{align*} W & = \int_{\mathcal C} \vec F \cdot d\vec r = - \int_{\mathcal C} \vec \nabla V \cdot d\vec r = V(\vec r(t_0)) - V(\vec r(t_1)) \\ & = \left[ \frac{1}{2} m v^2 \right]_{t_0}^{t_1}. \end{align*}

Dengan demikian, 12mv2(t0)+V(r(t0))=12mv2(t1)+V(r(t1)),\displaystyle \frac{1}{2} m v^2(t_0) + V \left(\vec r(t_0) \right) = \frac{1}{2} m v^2(t_1) + V \left(\vec r(t_1) \right), yaitu, energi total, yang merupakan jumlah energi kinetik dan energi potensial, bersifat kekal.

F1y=F2x,F2z=F3y,F3x=F1z,\displaystyle \frac{\partial F_1}{\partial y} = \frac{\partial F_2}{\partial x}, \qquad \frac{\partial F_2}{\partial z} = \frac{\partial F_3}{\partial y}, \qquad \frac{\partial F_3}{\partial x} = \frac{\partial F_1}{\partial z},

yaitu, curlF=0\text{curl} \vec F = \vec 0.

Rotasi (curl)

Rotasi suatu medan vektor F=F1ı+F2ȷ+F3k\vec F = F_1 \vec \imath+ F_2 \vec \jmath + F_3 \vec k yang memiliki turunan parsial kontinu adalah vektor yang diberikan oleh

curlF=×F=|ıȷkxyzF1F2F3|=(F3yF2z)ı+(F1zF3x)ȷ+(F2xF1y)k.\begin{align*} \text{curl} \vec F & = \vec \nabla \times \vec F = \left \vert \begin{array}{ccc} \vec \imath & \vec \jmath & \vec k \\ \frac{\partial}{\partial x} & \frac{\partial}{\partial y} & \frac{\partial}{\partial z} \\ F_1 & F_2 & F_3 \end{array} \right \vert \\ & = \left( \frac{\partial F_3}{\partial y} - \frac{\partial F_2}{\partial z} \right) \vec \imath+ \left( \frac{\partial F_1}{\partial z} - \frac{\partial F_3} {\partial x}\right) \vec \jmath + \left( \frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y}\right) \vec k. \end{align*}

Rotasi memiliki sifat-sifat berikut.

curl gradf=×f=0.\text{curl grad} f = \vec \nabla \times \vec \nabla f = \vec 0.

Teorema Stokes dan Teorema Green

Fluks medan vektor F\vec F melalui permukaan berorientasi SS adalah SFdS\iint_S \vec F \cdot d\vec S.

SFdS=RF(x,y,f(x,y))[fxıfyȷ+k]dxdy,\displaystyle \iint_S \vec F \cdot d\vec S = \iint_R \vec F\left(x,y,f(x,y)\right)\cdot \left[-\frac{\partial f}{\partial x} \vec \imath - \frac{\partial f}{\partial y} \vec \jmath + \vec k \right] dx dy,

SFdS=RF(r(s,t))(rs×rt)dsdt.\displaystyle \int_S \vec F \cdot d\vec S = \int_R \vec F\left(\vec r(s,t)\right) \cdot \left(\frac{\partial \vec r}{\partial s} \times \frac{\partial \vec r} {\partial t} \right) ds dt.

Teorema Stokes

Teorema Stokes menghubungkan sirkulasi suatu medan vektor mulus di sepanjang kurva tertutup 𝒞\mathcal C yang mulus sepotong demi sepotong dengan fluks curlF\text{curl} \vec F melalui sebarang permukaan mulus SS yang batasnya sama dengan 𝒞\mathcal C. Teorema ini berbunyi:

𝒞Fdr=ScurlFdA,\displaystyle \int_{\mathcal C} \vec F \cdot d\vec r = \int_S \text{curl}\vec F \cdot d\vec A,

di mana orientasi 𝒞\mathcal C ditentukan dari orientasi SS (atau sebaliknya) dengan kaidah tangan kanan. Perhatikan bahwa Teorema Stokes hanya berlaku jika F\vec F didefinisikan di setiap titik pada SS dan 𝒞\mathcal C.

Teorema Green

Teorema Green adalah versi planar dari Teorema Stokes dan berbunyi

𝒞Fdr=R(F2xF1y)dxdy,\displaystyle \int_{\mathcal C} \vec F \cdot d\vec r = \int_R \left(\frac{\partial F_2}{\partial x} - \frac{\partial F_1}{\partial y} \right) \ dx\ dy,

di mana F(x,y)=F1(x,y)ı+F2(x,y)ȷ\vec F(x,y) = F_1(x,y) \vec \imath + F_2(x,y) \vec \jmath adalah medan vektor mulus yang didefinisikan di setiap titik 𝒞\mathcal C serta di bagian dalam 𝒞\mathcal C, dan 𝒞\mathcal C adalah kurva tertutup sederhana yang mulus sepotong demi sepotong serta diorientasikan sedemikian rupa sehingga bagian dalamnya, RR, berada di sebelah kiri ketika kita bergerak di sepanjang 𝒞\mathcal C.

Teorema divergensi

Divergensi suatu medan vektor F=F1ı+F2ȷ+F3k\vec F = F_1 \vec \imath+ F_2 \vec \jmath + F_3 \vec k yang memiliki turunan parsial kontinu diberikan oleh

divF=F=F1x+F2y+F3z.\text{div} \vec F = \vec \nabla \cdot \vec F = \frac{\partial F_1}{\partial x} + \frac{\partial F_2} {\partial y} + \frac{\partial F_3}{\partial z}.

Perhatikan bahwa besaran ini adalah besaran skalar. Divergensi memiliki sifat-sifat berikut.

divF(P)=limε0Bε(P)FndSvol(Bε(P)).\text{div} \vec F(P) = \displaystyle \lim_{\varepsilon \rightarrow 0} \frac{\iint_{\partial B_\varepsilon(P)} \vec F \cdot \vec n\, dS}{\operatorname{vol}(B_\varepsilon(P))}.

div curlF=(×F)=0.\text{div curl} \vec F = \vec \nabla \cdot (\vec \nabla \times \vec F) = 0.

Menurut teorema divergensi, fluks suatu medan vektor mulus melalui permukaan tertutup SS yang mulus sepotong demi sepotong dan diorientasikan ke luar sama dengan integral volume divergensinya di seluruh daerah WW yang dibatasi oleh SS; yaitu,

SFdA=WdivFdV,\displaystyle \int_S \vec F \cdot d \vec A \ = \int_W \text{div}\vec F \ dV,

dengan asumsi bahwa F\vec F didefinisikan di setiap titik di dalam WW dan pada SS.

Contoh: persamaan kontinuitas

Perhatikan aliran v(r)\vec v(\vec r) suatu fluida dengan kerapatan ρ(r,t)\rho(\vec r,t). Misalkan VV adalah suatu daerah tetap dalam domain fluida. Massa fluida di dalam VV diberikan oleh m=VρdV\displaystyle m = \int_V \rho dV, dan laju perubahannya adalah

dmdt=VρtdV,\displaystyle \frac{d m}{d t} = \int_V \frac{\partial \rho}{\partial t}\ dV,

jika kita mengasumsikan bahwa ρ\rho mulus. Di sisi lain, laju perubahan ini sama dengan lawan dari fluks massa ρv\rho \vec v melalui batas SS dari VV, dan berdasarkan teorema divergensi,

dmdt=SρvdA=Vdiv(ρv)dV.\displaystyle \frac{d m}{d t} = - \int_S \rho\ \vec v \cdot d\vec A = - \int_V \text{div}(\rho\ \vec v)\ dV.

Dengan membandingkan kedua ungkapan dm/dt,dm/dt, yang sama untuk setiap daerah VV, kita memperoleh persamaan kontinuitas, yang berbunyi

ρt+div(ρv)=0.\displaystyle \frac{\partial \rho}{\partial t} + \text{div}(\rho\ \vec v) = 0.