Bab 1

Analisis Vektor

Kalkulus vektor tiga dimensi — operator diferensial nabla, teorema integral fundamental, dan sistem koordinat kurvilinier yang menjadi bahasa utama elektrodinamika klasik.

1.1 Vektor dalam Tiga Dimensi

Definisi 1.1 — Vektor

Vektor $\mathbf{A}$ dalam ruang tiga dimensi adalah besaran yang memiliki besaran (panjang) dan arah, direpresentasikan sebagai tripel terurut komponen terhadap basis ortonormal $\{\hat{\mathbf{x}},\,\hat{\mathbf{y}},\,\hat{\mathbf{z}}\}$:

$$ \mathbf{A} = A_x\,\hat{\mathbf{x}} + A_y\,\hat{\mathbf{y}} + A_z\,\hat{\mathbf{z}} $$

Notasi singkat: $\mathbf{A} = (A_x, A_y, A_z)$. Komponen-komponen dapat diperoleh melalui proyeksi: $A_x = \mathbf{A}\cdot\hat{\mathbf{x}}$, $A_y = \mathbf{A}\cdot\hat{\mathbf{y}}$, $A_z = \mathbf{A}\cdot\hat{\mathbf{z}}$.

Besar (norma) vektor didefinisikan melalui teorema Pythagoras generalisasi:

$$ |\mathbf{A}| = A = \sqrt{A_x^2 + A_y^2 + A_z^2} $$

Vektor satuan searah $\mathbf{A}$ diperoleh dengan membagi vektor oleh besarnya:

$$ \hat{\mathbf{A}} = \frac{\mathbf{A}}{|\mathbf{A}|} \notag \qquad\text{sehingga}\qquad |\hat{\mathbf{A}}| = 1 $$

1.1.1 Operasi Aljabar Vektor

Penjumlahan dan pengurangan vektor dilakukan komponen per komponen. Jika $\mathbf{A} = (A_x, A_y, A_z)$ dan $\mathbf{B} = (B_x, B_y, B_z)$:

$$ \mathbf{A} \pm \mathbf{B} = (A_x \pm B_x)\,\hat{\mathbf{x}} + (A_y \pm B_y)\,\hat{\mathbf{y}} + (A_z \pm B_z)\,\hat{\mathbf{z}} $$

Perkalian vektor dengan skalar $\alpha$:

$$ \alpha\,\mathbf{A} = \alpha A_x\,\hat{\mathbf{x}} + \alpha A_y\,\hat{\mathbf{y}} + \alpha A_z\,\hat{\mathbf{z}} \notag \qquad\text{dan}\qquad |\alpha\,\mathbf{A}| = |\alpha|\,|\mathbf{A}| $$

Geometri penjumlahan vektor memenuhi Hukum Jajaran Genjang: vektor $\mathbf{A}+\mathbf{B}$ merupakan diagonal jajaran genjang yang dibentuk oleh $\mathbf{A}$ dan $\mathbf{B}$ ketika keduanya ditempatkan ekor-ke-ekor.

z y x A A_z A_y A_x
Gambar 1.1 — Dekomposisi vektor $\mathbf{A}$ ke dalam komponen $A_x$, $A_y$, $A_z$ pada basis kartesius. Garis putus-putus menunjukkan proyeksi paralel terhadap sumbu.
Contoh 1.1

Diberikan $\mathbf{A} = 2\hat{\mathbf{x}} - \hat{\mathbf{y}} + 3\hat{\mathbf{z}}$. Hitung besar vektor dan vektor satuan $\hat{\mathbf{A}}$.

Besaran: $|\mathbf{A}| = \sqrt{2^2 + (-1)^2 + 3^2} = \sqrt{4+1+9} = \sqrt{14} \approx 3.742 \notag$

Vektor satuan: $\hat{\mathbf{A}} = \dfrac{2\hat{\mathbf{x}} - \hat{\mathbf{y}} + 3\hat{\mathbf{z}}}{\sqrt{14}} \approx 0.535\,\hat{\mathbf{x}} - 0.267\,\hat{\mathbf{y}} + 0.802\,\hat{\mathbf{z}} \notag$

Contoh 1.2 — Penjumlahan Vektor dan Jarak

Dua titik dalam ruang memiliki vektor posisi $\mathbf{r}_1 = \hat{\mathbf{x}} + 2\hat{\mathbf{y}} - \hat{\mathbf{z}}$ dan $\mathbf{r}_2 = 3\hat{\mathbf{x}} - \hat{\mathbf{y}} + 4\hat{\mathbf{z}}$. Hitung vektor perpindahan dari $P_1$ ke $P_2$ dan jarak antara kedua titik.

Vektor perpindahan: $\Delta\mathbf{r} = \mathbf{r}_2 - \mathbf{r}_1 = (3-1)\hat{\mathbf{x}} + (-1-2)\hat{\mathbf{y}} + (4-(-1))\hat{\mathbf{z}} = 2\hat{\mathbf{x}} - 3\hat{\mathbf{y}} + 5\hat{\mathbf{z}} \notag$

Jarak: $|\Delta\mathbf{r}| = \sqrt{4 + 9 + 25} = \sqrt{38} \approx 6.164 \notag$

Contoh 1.3 — Vektor Posisi dalam Koordinat Bola

Tunjukkan bahwa vektor satuan radial dalam koordinat bola dapat ditulis sebagai $\hat{\mathbf{r}} = \sin\theta\cos\phi\,\hat{\mathbf{x}} + \sin\theta\sin\phi\,\hat{\mathbf{y}} + \cos\theta\,\hat{\mathbf{z}}$.

Vektor posisi dalam Kartesius: $\mathbf{r} = x\hat{\mathbf{x}} + y\hat{\mathbf{y}} + z\hat{\mathbf{z}} = r\sin\theta\cos\phi\,\hat{\mathbf{x}} + r\sin\theta\sin\phi\,\hat{\mathbf{y}} + r\cos\theta\,\hat{\mathbf{z}} \notag$

Karena $\hat{\mathbf{r}} = \mathbf{r}/r$, maka: $\hat{\mathbf{r}} = \sin\theta\cos\phi\,\hat{\mathbf{x}} + \sin\theta\sin\phi\,\hat{\mathbf{y}} + \cos\theta\,\hat{\mathbf{z}} \notag$

Verifikasi: $|\hat{\mathbf{r}}| = \sqrt{\sin^2\theta\cos^2\phi + \sin^2\theta\sin^2\phi + \cos^2\theta} = \sqrt{\sin^2\theta(\cos^2\phi+\sin^2\phi) + \cos^2\theta} = \sqrt{\sin^2\theta + \cos^2\theta} = 1 \;\checkmark \notag$

Vektor Posisi dan Notasi Penting

Dalam elektrodinamika, vektor posisi $\mathbf{r} = x\hat{\mathbf{x}} + y\hat{\mathbf{y}} + z\hat{\mathbf{z}}$ dan besarannya $r = |\mathbf{r}| = \sqrt{x^2+y^2+z^2}$ muncul terus-menerus. Vektor pemisah dari titik sumber $\mathbf{r}'$ ke titik pengamatan $\mathbf{r}$ dinotasikan $\boldsymbol{\mathscr{r}} = \mathbf{r} - \mathbf{r}'$ dengan besar $\mathscr{r} = |\mathbf{r}-\mathbf{r}'|$. Konvensi ini sangat penting dalam rumus integral Coulomb dan hukum Biot–Savart.

1.2 Perkalian Titik dan Perkalian Silang

1.2.1 Produk Titik (Dot Product)

Definisi 1.2 — Produk Titik

Produk titik (hasil kali dalam) dua vektor didefinisikan secara aljabar sebagai:

$$ \mathbf{A} \cdot \mathbf{B} = A_x B_x + A_y B_y + A_z B_z $$

Secara geometris, produk titik memberikan proyeksi skalar $\mathbf{A}$ pada $\mathbf{B}$ dikalikan besar $\mathbf{B}$:

$$ \mathbf{A} \cdot \mathbf{B} = |\mathbf{A}|\,|\mathbf{B}|\,\cos\theta $$

di mana $\theta$ adalah sudut terkecil antara kedua vektor ($0 \leq \theta \leq \pi$).

Komponen $\mathbf{A}$ searah $\mathbf{B}$ (proyeksi skalar) dan vektor proyeksi diberikan oleh:

$$ \text{comp}_{\mathbf{B}} \mathbf{A} = \frac{\mathbf{A} \cdot \mathbf{B}}{|\mathbf{B}|} = |\mathbf{A}|\cos\theta \notag \qquad\qquad \text{proj}_{\mathbf{B}} \mathbf{A} = \frac{\mathbf{A} \cdot \mathbf{B}}{|\mathbf{B}|^2}\,\mathbf{B} $$
B A θ |A|cos θ
Gambar 1.2 — Interpretasi geometris produk titik: proyeksi skalar $\mathbf{A}$ pada $\mathbf{B}$ bernilai $|\mathbf{A}|\cos\theta$. Jika $\theta = 90°$, proyeksi nol — kedua vektor ortogonal.

Sifat-sifat penting produk titik: komutatif ($\mathbf{A}\cdot\mathbf{B} = \mathbf{B}\cdot\mathbf{A}$), distributif terhadap penjumlahan, dan $\mathbf{A}\cdot\mathbf{A} = |\mathbf{A}|^2$. Dua vektor ortogonal jika dan hanya jika $\mathbf{A}\cdot\mathbf{B} = 0$.

Produk titik basis ortonormal:

$$ \hat{\mathbf{x}}\cdot\hat{\mathbf{x}} = \hat{\mathbf{y}}\cdot\hat{\mathbf{y}} = \hat{\mathbf{z}}\cdot\hat{\mathbf{z}} = 1 \notag \qquad\qquad \hat{\mathbf{x}}\cdot\hat{\mathbf{y}} = \hat{\mathbf{y}}\cdot\hat{\mathbf{z}} = \hat{\mathbf{z}}\cdot\hat{\mathbf{x}} = 0 $$
Contoh 1.4 — Menentukan Sudut antara Dua Vektor

Diberikan $\mathbf{A} = \hat{\mathbf{x}} + \hat{\mathbf{y}} + \hat{\mathbf{z}}$ dan $\mathbf{B} = 2\hat{\mathbf{x}} - \hat{\mathbf{z}}$. Tentukan sudut $\theta$ antara keduanya.

Produk titik: $\mathbf{A}\cdot\mathbf{B} = (1)(2) + (1)(0) + (1)(-1) = 1 \notag$

Besaran: $|\mathbf{A}| = \sqrt{1+1+1} = \sqrt{3}$, $\quad|\mathbf{B}| = \sqrt{4+0+1} = \sqrt{5} \notag$

Sudut: $\cos\theta = \dfrac{\mathbf{A}\cdot\mathbf{B}}{|\mathbf{A}|\,|\mathbf{B}|} = \dfrac{1}{\sqrt{15}} \approx 0.2582 \notag$

Sehingga $\theta = \arccos(1/\sqrt{15}) \approx 75.0° \notag$

Contoh 1.5 — Proyeksi Vektor dan Komponen Sejajar/Tegak Lurus

Dengan $\mathbf{A} = 3\hat{\mathbf{x}} + 2\hat{\mathbf{y}} - \hat{\mathbf{z}}$ dan $\mathbf{B} = \hat{\mathbf{x}} + \hat{\mathbf{y}} + 2\hat{\mathbf{z}}$, tentukan komponen $\mathbf{A}$ yang sejajar dan tegak lurus terhadap $\mathbf{B}$.

Produk titik: $\mathbf{A}\cdot\mathbf{B} = 3 + 2 - 2 = 3$, $\quad|\mathbf{B}|^2 = 1+1+4 = 6 \notag$

Komponen sejajar: $\mathbf{A}_{\parallel} = \text{proj}_{\mathbf{B}}\mathbf{A} = \dfrac{3}{6}\,\mathbf{B} = \dfrac{1}{2}\hat{\mathbf{x}} + \dfrac{1}{2}\hat{\mathbf{y}} + \hat{\mathbf{z}} \notag$

Komponen tegak lurus: $\mathbf{A}_{\perp} = \mathbf{A} - \mathbf{A}_{\parallel} = \dfrac{5}{2}\hat{\mathbf{x}} + \dfrac{3}{2}\hat{\mathbf{y}} - 2\hat{\mathbf{z}} \notag$

Verifikasi: $\mathbf{A}_{\parallel}\cdot\mathbf{A}_{\perp} = \frac{1}{2}\cdot\frac{5}{2} + \frac{1}{2}\cdot\frac{3}{2} + (1)(-2) = \frac{5}{4} + \frac{3}{4} - 2 = 0 \;\checkmark \notag$

Dan $\mathbf{A}_{\parallel} + \mathbf{A}_{\perp} = \mathbf{A} \;\checkmark \notag$

1.2.2 Produk Silang (Cross Product)

Definisi 1.3 — Produk Silang

Produk silang $\mathbf{A}\times\mathbf{B}$ menghasilkan vektor yang tegak lurus terhadap bidang yang dibentuk $\mathbf{A}$ dan $\mathbf{B}$, dengan arah ditentukan oleh aturan tangan kanan:

$$ \mathbf{A} \times \mathbf{B} = \begin{vmatrix} \hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ A_x & A_y & A_z \\ B_x & B_y & B_z \end{vmatrix} $$

Besarnya sama dengan luas jajar genjang yang dibentuk oleh kedua vektor:

$$ |\mathbf{A} \times \mathbf{B}| = |\mathbf{A}|\,|\mathbf{B}|\,\sin\theta $$
Perhatian

Produk silang tidak komutatif: $\mathbf{A}\times\mathbf{B} = -\mathbf{B}\times\mathbf{A}$. Arah hasil silang mengikuti aturan tangan kanan — putar dari $\mathbf{A}$ ke $\mathbf{B}$ melalui sudut terkecil $\theta$, maka ibu jari menunjuk arah $\mathbf{A}\times\mathbf{B}$.

Produk silang basis ortonormal (siklus $\hat{\mathbf{x}} \to \hat{\mathbf{y}} \to \hat{\mathbf{z}} \to \hat{\mathbf{x}}$):

$$ \hat{\mathbf{x}}\times\hat{\mathbf{y}} = \hat{\mathbf{z}}, \quad \hat{\mathbf{y}}\times\hat{\mathbf{z}} = \hat{\mathbf{x}}, \quad \hat{\mathbf{z}}\times\hat{\mathbf{x}} = \hat{\mathbf{y}} \notag \qquad \hat{\mathbf{x}}\times\hat{\mathbf{x}} = \hat{\mathbf{y}}\times\hat{\mathbf{y}} = \hat{\mathbf{z}}\times\hat{\mathbf{z}} = \mathbf{0} $$
A B A×B θ |A||B|sin θ
Gambar 1.3 — Produk silang $\mathbf{A}\times\mathbf{B}$ tegak lurus terhadap bidang $\mathbf{A}$–$\mathbf{B}$, dengan besar sama dengan luas jajar genjang. Arah ditentukan oleh aturan tangan kanan.
Contoh 1.6 — Produk Silang dan Luas Segitiga

Tiga titik membentuk segitiga: $P_1 = (0,0,0)$, $P_2 = (2,1,0)$, $P_3 = (1,3,2)$. Hitung luas segitiga ini.

Vektor sisi: $\mathbf{a} = P_2 - P_1 = 2\hat{\mathbf{x}} + \hat{\mathbf{y}}$, $\quad\mathbf{b} = P_3 - P_1 = \hat{\mathbf{x}} + 3\hat{\mathbf{y}} + 2\hat{\mathbf{z}} \notag$

Produk silang: $\mathbf{a}\times\mathbf{b} = \begin{vmatrix}\hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ 2 & 1 & 0 \\ 1 & 3 & 2\end{vmatrix} = (2-0)\hat{\mathbf{x}} - (4-0)\hat{\mathbf{y}} + (6-1)\hat{\mathbf{z}} = 2\hat{\mathbf{x}} - 4\hat{\mathbf{y}} + 5\hat{\mathbf{z}} \notag$

$|\mathbf{a}\times\mathbf{b}| = \sqrt{4+16+25} = \sqrt{45} = 3\sqrt{5} \notag$

Luas segitiga $= \dfrac{1}{2}|\mathbf{a}\times\mathbf{b}| = \dfrac{3\sqrt{5}}{2} \approx 3.354 \notag$

Contoh 1.7 — Verifikasi Ortogonalitas Hasil Silang

Dengan $\mathbf{A} = \hat{\mathbf{x}} + 2\hat{\mathbf{z}}$ dan $\mathbf{B} = 3\hat{\mathbf{x}} - \hat{\mathbf{y}} + \hat{\mathbf{z}}$, hitung $\mathbf{A}\times\mathbf{B}$ dan verifikasi bahwa hasilnya tegak lurus terhadap $\mathbf{A}$ dan $\mathbf{B}$.

$\mathbf{A}\times\mathbf{B} = \begin{vmatrix}\hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ 1 & 0 & 2 \\ 3 & -1 & 1\end{vmatrix} = (0\cdot1 - 2\cdot(-1))\hat{\mathbf{x}} - (1\cdot1 - 2\cdot3)\hat{\mathbf{y}} + (1\cdot(-1) - 0\cdot3)\hat{\mathbf{z}} = 2\hat{\mathbf{x}} + 5\hat{\mathbf{y}} - \hat{\mathbf{z}} \notag$

Verifikasi: $\mathbf{A}\cdot(\mathbf{A}\times\mathbf{B}) = (1)(2) + (0)(5) + (2)(-1) = 2 - 2 = 0 \;\checkmark \notag$

$\mathbf{B}\cdot(\mathbf{A}\times\mathbf{B}) = (3)(2) + (-1)(5) + (1)(-1) = 6 - 5 - 1 = 0 \;\checkmark \notag$

1.2.3 Produk Tripel

Produk tripel skalar (triple product) menghasilkan skalar yang besarnya sama dengan volume parallelepiped:

$$ \mathbf{A} \cdot (\mathbf{B} \times \mathbf{C}) = \begin{vmatrix} A_x & A_y & A_z \\ B_x & B_y & B_z \\ C_x & C_y & C_z \end{vmatrix} $$
Sifat Produk Tripel Skalar

Produk tripel skalar bersifat invarian terhadap perputaran titik dan silang: $\mathbf{A}\cdot(\mathbf{B}\times\mathbf{C}) = \mathbf{B}\cdot(\mathbf{C}\times\mathbf{A}) = \mathbf{C}\cdot(\mathbf{A}\times\mathbf{B})$. Nilainya positif jika $\{\mathbf{A},\mathbf{B},\mathbf{C}\}$ membentuk basis tangan kanan, negatif jika tangan kiri, dan nol jika ketiga vektor koplanar.

Produk tripel vektor memenuhi identitas BAC-CAB yang sangat berguna dalam derivasi elektrodinamika:

$$ \mathbf{A} \times (\mathbf{B} \times \mathbf{C}) = \mathbf{B}(\mathbf{A}\cdot\mathbf{C}) - \mathbf{C}(\mathbf{A}\cdot\mathbf{B}) $$
Contoh 1.8 — Volume Parallelepiped

Hitung volume parallelepiped yang dibentuk oleh $\mathbf{A} = \hat{\mathbf{x}} + \hat{\mathbf{y}}$, $\mathbf{B} = \hat{\mathbf{y}} + \hat{\mathbf{z}}$, $\mathbf{C} = \hat{\mathbf{z}} + \hat{\mathbf{x}}$.

$V = |\mathbf{A}\cdot(\mathbf{B}\times\mathbf{C})| = \left|\begin{vmatrix}1 & 1 & 0 \\ 0 & 1 & 1 \\ 1 & 0 & 1\end{vmatrix}\right| = |1(1\cdot1 - 1\cdot0) - 1(0\cdot1 - 1\cdot1) + 0| = |1 + 1| = 2 \notag$

Contoh 1.9 — Verifikasi Identitas BAC-CAB

Verifikasi identitas BAC-CAB untuk $\mathbf{A} = \hat{\mathbf{x}} + 2\hat{\mathbf{y}}$, $\mathbf{B} = 2\hat{\mathbf{x}} - \hat{\mathbf{z}}$, $\mathbf{C} = \hat{\mathbf{y}} + 3\hat{\mathbf{z}}$.

Ruang kiri — $\mathbf{A}\times(\mathbf{B}\times\mathbf{C})$:

$\mathbf{B}\times\mathbf{C} = \begin{vmatrix}\hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ 2 & 0 & -1 \\ 0 & 1 & 3\end{vmatrix} = (0\cdot3 - (-1)\cdot1)\hat{\mathbf{x}} - (2\cdot3 - (-1)\cdot0)\hat{\mathbf{y}} + (2\cdot1 - 0\cdot0)\hat{\mathbf{z}} = \hat{\mathbf{x}} - 6\hat{\mathbf{y}} + 2\hat{\mathbf{z}} \notag$

$\mathbf{A}\times(\mathbf{B}\times\mathbf{C}) = \begin{vmatrix}\hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ 1 & 2 & 0 \\ 1 & -6 & 2\end{vmatrix} = 4\hat{\mathbf{x}} - 2\hat{\mathbf{y}} - 8\hat{\mathbf{z}} \notag$

Ruang kanan — $\mathbf{B}(\mathbf{A}\cdot\mathbf{C}) - \mathbf{C}(\mathbf{A}\cdot\mathbf{B})$:

$\mathbf{A}\cdot\mathbf{C} = 0 + 2 + 0 = 2$, $\quad\mathbf{A}\cdot\mathbf{B} = 2 + 0 + 0 = 2 \notag$

$\mathbf{B}(2) - \mathbf{C}(2) = 4\hat{\mathbf{x}} - 2\hat{\mathbf{z}} - 2\hat{\mathbf{y}} - 6\hat{\mathbf{z}} = 4\hat{\mathbf{x}} - 2\hat{\mathbf{y}} - 8\hat{\mathbf{z}} \notag$

Kedua ruas sama: $4\hat{\mathbf{x}} - 2\hat{\mathbf{y}} - 8\hat{\mathbf{z}} \;\checkmark \notag$

Contoh 1.10 — Koplanaritas Vektor

Tentukan apakah $\mathbf{A} = 2\hat{\mathbf{x}} - \hat{\mathbf{y}} + \hat{\mathbf{z}}$, $\mathbf{B} = \hat{\mathbf{x}} + \hat{\mathbf{y}} - 2\hat{\mathbf{z}}$, dan $\mathbf{C} = 3\hat{\mathbf{x}} - 3\hat{\mathbf{y}} + 5\hat{\mathbf{z}}$ koplanar.

Tiga vektor koplanar jika dan hanya jika produk tripel skalar bernilai nol:

$\mathbf{A}\cdot(\mathbf{B}\times\mathbf{C}) = \begin{vmatrix}2 & -1 & 1 \\ 1 & 1 & -2 \\ 3 & -3 & 5\end{vmatrix} \notag$

$= 2(5-6) - (-1)(5-(-6)) + 1(-3-3) = 2(-1) + 1(11) + 1(-6) = -2+11-6 = 3 \neq 0 \notag$

Karena $\mathbf{A}\cdot(\mathbf{B}\times\mathbf{C}) = 3 \neq 0$, ketiga vektor tidak koplanar — mereka membentuk basis linear-independen di $\mathbb{R}^3$. $\notag$

1.2.4 Ringkasan Identitas Produk Vektor

Tabel Identitas Produk Vektor
Identitas Ekspresi
Komutatif titik $\mathbf{A}\cdot\mathbf{B} = \mathbf{B}\cdot\mathbf{A}$
Antikomutatif silang $\mathbf{A}\times\mathbf{B} = -\mathbf{B}\times\mathbf{A}$
Distributif titik $\mathbf{A}\cdot(\mathbf{B}+\mathbf{C}) = \mathbf{A}\cdot\mathbf{B} + \mathbf{A}\cdot\mathbf{C}$
Distributif silang $\mathbf{A}\times(\mathbf{B}+\mathbf{C}) = \mathbf{A}\times\mathbf{B} + \mathbf{A}\times\mathbf{C}$
Tripel skalar $\mathbf{A}\cdot(\mathbf{B}\times\mathbf{C}) = \mathbf{B}\cdot(\mathbf{C}\times\mathbf{A}) = \mathbf{C}\cdot(\mathbf{A}\times\mathbf{B})$
BAC-CAB $\mathbf{A}\times(\mathbf{B}\times\mathbf{C}) = \mathbf{B}(\mathbf{A}\cdot\mathbf{C}) - \mathbf{C}(\mathbf{A}\cdot\mathbf{B})$
Lagrange $|\mathbf{A}\times\mathbf{B}|^2 = |\mathbf{A}|^2|\mathbf{B}|^2 - (\mathbf{A}\cdot\mathbf{B})^2$
Jacobi $\mathbf{A}\times(\mathbf{B}\times\mathbf{C}) + \mathbf{B}\times(\mathbf{C}\times\mathbf{A}) + \mathbf{C}\times(\mathbf{A}\times\mathbf{B}) = \mathbf{0}$

Identitas Lagrange merupakan konsekuensi langsung dari $|\mathbf{A}\times\mathbf{B}| = |\mathbf{A}||\mathbf{B}|\sin\theta$ dan $\mathbf{A}\cdot\mathbf{B} = |\mathbf{A}||\mathbf{B}|\cos\theta$ dikombinasikan dengan $\sin^2\theta + \cos^2\theta = 1$. Identitas Jacobi berguna dalam pembuktian beberapa teorema tentang medan vektor.

1.3 Operator Diferensial: Gradient

Operator diferensial $\nabla$ (dibaca nabla) adalah vektor operator yang menjadi tulang punggung seluruh elektrodinamika:

$$ \nabla = \hat{\mathbf{x}}\frac{\partial}{\partial x} + \hat{\mathbf{y}}\frac{\partial}{\partial y} + \hat{\mathbf{z}}\frac{\partial}{\partial z} $$
Cara Membaca Operator Nabla

$\nabla$ bukan vektor biasa — ia adalah operator diferensial. Dalam ekspresi, $\nabla$ harus selalu "beroperasi pada" sesuatu yang ada di sebelah kanannya. Aturan praktis: treat $\nabla$ seperti vektor dalam aturan aljabar, tetapi selalu terapkan turunan parsial pada fungsi di sebelah kanannya. Jangan pernah menulis $\mathbf{A}\nabla$ tanpa konteks yang jelas.

Definisi 1.4 — Gradient

Gradient dari fungsi skalar $T(x,y,z)$ adalah vektor yang komponennya adalah turunan parsial $T$ terhadap masing-masing koordinat:

$$ \nabla T = \frac{\partial T}{\partial x}\,\hat{\mathbf{x}} + \frac{\partial T}{\partial y}\,\hat{\mathbf{y}} + \frac{\partial T}{\partial z}\,\hat{\mathbf{z}} $$
Teorema 1.1 — Sifat Fundamental Gradient

Gradient $\nabla T$ menunjuk ke arah kenaikan paling curam dari $T$, dan besarnya $|\nabla T|$ adalah laju kenaikan maksimum tersebut.

Turunan arah dari $T$ searah vektor satuan $\hat{\mathbf{u}}$ diberikan oleh:

$$ \frac{dT}{ds}\bigg|_{\hat{\mathbf{u}}} = \nabla T \cdot \hat{\mathbf{u}} = |\nabla T|\,\cos\alpha $$

di mana $\alpha$ adalah sudut antara $\nabla T$ dan $\hat{\mathbf{u}}$. Turunan arah maksimum terjadi ketika $\hat{\mathbf{u}}$ searah $\nabla T$ ($\alpha = 0$), dan nol ketika $\hat{\mathbf{u}}$ tegak lurus terhadap $\nabla T$ — yaitu, sepanjang garis kontur $T = \text{konstan}$.

T=3 T=2 T=1 ∇T tegak lurus kontur, menunjuk ke T lebih besar
Gambar 1.4 — Gradient $\nabla T$ selalu tegak lurus terhadap garis kontur $T = \text{konstan}$ dan menunjuk ke arah kenaikan $T$ paling curam.
Contoh 1.11 — Gradient dan Turunan Arah

Diberikan $T(x,y,z) = 3x^2y - y^2z$. Hitung $\nabla T$ dan turunan arah di titik $P = (1, -1, 2)$ searah vektor $\mathbf{v} = 2\hat{\mathbf{x}} + \hat{\mathbf{y}} - 2\hat{\mathbf{z}}$.

Gradient: $\nabla T = 6xy\,\hat{\mathbf{x}} + (3x^2 - 2yz)\,\hat{\mathbf{y}} - y^2\,\hat{\mathbf{z}} \notag$

Di titik $P$: $\nabla T\big|_P = 6(1)(-1)\hat{\mathbf{x}} + (3 - 2(-1)(2))\hat{\mathbf{y}} - (-1)^2\hat{\mathbf{z}} = -6\hat{\mathbf{x}} + 7\hat{\mathbf{y}} - \hat{\mathbf{z}} \notag$

Vektor satuan: $\hat{\mathbf{v}} = \dfrac{2\hat{\mathbf{x}} + \hat{\mathbf{y}} - 2\hat{\mathbf{z}}}{\sqrt{4+1+4}} = \dfrac{1}{3}(2\hat{\mathbf{x}} + \hat{\mathbf{y}} - 2\hat{\mathbf{z}}) \notag$

Turunan arah: $\nabla T\big|_P \cdot \hat{\mathbf{v}} = \dfrac{1}{3}[(-6)(2) + (7)(1) + (-1)(-2)] = \dfrac{1}{3}(-12+7+2) = -1 \notag$

Interpretasi: di titik $P$, fungsi $T$ menurun dengan laju 1 satuan per satuan panjang searah $\mathbf{v}$. Laju kenaikan maksimum adalah $|\nabla T|_P = \sqrt{36+49+1} = \sqrt{86} \approx 9.27$. $\notag$

Contoh 1.12 — Gradient Fungsi Radial

Hitung $\nabla r$ di mana $r = \sqrt{x^2+y^2+z^2}$. Ini adalah salah satu hasil paling penting dalam elektrodinamika.

$\dfrac{\partial r}{\partial x} = \dfrac{x}{\sqrt{x^2+y^2+z^2}} = \dfrac{x}{r}$, demikian pula untuk $y$ dan $z$. $\notag$

Sehingga: $\nabla r = \dfrac{x}{r}\,\hat{\mathbf{x}} + \dfrac{y}{r}\,\hat{\mathbf{y}} + \dfrac{z}{r}\,\hat{\mathbf{z}} = \dfrac{\mathbf{r}}{r} = \hat{\mathbf{r}} \notag$

Lebih umum, untuk $f(r)$ yang hanya bergantung pada $r$:

$$ \nabla f(r) = \frac{df}{dr}\,\nabla r = \frac{df}{dr}\,\hat{\mathbf{r}} = f'(r)\,\hat{\mathbf{r}} $$

Contoh: $\nabla\!\left(\dfrac{1}{r}\right) = -\dfrac{1}{r^2}\,\hat{\mathbf{r}} = -\dfrac{\mathbf{r}}{r^3}$ — ini muncul langsung dalam medan listrik muatan titik. $\notag$

Contoh 1.13 — Perbedaan Dua Titik

Dalam elektrodinamika, sering kali perlu menghitung $\nabla\!\left(\dfrac{1}{|\mathbf{r}-\mathbf{r}'|}\right)$ di mana $\mathbf{r}'$ adalah vektor posisi sumber (konstan terhadap turunan). Tunjukkan bahwa:

$$ \nabla\!\left(\frac{1}{|\mathbf{r}-\mathbf{r}'|}\right) = -\frac{\mathbf{r}-\mathbf{r}'}{|\mathbf{r}-\mathbf{r}'|^3} $$

Bukti. Misalkan $\boldsymbol{\mathscr{r}} = \mathbf{r}-\mathbf{r}' = (x-x')\hat{\mathbf{x}} + (y-y')\hat{\mathbf{y}} + (z-z')\hat{\mathbf{z}}$ dan $\mathscr{r} = |\boldsymbol{\mathscr{r}}|$. Karena $\nabla$ hanya beroperasi pada $\mathbf{r}$ (bukan $\mathbf{r}'$), kita peroleh: $\notag$

$\dfrac{\partial \mathscr{r}}{\partial x} = \dfrac{x-x'}{\mathscr{r}}$, sehingga $\nabla\mathscr{r} = \dfrac{\boldsymbol{\mathscr{r}}}{\mathscr{r}}$. Maka: $\notag$

$\nabla\!\left(\dfrac{1}{\mathscr{r}}\right) = -\dfrac{1}{\mathscr{r}^2}\nabla\mathscr{r} = -\dfrac{1}{\mathscr{r}^2}\cdot\dfrac{\boldsymbol{\mathscr{r}}}{\mathscr{r}} = -\dfrac{\boldsymbol{\mathscr{r}}}{\mathscr{r}^3} \;\blacksquare \notag$

Aturan Produk untuk Gradient

Karena $\nabla$ adalah operator diferensial, aturan produk mirip diferensiasi biasa berlaku, tetapi urutan penting:

$$ \nabla(fg) = f\,\nabla g + g\,\nabla f \qquad\qquad \nabla\!\left(\frac{f}{g}\right) = \frac{g\,\nabla f - f\,\nabla g}{g^2} $$

di mana $f$ dan $g$ adalah fungsi skalar. Perhatikan bahwa ini identik dengan aturan rantai biasa — $\nabla$ berperan seperti $d/dx$.

1.4 Divergensi dan Curl

Dua operasi turunan orde pertama yang paling penting dalam elektrodinamika diperoleh dengan menerapkan $\nabla$ pada vektor melalui produk titik dan produk silang.

Definisi 1.5 — Divergensi

Divergensi dari medan vektor $\mathbf{F}$ adalah skalar yang mengukur keluaran bersih (net outflow) per satuan volume di suatu titik:

$$ \nabla \cdot \mathbf{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z} $$

Secara intuitif, $\nabla \cdot \mathbf{F} > 0$ di titik sumber (source), $\nabla \cdot \mathbf{F} < 0$ di titik cekungan (sink), dan $\nabla \cdot \mathbf{F} = 0$ jika medan tidak memiliki sumber lokal (solenoidal).

Definisi 1.6 — Curl (Rotasi)

Curl dari medan vektor $\mathbf{F}$ adalah vektor yang mengukur kecenderungan rotasi medan di sekitar suatu titik:

$$ \nabla \times \mathbf{F} = \begin{vmatrix} \hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\[2pt] \dfrac{\partial}{\partial x} & \dfrac{\partial}{\partial y} & \dfrac{\partial}{\partial z} \\[6pt] F_x & F_y & F_z \end{vmatrix} $$

Secara eksplisit, komponen-komponen curl adalah:

$$ (\nabla \times \mathbf{F})_x = \frac{\partial F_z}{\partial y} - \frac{\partial F_y}{\partial z}, \quad (\nabla \times \mathbf{F})_y = \frac{\partial F_x}{\partial z} - \frac{\partial F_z}{\partial x}, \quad (\nabla \times \mathbf{F})_z = \frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y} \notag $$
Analogi Fisik

Bayangkan $\mathbf{F}$ sebagai medan kecepatan fluida. Divergensi mengukur apakah fluida "mengembang" atau "mengempis" secara lokal — seperti memasukkan pompa (sumber) atau drain (cekungan). Curl mengukur apakah fluida "berputar" di sekitar titik tersebut — seperti meletakkan bilah kipas kecil dan mengamati apakah ia berputar.

Contoh 1.14 — Divergensi dan Curl Medan Polinomial

Hitung divergensi dan curl dari $\mathbf{F} = x^2 y\,\hat{\mathbf{x}} + y^2 z\,\hat{\mathbf{y}} + z^2 x\,\hat{\mathbf{z}}$.

Divergensi: $\nabla\cdot\mathbf{F} = \dfrac{\partial}{\partial x}(x^2 y) + \dfrac{\partial}{\partial y}(y^2 z) + \dfrac{\partial}{\partial z}(z^2 x) = 2xy + 2yz + 2zx \notag$

Curl: $\nabla\times\mathbf{F} = \left(\dfrac{\partial(z^2 x)}{\partial y} - \dfrac{\partial(y^2 z)}{\partial z}\right)\hat{\mathbf{x}} - \left(\dfrac{\partial(z^2 x)}{\partial x} - \dfrac{\partial(x^2 y)}{\partial z}\right)\hat{\mathbf{y}} + \left(\dfrac{\partial(y^2 z)}{\partial x} - \dfrac{\partial(x^2 y)}{\partial y}\right)\hat{\mathbf{z}} \notag$

$= (0 - y^2)\hat{\mathbf{x}} - (z^2 - 0)\hat{\mathbf{y}} + (0 - x^2)\hat{\mathbf{z}} = -y^2\,\hat{\mathbf{x}} - z^2\,\hat{\mathbf{y}} - x^2\,\hat{\mathbf{z}} \notag$

Contoh 1.15 — Medan Radial $1/r^2$ (Hukum Coulomb)

Hitung divergensi dari $\mathbf{F} = \dfrac{\hat{\mathbf{r}}}{r^2} = \dfrac{\mathbf{r}}{r^3}$ untuk $r \neq 0$. Medan ini persis berbentuk medan listrik muatan titik.

Dalam koordinat Kartesius, $\mathbf{F} = \dfrac{x}{(x^2+y^2+z^2)^{3/2}}\hat{\mathbf{x}} + \dfrac{y}{(x^2+y^2+z^2)^{3/2}}\hat{\mathbf{y}} + \dfrac{z}{(x^2+y^2+z^2)^{3/2}}\hat{\mathbf{z}}$. $\notag$

Hitung $\dfrac{\partial F_x}{\partial x} = \dfrac{(x^2+y^2+z^2)^{3/2} - x\cdot\frac{3}{2}(x^2+y^2+z^2)^{1/2}\cdot 2x}{(x^2+y^2+z^2)^3} = \dfrac{r^3 - 3x^2 r}{r^6} = \dfrac{r^2 - 3x^2}{r^5} \notag$

Dengan cara yang sama: $\dfrac{\partial F_y}{\partial y} = \dfrac{r^2-3y^2}{r^5}$, $\quad\dfrac{\partial F_z}{\partial z} = \dfrac{r^2-3z^2}{r^5} \notag$

$\nabla\cdot\mathbf{F} = \dfrac{3r^2 - 3(x^2+y^2+z^2)}{r^5} = \dfrac{3r^2 - 3r^2}{r^5} = 0 \quad (r \neq 0) \notag$

Medan $\hat{\mathbf{r}}/r^2$ adalah solenoidal di mana-mana kecuali di asal. Di $r=0$, divergensi "meledak" — hasilnya melibatkan fungsi delta Dirac: $\nabla\cdot(\hat{\mathbf{r}}/r^2) = 4\pi\,\delta^3(\mathbf{r})$. Ini adalah inti dari hukum Gauss untuk muatan titik. $\notag$

Contoh 1.16 — Curl Medan Radial

Tunjukkan bahwa untuk setiap fungsi skalar $f(r)$, medan radial $\mathbf{F} = f(r)\,\hat{\mathbf{r}}$ memiliki curl nol.

Dalam koordinat bola (lihat Bagian 1.6), curl dari medan yang hanya memiliki komponen radial $F_r = f(r)$:

$(\nabla\times\mathbf{F})_\theta = \dfrac{1}{r\sin\theta}\dfrac{\partial F_r}{\partial \phi} - \dfrac{1}{r}\dfrac{\partial(r F_\phi)}{\partial r} = 0 - 0 = 0 \notag$

$(\nabla\times\mathbf{F})_\phi = \dfrac{1}{r}\left[\dfrac{\partial(r F_\theta)}{\partial r} - \dfrac{\partial F_r}{\partial \theta}\right] = 0 - \dfrac{f'(r)}{r}\dfrac{\partial r}{\partial \theta} = 0 \notag$

Karena $\partial r/\partial \theta = 0$, komponen ini juga nol. Jadi $\nabla\times(f(r)\,\hat{\mathbf{r}}) = \mathbf{0}$ di mana-mana. $\notag$

Konsekuensi fisik: Medan listrik statis yang hanya memiliki komponen radial (seperti medan muatan titik atau kulit bola) selalu konservatif — dapat ditulis sebagai gradient dari potensial skalar. $\notag$

Contoh 1.17 — Curl Medan Rotasi (Kawat Berarus)

Hitung curl dari $\mathbf{F} = \dfrac{-y\,\hat{\mathbf{x}} + x\,\hat{\mathbf{y}}}{x^2+y^2}$. Medan ini mewakili medan magnetik dari kawat tak hingga berarus sepanjang sumbu $z$.

$F_x = \dfrac{-y}{x^2+y^2}$, $\quad F_y = \dfrac{x}{x^2+y^2}$, $\quad F_z = 0 \notag$

$\dfrac{\partial F_z}{\partial y} - \dfrac{\partial F_y}{\partial z} = 0 - 0 = 0 \notag$

$\dfrac{\partial F_x}{\partial z} - \dfrac{\partial F_z}{\partial x} = 0 - 0 = 0 \notag$

$\dfrac{\partial F_y}{\partial x} - \dfrac{\partial F_x}{\partial y} = \dfrac{(x^2+y^2) - x\cdot 2x}{(x^2+y^2)^2} - \dfrac{-(x^2+y^2) + y\cdot 2y}{(x^2+y^2)^2} = \dfrac{y^2-x^2}{(x^2+y^2)^2} + \dfrac{x^2-y^2}{(x^2+y^2)^2} = 0 \quad (x^2+y^2 \neq 0) \notag$

Jadi $\nabla\times\mathbf{F} = \mathbf{0}$ di mana-mana kecuali di sumbu $z$ ($x=y=0$). Di sumbu, curl "meledak" dan menghasilkan fungsi delta: $\nabla\times\mathbf{F} = 2\pi\,\delta(x)\delta(y)\,\hat{\mathbf{z}}$. Ini adalah versi 2D dari hukum Ampère untuk kawat tak hingga. $\notag$

Identitas Orde Kedua

Dengan menerapkan $\nabla$ dua kali, diperoleh dua identitas fundamental yang merupakan konsekuensi langsung dari kesetaraan turunan parsial campuran ($\partial^2 f/\partial x\partial y = \partial^2 f/\partial y\partial x$):

$$ \nabla \cdot (\nabla \times \mathbf{F}) = 0 \qquad\text{(divergensi curl selalu nol)} $$
$$ \nabla \times (\nabla T) = \mathbf{0} \qquad\text{(curl gradient selalu nol)} $$

Identitas pertama menyatakan bahwa medan curl selalu solenoidal. Identitas kedua menyatakan bahwa gradient selalu irasional. Kedua identitas ini menjadi dasar pengenalan potensial listrik dan potensial vektor magnetik.

Contoh 1.18 — Verifikasi Identitas Orde Kedua

Verifikasi $\nabla\times(\nabla T) = \mathbf{0}$ untuk $T = x^2 y + y^2 z + z^2 x$.

Gradient: $\nabla T = (2xy+z^2)\hat{\mathbf{x}} + (x^2+2yz)\hat{\mathbf{y}} + (y^2+2zx)\hat{\mathbf{z}} \notag$

Komponen-$x$ dari curl: $\dfrac{\partial(y^2+2zx)}{\partial y} - \dfrac{\partial(x^2+2yz)}{\partial z} = 2y - 2y = 0 \notag$

Komponen-$y$ dari curl: $\dfrac{\partial(2xy+z^2)}{\partial z} - \dfrac{\partial(y^2+2zx)}{\partial x} = 2z - 2z = 0 \notag$

Komponen-$z$ dari curl: $\dfrac{\partial(x^2+2yz)}{\partial x} - \dfrac{\partial(2xy+z^2)}{\partial y} = 2x - 2x = 0 \notag$

Konfirmasi: $\nabla\times(\nabla T) = \mathbf{0} \;\checkmark \notag$

Operator Laplacian didefinisikan sebagai divergensi dari gradient:

$$ \nabla^2 T = \nabla \cdot (\nabla T) = \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} $$
Identitas Produk untuk Divergensi dan Curl

Aturan produk untuk $\nabla\cdot$ dan $\nabla\times$ memerlukan kehati-hatian karena $\nabla$ adalah operator:

$$ \nabla \cdot (f\,\mathbf{A}) = (\nabla f)\cdot\mathbf{A} + f\,(\nabla\cdot\mathbf{A}) $$
$$ \nabla \times (f\,\mathbf{A}) = (\nabla f)\times\mathbf{A} + f\,(\nabla\times\mathbf{A}) $$
$$ \nabla\cdot(\mathbf{A}\times\mathbf{B}) = \mathbf{B}\cdot(\nabla\times\mathbf{A}) - \mathbf{A}\cdot(\nabla\times\mathbf{B}) $$
$$ \nabla\times(\mathbf{A}\times\mathbf{B}) = (\mathbf{B}\cdot\nabla)\mathbf{A} - (\mathbf{A}\cdot\nabla)\mathbf{B} + \mathbf{A}(\nabla\cdot\mathbf{B}) - \mathbf{B}(\nabla\cdot\mathbf{A}) $$

Identitas terakhir sering muncul dalam derivasi hukum Ampère–Maxwell dan konsep tekanan Maxwell.

Contoh 1.19 — Aturan Produk Divergensi

Verifikasi $\nabla\cdot(r^n\,\hat{\mathbf{r}}) = (n+2)\,r^{n-1}$ menggunakan aturan produk.

Misalkan $f = r^n$ dan $\mathbf{A} = \hat{\mathbf{r}}$. Maka: $\notag$

$\nabla\cdot(r^n\,\hat{\mathbf{r}}) = (\nabla r^n)\cdot\hat{\mathbf{r}} + r^n(\nabla\cdot\hat{\mathbf{r}}) \notag$

Karena $\nabla r^n = nr^{n-1}\nabla r = nr^{n-1}\hat{\mathbf{r}}$, komponen pertama: $(nr^{n-1}\hat{\mathbf{r}})\cdot\hat{\mathbf{r}} = nr^{n-1} \notag$

Dalam koordinat bola, $\nabla\cdot\hat{\mathbf{r}} = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}(r^2\cdot 1) = \dfrac{2}{r}$. Jadi komponen kedua: $r^n \cdot \dfrac{2}{r} = 2r^{n-1} \notag$

Total: $\nabla\cdot(r^n\,\hat{\mathbf{r}}) = nr^{n-1} + 2r^{n-1} = (n+2)r^{n-1} \;\checkmark \notag$

Kasus khusus: Untuk $n=-2$: $\nabla\cdot(\hat{\mathbf{r}}/r^2) = 0$ untuk $r\neq 0$, konsisten dengan Contoh 1.15. $\notag$

Contoh 1.20 — Laplacian $1/r$

Hitung $\nabla^2(1/r)$ untuk $r\neq 0$.

Dalam koordinat bola, untuk fungsi yang hanya bergantung pada $r$: $\notag$

$\nabla^2\!\left(\dfrac{1}{r}\right) = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}\!\left(r^2\dfrac{\partial}{\partial r}\dfrac{1}{r}\right) = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}\!\left(r^2\cdot\dfrac{-1}{r^2}\right) = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}(-1) = 0 \quad (r\neq 0) \notag$

Namun, di $r=0$, Laplacian ini "meledak". Hasil lengkapnya: $\notag$

$$ \nabla^2\!\left(\frac{1}{r}\right) = -4\pi\,\delta^3(\mathbf{r}) $$

Ini adalah identitas Laplace paling penting dalam elektrodinamika — ia menyatakan bahwa $1/r$ adalah fungsi Green dari operator Laplacian di tiga dimensi, dan merupakan dasar solusi umum persamaan Poisson. $\notag$

1.5 Teorema Integral: Gauss, Stokes, dan Green

Tiga teorema integral fundamental menghubungkan operasi diferensial lokal (divergensi, curl) dengan integral global atas permukaan dan kurva tertutup.

1.5.1 Teorema Divergensi (Gauss)

Teorema 1.2 — Teorema Divergensi Gauss

Integral fluks vektor $\mathbf{F}$ menembus permukaan tertutup $\mathcal{S}$ sama dengan integral volume divergensi $\mathbf{F}$ di dalam volume $\mathcal{V}$ yang dibatasi $\mathcal{S}$:

$$ \oint_{\mathcal{S}} \mathbf{F} \cdot d\mathbf{a} = \int_{\mathcal{V}} (\nabla \cdot \mathbf{F})\, d\tau $$

di mana $d\mathbf{a} = \hat{\mathbf{n}}\,da$ dengan $\hat{\mathbf{n}}$ adalah vektor normal satuan mengarah ke luar permukaan tertutup.

V S Fluks keluar melalui permukaan tertutup S
Gambar 1.5 — Teorema Gauss: fluks total $\mathbf{F}$ yang keluar melalui permukaan tertutup $\mathcal{S}$ sama dengan integral volume $\nabla\cdot\mathbf{F}$ di dalamnya.
Contoh 1.21 — Verifikasi Teorema Gauss (Kubus)

Verifikasi teorema divergensi untuk $\mathbf{F} = x^2\hat{\mathbf{x}} + y^2\hat{\mathbf{y}} + z^2\hat{\mathbf{z}}$ pada kubus $0 \leq x,y,z \leq 1$.

Ruang kanan (integral volume): $\nabla\cdot\mathbf{F} = 2x + 2y + 2z \notag$

$\int_0^1\!\!\int_0^1\!\!\int_0^1 (2x+2y+2z)\,dx\,dy\,dz = \int_0^1\!\!\int_0^1 [x^2+2xy+2xz]_0^1\,dy\,dz = \int_0^1\!\!\int_0^1 (1+2y+2z)\,dy\,dz \notag$

$= \int_0^1 [y+y^2+2yz]_0^1\,dz = \int_0^1 (2+2z)\,dz = [2z+z^2]_0^1 = 3 \notag$

Ruang kiri (integral permukaan) — oleh simetri, cukup hitung satu sisi dan kalikan: $\notag$

Sisi $x=1$: $\int_0^1\!\!\int_0^1 (1)\,dy\,dz = 1$. Sisi $x=0$: $\int_0^1\!\!\int_0^1 (0)\,dy\,dz = 0$. $\notag$

Dengan cara sama untuk sisi-sisi $y$ dan $z$, total fluks $= (1+0) + (1+0) + (1+0) = 3 \;\checkmark \notag$

Contoh 1.22 — Fluks Melalui Bola (Hukum Gauss untuk Muatan Titik)

Hitung fluks $\mathbf{F} = \hat{\mathbf{r}}/r^2$ melalui permukaan bola $r = R$.

Pada permukaan bola, $d\mathbf{a} = \hat{\mathbf{r}}\,da$ dan $r = R$, sehingga: $\notag$

$\oint_{r=R} \mathbf{F}\cdot d\mathbf{a} = \oint_{r=R} \frac{\hat{\mathbf{r}}}{R^2}\cdot\hat{\mathbf{r}}\,da = \frac{1}{R^2}\oint da = \frac{1}{R^2}\cdot 4\pi R^2 = 4\pi \notag$

Penting: Hasilnya $4\pi$ independen dari $R$! Ini bukan kebetulan — ia secara langsung mengarah ke hukum Gauss $\oint \mathbf{E}\cdot d\mathbf{a} = Q_{\text{enc}}/\varepsilon_0$. Dengan teorema divergensi, karena $\nabla\cdot(\hat{\mathbf{r}}/r^2) = 0$ untuk $r\neq 0$, seluruh kontribusi fluks berasal dari singularitas di $r=0$, yaitu $\nabla\cdot(\hat{\mathbf{r}}/r^2) = 4\pi\delta^3(\mathbf{r})$. $\notag$

Contoh 1.23 — Volume Melalui Teorema Divergensi

Gunakan teorema divergensi untuk menghitung volume bola radius $R$ dengan memilih $\mathbf{F} = \mathbf{r}/3$.

$\nabla\cdot\mathbf{F} = \nabla\cdot\dfrac{\mathbf{r}}{3} = \dfrac{1}{3}\left(\dfrac{\partial x}{\partial x} + \dfrac{\partial y}{\partial y} + \dfrac{\partial z}{\partial z}\right) = \dfrac{3}{3} = 1 \notag$

Oleh teorema Gauss: $V = \int_{\mathcal{V}} 1\,d\tau = \int_{\mathcal{V}} (\nabla\cdot\mathbf{F})\,d\tau = \oint_{\mathcal{S}} \mathbf{F}\cdot d\mathbf{a} = \oint \frac{\mathbf{r}}{3}\cdot\hat{\mathbf{r}}\,da \notag$

Pada permukaan bola, $\mathbf{r} = R\hat{\mathbf{r}}$, sehingga $\mathbf{F}\cdot d\mathbf{a} = \dfrac{R}{3}\,da$. Maka: $\notag$

$V = \dfrac{R}{3}\oint da = \dfrac{R}{3}\cdot 4\pi R^2 = \dfrac{4\pi R^3}{3} \;\checkmark \notag$

Trik ini — memilih $\mathbf{F}$ sehingga $\nabla\cdot\mathbf{F}=1$ — dapat digunakan untuk menghitung volume region apa pun melalui integral permukaan. $\notag$

1.5.2 Teorema Stokes

Teorema 1.3 — Teorema Stokes

Integral garis $\mathbf{F}$ sepanjang kurva tertutup $\mathcal{C}$ sama dengan integral curl $\mathbf{F}$ melalui setiap permukaan terbuka $\mathcal{S}$ yang dibatasi oleh $\mathcal{C}$:

$$ \oint_{\mathcal{C}} \mathbf{F} \cdot d\mathbf{l} = \int_{\mathcal{S}} (\nabla \times \mathbf{F}) \cdot d\mathbf{a} $$

Arah integrasi pada $\mathcal{C}$ dan arah normal $\hat{\mathbf{n}}$ pada $\mathcal{S}$ memenuhi aturan tangan kanan: jika jari-jari tangan kanan mengikuti arah integrasi pada $\mathcal{C}$, ibu jari menunjuk arah $\hat{\mathbf{n}}$.

S C dl Arah C dan n̂ memenuhi aturan tangan kanan
Gambar 1.6 — Teorema Stokes: integral garis sepanjang batas $\mathcal{C}$ sama dengan integral curl melalui permukaan terbuka $\mathcal{S}$ yang dibatasi $\mathcal{C}$. Perhatikan bahwa $\mathcal{S}$ dapat berbentuk apa pun selama batasnya $\mathcal{C}$.
Contoh 1.24 — Verifikasi Teorema Stokes (Lingkaran)

Verifikasi teorema Stokes untuk $\mathbf{F} = -y\hat{\mathbf{x}} + x\hat{\mathbf{y}} + z\hat{\mathbf{z}}$ pada lingkaran $x^2+y^2=1$, $z=0$.

Ruang kiri (integral garis): Parametrisasi $\mathbf{r}(\phi) = \cos\phi\,\hat{\mathbf{x}} + \sin\phi\,\hat{\mathbf{y}}$, $d\mathbf{l} = (-\sin\phi\,\hat{\mathbf{x}}+\cos\phi\,\hat{\mathbf{y}})\,d\phi$. $\notag$

$\mathbf{F}\cdot d\mathbf{l} = (-\sin\phi)(-\sin\phi) + (\cos\phi)(\cos\phi) = \sin^2\phi + \cos^2\phi = 1 \notag$

$\oint_{\mathcal{C}} \mathbf{F}\cdot d\mathbf{l} = \int_0^{2\pi} 1\,d\phi = 2\pi \notag$

Ruang kanan (integral permukaan): $\nabla\times\mathbf{F} = \begin{vmatrix}\hat{\mathbf{x}} & \hat{\mathbf{y}} & \hat{\mathbf{z}} \\ \partial_x & \partial_y & \partial_z \\ -y & x & z\end{vmatrix} = 2\hat{\mathbf{z}} \notag$

Pada permukaan $z=0$ di dalam lingkaran: $d\mathbf{a} = \hat{\mathbf{z}}\,da$, sehingga $(\nabla\times\mathbf{F})\cdot d\mathbf{a} = 2\,da$. $\notag$

$\int_{\mathcal{S}} 2\,da = 2\cdot\pi(1)^2 = 2\pi \;\checkmark \notag$

Contoh 1.25 — Stokes dengan Permukaan Berbeda

Untuk $\mathbf{F}$ dan kurva $\mathcal{C}$ yang sama dengan Contoh 1.24, hitung integral permukaan menggunakan setengah bola $x^2+y^2+z^2=1$, $z\geq 0$ sebagai pengganti cakram datar.

Curl tetap $\nabla\times\mathbf{F} = 2\hat{\mathbf{z}}$. Pada permukaan bola, $d\mathbf{a} = \hat{\mathbf{r}}\,R^2\sin\theta\,d\theta\,d\phi$ dengan $R=1$. $\notag$

$(\nabla\times\mathbf{F})\cdot d\mathbf{a} = 2\hat{\mathbf{z}}\cdot\hat{\mathbf{r}}\sin\theta\,d\theta\,d\phi = 2\cos\theta\,\sin\theta\,d\theta\,d\phi \notag$

karena $\hat{\mathbf{z}}\cdot\hat{\mathbf{r}} = \cos\theta$. Maka: $\notag$

$\int_0^{2\pi}\!\!\int_0^{\pi/2} 2\cos\theta\sin\theta\,d\theta\,d\phi = 2\pi \int_0^{\pi/2}\sin(2\theta)\,d\theta = 2\pi\left[-\frac{1}{2}\cos(2\theta)\right]_0^{\pi/2} = 2\pi(1) = 2\pi \;\checkmark \notag$

Hasilnya identik dengan cakram datar — ini menunjukkan keindahan teorema Stokes: kita bebas memilih permukaan mana pun yang dibatasi oleh $\mathcal{C}$. $\notag$

Contoh 1.26 — Gaya EM pada Muatan Bergerak (Stokes dalam Fisika)

Dalam elektrodinamika, teorema Stokes digunakan untuk menunjukkan bahwa GGL induksi $\mathcal{E} = \oint \mathbf{E}\cdot d\mathbf{l} = -\dfrac{d\Phi_B}{dt}$ setara dengan $\nabla\times\mathbf{E} = -\dfrac{\partial\mathbf{B}}{\partial t}$. Sebagai ilustrasi, jika $\mathbf{B} = B_0 t\,\hat{\mathbf{z}}$ dengan $B_0$ konstan, hitung GGL pada lingkaran radius $a$ di bidang $xy$.

$\mathcal{E} = \int_{\mathcal{S}}(\nabla\times\mathbf{E})\cdot d\mathbf{a} = -\int_{\mathcal{S}}\frac{\partial\mathbf{B}}{\partial t}\cdot d\mathbf{a} = -\int_{\mathcal{S}} B_0\,\hat{\mathbf{z}}\cdot\hat{\mathbf{z}}\,da = -B_0\cdot\pi a^2 \notag$

Jadi $\mathcal{E} = -\pi a^2 B_0$, negatif menandakan arah GGL melawan perubahan fluks (hukum Lenz). $\notag$

1.5.3 Teorema Green (Bidang)

Kasus Khusus 2D

Teorema Green adalah bentuk khusus teorema Stokes di bidang $xy$. Untuk medan $\mathbf{F} = F_x\hat{\mathbf{x}} + F_y\hat{\mathbf{y}}$ dan daerah $\mathcal{R}$ dengan batas $\mathcal{C}$:

$$ \oint_{\mathcal{C}} (F_x\,dx + F_y\,dy) = \int_{\mathcal{R}} \left(\frac{\partial F_y}{\partial x} - \frac{\partial F_x}{\partial y}\right) dx\,dy $$

Perhatikan bahwa integran di ruas kanan adalah tepat komponen-$z$ dari $\nabla\times\mathbf{F}$, konsisten dengan interpretasi teorema Green sebagai kasus 2D dari Stokes.

Contoh 1.27 — Teorema Green untuk Luas

Gunakan teorema Green untuk menghitung luas elips $\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1$.

Pilih $F_x = -y/2$, $F_y = x/2$. Maka $\dfrac{\partial F_y}{\partial x} - \dfrac{\partial F_x}{\partial y} = \dfrac{1}{2} + \dfrac{1}{2} = 1$. $\notag$

Luas $= \displaystyle\int_{\mathcal{R}} 1\,dx\,dy = \oint_{\mathcal{C}} \left(-\frac{y}{2}\,dx + \frac{x}{2}\,dy\right) \notag$

Parametrisasi elips: $x = a\cos t$, $y = b\sin t$, $0\leq t\leq 2\pi$: $\notag$

$A = \int_0^{2\pi}\left[-\frac{b\sin t}{2}(-a\sin t) + \frac{a\cos t}{2}(b\cos t)\right]dt = \frac{ab}{2}\int_0^{2\pi}(\sin^2 t + \cos^2 t)\,dt = \frac{ab}{2}\cdot 2\pi = \pi ab \;\checkmark \notag$

1.5.4 Identitas Green

Identitas Green Pertama dan Kedua

Dengan menerapkan teorema divergensi pada $\nabla\cdot(f\,\nabla g)$, diperoleh:

Identitas Green pertama:

$$ \int_{\mathcal{V}} \left[(\nabla f)\cdot(\nabla g) + f\,\nabla^2 g\right]d\tau = \oint_{\mathcal{S}} f\,\frac{\partial g}{\partial n}\,da $$

Identitas Green kedua: (diperoleh dengan menukar $f$ dan $g$ lalu mengurangkan)

$$ \int_{\mathcal{V}} \left(f\,\nabla^2 g - g\,\nabla^2 f\right)d\tau = \oint_{\mathcal{S}} \left(f\,\frac{\partial g}{\partial n} - g\,\frac{\partial f}{\partial n}\right)da $$

di mana $\partial g/\partial n = \nabla g\cdot\hat{\mathbf{n}}$ adalah turunan normal ke luar. Identitas Green kedua adalah alat utama untuk membuktikan keunikan solusi persamaan Laplace dan Poisson — yaitu, membuktikan bahwa persamaan $\nabla^2 V = -\rho/\varepsilon_0$ dengan syarat batas yang diberikan memiliki tepat satu solusi.

1.6 Sistem Koordinat Kurvilinier

Dalam elektrodinamika, banyak masalah memiliki simetri yang lebih alami diekspresikan dalam koordinat silinder atau bola daripada koordinat Kartesius. Sistem koordinat kurvilinier umum $(u_1, u_2, u_3)$ didefinisikan oleh fungsi transformasi $\mathbf{r}(u_1, u_2, u_3)$.

Definisi 1.7 — Faktor Skala

Dalam koordinat kurvilinier $(u_1, u_2, u_3)$, vektor basis lokal didefinisikan sebagai $\mathbf{e}_i = \partial\mathbf{r}/\partial u_i$. Faktor skala adalah besar vektor basis: $h_i = |\mathbf{e}_i|$. Vektor satuan lokal: $\hat{\mathbf{u}}_i = \mathbf{e}_i/h_i$. Elemen panjang busur:

$$ dl^2 = h_1^2\,du_1^2 + h_2^2\,du_2^2 + h_3^2\,du_3^2 $$

1.6.1 Koordinat Silinder

z s P s z φ ŝ φ̂ O Koordinat silinder (s, φ, z)
Gambar 1.7 — Koordinat silinder $(s, \phi, z)$: $s$ adalah jarak dari sumbu $z$, $\phi$ adalah sudut azimuth dari sumbu $+x$, dan $z$ adalah tinggi. Vektor satuan $\hat{\mathbf{s}}$, $\hat{\boldsymbol{\phi}}$, $\hat{\mathbf{z}}$ membentuk basis ortonormal kanan di setiap titik.
Koordinat Transformasi Faktor Skala
$s$ (sumbu) $x = s\cos\phi,\; y = s\sin\phi$ $h_s = 1$
$\phi$ (azimuth) $\phi = \arctan(y/x)$ $h_\phi = s$
$z$ (vertikal) $z = z$ $h_z = 1$
Konvensi Notasi

Dalam buku ini, koordinat silinder ditulis $(s, \phi, z)$ — bukan $(r, \phi, z)$ — untuk menghindari kebingungan dengan koordinat radial bola $r$. Simbol $r$ selalu dicadangkan untuk jarak dari titik asal dalam koordinat bola.

Operator Diferensial dalam Koordinat Silinder

Dengan mensubstitusikan $h_s=1$, $h_\phi=s$, $h_z=1$ ke rumus umum kurvilinier:

Gradient:

$$ \nabla T = \frac{\partial T}{\partial s}\,\hat{\mathbf{s}} + \frac{1}{s}\frac{\partial T}{\partial \phi}\,\hat{\boldsymbol{\phi}} + \frac{\partial T}{\partial z}\,\hat{\mathbf{z}} $$

Divergensi:

$$ \nabla \cdot \mathbf{F} = \frac{1}{s}\frac{\partial}{\partial s}(s F_s) + \frac{1}{s}\frac{\partial F_\phi}{\partial \phi} + \frac{\partial F_z}{\partial z} $$

Curl:

$$ \nabla \times \mathbf{F} = \left(\frac{1}{s}\frac{\partial F_z}{\partial \phi} - \frac{\partial F_\phi}{\partial z}\right)\hat{\mathbf{s}} + \left(\frac{\partial F_s}{\partial z} - \frac{\partial F_z}{\partial s}\right)\hat{\boldsymbol{\phi}} + \frac{1}{s}\left[\frac{\partial}{\partial s}(s F_\phi) - \frac{\partial F_s}{\partial \phi}\right]\hat{\mathbf{z}} $$

Laplacian:

$$ \nabla^2 T = \frac{1}{s}\frac{\partial}{\partial s}\!\left(s\frac{\partial T}{\partial s}\right) + \frac{1}{s^2}\frac{\partial^2 T}{\partial \phi^2} + \frac{\partial^2 T}{\partial z^2} $$
Contoh 1.28 — Divergensi dalam Koordinat Silinder

Hitung $\nabla\cdot\mathbf{F}$ untuk $\mathbf{F} = s\,\hat{\mathbf{s}}$ (medan radial silinder) dalam koordinat silinder, lalu verifikasi dalam Kartesius.

Silinder: $\nabla\cdot(s\,\hat{\mathbf{s}}) = \dfrac{1}{s}\dfrac{\partial}{\partial s}(s\cdot s) + 0 + 0 = \dfrac{1}{s}\cdot 2s = 2 \notag$

Verifikasi Kartesius: $\mathbf{F} = s\,\hat{\mathbf{s}} = s\cos\phi\,\hat{\mathbf{x}} + s\sin\phi\,\hat{\mathbf{y}} = x\,\hat{\mathbf{x}} + y\,\hat{\mathbf{y}} \notag$

$\nabla\cdot\mathbf{F} = \dfrac{\partial x}{\partial x} + \dfrac{\partial y}{\partial y} + 0 = 1 + 1 = 2 \;\checkmark \notag$

Perhatikan bahwa faktor $1/s$ dalam rumus silinder tepat mengkompensasi fakta bahwa area cincin bertambah dengan $s$. $\notag$

Contoh 1.29 — Curl dalam Koordinat Silinder (Solenoid)

Medan magnetik di dalam solenoid ideal tak hingga: $\mathbf{B} = B_0\,\hat{\mathbf{z}}$ (konstan). Hitung $\nabla\times\mathbf{B}$.

$\nabla\times\mathbf{B} = \left(\dfrac{1}{s}\dfrac{\partial B_0}{\partial \phi} - 0\right)\hat{\mathbf{s}} + \left(0 - \dfrac{\partial B_0}{\partial s}\right)\hat{\boldsymbol{\phi}} + \dfrac{1}{s}\left[\dfrac{\partial}{\partial s}(s\cdot 0) - 0\right]\hat{\mathbf{z}} = \mathbf{0} \notag$

Konsisten dengan hukum Ampère: di dalam solenoid (jauh dari ujung), tidak ada arus yang mengalir melalui permukaan yang dibungkus, sehingga $\nabla\times\mathbf{B} = \mu_0\mathbf{J} = \mathbf{0}$. $\notag$

1.6.2 Koordinat Bola

z x P r s = r sin θ r cos θ θ φ O
Gambar 1.8 — Koordinat bola $(r, \theta, \phi)$: $r$ adalah jarak dari asal, $\theta$ adalah sudut polar dari sumbu $+z$ ($0 \leq \theta \leq \pi$), dan $\phi$ adalah sudut azimuth dari sumbu $+x$ ($0 \leq \phi < 2\pi$).
Koordinat Transformasi Faktor Skala
$r$ (radial) $x = r\sin\theta\cos\phi,\; y = r\sin\theta\sin\phi,\; z = r\cos\theta$ $h_r = 1$
$\theta$ (polar) $\theta = \arccos(z/r)$ $h_\theta = r$
$\phi$ (azimuth) $\phi = \arctan(y/x)$ $h_\phi = r\sin\theta$
Operator Diferensial dalam Koordinat Bola

Dengan mensubstitusikan $h_r=1$, $h_\theta=r$, $h_\phi=r\sin\theta$:

Gradient:

$$ \nabla T = \frac{\partial T}{\partial r}\,\hat{\mathbf{r}} + \frac{1}{r}\frac{\partial T}{\partial \theta}\,\hat{\boldsymbol{\theta}} + \frac{1}{r\sin\theta}\frac{\partial T}{\partial \phi}\,\hat{\boldsymbol{\phi}} $$

Divergensi:

$$ \nabla \cdot \mathbf{F} = \frac{1}{r^2}\frac{\partial}{\partial r}(r^2 F_r) + \frac{1}{r\sin\theta}\frac{\partial}{\partial \theta}(\sin\theta\, F_\theta) + \frac{1}{r\sin\theta}\frac{\partial F_\phi}{\partial \phi} $$

Curl:

$$ \nabla \times \mathbf{F} = \frac{1}{r\sin\theta}\left[\frac{\partial}{\partial \theta}(\sin\theta\, F_\phi) - \frac{\partial F_\theta}{\partial \phi}\right]\hat{\mathbf{r}} + \frac{1}{r}\left[\frac{1}{\sin\theta}\frac{\partial F_r}{\partial \phi} - \frac{\partial}{\partial r}(r F_\phi)\right]\hat{\boldsymbol{\theta}} + \frac{1}{r}\left[\frac{\partial}{\partial r}(r F_\theta) - \frac{\partial F_r}{\partial \theta}\right]\hat{\boldsymbol{\phi}} $$

Laplacian:

$$ \nabla^2 T = \frac{1}{r^2}\frac{\partial}{\partial r}\!\left(r^2\frac{\partial T}{\partial r}\right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial \theta}\!\left(\sin\theta\,\frac{\partial T}{\partial \theta}\right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2 T}{\partial \phi^2} $$
Contoh 1.30 — Konversi Koordinat Kartesius ke Bola

Tentukan koordinat bola $(r, \theta, \phi)$ dari titik $P$ dengan koordinat Kartesius $(1, 1, \sqrt{2})$.

$r = \sqrt{1^2+1^2+(\sqrt{2})^2} = \sqrt{1+1+2} = \sqrt{4} = 2 \notag$

$\theta = \arccos\!\left(\dfrac{z}{r}\right) = \arccos\!\left(\dfrac{\sqrt{2}}{2}\right) = \dfrac{\pi}{4} = 45° \notag$

$\phi = \arctan\!\left(\dfrac{y}{x}\right) = \arctan(1) = \dfrac{\pi}{4} = 45° \notag$

Jadi $P$ memiliki koordinat bola $(r, \theta, \phi) = (2, \pi/4, \pi/4)$. $\notag$

Contoh 1.31 — Konversi Koordinat Kartesius ke Silinder

Tentukan koordinat silinder $(s, \phi, z)$ dari titik yang sama $P = (1, 1, \sqrt{2})$.

$s = \sqrt{x^2+y^2} = \sqrt{1+1} = \sqrt{2} \notag$

$\phi = \arctan\!\left(\dfrac{y}{x}\right) = \arctan(1) = \dfrac{\pi}{4} \notag$

$z = \sqrt{2} \notag$

Jadi $P$ memiliki koordinat silinder $(s, \phi, z) = (\sqrt{2}, \pi/4, \sqrt{2})$. $\notag$

Contoh 1.32 — Divergensi dalam Koordinat Bola

Hitung $\nabla\cdot\mathbf{F}$ untuk $\mathbf{F} = r^2\cos\theta\,\hat{\mathbf{r}}$ dalam koordinat bola.

$\nabla\cdot\mathbf{F} = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}(r^2 \cdot r^2\cos\theta) + 0 + 0 = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}(r^4\cos\theta) = \dfrac{1}{r^2}\cdot 4r^3\cos\theta = 4r\cos\theta \notag$

Contoh 1.33 — Curl Medan Dipol Magnetik

Medan vektor potensial dari dipol magnetik $\mathbf{m} = m\hat{\mathbf{z}}$ adalah $\mathbf{A} = \dfrac{\mu_0}{4\pi}\dfrac{m\sin\theta}{r^2}\,\hat{\boldsymbol{\phi}}$. Hitung $\mathbf{B} = \nabla\times\mathbf{A}$.

Dengan $A_\phi = \dfrac{\mu_0 m\sin\theta}{4\pi r^2}$ dan $A_r = A_\theta = 0$: $\notag$

Komponen-$r$: $\dfrac{1}{r\sin\theta}\left[\dfrac{\partial}{\partial\theta}\!\left(\sin\theta\cdot\dfrac{\mu_0 m\sin\theta}{4\pi r^2}\right) - 0\right] = \dfrac{\mu_0 m}{4\pi r^3\sin\theta}\dfrac{\partial}{\partial\theta}(\sin^2\theta) = \dfrac{\mu_0 m}{4\pi r^3\sin\theta}\cdot 2\sin\theta\cos\theta = \dfrac{\mu_0 m}{4\pi}\dfrac{2\cos\theta}{r^3} \notag$

Komponen-$\theta$: $\dfrac{1}{r}\left[0 - \dfrac{\partial}{\partial r}\!\left(r\cdot\dfrac{\mu_0 m\sin\theta}{4\pi r^2}\right)\right] = \dfrac{1}{r}\left[-\dfrac{\mu_0 m\sin\theta}{4\pi}\dfrac{\partial}{\partial r}\!\left(\dfrac{1}{r}\right)\right] = \dfrac{1}{r}\left[-\dfrac{\mu_0 m\sin\theta}{4\pi}\!\left(\dfrac{-1}{r^2}\right)\right] = \dfrac{\mu_0 m}{4\pi}\dfrac{\sin\theta}{r^3} \notag$

Komponen-$\phi$: $0$ (karena $A_r$ dan $A_\theta$ nol, dan $A_\phi$ tidak bergantung pada $\phi$). $\notag$

Sehingga: $\notag$

$$ \mathbf{B} = \frac{\mu_0 m}{4\pi r^3}\left(2\cos\theta\,\hat{\mathbf{r}} + \sin\theta\,\hat{\boldsymbol{\theta}}\right) $$

Ini adalah ekspresi standar medan dipol magnetik yang akan kita temukan kembali di Bab 5 melalui multipol ekspansi. $\notag$

Contoh 1.34 — Laplacian dalam Koordinat Bola (Harmonik Sferis)

Hitung $\nabla^2(r\cos\theta)$ dan $\nabla^2(r^2\cos^2\theta)$. Fungsi-fungsi ini berkaitan dengan harmonik sferis $Y_1^0$ dan $Y_2^0$.

Pertama: $T = r\cos\theta = z$. Karena $z$ adalah koordinat Kartesius, $\nabla^2 z = 0$. Verifikasi dalam koordinat bola: $\notag$

$\nabla^2(r\cos\theta) = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}\!\left(r^2\cos\theta\right) + \dfrac{1}{r^2\sin\theta}\dfrac{\partial}{\partial\theta}\!\left(\sin\theta\cdot(-r\sin\theta)\right) + 0 \notag$

$= \dfrac{2\cos\theta}{r} + \dfrac{1}{r\sin\theta}\dfrac{\partial}{\partial\theta}(-r\sin^2\theta) = \dfrac{2\cos\theta}{r} + \dfrac{-r}{r\sin\theta}\cdot 2\sin\theta\cos\theta = \dfrac{2\cos\theta}{r} - \dfrac{2\cos\theta}{r} = 0 \;\checkmark \notag$

Kedua: $T = r^2\cos^2\theta$. Langsung dalam koordinat bola: $\notag$

$\dfrac{1}{r^2}\dfrac{\partial}{\partial r}\!\left(r^2\cdot 2r\cos^2\theta\right) = \dfrac{1}{r^2}\dfrac{\partial}{\partial r}(2r^3\cos^2\theta) = 6\cos^2\theta \notag$

$\dfrac{1}{r^2\sin\theta}\dfrac{\partial}{\partial\theta}\!\left(\sin\theta\cdot r^2\cdot 2\cos\theta\cdot(-\sin\theta)\right) = \dfrac{1}{\sin\theta}\dfrac{\partial}{\partial\theta}(-2\cos\theta\sin^2\theta) \notag$

$= \dfrac{1}{\sin\theta}(-2)[-\sin\theta\sin^2\theta + \cos\theta\cdot 2\sin\theta\cos\theta] = \dfrac{1}{\sin\theta}(-2)\sin\theta(-\sin^2\theta+2\cos^2\theta) \notag$

$= 2(3\cos^2\theta - 1) \notag$

$\nabla^2(r^2\cos^2\theta) = 6\cos^2\theta + 2(3\cos^2\theta-1) = 6\cos^2\theta + 6\cos^2\theta - 2 = 12\cos^2\theta - 2 = 2(6\cos^2\theta-1) \notag$

Hasil ini bukan nol — artinya $r^2\cos^2\theta$ bukan solusi persamaan Laplace (bukan harmonik sferis murni). Namun $r^2(3\cos^2\theta-1)$ adalah harmonik sferis $Y_2^0$ yang memenuhi $\nabla^2 Y = 0$. $\notag$

1.6.3 Grad, Div, dan Curl dalam Koordinat Kurvilinier Umum

Dalam koordinat kurvilinier umum dengan faktor skala $h_1, h_2, h_3$, operator diferensial mengambil bentuk:

Operator Diferensial dalam Koordinat Kurvilinier

Gradient:

$$ \nabla T = \hat{\mathbf{u}}_1\,\frac{1}{h_1}\frac{\partial T}{\partial u_1} + \hat{\mathbf{u}}_2\,\frac{1}{h_2}\frac{\partial T}{\partial u_2} + \hat{\mathbf{u}}_3\,\frac{1}{h_3}\frac{\partial T}{\partial u_3} $$

Divergensi:

$$ \nabla \cdot \mathbf{F} = \frac{1}{h_1 h_2 h_3}\left[ \frac{\partial}{\partial u_1}(h_2 h_3 F_1) + \frac{\partial}{\partial u_2}(h_1 h_3 F_2) + \frac{\partial}{\partial u_3}(h_1 h_2 F_3) \right] $$

Curl:

$$ \nabla \times \mathbf{F} = \frac{1}{h_1 h_2 h_3}\begin{vmatrix} h_1\hat{\mathbf{u}}_1 & h_2\hat{\mathbf{u}}_2 & h_3\hat{\mathbf{u}}_3 \\[4pt] \dfrac{\partial}{\partial u_1} & \dfrac{\partial}{\partial u_2} & \dfrac{\partial}{\partial u_3} \\[6pt] h_1 F_1 & h_2 F_2 & h_3 F_3 \end{vmatrix} $$

Laplacian:

$$ \nabla^2 T = \frac{1}{h_1 h_2 h_3}\left[ \frac{\partial}{\partial u_1}\!\left(\frac{h_2 h_3}{h_1}\frac{\partial T}{\partial u_1}\right) + \frac{\partial}{\partial u_2}\!\left(\frac{h_1 h_3}{h_2}\frac{\partial T}{\partial u_2}\right) + \frac{\partial}{\partial u_3}\!\left(\frac{h_1 h_2}{h_3}\frac{\partial T}{\partial u_3}\right) \right] $$
Kapan Menggunakan Koordinat Apa?

Kartesius $(x,y,z)$: masalah tanpa simetri tertentu, atau bidang-bidang datar. Contoh: kapasitor pelat sejajar, gelombang plan.
Silinder $(s,\phi,z)$: masalah dengan simetri aksial — kawat tak hingga, silinder, solenoid, kabel koaksial.
Bola $(r,\theta,\phi)$: masalah dengan simetri bola — muatan titik, kulit bola, distribusi radial, dipol.

Pilihan koordinat yang tepat sering kali menyederhanakan masalah secara drastis — persamaan yang tampak tidak terpecahkan dalam Kartesius dapat menjadi persamaan ODE biasa dalam koordinat yang sesuai.

Tabel Perbandingan Faktor Skala
Sistem $(h_1, h_2, h_3)$ $dl^2$ $d\tau$
Kartesius $(x,y,z)$ $(1, 1, 1)$ $dx^2+dy^2+dz^2$ $dx\,dy\,dz$
Silinder $(s,\phi,z)$ $(1, s, 1)$ $ds^2+s^2 d\phi^2+dz^2$ $s\,ds\,d\phi\,dz$
Bola $(r,\theta,\phi)$ $(1, r, r\sin\theta)$ $dr^2+r^2 d\theta^2+r^2\sin^2\!\theta\,d\phi^2$ $r^2\sin\theta\,dr\,d\theta\,d\phi$

Elemen volume $d\tau = h_1 h_2 h_3\,du_1\,du_2\,du_3$ dan elemen luas permukaan yang tegak lurus terhadap $\hat{\mathbf{u}}_i$ adalah $da_i = h_j h_k\,du_j\,du_k$ (dengan $i,j,k$ berbeda dan siklis).

1.7 Transformasi Vektor Satuan antar Koordinat

1.8 Tabel Ringkasan Identitas Kalkulus Vektor

1.9 Fungsi Delta Dirac dalam Tiga Dimensi

1.10 Tensor Dyadik (Pengantar)