Dokumen ini disusun untuk menjembatani konsep kualitatif menuju pemodelan komputasional formal, berfokus pada dua pilar utama: Efek Pusaran Ruang-Waktu (Relativitas Umum) dan Mekanika Pelepasan Massa Sentrifugal (Dinamika Fluida).
1. FORMULASI MATEMATIS UTAMA
A. Geometri Ruang-Waktu dan Dragging Metric (Efek Lense-Thirring)
Untuk mengukur sejauh mana putaran Matahari purba menyeret ruang-waktu di sekitarnya sehingga membentuk "saluran pusaran", kita mengadopsi Metrik Kerr untuk benda berputar dengan momentum sudut $J$:
$$ds^2 = -\left(1 - \frac{r_s r}{\rho^2}\right) c^2 dt^2 - \frac{2 r_s r a \sin^2\theta}{\rho^2} c dt d\phi + \frac{\rho^2}{\Delta} dr^2 + \rho^2 d\theta^2 + \left(r^2 + a^2 + \frac{r_s r a^2 \sin^2\theta}{\rho^2}\right) \sin^2\theta d\phi^2$$
Di mana:
$r_s = \frac{2GM}{c^2}$ (Jari-jari Schwarzschild Matahari)
$a = \frac{J}{Mc}$ (Parameter spin per satuan massa)
$\rho^2 = r^2 + a^2 \cos^2\theta$
$\Delta = r^2 - r_s r + a^2$
Kecepatan Sudut Pusaran Ruang-Waktu ($\Omega_{space}$):
Besarnya "seretan" ruang-waktu yang memaksa materi mengalir mengikuti alur lintasan dirumuskan sebagai:
$$\Omega_{space}(r, \theta) = -\frac{g_{t\phi}}{g_{\phi\phi}} = \frac{r_s r a c}{\left(r^2 + a^2\right)\rho^2 + r_s r a^2 \sin^2\theta}$$
Arti Fisika: Pada ekliptika ($\theta = \pi/2$) dan jarak dekat ($r \to \text{Merkurius}$), $\Omega_{space}$ memberikan percepatan sudut tambahan pada serpihan materi, menciptakan alur fluida yang memaksanya bergerak searah pusaran rotasi Induk.
B. Batas Keseimbangan Sentrifugal dan Pelepasan Massa
Materi di permukaan/lapisan outer-shell Matahari purba ($r = R_\odot$) akan terlempar (ejected) jika gaya sentrifugal rotasi melampaui penarikan gravitasi interior:
$$\frac{v_{\phi}^2}{R_\odot} \ge \frac{G M(r)}{R_\odot^2} \implies \omega_{rot}^2 R_\odot^3 \ge G M(r)$$
Jika Matahari purba berputar dengan kecepatan sudut kritis $\omega_{crit} = \sqrt{\frac{GM}{R_\odot^3}}$, terjadi disintegrasi massal pada zona khatulistiwa.
Persamaan Diferensial Laju Kehilangan Massa ($\dot{M}$):
$$\frac{dM}{dt} = - \int_{\text{khatulistiwa}} \rho_{shell}(v_{\phi} - v_{escape}) A_{eff} \, d\theta$$
C. Pemodelan Pusaran Balik (Back-Eddy) Venus
Untuk memodelkan rotasi terbalik Venus, kita menggunakan persamaan Navier-Stokes untuk fluida ruang-waktu/debu terkompresi pada koordinat silinder ($r, \phi$).
Kecepatan aliran pusaran lokal $u_\phi(r)$ di antara dua saluran arus utama (Merkurius pada $r_1$ dan Bumi pada $r_2$) mematuhi kondisi sirkulasi Vortisitas ($\mathbf{\omega} = \nabla \times \mathbf{u}$):
$$\frac{\partial \mathbf{\omega}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{\omega} = \nu \nabla^2 \mathbf{\omega}$$
Ketika dua saluran $r_1$ dan $r_2$ memiliki kecepatan arus $u_1$ dan $u_2$ yang tinggi, wilayah batas di antaranya ($r_1 < r_{Venus} < r_2$) mengalami tegangan geser invers (shear-stress inversion):
$$\left( \frac{\partial u_\phi}{\partial r} \right)_{inter-channel} < 0$$
Syarat ini secara otomatis memicu pembentukan pusaran sekunder berarah negatif ($\omega_z < 0$), yang memaksa calon gumpalan materi Venus berputar berlawanan arah (retrograde).
2. SPESIFIKASI ARSITEKTUR SIMULASI FISIKA
Untuk membuktikan formulasi di atas secara numerik, berikut arsitektur simulasi yang dirancang:
A. Parameter Input & Variabel
| Parameter | Nilai Asumsi / Rentang Awal | Deskripsi |
| Massa Induk ($M_0$) | $1,0 \text{ hingga } 1,05 \, M_\odot$ | Massa Matahari Purba sebelum lontaran |
| Profil Spin ($a$) | $0,85 \text{ hingga } 0,98$ | Rotasi super-kencang mendekati batas ekstrem |
| Viskositas Medium ($\nu$) | $10^{12} - 10^{15} \, \text{m}^2/\text{s}$ | Viskositas magnetohidrodinamik plasma purba |
| Grid Koordinat | $2 \text{D} / 3\text{D}$ Mesh Relativistik | $r \in [0.1, 2.0 \, \text{AU}]$, $\phi \in [0, 2\pi]$ |
B. Logika Algoritma Simulasi (Step-by-Step)
Inisialisasi Lapangan Gravitasi: Set metrik ruang-waktu Kerr berputar cepat di pusat grid.
Pemicuan Disintegrasi: Naikkan kecepatan putar $M_0$ hingga melampaui $\omega_{crit}$, lalu lepaskan paket partikel ($N = 10^6$ partikel) dari zona khatulistiwa.
Penyaluran Arus (Vortex Tracing): Hitung lintasan partikel menggunakan Persamaan Geodesik terganggu:
$$\frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha\beta} \frac{dx^\alpha}{d\tau} \frac{dx^\beta}{d\tau} = F_{\text{viskos}}^{\mu}$$
Pemeriksaan Batas Roche (Zona Merkurius):
Jika $r < r_{\text{Roche}}$, hancurkan gumpalan, kembalikan material ke $M_0$.
Jika $r \ge 0.39 \, \text{AU}$, izinkan partikel mengalami akresi gravitasi lokal.
Kalkulasi Vortisitas Lokal (Zona Venus): Ukur momentum sudut total ($\mathbf{L}_{local}$) gumpalan materi yang terbentuk di antara saluran $r_1$ dan $r_2$.
3. HASIL YANG DIHARAPKAN UNTUK VALIDASI
Kurva Kepadatan vs Jarak: Terbukti secara numerik bahwa materi padat ber-kepadatan tinggi ($\rho > 5 \text{ g/cm}^3$) mengendap secara stabil pada $r \ge 0,39 \text{ AU}$.
Vektor Vortisitas Venus: Grafik simulasi menunjukkan tanda minus ($\mathbf{L}_z < 0$) khusus pada gumpalan materi yang terbentuk pada fase akhir di zona $r \approx 0,72 \text{ AU}$.
