--- title: "GMDIS-06｜多介質耦合動力學：固體—空氣—聲場—熱場、界面功率與能量一致性" english_title: "Multi-Medium Coupled Dynamics: Solid–Air–Acoustic–Thermal Interaction, Interface Power, and Energy Consistency" series: "GMDIS — Global Multiscale Dynamics and Informational Sufficiency Theory" paper: "06 / 14" version: "v0.1" date: "2026-10-16" locale: "zh-TW" author: "Neo.K" research_collaboration: "Aletheia / GPT-6" institution: "EveMissLab / 一言諾科技有限公司" status: "Detailed Research Manuscript — Model-specific analytical results; physical validation pending" license_note: "研究原稿；授權方式以權利人後續決定為準" canonical_format: "UTF-8 Markdown; mathematics enclosed only by $...$ or $$...$$" predecessors: ["GMDIS-01", "GMDIS-02", "GMDIS-03", "GMDIS-04", "GMDIS-05"] successor: "GMDIS-07 — 持續作用與動態穩定結構" --- # GMDIS-06｜多介質耦合動力學 ## 固體—空氣—聲場—熱場、界面功率與能量一致性 **Multi-Medium Coupled Dynamics: Solid–Air–Acoustic–Thermal Interaction, Interface Power, and Energy Consistency** **系列：** 全域多尺度動力學與資訊充分性理論（GMDIS） **篇次：** 06／14｜**版本：** v0.1｜**日期：** 2026-10-16（Asia/Taipei） **作者：** Neo.K｜**AI 協力：** Aletheia（GPT-6）｜**機構：** EveMissLab／一言諾科技有限公司 **機構：** EveMissLab / 一言諾科技有限公司 **日期：** 2026-10-16 **版本：** v0.1 --- ## 摘要 本研究承接 GMDIS-01 至 GMDIS-05，建立一套可由具體力學方程、界面條件與能量恆等式核對的多介質耦合動力學研究層。前篇 GMDIS-05 已以固定弦與集中振盪器證明一個線性場—集中系統的適定性與能量估計；本篇將問題推進至真實節拍器世界可能涉及的彈性固體、可壓縮空氣、聲學壓力、流體動量、熱交換、擒縱驅動及接觸事件。此處「多介質」不是將所有影響統稱為勢能，更不是將 SIO 資訊纖維宣稱為物理微觀實體，而是明定各物理子域的狀態型別、幾何域、耦合界面、守恆／耗散帳本與可觀察量。 本篇的直接數學核心是**一維聲腔—集中機械振子之功率守恆耦合模型**。對聲壓 $p$ 、聲學質點速度 $v$ 、活塞／等效節拍器坐標 $q$，在小振幅固定參考域假設下，以一維線性聲學方程與邊界速度匹配構造介面。證得理想耦合時機械與聲學的邊界功率精確抵銷；加入機械阻尼、聲學體積阻尼、集總熱容及熱外洩後，得到封閉且單位一致的總能量帳本。進一步建立能量相容的離散聲腔與集中振子的有限維矩陣模型，給出對稱的離散交換算子及隱式中點時間積分，並在附錄提供可獨立執行的 NumPy 測試。這些結果是**指定模型下的解析推導與數值測試**，不是 100 或 1,000 台實物節拍器的已證實結果。 對較完整的固體—流體—熱系統，本文區分聲學線性化近似、可壓縮 Navier–Stokes–Fourier、彈性／黏彈性固體、移動介面 FSI、熱黏性邊界層與非光滑擒縱接觸，並說明為何不能同時把全流體總能量與其聲學擾動能當成兩份獨立儲能。本篇提出介面功率不變量、能量殘差、跨尺度擾動敏感度、參數無因次群及逐層增模／消融契約，作為未來 GMDIS-07 穩定性、GMDIS-08 跨尺度記憶、GMDIS-11 反事實及 GMDIS-12／13 大規模基準測試的物理骨架。 **核心限制：** 強作用不必定支配長期結果；弱共振、記憶、幾何、摩擦與熱邊界都可能改變可行未來。GRST 的 Shi 仍是「規則誘導之未來空間非均勻結構」，不是焦耳、牛頓或物理勢能。SIO 描述的合法自由、構形約束與非局部陰影，亦須與有限速度的物理信號傳播區分。多場建模是否真的改善同步狀態、恢復率及反事實預測，必須接受既有物理模型的強基線與資料檢驗。 **關鍵詞：** GMDIS；多介質耦合；彈性介質；流體—固體交互作用；聲學；能量流；熱力學；界面功率；有限元／有限體積；離散功率埠；節拍器；GRST；SIO；資訊充分性；72 格分類。 --- # 0. 本篇在系列中的責任與證據邊界 ## 0.1 直接承接的五項契約 - **GMDIS-01：** 從全域物理世界生成可能的集體運動，不預設「同步」是唯一結局；世界模型應能生成反相、分群、多穩態、混沌與停止等候選結果。 - **GMDIS-02：** 不把相位或吸引子標籤當作完整動力學；觀測充分性必須相對於指定未來任務、干預集合、時間窗與容許誤差。 - **GMDIS-03：** SIO 的**硬約束自由度陰影**與由物理作用導出的**動態影響陰影**分開；真實作用具有邊界、通道與傳播尺度。 - **GMDIS-04：** 牛頓力、應力、動量、功率、能量、耗散、熱與廣義張力各自擁有不同物理型別；GRST 勢不可當作新能量或新基本力。 - **GMDIS-05：** 要先指定能量空間、跡與界面算子；弱形式、適定性及數值近似不可脫離模型假設直接宣稱。 本篇特意不再複寫上述系列的整部母論；新增責任是把「空氣、聲波、固體及熱」寫成可檢查的物理子系統，以及揭露它們跨界的動力後果。 ## 0.2 證據標籤 | 標籤 | 意義 | 本文使用規則 | |---|---|---| | D-06 | 定義／表示選擇 | 只在明定型別與語境內成立，不冒充經驗事實 | | M-06 | 數學命題 | 在列出的方程、參數、正則性與邊界條件下推導 | | C-06 | 構造性反例 | 顯示某種一般性推論不成立，不代表任何真實裝置必然如此 | | N-06 | 數值一致性測試 | 可以執行並重算，不等於外部實物驗證 | | H-06 | 新增的可反駁研究假說 | 需要基線、消融、資料與判敗標準 | | L-06 | 學術既有方法 | 必須標示作者、年代及其適用的問題域 | | Future | 後續研究義務 | 不得由本篇模型自動推導成已證結論 | ## 0.3 最低成功與明確失敗 本篇最低成功不是讓所有場「統一」，而是完成：(i) 一個跨介質能量正確的可執行最小模型；(ii) 一份從線性聲學到真實流體—固體的正確模型提升路線；(iii) 能讓物理實驗判定何時可以忽略聲音、空氣與熱、何時不應忽略的消融準則。 若最後發現共同底座的機械作用已能完整且較簡潔地解釋所研究節拍器裝置，而聲、流、熱額外維度沒有預測增益，則相應多介質新增模型必須降級為非必要修飾，不能硬保留「萬象張力」之名詞。 --- # 1. 世界的分域與型別：什麼才算一個多介質物理模型？ ## 1.1 D-06-01：分域與參考構形 令世界的物理研究區域含： $$ \Omega=\Omega_s\cup\Omega_f\cup\Gamma_{sf}, $$ 其中 $\Omega_s$ 是固體（含底板、支架等）， $\Omega_f$ 是空氣或其他流體， $\Gamma_{sf}$ 是共同物理界面。有限個節拍器內部坐標採 $q\in\mathbb R^N$ ；必要時細分擒縱機構的連桿與離散模式 $\xi\in\Xi$。 上述集合只是**瞬時物理分域**；如果界面位移不再可忽略，必須寫 $\Omega_f(t)$ 、 $\Omega_s(t)$ 、 $\Gamma_{sf}(t)$ 並配備參考構形到實際構形的映射。不能把移動壁面問題的微分方程直接放在固定域上，卻沒有小振幅線性化理由。 **建模型別：** | 物理分量 | 變數 | SI 單位 | 典型儲能或守恆角色 | 模型依賴 | |---|---|---|---|---| | 集中機構 | $q,\dot q$ | m, m/s；角座標則 rad | 機械動能、位能、外加驅動 | 質量、幾何、擒縱契約 | | 固體位移／速度 | $u_s,\dot u_s$ | m, m/s | 動量、彈性或黏彈性能 | 材料及邊界 | | 固體應力 | $\sigma_s$ | Pa | 牽引與功率通道 | 構成關係 | | 空氣密度／流速 | $\rho_f,v_f$ | kg/m³, m/s | 質量及動量輸運 | 可壓／不可壓、黏性 | | 聲壓擾動 | $p$ | Pa | *線性化擾動*能量 | 小擾動、背景靜止或明示線性化 | | 聲學速度擾動 | $v_a$ | m/s | *線性化擾動*動能 | 不等於完整流場必有額外能量 | | 溫度與熱通量 | $T,j_h$ | K, W/m² | 內能、傳熱與熵 | 熱容、導熱率、本構 | | 接觸／擒縱 | $\xi,J$ | 標記、N·s | 切換、衝量、可能耗散／補能 | 接觸規則與重設映射 | 表中的 $\rho_f$ 是**質量密度**，不能與 GRST 的任何資訊密度混寫。聲學 $p$ 是對背景壓力的微小偏差，不是熱力學總壓力。 ## 1.2 D-06-02：模型級別，而非所有變數永遠都要同時計算 我們定義下列有序的**近似家族**（非宇宙本體層級）： | 層 | 模型標記 | 保留的作用 | 明確捨去的作用 | |---|---|---|---| | M0 | 集中機械 | 擺、底座、彈簧、摩擦、擒縱 | 空間聲場與完整熱場 | | M1 | 機械＋理想線性聲學 | 聲壓、流體速度、壁面功率 | 黏性熱邊界層與大振幅流動 | | M2 | M1＋有效阻尼＋集總熱容 | 聲／機械耗散去向及環境熱交換 | 溫度空間梯度及熱回饋的完整本構 | | M3 | 黏性可壓縮流體＋熱傳 | 黏性、傳熱、非線性對流、可能流動噪聲 | 若另有其他精細機構仍可近似 | | M4 | 真實移動 FSI＋黏彈性＋擒縱接觸 | 幾何非線性、材料與接觸事件 | 未實證的更細物理保持未定 | **關鍵：M1 的聲學場是從 M3 在適當條件下抽出的擾動近似，不是要把兩套方程疊到一起再加總能量。** 熱聲耦合升級時，要重新界定哪些量是背景、哪些是擾動、哪些是內能，以避免雙重計帳。 ## 1.3 空間尺度與傳播延遲 材料與聲速會建立有限傳播速度的有效尺度。典型聲傳時間 $t_a=L/c$ ；結構彈性波有自己的傳播速度（取決於材料及模態），熱傳 Fourier 方程在數學上屬拋物型且具有理想化的瞬時非零影響尾；因此不能把「熱擴散方程的支撐」與「實驗中瞬時可檢測的訊號」混同，更不能由 SIO 的關係陰影推出超光速或無延遲的物理作用。 對有限時窗 $[0,T]$ 研究動態影響範圍時，應同時標記：(i) 實際因果傳播機制；(ii) 有效可測門檻；(iii) 觀察裝置帶寬與噪聲；(iv) 所採連續近似。 --- # 2. 最低可驗證世界：一維活塞—聲腔—機械振子 本篇主張先由一個**確切能量相容的場—集中模型**出發，再談更難的三維世界。它具有真實空氣中的聲學傳播，但不假裝等於完整擒縱節拍器。 ## 2.1 D-06-03：假設 1. 流體是均勻、靜止背景空氣；平衡密度 $\rho_0>0$ 、小訊號聲速 $c>0$ 、體積彈性模數 $K_a=\rho_0c^2$，均視作常數。 2. 聲腔初始位置為 $x\in(0,L)$，截面積 $A>0$，振幅小於特徵長度；對聲場使用**固定參考區域的一階線性化**。 3. $x=0$ 處活塞的坐標 $q$ 朝 $+x$ 指向聲腔內部，活塞正速度 $\dot q$ 使空氣流向 $+x$ ； $x=L$ 為不動的理想剛性壁。 4. 活塞等效質量 $m>0$ 、線性彈簧剛度 $k>0$ 、機械黏性阻尼 $\beta\ge0$ 、外力 $F(t)$ （均沿 $+x$ 定正）。 5. 聲學速度 $v(x,t)$ 、聲壓擾動 $p(x,t)$ 取小訊號一階系統；無流動穿透剛性壁，不含隨機或量子作用。 6. 不將小振幅的固定域近似外推至活塞大位移、湍流、分離流或聲學強非線性。 ## 2.2 D-06-04：聲場方程與機械方程 聲學一階方程： $$ \partial_t p+K_a\partial_xv=0, $$ $$ \rho_0\partial_tv+\partial_xp=0. $$ 聲學邊界速度： $$ v(0,t)=\dot q(t),\qquad v(L,t)=0. $$ 集中機械動力學： $$ m\ddot q+\beta\dot q+kq=F(t)-A p(0,t). $$ **正負號不是裝飾。** 壓力若為正，空氣將活塞往 $-x$ 推，故集中方程中為 $-Ap(0,t)$。對 $x=0$ 的流體域，其外法向指向 $-x$，使流體透過活塞收到的邊界功率為 $+Ap(0)\dot q$，兩者正好相反。 量綱檢查： $\rho_0\partial_tv$ 與 $\partial_xp$ 均為 $\mathrm{Pa/m}$ ； $K_a\partial_x v$ 與 $\partial_tp$ 均為 $\mathrm{Pa/s}$ ； $A p$ 為 N。完整模型不包含「規則勢等於 $p$ 」的任何型別等價。 ## 2.3 D-06-05：聲學擾動能與機械儲能 定義聲學能量： $$ E_a(t)=A\int_0^L\left(\frac{p(x,t)^2}{2K_a}+\frac{\rho_0v(x,t)^2}{2}\right)dx. $$ 定義機械能： $$ E_m(t)=\frac m2\dot q(t)^2+\frac k2q(t)^2. $$ 第一個 $p^2$ 項代表對背景態的二次小擾動壓縮能；它**不是**空氣完整熱力學內能的第二份額外儲能。只有在 M1 模型內，它才能被當作獨立的聲學總能成分。 ## 2.4 M-06-01：聲場局部能量守恆 **命題。** 若 $p,v$ 足夠光滑並滿足 2.2 的無耗散聲學方程，則： $$ \partial_t\left(\frac{p^2}{2K_a}+\frac{\rho_0v^2}{2}\right)+\partial_x(pv)=0. $$ **證明。** 第一條方程乘以 $p/K_a$，第二條乘以 $v$，兩式相加： $$ \frac p{K_a}\partial_tp+\rho_0v\partial_tv+p\partial_xv+v\partial_xp=0. $$ 前二項是能量密度的時間導數；後二項恰為 $\partial_x(pv)$。證畢。 對固定參考域積分、利用 $v(L)=0$ 及 $v(0)=\dot q$： $$ \frac{dE_a}{dt}=Ap(0,t)\dot q(t). $$ 注意正能流不是「所有空氣都做正功」，它的正負依瞬時壓力與速度乘積而變；一週期內亦可能正負交換。 ## 2.5 M-06-02：理想介面總能量精確對帳 機械方程乘以 $\dot q$，得： $$ \frac{dE_m}{dt}=F\dot q-\beta\dot q^2-Ap(0)\dot q. $$ 因此： $$ \boxed{\frac{d}{dt}(E_m+E_a)=F\dot q-\beta\dot q^2.} $$ **證明。** 由 M-06-01 的 $+Ap\dot q$ 與集中系統的 $-Ap\dot q$ 完全相消，剩餘外力輸入與機械阻尼耗散。證畢。 這是本篇第一個強制不變量：**界面功率在整體能量帳本內不可被重算成創生項或損失項。** 若計算所得的聲學子域與機械子域同時獲得 $+Ap\dot q$，通常是法向或耦合符號出錯；若兩邊都刪去，又無法追蹤各自的介質能量交換。 ## 2.6 原理的邊界 - 上述是在光滑解或可合法處理邊界跡的適定弱解下成立的能量恆等式；**未**在此證明任意低正則初態下的存在唯一性。 - 若外部聲學輻射或開放邊界存在， $x=L$ 的零聲速壁條件不再成立，右端必須增加對外聲功率。 - 若活塞位移很大、流體具有背景流、壓力不再可線性化，應採真實移動控制體積的能量式，不能繼續無條件使用 $E_a$。 - 即使 $F=\beta=0$ 而 $E_m+E_a$ 守恆，也**不**推出兩個以上振盪器會同步。 --- # 3. 加入聲學耗散與熱：耗散並不等於能量消失 ## 3.1 D-06-06：有效聲學阻尼 以常數 $r\ge0$ 表示**理想化的分布式線性拖曳**，將動量式改為： $$ \rho_0\partial_tv+\partial_xp+r v=0. $$ $r$ 的單位為 $\mathrm{kg/(m^3s)}$，使 $rv$ 與壓力梯度具有同一單位。它是一個可重算的**有效模型參數**，不是由完整 Navier–Stokes 黏性本構直接推導出的普遍「空氣阻力常數」。特別是狹縫黏性、熱邊界層及聲輻射通常有頻率依賴與空間結構，需要進一步升級模型。 利用和 M-06-01 相同的乘法： $$ \frac{dE_a}{dt}=A p(0)\dot q-Ar\int_0^L v^2\,dx. $$ 聲學機械擾動能被移出，但在完整熱力學世界中，這些能量可能成為內能或穿過外部邊界；不能說「弱的空氣勢能被大的節拍器勢能吃掉」。 ## 3.2 M-06-03：帶分布阻尼的機械—聲學帳本 若 $F$ 為外力，得到： $$ \frac{d}{dt}(E_m+E_a) =F\dot q-\underbrace{\beta\dot q^2}_{D_m} -\underbrace{Ar\int_0^L v^2\,dx}_{D_a}. $$ 其中 $D_m,D_a\ge0$，其型別皆為 W。聲場吸收機械能有時是瞬時能量儲存，並非每個時刻都為耗散；真正的非負耗散由所選阻尼／熱黏性模型明定。 ## 3.3 D-06-07：集總熱儲能 令熱容 $C_T>0$ （J/K）、相對環境溫度 $\vartheta(t)=T(t)-T_\infty$ （K）、環境散熱係數 $h\ge0$ （W/K）。令 $Q_T(t)$ 為額外輸入熱功率（W）。**在指定的能量劃界內**，若 $D_m$ 與 $D_a$ 最終全部進入被模型計入的熱節點，可以建立： $$ C_T\dot\vartheta=D_m+D_a-h\vartheta+Q_T. $$ 定義相對參考溫度的熱內能差： $$ E_T=C_T\vartheta. $$ 則： $$ \boxed{\frac{d}{dt}(E_m+E_a+E_T)=F\dot q+Q_T-h\vartheta.} $$ 這是 M-06-04（模型內第一定律）。其中 $E_T$ 可以為負，因為它是**相對參考態的內能差**，不是要求絕對內能必須大於零的候選儲能泛函。若要討論熱力學熵，還需要絕對溫度 $T>0$ 及適當的熱傳本構。 **兩項限制：** (a) 若 $r$ 是「聲波向外輻射」而非被內部吸收，對外功率必須另列，不能全部灌入 $C_T$ ；(b) 若 $\rho_0,c,r,\beta$ 隨溫度改變，此簡化方程的雙向熱—聲／熱—機械作用需要額外推導，不可直接把它稱作已封閉的完整熱聲理論。 ## 3.4 熱被移出不等於熱量消失 設 $\vartheta>0$ 時由系統向環境放熱 $h\vartheta$。若將環境也算入閉合世界，熱外洩應以環境能量增加記帳；在開放子系統中則是負的邊界通量。研究邊界擴張時，**同一個 $h\vartheta$ 應由「外流」重新改寫為「內部轉移」，不能同時在兩邊都加一次**。 一旦定義外界，GRST 世界線的比較條件就必須包括這個熱邊界。不同的系統劃界會改變能量子帳本，但不應無理由地改變同一封閉物理世界的總記帳。 --- # 4. 從一維聲學走向連續彈性固體 ## 4.1 D-06-08：固體場的最低線性模型 在參考固體域 $\Omega_s$，令位移 $u_s(X,t)\in\mathbb R^d$ 、速度 $v_s=\partial_tu_s$ 、密度 $\rho_s>0$。小應變為： $$ \varepsilon(u_s)=\frac12\big(\nabla u_s+(\nabla u_s)^\top\big). $$ 對彈性張量 $\mathbb C$ 和黏性張量 $\mathbb D$，最低線性黏彈模型採： $$ \sigma_s=\mathbb C:\varepsilon(u_s)+\mathbb D:\varepsilon(v_s), $$ $$ \rho_s\partial_{tt}u_s=\nabla\cdot\sigma_s+f_s. $$ 假定 $\mathbb C$ 滿足適當對稱與正定， $\mathbb D$ 為半正定耗散算子，且位移邊界約束足以排除未固定的零能剛體模態（或另外用商空間表示）。不得因為「幾何是資訊的顯影」就省略真實物理幾何及材料張量。 ## 4.2 D-06-09：固體能量 定義： $$ E_s=\frac12\int_{\Omega_s}\rho_s|v_s|^2\,dV +\frac12\int_{\Omega_s}\varepsilon(u_s):\mathbb C:\varepsilon(u_s)\,dV. $$ 在足夠光滑的情況下，對運動式乘 $v_s$ 並積分可得： $$ \frac{dE_s}{dt} =\int_{\Omega_s}f_s\cdot v_s\,dV +\int_{\partial\Omega_s}(\sigma_s n_s)\cdot v_s\,dS -\int_{\Omega_s}\varepsilon(v_s):\mathbb D:\varepsilon(v_s)\,dV. $$ **M-06-05。** 上式在小應變、線性黏彈與適當邊界條件下成立。最後一項非負，表示機械耗散功率。它是否全部轉為模型中的固體內熱，取決於所選熱力學閉合及能量劃界。 ## 4.3 固體內應力與聲壓不是同一個量 聲場的 $p$ 是標量壓力擾動；固體的 $\sigma_s$ 是二階 Cauchy 應力張量，除純靜水壓特殊情況外還包含剪應力。即使兩者單位同為 Pa，也沒有普遍等號。跨介面需使用**法向牽引力**及**界面運動速度**，而不是把 $p$ 與 $\sigma_s$ 用未定義的代數直接相加。 ## 4.4 二維／三維幾何何時真正進入因果鏈？ 對只允許同向平移的一塊理想剛體底座，若各節拍器質量、方向和耦合相同，單純交換它們在底座上畫出的座標不必改變模型。若底座允許轉動，力矩 $\tau=r\times F$ 使施力位置成為動力變數；若允許彎曲與扭轉，載荷投影到不同結構模態，其空間耦合核會不同。 這是 GMDIS-03 的**幾何零效應基線**，本篇繼續強制保留。不能先假定 10×10 的排列具有物理影響，再用已出現的同步畫面逆推出「位置規則勢」。 --- # 5. 從聲學近似到真正可壓縮、黏性且傳熱的空氣 ## 5.1 為什麼理想聲學方程需要物理上游 理想聲學是對某類穩定背景下**小擾動**的有效方程，具有清楚的波動速度和近似能量守恆式；但它沒有涵蓋流體慣性對流、聲致流、黏性邊界層、紊流、熱傳與非線性聲波。在實際節拍器周圍，是否需要更高階模型由聲學馬赫數、幾何尺度、邊界層厚度、振幅及驗證資料決定，而不是由理論偏好決定。 ## 5.2 D-06-10：可壓縮黏性流體的最低守恆方程 在流體域（固定 Eulerian 空間表示下），選密度 $\rho$ 、速度 $v_f$ 、總壓力 $P$ 、單位質量內能 $e$ 、溫度 $T$。質量守恆： $$ \partial_t\rho+\nabla\cdot(\rho v_f)=0. $$ 動量平衡： $$ \partial_t(\rho v_f)+\nabla\cdot(\rho v_f\otimes v_f) =\nabla\cdot\sigma_f+\rho f_f. $$ Cauchy 應力分解： $$ \sigma_f=-P I+\tau_f, $$ $$ \tau_f=2\mu D(v_f)+\lambda(\nabla\cdot v_f)I, $$ $$ D(v_f)=\frac12\big(\nabla v_f+(\nabla v_f)^\top\big). $$ $\mu\ge0$ 為動力黏度（Pa·s）， $\lambda$ 為第二黏性係數；為保證黏性功率非負，三維等向 Newtonian 本構需滿足 $\lambda+2\mu/3\ge0$。如採更一般材料，條件應寫在 $\tau_f:D(v_f)\ge0$ 的耗散性上。 內能平衡在足夠光滑情況下為： $$ \rho\frac{De}{Dt} =-P\nabla\cdot v_f+\tau_f:\nabla v_f -\nabla\cdot j_h+\rho r_h, $$ 其中 $D/Dt=\partial_t+v_f\cdot\nabla$ 為物質導數， $r_h$ 是單位質量的外部熱源率。配合 Fourier 導熱： $$ j_h=-\kappa_T\nabla T, $$ 以及物態方程 $P=P(\rho,T)$ 、熱容關係 $e=e(\rho,T)$，才是一組候選閉合系統。不同的狀態方程、非牛頓流體或熱聲效應需分別建模。 **單位界線：** 此處 $P$ 是總熱力學壓力，不等於第二節以 $p$ 表示的擾動壓力。將完整流體的總能量 $\int\rho(e+|v_f|^2/2)\,dV$ 與由其線性化提取的 $E_a$ 無條件相加，會造成雙重計數。 ## 5.3 第二定律的必需條件 取 Fourier 導熱與上述 Newtonian 黏性，在合適局部平衡熱力學假設下，體積熵生成候選形式為： $$ \sigma_{\mathrm{ent}} =\frac{\tau_f:D(v_f)}{T} +\frac{\kappa_T|\nabla T|^2}{T^2}\ge0, $$ 其中 $T>0$ 為絕對溫度，且本構滿足正耗散條件。若熱源與開放質量流存在，需在完整熵平衡中加上相應通量和外源，不能單憑上式涵蓋所有非平衡過程。 這個第二定律條件與 M-06-02 的「擾動機械能非增」不同。即使某個子系統機械能增加，也可能因外界驅動而保持完全合法的熱力學帳本；同樣，總能量守恆亦不表示熵不增加。 ## 5.4 聲學如何從流體方程近似得到？ 令： $$ \rho=\rho_0+\rho',\qquad P=P_0+p,\qquad v_f=v_a, $$ 背景靜止，擾動微小，且取有效等熵小訊號關係： $$ p=c^2\rho'. $$ 在忽略二階擾動、黏性、熱擴散與外力之候選近似下，連續性方程轉為： $$ \partial_t\rho'+\rho_0\nabla\cdot v_a=0, $$ 動量式轉為： $$ \rho_0\partial_t v_a+\nabla p=0. $$ 合併得到第二節的 $\partial_tp+\rho_0c^2\nabla\cdot v_a=0$。此步表明**M1 是 M3 在限制條件下的一種降階候選，而非與 M3 平行存在的額外物質。** ## 5.5 黏性與熱邊界層為何可能比體積平均更重要？ 接近固體壁面，無滑移速度與溫度邊界可能形成很薄的黏性／熱邊界層。若驅動角頻率 $\omega>0$，黏性運動學係數 $\nu=\mu/\rho_0$，熱擴散率 $\alpha_T=\kappa_T/(\rho_0 c_p)$，常見線性尺度為： $$ \delta_v=\sqrt{\frac{2\nu}{\omega}},\qquad \delta_T=\sqrt{\frac{2\alpha_T}{\omega}}. $$ 這些**長度尺度**不是「新的勢能」；其意義是判定是否需解析近壁面梯度、改用有效阻抗邊界，或測量其對聲功率耗散與局部流動的修正。若 $\delta_v/L$ 很小，未必代表黏性效應對高品質因子共振腔可忽略；體積平均小不等於長時間累積小。Bach 與 Bruus（2018）及後續熱黏聲學研究提供成熟對照，[E-06-04]、[E-06-05]。 ## 5.6 不可壓縮近似不是「音波仍照原樣存在」 若從完整流體改採嚴格不可壓縮約束 $\nabla\cdot v_f=0$，壓力通常成為維持約束的乘子。此極限可能消去可壓縮聲波的獨立自由度；因此不能同時宣稱「嚴格不可壓縮」又使用第二節 $p^2/(2K_a)$ 當成同一流體的完整聲學儲能。若要混合不可壓縮流體區與可壓縮聲區，必須建立合法的多域耦合模型。 --- # 6. 固體—流體介面：作用真正在哪裡交換？ ## 6.1 D-06-11：介面的方向與運動契約 令 $n_f$ 為**流體區向外**單位法向量， $n_s$ 為**固體區向外**單位法向量，在同一光滑介面滿足 $n_s=-n_f$。若黏性流體與固體無滑移接觸，移動介面上的速度匹配為： $$ v_f\big|_{\Gamma_{sf}(t)}=\partial_tu_s\big|_{\Gamma_{sf}(t)}. $$ 若採無黏的理想聲學耦合，通常只需法向速度匹配： $$ v_a\cdot n_f=\partial_tu_s\cdot n_f. $$ 此時流體壓力為法向牽引力，切向接觸法則另論。全黏性應力平衡的理想零質量介面為： $$ \sigma_f n_f+\sigma_s n_s=0. $$ 此式的物理意義是同一內部介面的雙向牽引互為相反數，並不是說材料的應變、流速與溫度都相等。若允許界面滑移、滲透、表面張力、膜質量或黏著層，必須改寫這組契約。 ## 6.2 M-06-06：介面功率相消 在同一物理介面上，若速度連續且牽引互為相反數，則： $$ \int_{\Gamma_{sf}}(\sigma_f n_f)\cdot v_f\,dS +\int_{\Gamma_{sf}}(\sigma_s n_s)\cdot v_s\,dS=0. $$ **證明。** 在介面上 $v_f=v_s=v_\Gamma$ ；代入 $\sigma_f n_f=-\sigma_s n_s$，被積函數逐點相消。證畢。 這不意味兩側實際「沒有相互作用」。它只說在**將兩側都列為內部子系統、且沒有額外儲能或耗散的理想介面**上，力學功率的正負帳目必須合計為零。 ## 6.3 一個具儲能的界面不能仍套用零總功率 若界面有質量 $m_I$ 、彈性儲能 $V_I$ 、黏性耗散 $D_I\ge0$，則兩側輸入界面之淨功率不必直接為零，而應成為： $$ P_{s\to I}+P_{f\to I}=\frac{d}{dt}(K_I+V_I)+D_I. $$ 若介面還接有控制器或擒縱器，須加其能量埠。這個擴充是防止數值模型出現假耗散或假驅動的最低條件。將理想介面突然改為彈性膠層，但不增加能量節點，會使守恆帳本必然缺項。 ## 6.4 移動邊界需要 Reynolds 傳輸或 ALE 形式 如果 $\Omega_f(t)$ 真正移動，對任意能量密度 $e_f$ 的總時間導數需區分局部時間變化與邊界輸運。令界面／計算網格速度 $w_\Gamma$，在正規條件下 Reynolds 傳輸定理給出： $$ \frac{d}{dt}\int_{\Omega_f(t)}e_f\,dV =\int_{\Omega_f(t)}\partial_te_f\,dV +\int_{\partial\Omega_f(t)}e_f(w_\Gamma\cdot n_f)\,dS. $$ 在 ALE 網格表示中，對流相對速度是 $v_f-w_{\mathrm{mesh}}$，而非永遠 $v_f$。若忽略幾何守恆律（geometric conservation law）、移動控制體積與網格速度，甚至常數物理場也可能出現假的能量變化。既有 ALE 文獻 [E-06-06]、[E-06-07] 是必須超越的基線，而非 GMDIS 自行發明的算法。 ## 6.5 附加質量（added mass）：空氣不能只有「力很小」一個判準 流體在固體加速時可被迫共同運動，使固體顯示出額外慣性負擔；其效應與流體密度、構形、邊界條件及頻率有關。Causin、Gerbeau、Nobile [E-06-08] 已指出，某些 FSI 分割算法因 added-mass coupling 出現數值不穩定，哪怕連續物理模型本身並非失穩。 因此將聲固模型分成「先算空氣，再算底座」的時間步流程，不能直接被認為物理上正確。必須同時檢查時間步依賴、介面殘差、同步迭代及能量帳本。 --- # 7. 最低離散多介質模型：聲腔體積單元＋機械活塞 本節是可復現模型的工程核心。為避免高維移動域 PDE 的適定性問題與大規模數值工程在同一步驟混雜，採用**固定參考域、一維線性聲場的功率相容半離散**，並驗證它的能量恆等式。此模型不聲稱提供三維 Navier–Stokes—熱彈性—非光滑接觸的存在定理。 ## 7.1 D-06-12：離散狀態、能量與耦合矩陣 取 $n$ 個聲壓自由度 $p\in\mathbb R^n$， $n-1$ 個內部面上定義的速度自由度 $v\in\mathbb R^{n-1}$，以及集中機械坐標 $q,\dot q\in\mathbb R$。 令 $C_a\in\mathbb R^{n\times n}$ 、 $M_a\in\mathbb R^{(n-1)\times(n-1)}$ 均對稱正定， $G\in\mathbb R^{(n-1)\times n}$ 為離散梯度／面通量耦合矩陣， $R_a\succeq0$ 為聲學速度阻尼矩陣， $b\in\mathbb R^n$ 為活塞聲學功率埠。 在均勻截面 $A$ 、等距聲壓單元長度 $\Delta x=L/n$ 的一階有限體積模型中，可取： $$ C_a=\frac{A\Delta x}{K_a}I_n, \qquad M_a=\rho_0A\Delta x I_{n-1}, $$ $$ G_{j,j}=A,\qquad G_{j,j+1}=-A,\qquad 1\le j\le n-1, $$ $$ b=Ae_1. $$ 其中 $e_1$ 是第一個標準基向量。 $G$ 內部面流速向 $+x$ 為正；剛性末端沒有獨立流出埠。這個配置是能量相容的**有限體積範例**，不是聲腔所有可能網格或高階方法的唯一選擇。 ## 7.2 D-06-13：半離散耦合演化方程 令 $w=\dot q$。採： $$ \dot q=w, $$ $$ m\dot w=-kq-\beta w-b^\top p+F(t), $$ $$ C_a\dot p=-G^\top v+bw, $$ $$ M_a\dot v=Gp-R_av. $$ 其功率配對為 $(-b^\top p)w$ 與 $p^\top bw$ ；兩者恰抵。內部面則是 $-p^\top G^\top v$ 與 $v^\top Gp$，再次相消。這是離散版「**同一界面只有一份功率交換**」的最低代數標準。 ## 7.3 D-06-14：離散機械—聲學總能量 $$ E_h=\frac m2w^2+\frac k2q^2+\frac12p^\top C_ap+\frac12v^\top M_av. $$ $C_a$ 的量綱是 $\mathrm{m^3/Pa}$，故 $p^\top C_a p$ 是 J； $M_a$ 具有 kg 單位，故 $v^\top M_av$ 也是 J； $b$ 具有 $\mathrm{m^2}$ 單位，故 $b^\top p$ 是 N、 $p^\top b w$ 是 W。 ## 7.4 M-06-07：半離散能量恆等式 **命題。** 若 $m,k>0$ 、 $\beta\ge0$ 、 $C_a,M_a$ 對稱正定、 $R_a$ 對稱半正定，則任意古典解滿足： $$ \boxed{\dot E_h=Fw-\beta w^2-v^\top R_a v.} $$ **證明。** 分別將機械方程乘以 $w$ 、聲壓方程左乘 $p^\top$ 、聲學速度方程左乘 $v^\top$： $$ \dot E_h=Fw-\beta w^2 -\underbrace{w b^\top p}_{\mathrm{mech\ port}} +\underbrace{p^\top b w}_{\mathrm{acoustic\ port}} -\underbrace{p^\top G^\top v}_{\mathrm{pressure\ gradient}} +\underbrace{v^\top Gp}_{\mathrm{velocity\ flux}} -v^\top R_av. $$ 由純量轉置與矩陣乘法相消，得結果。證畢。 **推論。** 當 $F=\beta=R_a=0$ 時，離散能量守恆。當 $F=0$ 且耗散半正定時， $E_h$ 非增；但此式本身**不能**證明單位會同步、能量必趨零或每個模態均被阻尼。若 $n=1$，沒有聲學內部速度自由度，模型僅為氣體可壓縮容積元件，不能用來聲稱已模擬聲波傳播。 ## 7.5 M-06-08：有限維適定性（可實際證明，而非援引完整 PDE） 將 $z=(q,w,p,v)$ 寫成線性非自治輸入系統： $$ \dot z=L_hz+f_h(t). $$ 由 $C_a,M_a$ 可逆， $L_h$ 是有限矩陣。若 $F\in L^1_{\mathrm{loc}}([0,\infty))$，則對任意初值 $z(0)$，存在唯一絕對連續解： $$ z(t)=e^{tL_h}z(0)+\int_0^te^{(t-s)L_h}f_h(s)\,ds. $$ 本命題是標準線性 ODE 結論。**它沒有證明 $h\to0$ 網格收斂到前述移動邊界連續 FSI 解，也不涵蓋帶碰撞／湍流的非線性模型。** 連續波問題的逼近、穩定性與收斂率需額外處理。 ## 7.6 D-06-15：離散熱節點與第一定律 令 $\vartheta$ 是相對環境溫度，並明定 $\beta w^2$ 與 $v^\top R_av$ 被內部熱節點吸收，則： $$ C_T\dot\vartheta=\beta w^2+v^\top R_av-h_T\vartheta+Q_T. $$ 設總相對參考能量 $E_h^\star=E_h+C_T\vartheta$，則： $$ \boxed{\dot E_h^\star=Fw+Q_T-h_T\vartheta.} $$ 若熱能部分流向未建模的接觸物、底板、聲外場，應按實際比例與熱源位置改寫；上式的獨立正確性來自**本篇明定的閉合假設**，不能把它套到所有裝置上。 ## 7.7 觀察不足的明確構造 取兩組初態： $q(0)=w(0)=0$ 、其餘機械參數相同，卻有 $p^{(1)}(0)=0$ 、 $p^{(2)}(0)=P e_1$，並令兩組 $v(0)=0$。初始節拍器的位置、速度、由這兩者推出的簡化相位都一樣，但立即加速度不同： $$ \ddot q^{(2)}(0)-\ddot q^{(1)}(0)=-\frac{AP}{m}. $$ 這個差異沒有神祕新作用：它只是**未納入機械觀察的空氣壓力狀態**。因此 GMDIS-02 的「投影不足」在本篇獲得了具有單位與功率路徑的實例；並且已由既有聲固耦合力學充分描述，不能自動聲稱 GRST 具有獨立的物理預測增益。 --- # 8. 時間離散化：連續帳本如何不被數值格式破壞？ ## 8.1 D-06-16：隱式中點法 對第 7 節線性模型，取步長 $\Delta t>0$： $$ \frac{z_{n+1}-z_n}{\Delta t}=L_hz_{n+1/2}+f_{n+1/2}, $$ $$ z_{n+1/2}=\frac12(z_n+z_{n+1}). $$ 若 $F_{n+1/2}$ 是代表該步外力的指定中點值，對二次能量的增量有精確代數恆等式： $$ E_{h,n+1}-E_{h,n} =\Delta t\left(F_{n+1/2}w_{n+1/2}-\beta w_{n+1/2}^2 -v_{n+1/2}^\top R_av_{n+1/2}\right). $$ **M-06-09。** 上式在精確線性代數求解與二次能量假設下成立；由 $x_{n+1}^\top Hx_{n+1}-x_n^\top Hx_n=(x_{n+1}+x_n)^\top H(x_{n+1}-x_n)$ 的恆等式即得。這不是聲稱任何 Runge–Kutta 或任意時間步都在所有非線性熱流問題中守恆。 ## 8.2 含熱節點的同一時間階段 若取中點散熱及中點機械／聲學耗散，熱節點更新為： $$ C_T\frac{\vartheta_{n+1}-\vartheta_n}{\Delta t} =\beta w_{n+1/2}^2+v_{n+1/2}^\top R_av_{n+1/2} -h_T\vartheta_{n+1/2}+Q_{T,n+1/2}. $$ 因而： $$ E_{h,n+1}^\star-E_{h,n}^\star =\Delta t\left(F_{n+1/2}w_{n+1/2} +Q_{T,n+1/2}-h_T\vartheta_{n+1/2}\right). $$ 這個離散帳本是一項**檢查數值物理是否被悄悄扭曲**的具體手段。若僅輸出漂亮的聲壓波形與節拍器同步圖，卻沒有能量殘差，便無法判定數值方法是否把界面交換變成了假能量源。 ## 8.3 為什麼能量正確還不足以說模擬準確？ 隱式中點法可在二次線性能量模型保持嚴格守恆／耗散帳本，但時間步過大仍會造成相位誤差、波速色散、共振峰偏移及無法解析快事件。必須再測量： - 時間步細化下的相位、振幅、模式與功率收斂； - 空間網格細化下的聲學本徵頻率及壓力場收斂； - 快慢時間尺度的可解析性； - 接觸事件定位及壓力波反射； - 是否存在時間步改變就消失的假同步狀態。 **能量一致性是必要品質閘門，不是整套多場物理模擬的充分性定理。** ## 8.4 分割耦合與單體耦合不可混為一談 - **單體（monolithic）法：** 一個整體系統同時求機械與聲場未知量，易於維持一致介面，但矩陣大、預條件要求高。 - **強耦合分割（iterated partitioned）法：** 交替求流體與固體直到介面力、速度和能量殘差滿足閘門；仍須處理 added mass。 - **弱耦合分割（loosely coupled）法：** 先後各算一次即可，成本低但可能不穩定甚至有假能量，不能對高耦合區域直接推廣。 任何方法都須提供同樣的輸入、界面條件、離散誤差規格與比較成本。只因「全域計算」在概念上採耦合模型，不意味某實作必然優於成熟 monolithic FSI 或多物理求解器。 --- # 9. 解析案例：為什麼空氣、聲音與幾何都不能僅按「力大小」排序？ ## 9.1 C-06-01：活塞固定時的聲腔正常模態 先將活塞固定，令 $\dot q=0$，兩端都是剛性壁，取理想無耗散聲場。對任意正整數 $n$： $$ k_n=\frac{n\pi}{L},\qquad \omega_n=ck_n, $$ 可設： $$ p_n(x,t)=P_n\cos(k_nx)\cos(\omega_nt), $$ $$ v_n(x,t)=\frac{P_n}{\rho_0c}\sin(k_nx)\sin(\omega_nt). $$ 代入 $\partial_tp+K_a\partial_xv=0$ 、 $\rho_0\partial_tv+\partial_xp=0$，恰好抵銷；又因 $\sin(k_nL)=0$，兩端速度為零。因此理想封閉聲腔的非零聲學本徵頻率為： $$ f_n=\frac{nc}{2L}. $$ **M-06-10。** 這是**兩端固定壁、常係數無耗散聲場**的結果；若一端允許活塞移動並與彈簧／底板耦合，原本的獨立聲腔頻率通常改變，不能將 $f_n$ 當成耦合系統的精確本徵頻率。 此例強調：聲學場即使沒有讓活塞當下移動，也可能含非零的壓縮／流體動能；同一瞬間觀察到節拍器的相位，無法讀出全部的聲學正常模態。 ## 9.2 C-06-02：低頻聲腔在活塞上表現為附加彈簧 在 $\omega L/c\ll1$ 且無顯著熱交換、活塞小振幅、近乎均勻等熵壓縮的條件下，聲腔容積為 $V=AL$，活塞往內移動 $q$ 造成一階容積縮小 $\Delta V\approx -Aq$。因此： $$ p\approx K_a\frac{Aq}{AL}=\frac{K_a}{L}q. $$ 活塞受到空氣恢復力 $-Ap$，相當於增加： $$ k_{\mathrm{gas}}=\frac{AK_a}{L}. $$ 故低頻近似的機械有效剛度為 $k_{\mathrm{eff}}=k+k_{\mathrm{gas}}$，不是因為「空氣力被吸收」而消失。 **重要限制：** 若空氣與壁面進行顯著熱交換，壓縮接近等溫，應使用相應體積模量，不能仍用絕熱聲速所定義的 $K_a$。若活塞／底座的幾何位置或作用有效面積改變， $k_{\mathrm{gas}}$ 也可能改變。 ## 9.3 C-06-03：同樣的活塞相位，隱藏聲壓使瞬時加速度不同 在第二節的線性模型，若兩個世界在 $t=0$ 都有 $q(0)=0$ 、 $\dot q(0)=0$，但世界 A 的初始聲場全零，世界 B 取： $$ p_B(x,0)=P\cos(\pi x/L),\qquad v_B(x,0)=0, $$ 則初始聲學速度均符合 $v_B(0)=\dot q=0$ 與 $v_B(L)=0$。在光滑可代入的初始資料下： $$ \ddot q_B(0)-\ddot q_A(0)=-\frac{AP}{m}. $$ 這是嚴格的**機械觀察投影不足反例**。它既不需新基本力，也不需要量子效應；原因就是被投影省略的聲壓已經是完整狀態的一部分。 但若實驗對兩組世界增加了聲學感測器，或用足夠長的延遲觀測重建內部聲場，上述不可區分性可能被打破。故「相位不足」必須相對於指定觀察器／歷史長度／干預域成立，不能被升級為對一切相位法的全稱否定。 ## 9.4 C-06-04：聲學輸入阻抗不是純粹耗散 考慮與 2 節同一理想聲腔，但以固定角頻率 $\omega$ 的複數振幅求解。採 $e^{i\omega t}$ 時間慣例， $k_a=\omega/c$，閉端 $v(L)=0$。由一維聲場可得活塞端的壓力—速度比： $$ Z_{\mathrm{in}}(\omega)=\frac{\widehat p(0)}{\widehat v(0)} =-i\rho_0c\cot(k_aL). $$ 此處 $Z_{\mathrm{in}}$ 單位為 $\mathrm{Pa\,s/m}$，若要換成機械力—速度阻抗，需乘以面積 $A$ 並留意作用反向。理想模型中 $Z_{\mathrm{in}}$ 是純虛數；因此聲腔在此近似下主要形成可逆的能量儲存／反應性載荷，而非平均耗散。引入熱黏性或開放輻射時才出現相應的正實耗散部分（符號取決於端口方向慣例）。 當 $\omega L/c$ 很小， $\cot(k_aL)\approx1/(k_aL)$，將 $\widehat v=i\omega\widehat q$ 代入，得到 9.2 的 $\widehat p\approx (K_a/L)\widehat q$。這展示**同一空氣介質隨頻率與邊界改變，可以呈現恢復剛度、聲學記憶或共振響應**，不能用單一「空氣勢能較弱」判斷其長期地位。 在極理想閉合腔的共振條件附近，阻抗函數可能出現奇異行為；真實有限阻尼、活塞耦合、非線性與裝置有限剛度會改變此結果。不可把理想奇點直接當成可測的無限聲力。 ## 9.5 C-06-05：很強的恆定載荷也可能零平均輸入功率 設集中機械系統承受恆定外力 $F_0$，在阻尼下趨向靜態平衡 $q_*=F_0/k$ （此例暫不耦合聲腔）。於平衡態 $\dot q=0$，所以外力瞬時功率 $F_0\dot q=0$，即使 $|F_0|$ 很大。與此對照，弱而頻率適當的週期力能持續對振動模式做平均正功；其可否被放大取決於相位、共振品質因子、阻尼和輸入頻譜。 這證偽「持續較強的牛頓力必然壓倒較弱作用」的無條件版本，但並不否定在**耦合條件與耗散特性已指定**後，某些持續驅動能夠主導長期模式。 ## 9.6 C-06-06：界面方法可以讓數值世界和物理世界給出相反穩定性 對同一個聲—固邊界問題，若在離散實作中把第 7 節的 $-b^\top p$ 誤寫成 $+b^\top p$，而聲學端仍保留 $+bw$，則總機械—聲學能量導數含有： $$ \dot E_h=Fw-\beta w^2-v^\top R_av+2wb^\top p. $$ 這個沒有一般固定符號的額外項完全來自**錯誤的界面代數**，而非新物理。當 $wb^\top p>0$ 時它偽造額外輸入，可能製造假振幅成長或假同步；當小於零則產生假耗散。因此 GMDIS 對「世界線」的預測，必須先通過界面功率檢驗，才有資格討論規則勢。 --- # 10. 節拍器不是活塞：從最低基準走向擒縱、碰撞與非線性 ## 10.1 D-06-17：多節拍器擺角與共享底座 令第 $i$ 個節拍器的擺角為 $\theta_i$ 、角速度 $\omega_i$ 、內部擒縱狀態 $\xi_i$ ；共同底座廣義坐標為 $Q$，可同時含平移與旋轉。其真實幾何位置 $r_i$ 進入關節動能、質心運動、接觸位置及底座力矩。理想一般形式可寫： $$ M(z)\ddot z+C(z,\dot z)\dot z+\nabla V(z) =F_{\mathrm{esc}}(z,\dot z,\xi)+F_{\mathrm{contact}}+F_{\mathrm{fluid}}+F_{\mathrm{ext}}, $$ 其中 $z=(Q,\theta_1,\ldots,\theta_N)$ ；此式是系統形式記法，不保證 $C$ 只有一個唯一矩陣定義，也不自動包含所有接觸非光滑項。每一個 $F$ 必須標明是真正的廣義力、約束乘子或外部功率埠。與空氣耦合時，必須把 $F_{\mathrm{fluid}}$ 連接到流場界面牽引，不能將它獨立生成再在流場帳本重算。 ## 10.2 擒縱做功的能源從哪裡來？ 真正會長期規律擺動的機械節拍器需要擒縱機構持續從發條、配重或其他能源釋放能量。把擒縱當作一個純粹「穩定力」而不計能源，會錯把非平衡定態寫成永動系統。 對每個 $i$，可定義擒縱輸入功率 $P_{\mathrm{esc},i}=\tau_{\mathrm{esc},i}\omega_i$。若能源儲量為 $E_{\mathrm{bat},i}$，無補充時需有： $$ \dot E_{\mathrm{bat},i}+P_{\mathrm{esc},i}+D_{\mathrm{drive},i}=0, $$ 在所選功率方向與模型假設下 $D_{\mathrm{drive},i}\ge0$。若擒縱器內部彈簧、撞擊及摩擦有可觀儲能，還需分別列出；不可以把它們全部藏入一個「非線性力」後仍聲稱總機械能守恆。 ## 10.3 非光滑事件與重設映射 碰撞、齒輪擒縱的卡合／釋放及單邊接觸可以造成速度瞬時跳變。令連續演化段標記為 $\xi$，事件面為 $h_\xi(z,\dot z)=0$，重設映射為： $$ (z^+,\dot z^+,\xi^+)=\mathcal R_\xi(z^-,\dot z^-,\xi^-). $$ 例如理想無外部供能的非彈性碰撞可以令機械動能降低，但降低部分要進入熱、聲、材料變形或未建模的耗散帳本。若撞擊由擒縱儲能驅動而增加機械能，則必須從能源儲量扣除對應做功。未交代這些條件，不應把光滑 $C_0$ 半群直接套入跨事件跳變。 ## 10.4 為什麼空氣中的聲波與底座振動需要分別測量？ 節拍器可能同時透過底座振動與空氣聲學作用於其他單位，但兩種通道在傳播速度、阻抗、衰減、頻率範圍、耦合位置及外界邊界上不同。已知有**聲學耦合電子節拍器**的實驗例子，不能從「一般機械節拍器的共同底座通常最重要」推論聲音永遠可略；反過來也不能因聲波可同步某種電子節拍器而認定聲學主導一般機械裝置。 必須設計**介質通道消融**：固定底板、設聲隔離、改變表面吸音／反射、固定溫度但調整空氣密度等，並以相同擒縱機構與初態比較。實驗順序需防止改變空氣壓力同時改變擺臂阻力，從而把不同物理因果通道混成一個效果。 ## 10.5 長時間穩定運動是跨能量埠平衡，不是張力全部消除 若進入統計穩定的週期或有界非平衡定態，可以在足夠長的觀察窗內考察平均功率帳本。例如單個耦合子系統在週期穩態： $$ \langle P_{\mathrm{esc}}+P_{\mathrm{ext}}\rangle =\langle D_{\mathrm{mechanical}}+D_{\mathrm{acoustic}}+P_{\mathrm{radiated}}\rangle, $$ 前提是該時窗中**總儲能平均變化為零**，並且各項通道沒有重複計算。此式不表示每個時刻各項功率都相等，也不表示不同節拍器一定同相，更不表示宇宙規則勢數值等於右側瓦特數。 --- # 11. 多尺度與無因次控制：決定何時需要更精細的世界 ## 11.1 D-06-18：最低時間與空間尺度 對特徵長度 $L_*$ 、典型速度 $U_*$ 、角頻率 $\omega_*$ 、聲速 $c$ 、運動黏度 $\nu$ 、熱擴散率 $\alpha_T$，以下組合可作為模型選型的**提示**，但不是萬能閾值： | 無因次群 | 定義 | 物理意義與警戒 | |---|---|---| | 聲學尺度 $k_aL_*$ | $\omega_*L_*/c$ | 波長與裝置大小；大時可能有多個聲學模態 | | Mach | $U_*/c$ | 壓縮性與小訊號假設的條件之一；小 Mach 不自動消除聲學共振 | | Reynolds | $\rho_0U_*L_*/\mu$ | 黏性與慣性比值；需指定速度與特徵長度 | | 黏性邊界層比 | $\delta_v/L_*$ | 近壁面尺度與空間解析度 | | 熱邊界層比 | $\delta_T/L_*$ | 熱交換與聲學耦合的近壁面尺度 | | Strouhal | $\omega_*L_*/U_*$ | 非定常相對於對流；若 $U_*=0$ 該定義不適用 | | 流／固質量比 | $M_{f,\mathrm{eff}}/M_s$ | 對附加質量、FSI 分割收斂的提示 | | Péclet | $\rho_0c_pU_*L_*/\kappa_T$ | 熱對流與熱擴散相對尺度 | 若 $\rho_0,c,\mu,\kappa_T$ 隨溫度變動，這些群也可能隨時間與位置改變；不能只用一次初始值作為全部世界線上的固定判斷。 ## 11.2 不把細尺度直接等同更好 解析每一個聲學黏性邊界層、材料晶格與空氣分子，計算成本可以極高，且未必提高指定任務的預測力。GMDIS-02 的「任務充分性」在此具有直接工程作用： 若只需預測 10 秒後是否仍鎖相，適度阻抗等效或模型降階可能已足夠；若需對改變封閉聲腔幾何的反事實做精準預測，原本的等效阻尼模型就可能不足。 **更精細不是學術上的自動勝利，充分且可重建才是研究目標。** ## 11.3 介面相容與模型粗粒化 從細尺度流體／固體模型萃取粗粒化聲學阻抗時，至少指定： - 空間平均、頻率帶與分析時間窗； - 哪些守恆量與可觀測量必須保留； - 耗散是否變成顯式阻尼或記憶核； - 不同模型是否使用相同外部功率埠； - 哪些外部干預要求重新估計縮減模型。 若 $M1\to M2$ 的縮減對指定輸出充分，不需要因為理論追求「萬有張力」而強迫啟用 M3。反之，M1 對熱邊界干預若失效，應承認缺失並提升模型。 --- # 12. GRST、SIO 與 24／72：多介質世界的語義接縫 ## 12.1 SIO：介質的佔據與自由度，不是聲學能量 在固定規則 $R$ 、物理容器 $\mathfrak C$ 下，可行物理配置為： $$ \Omega_{R,\mathfrak C}\subseteq\mathcal X, $$ 其中 $\mathcal X$ 包含各種場與集中坐標的乘積，但合法配置需滿足界面約束、有限能量、材料本構及守恆條件。SIO 關心某個物體、隔板、約束、聲學邊界或接觸狀態，如何改變整個合法配置與可行路徑的集合。 **不能**由「一片吸音材料的幾何支撐很小」直接推出「聲學能量被無限遠地瞬時消耗」；前者是容器約束資訊，後者需由明確 PDE、材料阻抗與有限傳播時間決定。 ## 12.2 兩種陰影在本篇的各自地位 **硬約束陰影：** 插入完全封閉隔板、改變邊界連通性或鎖定底板轉動自由度，可以改變合法配置集合，甚至在固定容器的路徑定義下立即改變全域可達關係。 **動態影響陰影：** 對單一節拍器施加微小脈衝後，在指定物理模型與時間窗中，其他單位的壓力、速度、恢復時間或音波響應出現變化。這須以脈衝響應、傳遞函數或受控敏感度衡量，且遵守物理傳播規則。 對一個觀察映射 $Y_i(U)$ 及局部干預 $\delta_j$，可定義有限時間敏感度（候選）為： $$ S_{ij}(T)=\sup_{0\le t\le T}\frac{\|Y_i(U_t^{\delta_j})-Y_i(U_t^0)\|}{\|\delta_j\|}, $$ 其中分母、範數與干預族需事先固定。 $S_{ij}(T)$ 在聲學範圍不代表超距作用；如果缺乏適當有界性，極限敏感度也未必存在。數值測量時要處理時窗、底層噪聲及解析度。 ## 12.3 GRST：勢是未來結構剖面，不是聲壓積分 本篇使用既有 GRST-12 的定義：全域規則勢是規則對可行未來空間誘導出的結構性非均勻。對已具有物理方程的世界，可以萃取候選剖面： $$ \Phi_R(U_0;T,\ell)= (\mathcal R_T,\,\mathcal A_T,\,\mathcal D_T,\,\mathcal C_T,\,\mathcal S_T), $$ 其中在先行定義之下，分別代表可達性剖面、擾動放大、累積耗散、跨介質能量交換與模式存續／恢復性。各分量可以是函數、向量或集合，不要求相同量綱或自動可合併成單一 scalar。 **沒有任何一項 $\Phi_R$ 可直接宣稱是新自然力**。若所提剖面能被 port-Hamiltonian 能量帳本、傳遞函數、可達 Gramian、Lyapunov／Floquet 特徵或成熟 FSI 模型完全取代，則應明記其詮釋或整合價值，而不是臨時創造不可證明的獨立物理實體。 ## 12.4 24／72 的原始四軸不改動 既有計算分類的四個軸：底空間 B（連續／離散）、更新組織 U（序列／選擇／並行／識別）、觀察 O（連續／離散／拒單測）、最低充分轉移律 L（確定／古典隨機／量子通道），仍是**計算與表示事件的可尋址分類**，而不是四個自然基本力學軸。 | 模型／數值事件 | 候選格 | 判定邊界 | |---|---|---| | 線性聲場的有限體積並行更新、連續聲壓輸出 | C-P-C-F | 連續物理近似在有限精度網格上求解；分類相對於選用表示 | | 擒縱事件檢測與離散狀態切換 | C-J-D-F | 若事件律確定且在指定連續狀態上檢測 | | 隨機熱擾動的古典場數值演化 | C-P-C-K | 僅當隨機核是指定最低充分轉移律 | | 全場觀測多指標共同重建失敗時 | C-P-X-F | X 必須先指定觀測族與容許表示類並提出不足證據 | | 模型辨識器依既有索引選擇求解器 | D-R-D-F | R 是計算／分析器的識別讀出，不是宣稱物理粒子會檢索 | **Q 類僅在必須保留量子相干或非交換量子通道時選用。** 由於聲波存在相位或擺鐘「相干」，並不能因此把古典節拍器標為量子 Q。對大型多層系統，使用格點路徑／超圖與轉移契約，不把整個世界永久塞入唯一一格。 ## 12.5 新增語義不是事後命名 本篇允許兩種分類導向的實際工程假說：(i) 24／72 標註可降低跨模型接口錯誤，(ii) 可提升適當模型選擇、粗粒化與重建效率。這些都屬 H-06，必須設置同樣規模但沒有 24／72 標註的成熟工作流當基線；「分類名字很完整」不是成功指標。 --- # 13. 從模型到實驗：多介質世界的反事實矩陣 ## 13.1 D-06-19：固定比較任務 一切實驗前先聲明「模型到底要預測什麼」。本系列建議至少區分五種任務： 1. **局部軌跡任務：** 預測指定節拍器的角位置、速度、擒縱事件時間與後續相位。相位誤差可用週期座標評估，但不能取代機械能、功率及場狀態。 2. **集體模式任務：** 區分同相、反相、群集、漂移、相位滑移、停止及多穩態；分類閾值須在觀察資料之前定義。 3. **穩定性任務：** 對已形成的運動施加標準化擾動，測量是否恢復、恢復時間、可能的轉態率與放大程度。 4. **介質干預任務：** 改變聲學壁面、流體性質或結構耦合後，預測其對集體行為的方向與大小。 5. **跨尺度重建任務：** 只保留某些觀測後，檢查能否恢復必要的場變數、能量帳本或後續干預反應。 以上不能混成一個「同步準確率」。對部分任務，相位—振幅、輸入輸出傳遞函數或既有狀態估測器可能完全足夠；此時 GMDIS 不能為了展示新框架而增加無必要的模型複雜度。 ## 13.2 D-06-20：四類干預保持分離 依 GRST-07 的雙層反事實，分成：(a) 初始狀態干預；(b) 模型參數與規則干預；(c) 幾何／邊界干預；(d) 觀測器干預。研究表須確實說明哪些物理條件不變。 | 干預標籤 | 具體操作 | 主要作用路徑 | 需固定的干擾因素 | |---|---|---|---| | I0 | 只改一個擺角／擺速 | 狀態干預 | 全部幾何、底座及空氣參數 | | I1 | 改底座整體質量但不改剛度 | 慣性路徑 | 幾何及阻尼改變須量測 | | I2 | 改底座可轉動／可形變自由度 | 結構模式 | 不可同時改節拍器數量與能源 | | I3 | 改節拍器空間排列 | 幾何力矩及模態耦合 | 必須先確認模型具有位置敏感自由度 | | I4 | 改聲學牆面阻抗 | 反射／吸收／聲功率 | 底板機械支撐與熱邊界 | | I5 | 改流體密度／壓力 | 流體慣性、聲速、阻尼 | 擒縱摩擦與發條輸出可能受連帶影響 | | I6 | 改溫度或熱交換條件 | 材料、黏度、聲速、熱傳 | 除熱效應外的擒縱或機械參數漂移 | | I7 | 改觀測裝置或頻率帶 | 觀測投影 | 不應把感測器改變誤稱規則改變 | | I8 | 改模型求解器／離散網格 | 計算表示 | 真實物理世界與外界條件保持固定 | 對 I5、I6 等多重影響干預，理想上應透過校準與附加控制實驗拆解因果路徑；若無法拆解，應回報「複合干預」而非單獨歸因於聲學或熱學。 ## 13.3 幾何零效應必須成為預先接受的負結果 設只有一維剛體平移自由度的共同底座，所有節拍器方向、質量及同底座的作用規則一致。只調換物體的平面排列但不改有效接點、方向與系統慣性，模型的集體動力可能完全不變。這種「幾何標籤更改而動態不變」是合理的**零效應控制組**。 若轉成可旋轉／可彎曲底座，同一排列改動可能影響力矩與模態。二者的差異必須能從物理自由度的增刪及施力空間分布推出，不得僅以漂亮的網格圖像當作證據。 ## 13.4 聲學介面消融與真正的音波耦合 一組最小比較可以是： - A：共享可動底座、正常空氣聲場； - B：相同共享底座，對聲場施加高吸收／減小反射的候選改變； - C：底座機械耦合降低，仍允許聲學耦合； - D：同時降低底座與聲學通道； - E：在指定條件下開放外部聲學驅動； - F：不改物理世界，僅從模擬刪去聲學模組作模型消融。 A–E 是實物干預設計候選，F 是數值／模型消融，**不可把 F 與真正移除物理聲場混為一談**。真實吸音罩會改變容積、空氣流通與可能的溫度分布；同時需要量測其副作用。 ## 13.5 外部實物並非本篇已完成 本篇沒有製造或測量任何 100 台、1000 台節拍器矩陣。將理論上的「可達性變化」「聲學影響範圍」「能量重新分配」直接當作實測結論，屬不合格證據。本篇實際完成的是線性化模型推導及附錄的有限維數值一致性核對。 --- # 14. 基準階梯與實驗可反駁性 ## 14.1 D-06-21：由 E0 至 E8 的研究階梯 | 基準 | 系統與比較 | 最低可驗證指標 | 成功或失敗的邊界 | |---|---|---|---| | E0 | 獨立振子；無底板、無聲學耦合 | 已知解、基本能量帳本 | 任何無來源的同步都視作模型錯誤 | | E1 | GMDIS-05 弦—集中振子 | 弦／振子能量、半群／Galerkin 參照 | 不得偷換成本篇新聲固定理 | | E2 | 本篇 1D 聲腔—活塞，理想無損 | 聲壓波速、介面功率、總能量 | 零耗散不能憑空失去能量 | | E3 | E2＋有效阻尼與集總熱 | 非負耗散、熱上升、散熱功率 | 損失必能找到模型內／外去向 | | E4 | 剛體底板＋多節拍器 | 質量矩陣與幾何零效應 | 只平移底座不應憑空產生位置力矩 | | E5 | 平面可旋轉、可彎曲底座 | 幾何模態與受控擾動 | 需對照成熟多體／FSI 基線 | | E6 | 聲固耦合＋可調聲學邊界 | 聲場與結構場的能量傳遞 | 必須避免假波長／假共振峰 | | E7 | 熱黏性與擒縱非光滑模型 | 熱流／熵、衝量、能源帳本 | 不將每次碰撞視為可微連續流 | | E8 | 2、10、100、1000 單位矩陣 | 預測、恢復、干預、成本與誤差 | 系統規模大不等於新理論勝利 | 由 E2/E3 到 E6/E7，**模型複雜度及未辨識參數也會同步增加**。若額外自由度不能改善指定任務或建立獨立的干預預測，應選擇更簡潔的有效模型，不必強求滿級。 ## 14.2 必須與哪些成熟方法競爭？ 1. **完整牛頓／Lagrange 多體模型：** 限於實驗所需的剛體與可形變自由度，不以相位模型作唯一反方。 2. **聲學傳遞函數、阻抗與模態模型：** 包括一維閉端聲腔、聲學有限元及頻率依賴邊界。 3. **熱黏性聲學與經典 FSI：** 具實際界面條件、阻抗、ALE／immersed boundary 或等效成熟方法。 4. **相位—振幅、Floquet、Lyapunov 與 Koopman 等長期動力工具：** 依待測任務採用，避免把「傳統方法」錯當只有 Kuramoto 一種。 5. **Port-Hamiltonian／GENERIC：** 作為能量交換、耗散和熱力學一致性的強基線，尤其重要。 6. **有記憶的縮減模型：** 對投影掉的聲學／結構自由度，先允許 Volterra、Mori–Zwanzig 或識別傳遞核等成熟方法補足。 若任何成熟方法在相同資料預算、物理邊界與驗證任務下完全達到要求，GMDIS 應紀錄其成功，而不是宣稱新理論自動更高階。 ## 14.3 H-06-01：跨介質模型的任務增益（待驗） **假說：** 對某些高反射聲腔、低阻尼底板或具有聲學結構共振的配置，顯式聲場／結構耦合比「僅振子＋常數等效阻尼」更能預測幾何、聲學邊界與短時擾動干預下的模式轉換。 **主要判敗條件：** 若在預先保留的參數組、初態與實驗干預上，受控基線持續達到相同或更好之預測誤差、校準與資料效率，而顯式聲場沒有額外外部泛化價值，則不得將此 H-06 標成成功。 ## 14.4 H-06-02：持續作用的條件性主導（待驗） **假說：** 若某模式與持續擒縱輸入有效耦合，且其在指定時窗中獲得的淨平均功率與擾動恢復機制足以彌補耗散，該模式較可能長期存續；但小幅弱共振通道可使穩定模式換位。 **判敗條件：** 若加入平均功率、頻譜匹配、損失、回饋與傳遞延遲等指標後，對模式存續沒有超出既有模型的解釋或預測增益，則此組「張力」描述只保留為整理語言。 ## 14.5 H-06-03：SIO 動態陰影是否有額外可用資訊（待驗） **假說：** 受控局部干預所得的跨介質敏感度剖面，在指定任務中可能比單一相位同步指標更能預測受擾恢復與分岔。 **判敗條件：** 若完整傳遞函數、Green's function、可觀測 Gramian 或成熟 FSI 量已完全吸收該剖面的相關資訊，需承認它不是新的獨立物理不變量；可保留為模型介面或視覺化資料結構。 ## 14.6 H-06-04：24／72 的模型配置價值（待驗） **假說：** 對需要在 ODE、PDE、事件、隨機及觀測模型間切換的長管線，固定四軸標註與介面契約可能降低模型型別錯誤，改善端到端成本或資料追溯性。 **判敗條件：** 若與同等規模的型別系統／多物理組件圖相比，沒有提高可重現性或降低錯誤，此標註即不獲得額外效能主張。 ## 14.7 不能事後改閾值的評估協議 對每個實驗先封存：時間窗、收斂容差、感測器頻寬、訓練／驗證切分、模型校準參數集合、干預可行域及負面結果保存策略；所有看見結果後新增的閾值、選擇或統計量標為**探索性**。多次比較若存在多重測試問題，須實施適當校正或明確限制結論。 --- # 15. 統一驗證矩陣：定律、演化、觀測與實作要分開報告 ## 15.1 D-06-22：不能只輸出最終同步影片 最低驗證包需要同時有： | 類型 | 必要輸出 | 通常失敗徵象 | |---|---|---| | 物理型別 | SI 單位、材料常數、方向約定 | 牛頓與焦耳直接相加、密度與資訊量混淆 | | 守恆 | 總能量、機械能、聲能、熱能、外界埠 | 界面交換被算兩次、漏掉聲輻射 | | 熱力 | 熱功率去向、正溫、熵生成 | 模擬冷卻卻找不到熱邊界 | | 數值 | 空間／時間收斂、穩定與能量殘差 | 時間步改變後同步模式消失 | | 觀測 | 採樣頻率、噪聲、解析度、缺失資料 | 觀察精度不足誤判鎖相 | | 干預 | 狀態／規則／邊界／觀察層分離 | 多個物理條件一起變動造成偽因果 | | 機制 | 聲阻抗、結構模態、延遲、能源補給 | 只能回報相位，不能說明作用來源 | | 比較 | 強基線、成本、外部泛化、負結果 | 模型越複雜自動被視為越正確 | ## 15.2 能量殘差不能只有一個無量綱數字 對時間步長 $\Delta t$，可定義離散機械聲學殘差： $$ \mathcal R_E^n=\left(E_{h,n+1}-E_{h,n}\right) -\Delta t\left(F_{n+1/2}w_{n+1/2} -\beta w_{n+1/2}^2-v_{n+1/2}^\top R_av_{n+1/2}\right). $$ 而含熱之總殘差可定義： $$ \mathcal R_{E^\star}^n= E_{h,n+1}^\star-E_{h,n}^\star -\Delta t\left(F_{n+1/2}w_{n+1/2}+Q_{T,n+1/2}-h_T\vartheta_{n+1/2}\right). $$ 兩者都具有 J 單位；若要比較不同能量規模，可使用指定非零基準尺度 $E_\mathrm{ref}$ 正規化，但仍需同時保留絕對殘差，避免分母很小時產生誤導。對全流體 moving-FSI，必須重新推導相應的流體能量、移動邊界功及熱流，**不能**直接把本式的數值零殘差當成真實多場模型的證明。 ## 15.3 驗證與驗證對象的分離 - **Code verification：** 程式有沒有正確求解**寫下的方程**？靠能量恆等式、解析解、manufactured solutions、網格與時間步收斂。 - **Model validation：** 所寫的方程與參數對**外部實際裝置**是否足夠準確？靠實物量測與未知條件泛化。 - **Conceptual claim validation：** GRST／SIO 是否在已有物理指標之外提供可重複的新增結構／預測能力？靠預註冊消融與競爭性基線。 這三個問題不能用同一個 PASS 蓋章。 --- # 16. 本篇完成的數值一致性實驗 N-06-01 ## 16.1 已執行與不曾執行的範圍 本篇附錄 C 提供 NumPy 可執行程式；它以第 7 節 $n=12$ 個聲壓元件、 $n-1$ 個面速度元件，耦合一個集中機械振子，使用隱式中點法進行 3000 個時間步，步長 $10^{-5}$ s。程式同時計算機械—聲學能量帳本及熱節點帳本。 在本稿撰寫時已實際執行的程式環境，是容器中的 Python／NumPy 浮點數值計算；**沒有**直接測量任何真實空氣聲腔、實體節拍器、三維固體或熱電偶資料，也沒有對整個非線性流體—固體方程作物理數值收斂證明。 ## 16.2 N-06-01：測試參數 | 參數 | 指定值 | 單位／備註 | |---|---:|---| | $n$ | 12 | 聲壓單元數 | | $\rho_0$ | 1.20 | kg/m³ | | $c$ | 343.0 | m/s | | $A$ | 0.004 | m² | | $L$ | 0.25 | m | | $m$ | 0.35 | kg | | $k$ | 13.0 | N/m | | $\beta$ | 0 或 0.04 | kg/s；無耗散與耗散兩組 | | $R_a$ | 零矩陣或對角線 0.00012 | 有效聲學速度阻尼；矩陣元素單位 kg/s | | $C_T$ | 8.0 | J/K | | $h_T$ | 0.06 | W/K | | $F,Q_T$ | 0 | 無外部施力／供熱 | | $q(0)$ | 0.002 | m；其餘機械聲學初態設為零 | | $\Delta t$ | $10^{-5}$ | s | | 步數 | 3000 | 總時窗 0.03 s | 初始 $q$ 位移不是宣稱活塞真的被物理準靜態地壓進封閉氣體卻毫無聲壓；它是一組用於線性初值問題的數學指定資料。真實實驗如何準備初態必須另外建模與測量。 ## 16.3 N-06-01：直接重算結果 | 指標 | 無耗散 ( $\beta=0,R_a=0$ ) | 帶耗散 ( $\beta=0.04,R_a\neq0$ ) | |---|---:|---:| | 初始聲學＋機械能 (J) | $2.6000000000000002\times10^{-5}$ | 同左 | | 末態聲學＋機械能 (J) | $2.5999999999972053\times10^{-5}$ | $2.599938054420463\times10^{-5}$ | | 最大單步機械／聲學帳本殘差 (J) | $4.743384504624082\times10^{-20}$ | $4.1521853108832586\times10^{-20}$ | | 最大單步含熱帳本殘差 (J) | $4.743384504624082\times10^{-20}$ | $3.914366226143245\times10^{-20}$ | | 熱節點末態相對升溫 (K) | 0 | $7.742528925811134\times10^{-11}$ | 以上是本篇**指定離散 ODE 的已執行浮點一致性**。殘差達到極小量，是因時間法對此二次能量的代數結構相容；它不證明大型真實物理模型沒有模型誤差、界面幾何誤差或材料誤差。若更換 BLAS、CPU、NumPy 版本或運算精度，最後幾位浮點尾數可以不同。 ## 16.4 N-06-02：隱藏聲學壓力的瞬時干預 在相同 $q(0)=w(0)=0$ 的兩個世界，僅令其中一組第一個聲壓單元多 $100$ Pa。按第 7 節機械式，瞬時加速度差為： $$ \Delta\ddot q=-\frac{A\cdot100}{m} =-\frac{0.004\cdot100}{0.35} \approx-1.142857142857143\ \mathrm{m/s^2}. $$ 此結果由解析式直接重算。它展示：當相位、機械位置與速度相同，**未觀察的介質狀態可以立即影響被觀測的局部加速度**。若未來研究發現傳統全狀態模型已能完整預測此現象，應將之認定為成熟耦合力學的成果，而不是新物理力的發現。 --- # 17. 學術邊界、判敗義務與向後續系列的交付 ## 17.1 不宣稱的成果 1. 本篇**未**證明空氣或聲音的能量永遠低於共同底板作用，亦未證明持續強力必然導向同步。 2. 本篇**未**建立一種普適的「萬有張力守恆律」；守恆與耗散必須由明確物理系統與封閉邊界定義。 3. 本篇**未**證明存在獨立於成熟多物理學的「規則勢物理場」。 4. 本篇**未**證明熱—黏—彈—接觸的完整非線性系統對所有初態全球適定。 5. 本篇**未**證明所有多場效應都無法被相位—振幅、延遲觀測、模態與記憶核模型準確重建。 6. 本篇**未**把有限聲場空間格點或 SIO 構成纖維宣稱為宇宙基本構成物。 7. 本篇**未**將 72 格視為 72 種自然基本力量。 8. 本篇**未**用數值能量殘差替代模型對真實世界的有效性驗證。 ## 17.2 四層研究真值必須分離 - **數學：** 2、3、4、6、7、8、9 節的限定模型推導若可重新演算，即為該模型條件下的可證結果。 - **數值：** 附錄 C 的代碼在已指定參數、時間法與算術下通過能量等式測試，屬計算實例。 - **自然科學：** 與真實物件、流場、聲壓、熱傳、可轉動底板及擒縱的對應，目前仍待實驗量測。 - **新穎性：** GRST/SIO/24-72 對物理預測或工程流程有多少增益，必須經盲測、基線競爭與失敗保存才可評價。 ## 17.3 GMDIS-07 的直接輸入 下一篇〈持續作用與動態穩定結構〉可以使用本篇的機械、聲學、熱學及能源埠；但必須另外建立穩定集合、擾動恢復、週期解、Floquet 譜、非光滑事件的回歸映射與多穩態分岔條件。**「有能量帳本」不等於「有吸引子」、「有週期解」或「必然同步」。** ## 17.4 GMDIS-08 至 GMDIS-14 的接口 | 後篇 | 本篇交付內容 | 尚需補證／驗證 | |---|---|---| | 08 跨尺度投影、記憶與自由度傳輸 | 壓力／固體場的隱藏狀態與阻抗例子 | 消去場後的精確／近似記憶核與誤差界 | | 09 24／72 物理語義映射 | 模型級別與計算事件的標註表 | 物理演化與識別計算的型別分離 | | 10 多不變量剖面 | 能量、傳播、恢復與觀測量候選 | 最小充分觀察族及 X 判定 | | 11 規則勢與反事實 | I0–I8 干預、界面／規則類型 | 識別與不可識別條件、態勢換位定量 | | 12 節拍器矩陣基準 | E0–E8 分階模型與物理接口 | 2、10、100、1000 單位可重現模擬／實測 | | 13 基線競爭與消融 | 全部強基線與判敗條款 | 預註冊外部資料、盲測、成本與負面結果 | | 14 理論統合 | 不可混淆的物理量與可操作結果 | 僅保留受證據支持的主張；必要時降維 | --- # 18. 結論：集體秩序必須能通過真正的跨介質物理帳本 從 100 到 1000 個節拍器的同步圖景，很容易把研究注意力集中在相位、平均序參量或長期吸引子上。GMDIS 的起點不同：我們希望知道整個世界在多種介質及不同尺度之間，究竟發生了哪些因果上合法的作用，使某些集體運動可以存續，另一些運動衰減、轉換、分岔或根本不出現。 本篇建立的最低物理答案，並不是宣稱存在新的「萬有勢能」：它是一條可以核對的**界面功率交換律與開放系統第一定律**。空氣壓力可以壓縮儲能、聲波可傳播與反射、底座可以彎曲及傳遞機械動量、流體黏性與材料摩擦可將機械擾動能轉成熱、擒縱可以持續注入已記帳的能源。這些作用強弱不同、時間尺度不同，而且可能在不同物理條件下交換主導地位；只以瞬時力值、聲學能量大小或某個同步指標排序，通常不能完成全部動力學解釋。 **真正的研究命題是：** 若以 GRST 討論規則誘導的未來非均勻性，以 SIO 討論容器與剩餘自由度，以 24／72 區分模型表示與轉移律，那麼每一次跨學科映射都應回到可驗證的多介質介面、可追蹤的功率路徑、能量／熵條件、外部反事實與失敗邊界。這不是否定既有相位、吸引子、FSI 或熱力學，而是迫使研究者回答：**我們究竟省略了哪一層世界，省略後對什麼任務仍充分，又在哪一種干預下不再充分？** 本篇提供了一個可以逐層擴張並接受失敗的物理基準。下一篇將在其上研究持續驅動、耗散與長期穩定性，但不預先把穩定與同步視為必然。 --- # 參考文獻與來源審計 ## A. 既有 GMDIS 系列原始檔（跨篇依賴） 1. **[GMDIS-01]** Neo.K × Aletheia（2026-10-16），《從同步現象到全域多尺度動力學：研究母論與可證偽邊界》，v0.1。相對路徑：[`GMDIS_01_From_Synchrony_to_Global_Multiscale_Dynamics_v0.1_2026-10-16.md`](./GMDIS_01_From_Synchrony_to_Global_Multiscale_Dynamics_v0.1_2026-10-16.md)。 2. **[GMDIS-02]** Neo.K × Aletheia（2026-10-16），《相位與吸引子的資訊充分性邊界：投影、反例與重建條件》，v0.1。相對路徑：[`GMDIS_02_Phase_Attractor_Informational_Sufficiency_Boundaries_v0.1_2026-10-16.md`](./GMDIS_02_Phase_Attractor_Informational_Sufficiency_Boundaries_v0.1_2026-10-16.md)。 3. **[GMDIS-03]** Neo.K × Aletheia（2026-10-16），《GRST × SIO：底空間、自由度陰影與未來勢賦形》，v0.1。相對路徑：[`GMDIS_03_GRST_SIO_Substrate_Freedom_Shadows_Future_Shi_Shaping_v0.1_2026-10-16.md`](./GMDIS_03_GRST_SIO_Substrate_Freedom_Shadows_Future_Shi_Shaping_v0.1_2026-10-16.md)。 4. **[GMDIS-04]** Neo.K × Aletheia（2026-10-16），《萬有張力的物理語義校準：從廣義影響到力、應力、功率與跨尺度能量帳本》，v0.1。相對路徑：[`GMDIS_04_Physical_Semantic_Calibration_of_OmniTension_v0.1_2026-10-16.md`](./GMDIS_04_Physical_Semantic_Calibration_of_OmniTension_v0.1_2026-10-16.md)。 5. **[GMDIS-05]** Neo.K × Aletheia（2026-10-16），《多場泛函分析底座：狀態積空間、耦合算子、弱解、適定性與能量估計》，v0.1。相對路徑：[`GMDIS_05_Multifield_Functional_Analytic_Foundation_v0.1_2026-10-16.md`](./GMDIS_05_Multifield_Functional_Analytic_Foundation_v0.1_2026-10-16.md)。 ## B. EveMissLab 既有理論來源（內部理論接縫，不代表自然科學結論已證實） - **[I-06-01]** Neo.K（2026-10-05），[《GRST-12：全域規則勢統合理論》](https://unboundedaxiom.org/p/lm-005022/)，v0.1。用於區分規則、可行自由、Shi 非均勻結構、世界線與觀察。 - **[I-06-02]** Neo.K（2026-10-04），[《GRST-03：自由度不是勢》](https://unboundedaxiom.org/p/lm-005013/)，v0.1。用於自由纖維及規則賦形接縫。 - **[I-06-03]** Neo.K × Aletheia（2026-10-02），[《空間資訊佔據、構成纖維與容器自由度：一個主客邊—ACTM 統合框架》](https://unboundedaxiom.org/p/lm-005071/)，v0.1。用於容器、約束與非局部自由度陰影的定義；**未**把資訊纖維當成物理細絲。 - **[I-06-04]** Neo.K × AI 協作，[《計算的二十四重範式》](https://unboundedaxiom.org/p/lm-001708/)，正式版候選稿 v4.0，2026-07-19。用於三軸的語境相對分類與拒單測定義。 - **[I-06-05]** Neo.K × Aletheia，[《從二十四重計算形態學到七十二格計算動力學》](https://unboundedaxiom.org/p/lm-001687/)，v0.1，2026-07-19。用於第四軸最低充分轉移律。 - **[I-06-06]** Neo.K，[《無限維規則論：存在、系統、力量與張力的統一本體論》](https://unboundedaxiom.org/p/lm-000762/)，2026-05-26。作為廣義有向影響的歷史來源；本篇不接受「所有力量皆守恆」作為物理定律。 - **[I-06-07]** Neo.K × AI 協作，《OmniTension：基於張力語法的萬象力學統一理論總綱》。作為歷史研究動機，不把其中尚未物理校準的張力算式直接當作牛頓、熱力學或場論的替代品。 - **[I-06-08]** Neo.K × Aletheia（2026-09-22），《全域耦合計算場論的形式核心》。作為計算模型組織的前置背景；與自然物理耦合實測不同。 ## C. 公開科學與工程文獻（E-06 編號） > **來源檢查方式：** 書目與 DOI／出版社或可信學術索引中的題名及摘要於 2026-10-16 線上核對；本文並非宣稱逐頁審閱所有付費全文。下列文獻是成熟研究基線與方法來源，不是用來佐證 GMDIS 的尚未實驗主張。投稿前應另行全文複核引用內容及頁碼。 **[E-06-01]** Pantaleone, J. (2002). “Synchronization of metronomes.” *American Journal of Physics*, **70**(10), 992–1000. DOI：[10.1119/1.1501118](https://doi.org/10.1119/1.1501118)。共同移動底座節拍器的傳統分析基線。 **[E-06-02]** Hoogesteger, M. M., Pogromsky, A. Y., & Nijmeijer, H. (2026). “Huygens' synchrony: Emergent modes of planar symmetrically coupled pendula.” *Chaos*, **36**(9), 093135. DOI：[10.1063/5.0331194](https://doi.org/10.1063/5.0331194)。對平面共同底座的排列幾何與平移／轉動模式提供實驗及多體模擬對照；不可由此推論所有排列均產生不同運動。 **[E-06-03]** Kada, H., Kondo, K., Sumino, S., & Tokuda, I. T. (2017). “Experimental study on synchronization of sound-coupled electronic metronomes.” *Nonlinear Theory and Its Applications, IEICE*, **8**(1), 15–24. DOI：[10.1587/nolta.8.15](https://doi.org/10.1587/nolta.8.15)。證明聲學訊號可成為某種電子節拍器裝置的實驗耦合通道；非一般機械底座的自動等價證據。 **[E-06-04]** Bach, J. S., & Bruus, H. (2018). “Theory of pressure acoustics with viscous boundary layers and streaming in curved elastic cavities.” *Journal of the Acoustical Society of America*, **144**(2), 766 起. DOI：[10.1121/1.5049579](https://doi.org/10.1121/1.5049579)。近壁黏性、聲致流與彈性腔體有效邊界的成熟處理。 **[E-06-05]** Joergensen, J. H., & Bruus, H. (2021). “Theory of pressure acoustics with thermoviscous boundary layers and streaming in elastic cavities.” *Journal of the Acoustical Society of America*, **149**(5), 3599 起. DOI：[10.1121/10.0005005](https://doi.org/10.1121/10.0005005)。聲學溫度場、熱黏性邊界層及聲致流的有效模型。 **[E-06-06]** Donea, J., Giuliani, S., & Halleux, J. P. (1982). “An arbitrary Lagrangian–Eulerian finite element method for transient dynamic fluid-structure interactions.” *Computer Methods in Applied Mechanics and Engineering*, **33**, 689–723. DOI：[10.1016/0045-7825(82)90128-1](https://doi.org/10.1016/0045-7825(82)90128-1)。移動界面及 ALE 傳統文獻。 **[E-06-07]** Donea, J., Huerta, A., Ponthot, J.-Ph., & Rodríguez-Ferran, A. (2004). “Arbitrary Lagrangian–Eulerian Methods.” In *Encyclopedia of Computational Mechanics*. Wiley. DOI：[10.1002/0470091355.ecm009](https://doi.org/10.1002/0470091355.ecm009)。ALE 綜述及數值實作邊界。 **[E-06-08]** Causin, P., Gerbeau, J. F., & Nobile, F. (2005). “Added-mass effect in the design of partitioned algorithms for fluid–structure problems.” *Computer Methods in Applied Mechanics and Engineering*, **194**(42–44), 4506–4527. DOI：[10.1016/j.cma.2004.12.005](https://doi.org/10.1016/j.cma.2004.12.005)。added-mass 耦合導致的分割演算法不穩定基線。 **[E-06-09]** van der Schaft, A. J., & Maschke, B. M. (2002). “Hamiltonian formulation of distributed-parameter systems with boundary energy flow.” *Journal of Geometry and Physics*, **42**(1–2), 166–194. DOI：[10.1016/S0393-0440(01)00083-3](https://doi.org/10.1016/S0393-0440(01)00083-3)。邊界能量流與分布參數 port-Hamiltonian 成熟架構。 **[E-06-10]** Grmela, M., & Öttinger, H. C. (1997). “Dynamics and thermodynamics of complex fluids. I. Development of a general formalism.” *Physical Review E*, **56**(6), 6620–6632. DOI：[10.1103/PhysRevE.56.6620](https://doi.org/10.1103/PhysRevE.56.6620)。GENERIC 非平衡可逆／不可逆統一形式。 **[E-06-11]** Öttinger, H. C., & Grmela, M. (1997). “Dynamics and thermodynamics of complex fluids. II. Illustrations of a general formalism.” *Physical Review E*, **56**(6), 6633–6655. DOI：[10.1103/PhysRevE.56.6633](https://doi.org/10.1103/PhysRevE.56.6633)。GENERIC 的領域範例。 **[E-06-12]** Kotyczka, P., & Lefèvre, L. (2019). “Discrete-time port-Hamiltonian systems: A definition based on symplectic integration.” *Systems & Control Letters*, **133**, 104530. DOI：[10.1016/j.sysconle.2019.104530](https://doi.org/10.1016/j.sysconle.2019.104530)。離散時間互連與能量一致性的強基線。 **[E-06-13]** Kuhl, E., Hulshoff, S., & de Borst, R. (2003). “An arbitrary Lagrangian Eulerian finite-element approach for fluid–structure interaction phenomena.” *International Journal for Numerical Methods in Engineering*, **57**, 117–142. DOI：[10.1002/nme.749](https://doi.org/10.1002/nme.749)。FSI 時空離散及活塞測試的傳統參照。 ## D. 本篇引用紀律 - 本篇引用既有論文只支援已發表模型、方法及其適用邊界；**不**將他人成果改寫成 GMDIS 首次發現。 - GMDIS-04 已明確校準早期 OmniTension 中「力量守恆」或「力就是張力梯度」等概念；本篇採受控繼承而不悄悄恢復該類未證物理等價。 - 若改用別的聲學／流體模型，必須更新界面能量式與參考文獻，不可以保留現稿的公式並假定照樣成立。 - 所列 DOI 和網路資料不代表外部學者支持 GRST／SIO 或本系列全部概念。 --- # 附錄 A｜符號、SI 單位與模型級別速查 | 符號 | 定義 | 典型單位 | 關鍵注意事項 | |---|---|---|---| | $L$ | 一維聲腔長度 | m | 真實移動邊界時不是永遠固定 | | $A$ | 活塞／聲腔截面積 | m² | 對多維彎曲面須改用界面積分 | | $q,w$ | 活塞位移／速度 | m, m/s | 真實擺角用 rad 與慣量、力矩替代 | | $m,k,\beta$ | 機械質量、剛度、阻尼 | kg, N/m, kg/s | 受控參數，不等於 GRST 勢 | | $\rho_0,c,K_a$ | 背景密度、聲速、體積模量 | kg/m³, m/s, Pa | $K_a=\rho_0c^2$ 依指定等熵小訊號模型 | | $p,v$ | 聲壓擾動、流體質點速度 | Pa, m/s | 非完整總壓力／非任何情況的全流場 | | $E_a,E_m$ | 聲擾動能、集中機械能 | J | 子系統能量劃界須一致 | | $D_m,D_a$ | 機械、聲學模型內耗散率 | W | 若熱／外輻射去向不同，分別記錄 | | $C_T,h_T$ | 集總熱容、熱外洩係數 | J/K, W/K | 不自動構成完整熱聲耦合 | | $\rho_f,v_f,P$ | 真實流體密度、速度、總壓力 | kg/m³, m/s, Pa | M3 模型的狀態，不與 M1 重複計能 | | $\sigma_s,\sigma_f$ | 固／流應力 | Pa | 透過法向牽引進介面功率 | | $C_a,M_a$ | 離散聲壓容積順應／速度慣性矩陣 | m³/Pa, kg | 有效離散能量的正定矩陣 | | $G,b,R_a$ | 聲學梯度／流量／阻尼耦合 | m², m², kg/s | $G$ 與 $-G^\top$ 互為功率配對 | | $\Gamma_{sf}$ | 固液／固氣共同界面 | m² 所在曲面 | 不能忽略移動界面與法向定義 | | $\mathfrak C,\mathcal F$ | SIO 容器與自由集合 | 非普遍 SI 型別 | 是結構概念，不是任意物理能量 | | $\Sigma_R,\sigma_R$ | GRST 全域勢與局部態勢 | 非普遍 SI 型別 | 不是壓力／新基本力／勢能 | **英譯提示：** 本篇使用 “multimedia coupled dynamics” 表示多種物理媒介共同作用；為減少與多媒體影音（multimedia）之語義歧義，對外摘要可改寫為 “multi-medium coupled dynamics” 或 “multiphysics coupling”，並保留中文正題不變。英文副標本次沿用既定研究系列用語，但在正式投稿前應依目標學科調整。 --- # 附錄 B｜三組推導的最小覆核程序 ## B.1 聲—固邊界功率 理想聲場： $$ \partial_t\left(\frac{p^2}{2K_a}+\frac{\rho_0v^2}{2}\right)=-\partial_x(pv). $$ 流體域 $x\in(0,L)$，外側剛性壁 $v(L)=0$，活塞端 $v(0)=w$。乘以面積積分： $$ \dot E_a=-A[pv]_{x=0}^{x=L}=Ap(0)w. $$ 機械端受到反向流體力 $-Ap(0)$，故 $\dot E_m=Fw-\beta w^2-Ap(0)w$。兩式之界面項相消。若方向約定反轉，兩邊符號必須**同時**反轉以維持總功率不變。 ## B.2 面通量矩陣的反對稱耦合 $C_a\dot p=-G^\top v+bw$， $M_a\dot v=Gp-R_av$。聲學能量時間導數為： $$ \dot E_a^{(h)} =p^\top(-G^\top v+bw)+v^\top(Gp-R_av) =p^\top bw-v^\top R_av. $$ 機械端 $m\dot w=-kq-\beta w-b^\top p+F$，故： $$ \dot E_m=Fw-\beta w^2-wb^\top p. $$ 兩者相加得 M-06-07。此推導與本篇任何宏觀「張力」名詞完全無關，是獨立的物理／線性代數檢核。 ## B.3 二次能量與隱式中點的代數恆等式 令 $E(z)=\frac12 z^\top H z$， $H=H^\top$。則： $$ E(z_{n+1})-E(z_n) =\left(\frac{z_{n+1}+z_n}{2}\right)^\top H(z_{n+1}-z_n). $$ 若 $\Delta z=\Delta t\,L_hz_{1/2}$，則： $$ \Delta E=\Delta t\,z_{1/2}^\top HL_hz_{1/2}. $$ 只要所定義 $L_h$ 的對稱能量部分恰好等於 $-\beta w^2-v^\top R_av$，便得到第 8 節的精確離散帳本。對一般非二次能量不能直接使用同一結論，必要時應採 discrete gradient 或其他結構保持算法並重新驗證。 --- # 附錄 C｜可獨立執行的 NumPy 數值驗證程式 本程式只依賴 Python 3 與 NumPy。它不需要下載物理實驗資料、不會調用網路服務，也不會存取任何私人檔案。預設重現第 16 節的兩組實驗。執行命令： ```bash python -m pip install numpy python gmdis06_validation.py ``` 以下程式碼為實際測試過的完整來源： ```python import numpy as np def validate(n=12, dt=1.0e-5, steps=3000, damp=True): # A power-conjugate finite-volume acoustic cavity + a piston oscillator. rho, c = 1.20, 343.0 area, length = 0.004, 0.25 m, k, beta = 0.35, 13.0, (0.04 if damp else 0.0) CT, hT = 8.0, 0.06 dx = length / n comp = np.eye(n) * (area * dx / (rho * c*c)) mass_v = np.eye(n - 1) * (rho * area * dx) drag = np.eye(n - 1) * (0.00012 if damp else 0.0) G = np.zeros((n - 1, n)) for j in range(n - 1): G[j, j] = area G[j, j + 1] = -area b = np.zeros(n) b[0] = area N = 2 + n + (n - 1) z = np.zeros(N) z[0] = 0.002 # initial piston offset, metres z[1] = 0.00 z[2] = 0.0 L = np.zeros((N, N)) L[0, 1] = 1.0 L[1, 0] = -k/m L[1, 1] = -beta/m L[1, 2:2+n] = -b/m L[2:2+n, 1] = np.linalg.solve(comp, b) L[2:2+n, 2+n:] = -np.linalg.solve(comp, G.T) L[2+n:, 2:2+n] = np.linalg.solve(mass_v, G) L[2+n:, 2+n:] = -np.linalg.solve(mass_v, drag) H = np.zeros((N, N)) H[0, 0], H[1, 1] = k, m H[2:2+n, 2:2+n] = comp H[2+n:, 2+n:] = mass_v assert np.allclose(L.T @ H + H @ L, -2 * np.diag([0.0, beta] + [0.0]*n + [drag[0,0] if n > 1 else 0.0]*(n-1)), rtol=1e-8, atol=1e-12) I = np.eye(N) P = np.linalg.solve(I - .5*dt*L, I + .5*dt*L) theta = 0.0 max_mechanical_ledger = 0.0 max_full_ledger = 0.0 energy0 = float(.5 * z @ H @ z) for _ in range(steps): z1 = P @ z mid = .5 * (z + z1) dloss = beta * mid[1]**2 + mid[2+n:] @ drag @ mid[2+n:] energy1 = float(.5 * z1 @ H @ z1) energy0_step = float(.5 * z @ H @ z) mech_res = (energy1 - energy0_step) + dt * dloss max_mechanical_ledger = max(max_mechanical_ledger, abs(mech_res)) theta1 = ((CT/dt - hT/2)*theta + dloss)/(CT/dt + hT/2) full_res = (energy1 + CT*theta1) - (energy0_step + CT*theta) + dt*hT*(theta1+theta)/2 max_full_ledger = max(max_full_ledger, abs(full_res)) z, theta = z1, theta1 return {"n":n,"dt":dt,"steps":steps, "damp":damp, "energy0_J":energy0,"energy_final_J":float(.5*z@H@z), "temperature_rise_K":theta, "max_step_mech_res_J":max_mechanical_ledger, "max_step_full_res_J":max_full_ledger, "peak_initial_accel_delta_m_s2":-area * 100.0/m} if __name__ == '__main__': for on in [False, True]: res=validate(n=12,dt=1e-5,steps=3000,damp=on) for k, v in res.items():print(f'{k}: {v}') print() ``` **重現結果的正確理解：** 在兩個測試中，約 $10^{-20}$ J 的逐步能量殘差是浮點算法對特定二次能量方程的數值一致性結果。聲學傳播、接觸、熱黏性及實物裝置的外部有效性仍需額外驗證。程式以單一線性熱節點接收所有有效阻尼，**不是**完整熱物理仿真。 --- # 附錄 D｜最小實驗記錄契約（研究預註冊草案） 以下 JSON 只作資料交換／記錄範例，**不是**聲稱 E4–E8 已執行。真實實驗可以另行擴張欄位，但不得刪除核心的單位、模型、外部作用與證據狀態。 ```json { "schema": "GMDIS-06-benchmark-v0.1", "paper": "GMDIS-06", "case_id": "E2-linear-acoustic-piston-demo", "study_type": "numerical_model_verification", "status": "executed_numerical_only", "world": { "reference_geometry": "1d_fixed_acoustic_cavity_small_signal", "fluid_medium": "air_linearized", "rho_kg_m3": 1.2, "sound_speed_m_s": 343.0, "area_m2": 0.004, "length_m": 0.25, "solid_model": "concentrated_linear_oscillator", "moving_mesh": false, "nonlinear_escapement": false, "real_physical_metronome": false }, "interfaces": { "power_contract": "piston_force_times_velocity_opposes_fluid_input", "velocity_continuity": "linearized_normal_velocity", "boundary_back_wall": "rigid_no_normal_flow" }, "representation": { "spatial_scheme": "pressure_cell_face_velocity_mimetic", "time_scheme": "implicit_midpoint", "pressure_cells": 12, "dt_seconds": 0.00001, "steps": 3000, "quantum_channel_required": false }, "inputs": { "initial_piston_displacement_m": 0.002, "force_input_N": 0.0, "thermal_external_input_W": 0.0 }, "observables": [ "q_m", "p_Pa", "v_m_s", "mechanical_acoustic_energy_J", "relative_thermal_temperature_K", "interfacial_power_W", "energy_ledger_residual_J" ], "claims": { "analytical_energy_ledger": "derived_model_conditional", "numerical_code_verification": "executed", "pde_spatial_convergence": "not_yet_demonstrated", "physical_device_validation": "not_tested", "GRST_independent_predictive_gain": "not_tested", "SIO_independent_predictive_gain": "not_tested" }, "required_comparators": [ "classical_acoustic_equations", "mechanical_ode_only", "port_hamiltonian_equivalent", "spatial_and_temporal_refinement" ], "open_questions": [ "measured_fluid_structure_interface", "acoustic_impedance_identification", "thermoviscous_losses", "real_escapement_energy_source", "synchronization_bifurcation" ] } ``` **正式存檔規則：** 每次實驗獨立保存原始數據、原始觀測時間戳、硬體配置、材料參數、程式版本、隨機種子（若有）、能量邊界及 SHA-256；不以結果是否漂亮作為是否保存的條件。對外公開時應移除個人身分、私密設備或商業內部敏感資訊。 --- # 附錄 E｜審查、負面結果及未來更新檢查表 ## E.1 物理恆等式 - [ ] 已分別標明 $p$ （聲壓擾動）與 $P$ （完整熱力學壓力）。 - [ ] 未把 $\rho$ （質量密度）當作資訊占據密度。 - [ ] 所有界面法向、速度匹配、牽引符號均有明確方向約定。 - [ ] 聲—固界面的一份功率只在各子域出現一正一負的對應項。 - [ ] 若界面有儲能或耗散，未繼續套用零儲能理想界面等式。 - [ ] 未同時計入 M3 完整流體內能與由它投影出的 M1 聲學擾動能。 - [ ] 熱損失去向有具體介質、外部埠或顯式未解析項。 - [ ] 外部擒縱能源、機械阻尼、碰撞能量跳變皆可逐項追蹤。 - [ ] 無耗散極限、無耦合極限及幾何零效應均已保留。 ## E.2 數學與數值 - [ ] 已區分光滑方程積分恆等式與低正則弱解適定性。 - [ ] 不將有限維 ODE 的存在唯一性冒充完整 3D FSI 的適定性。 - [ ] 若改變壁面位置與流體域，已加入相容的移動網格／控制體積項。 - [ ] 有獨立可重算的機械、聲學、熱及總能量殘差。 - [ ] 有時間步、空間網格、接口收斂及模態頻率的分離驗證。 - [ ] 熱力學第二定律未由「機械能非增」直接推定。 - [ ] 帶擒縱碰撞的非光滑事件，不被當作全時刻連續可微。 - [ ] 所有模擬數字都標記數值示例，而非實物測量。 ## E.3 理論接縫與科研誠實 - [ ] GRST 的 Shi 只作未來結構與規則條件性分析，未量綱偷換成焦耳或牛頓。 - [ ] SIO 的合法自由度陰影與有限時間動態影響陰影已分開。 - [ ] 24／72 並未被認定為 72 種基本物理定律。 - [ ] 古典聲相干與量子 Q 軸的意義沒有混淆。 - [ ] 強化後的既有相位—振幅、聲學阻抗、FSI、pH、GENERIC、記憶核模型具有公平基線地位。 - [ ] 若新增多場模型沒有任務增益，允許降階而不是補救性改名。 - [ ] 負面結果與計算失敗記錄保留；不倒填模型成功指標。 - [ ] 沒有把尚未完成的真實節拍器矩陣實驗報成 PASS。 ## E.4 文件與學術流程 - [ ] UTF-8 解碼無錯誤、無取代字元 U+FFFD 與不合法控制字元。 - [ ] 每組 ` $...$ ` 與 `$$...$$` 完整配對且花括號配對。 - [ ] 原始 Markdown 保留真實 LaTeX，不由聊天畫面反向複製渲染結果。 - [ ] 內附程式與文中數值輸出可獨立重算。 - [ ] 正式投稿前將所有重要參考文獻核對到原文與引用頁碼。 - [ ] 若未來更新物理方程，保留原稿及重新驗證結果。 --- ## 文件狀態聲明 **GMDIS-06 v0.1 已完成詳細研究原稿。** 本篇在清楚界定的小擾動聲腔—集中機械模型與功率相容有限維離散化之下，具有可直接推導並由 NumPy 重新驗證的能量恆等式；所提出的完整可壓縮流體—固體—熱—擒縱模型仍是研究路線，未獲一般適定性或實物實驗驗證。GRST／SIO 的新增自然科學價值仍需受公平的成熟理論基線、反事實與外部資料檢驗。 **下一篇：** GMDIS-07〈持續作用與動態穩定結構〉。 **Canonical source：** 本 UTF-8 Markdown 原始檔。全部正式數學分隔符採用 ` $...$ ` 或 `$$...$$`；聊天畫面不是數學原稿的唯一可信來源。