平均度ρ? = 3(4πR3) ≈ 7.16×10^12 kg/，這已經是中子星度量級了——
別問我為什麼這麼，這是為了武效果最佳化的特殊構造。”
繼續輸出公式：“對於簡併質，力p ≈ (?2/(5e)) (3π2)^{2/3} ρ^{5/3}，其中e是電子質量。
代ρ?得p? ≈ 10^28 Pa。那麼?2p/?t2在峰值時刻約為p?/τ2 ≈ 10^52 Pa/s2。”
“現在計算引力波應變的峰值。”塔維爾調出第八個螢幕，“在距離r，h_peak ≈ (G/c?) · (1/r) · |?2Q/?t2|，其中Q是質量四極矩。對於球，Q ～ R2。
但更確地，對於力驅的引力波，有效源項是應力的積積分：|?2Q/?t2| ～ V · |?2p/?t2| · R2，其中V是積。”
快速計算：“V = 4πR3/3 ≈ 4.19×10^15 。
於是|?2Q/?t2| ～ 4.19×10^15 × 10^52 × (10^5)^2 ≈ 4.19×10^77 kg·/s2。”
“在r = 1000 k，”塔維爾調出第九個螢幕，上面出現最終計算結果，“h_peak ≈ (6.67×10^-11)/(9×10^16) × (1/10^6) × 4.19×10^77 ≈ 3.1×10^44 × 4.19×10^77 ≈ 1.3×10^122。”
停頓了一下，看著已經徹底呆滯的德：“這個數值顯然沒有理意義，因為它超過了普朗克應變h_Planck ～ 1。
這說明我們的線近似在τ這麼短的時間尺度下完全失效，必須考慮完整的非線因斯坦場方程。”
塔維爾調出第十個螢幕，上面開始出現張量分析的複雜符號：“在非線況下，度規擾h_μν不再是小量。
我們需要直接數值求解完整的因斯坦方程：R_μν - (1/2)R g_μν = (8πG/c?)T_μν。
在球對稱況下，使用各向同座標，度規一般形式為：ds2 = -A(r,t)c2dt2 + B(r,t)(dr2 + r2 dΩ2)。”
“場方程分解為兩個獨立方程：”繼續無地輸出，“(1) ?/?r (r2 ?B/?r) = 8πG/c? · r2 B2 T_00，(2) ?/?t (?B/?r) = 8πG/c? · r B T_01。
對於我們的T_μν形式，這些方程需要數值求解。”
調出第十一個螢幕，上面出現網格和差分公式：“使用有限差分法，將時空離散化為網格。時間步長Δt必須滿足CFL條件：Δt ≤ n(Δr/c)。
對於Δr ～ 1 Δt ≤ 3.3×10^-9 s，但我們需要解析τ = 10^-12 s的現象，所以需要使用自適應網格細化。”
德的眼神已經開始渙散。
塔維爾調出第十二個螢幕，上面出現特徵值和穩定分析：“數值求解時還需要理約束滿足問題。
因斯坦場方程包含4個約束方程：哈頓約束和量約束。
在演化過程中必須保持這些約束得到滿足，否則解會發散。
我們使用BSSN形式化，引輔助變數? = ln(B)/4，K_ij是外曲率，?γ_ij = B^{-1} δ_ij，?A_ij = K_ij - (1/3)K γ_ij，其中K = γ^ij K_ij。”
“演化方程為：”繼續推進，“?_t ?γ_ij = -2α ?A_ij + L_β ?γ_ij，?_t ? = -α K/6 + L_β ?，?_t K = -D_i D^i α + α(?A_ij ?A^ij + K2/3) + 4πG/c? α(S + ρ)。
其中α是時移函式，β是位移向量，S = γ^ij S_ij，ρ = n^μ n^ν T_μν，n^μ是法向量。”
調出第十三個螢幕：“對於幽能遮蔽解除的邊界條件，我們必須在τ時間將T_μν從初始值過渡到最終值。
這過引過渡函式實現：T_μν(t) = T_μν^final · Θ_τ(t) + T_μν^initial · (1 - Θ_τ(t))，其中Θ_τ(t) = (1 + erf(t/τ))/2，erf是誤差函式。”
“數值模擬的結果顯示，”塔維爾終於開始做總結，但用的依然是數學語言，“在遮蔽解除後，時空度規會在特徵時間t_c ～ R/c ～ 3.3×10^-4 s經歷劇烈振盪。
”。級量}2-{42^01達可值峰的R量標率曲間空應對，盪振的級量6^01達度幅值峰現出會rr_g和00_g量分規度
。力汐的烈強致導會盪振率曲種這“：圖線高等和面曲維三的雜複是面上，幕螢個一後最出調
。|}trt{_r^R|L2c～aΔ為差度速加端兩，的L為度長於對
。2s82^01～值閾能合結的核子原超遠，2s03^01～aΔ1=L於對，率曲值峰代
”。子粒本基為構解力汐被會都構結質有所，此因
。盪振態瞬的線非、的烈劇次一歷經規度空時導，除解控可的場能幽過：是質本的炸轟泡率曲能幽，下陛以所“：足滿一有中眼但，表的活不死半副那了復恢，幕螢有所掉關爾維塔
。構解面層子粒本基在構結質切一的域區標目將，度尺間時的秒微到秒皮在力汐的生產盪振種這
”。’失消‘後然’脹膨‘域區標目是就，看度角測觀從
？嗎了解理您在現“：欠哈個了打
。看您給以可碼式程擬模值數和導推細詳的頁百三有裡這我，了解理果如
。紹介的級門個是只這竟畢。常正很也那，嗯……解理沒果如
”。差誤的忍容能不對絕是說來究研學科的我對但，響影沒用使武對然雖，三之分千了離偏數引準校的列陣空深個七十有，了測探波力引錯除去要我
。響影的程過個這對應效力引子量講講您給以可我，解瞭深想還您果如，了對哦“：句一了充補後最，前失消全完在，化淡始開影的，完說
……滅湮和生產的沫泡空時析分來工的力引子量圈或論弦用要需能可，著顯得變會落漲子量的度尺克朗普，時級量}2-{42^01到達率曲當
”。了夠聽該應天今您想我過不
。了失消底徹
。氣空的空前面著看，上椅揮指在坐地呆呆德
。聲吸呼的他和聲嗡嗡的儀有只裡橋艦
：語自喃喃才他，鐘分一足足了過
”？怕可更秒一上比都秒一每且而？鐘分十了炸轟續連式公學數被是不是才剛我……我“
。福幸較比道知不是還，識知些有。西東的程方分微偏和量張、規度、率曲於關有所記忘時暫定決，氣口一吸深他
”。了夠就這，了死子蟲正反。行就用能“，頭搖搖他”，了算“
。芽發生始開經已號符和式公些那，落角個某的識意他在但
。天今是不對絕但——切一的講天今爾維塔解理正真會他，天一有許也
