昨天那篇算了一整天的冷卻鏈，中間有三個數字是直接借來用的：牛肉的導熱係數 k = 0.48 W/m·K、熱擴散係數 α = 1.4 × 10⁻⁷ m²/s、水浴的表面熱傳係數 h = 95 W/m²·K。然後又搬出一條長得像咒語的東西——θ = A₁·exp(−λ₁²·Fo)，λ₁ 由 λ tan λ = Bi 決定——把它當黑盒子用了五次，算出 50 分鐘、147 分鐘、21 分鐘、124 分鐘、4.6 小時。
那條式子從哪裡來？A₁ 是什麼？為什麼平板、圓柱、球的公式不一樣，而且不一樣的地方偏偏是指數而不是係數？h = 95 這個數字有多可靠，換成蒸氣箱、烤箱、冰箱各是多少？還有一個更實際的問題：既然 h 越大越快，那買一台更強的循環機能快多少？
今天把這本帳補完。這是第三本帳——前兩本是致死帳（F 值積分）與穩定化帳（限時冷卻），這一本是熱怎麼進去、又怎麼出來的帳。它是前兩本的共同前提：所有時間—溫度標準都寫在「產品最冷點」上，而最冷點的溫度歷程完全由這本帳決定。最後再用它把 FSIS Appendix A 那條 50–130°F 六小時的隱藏限制實算一遍，並算一算再加熱那條 74°C／兩小時的規定，在物理上到底卡住了多厚的東西。
📖 學
一、那條式子的來歷：分離變數與特徵值問題
熱傳導方程本身很單純：
Baldwin 寫成的那個帶 β 的形式，其實就是把 ∇² 在對稱座標下展開：
β = 0 是直角座標（大平板，r 從中心平面量起）、β = 1 是柱座標（長圓柱）、β = 2 是球座標。這三個不是三條不同的方程，是同一條方程在三種座標系下的樣子。
要解它，先無因次化。定義三個無因次量：
θ 是「還沒走完的路占全程的比例」（一開始 θ = 1，完全平衡時 θ = 0）；η 是無因次位置；Fo 是無因次時間。代進去之後，方程裡的 α、R 全部消失，只剩：
邊界條件是牛頓冷卻定律：表面傳出去的熱等於內部傳到表面的熱，−k ∂T/∂r|_R = h(T(R) − T∞)。無因次化之後：
**整個問題只剩兩個參數：Fo 和 Bi。**這就是無因次化的威力——不管你煮的是牛肉還是馬鈴薯、放在水浴還是烤箱、厚 2 公分還是 20 公分，只要 Bi 一樣，θ 對 Fo 的曲線就完全一樣。
接下來分離變數：假設 θ(η, Fo) = X(η)·G(Fo)。代入之後兩邊分別只含一個變數，所以都等於同一個常數（記為 −λ²）：
指數項就是這樣冒出來的。時間那一半的解永遠是指數衰減，沒有第二種可能。而衰減速率 λ² 不是自由參數——它由空間那一半的邊界條件挑出來，只有特定的離散值可以用。這些值就是特徵值（eigenvalue）。
空間那一半，三種座標各給出不同的特徵函數與不同的特徵方程：
| 幾何 | 特徵函數 X(η) | 特徵方程 | λ₁ 的上界 |
|---|---|---|---|
| 大平板（β = 0） | cos(λη) | λ tan λ = Bi | π/2 = 1.5708 |
| 長圓柱（β = 1） | J₀(λη) | λ·J₁(λ)/J₀(λ) = Bi | 2.4048 |
| 球體（β = 2） | sin(λη)/(λη) | 1 − λ cot λ = Bi | π = 3.1416 |
J₀、J₁ 是第一類貝索函數。三條特徵方程都是超越方程，有無窮多個根 λ₁ < λ₂ < λ₃ < …，完整解是全部加起來：
係數 A_n 由初始條件（θ = 1 處處成立）做正交展開得到。平板的：
圓柱的：
球體的：
到這裡為止都還是完整解。昨天用的「一項級數解」是把 n ≥ 2 的項全部丟掉：
丟得掉的原因是 λ₂ ≫ λ₁。以 Bi = 2 的平板為例，λ₁ = 1.0769、λ₂ ≈ 3.64。第二項的衰減速率是 λ₂² = 13.3，第一項只有 1.16——第二項衰減得快了 11.5 倍。所以只要等一下下，第二項就沒了。「等一下下」的量化門檻就是 Fo > 0.2，這條線後面會單獨算。
二、為什麼形狀改變的是指數，不是係數
這是整件事最漂亮、也最容易被誤解的一點。
很多人的直覺是：球比較快、板比較慢，差別大概是個「幾倍」的常數。錯了。看一下實際的特徵值——用數值方法解上面那三條特徵方程，Bi = 2.0 時：
| 幾何 | λ₁ | λ₁² | A₁ |
|---|---|---|---|
| 平板 | 1.0769 | 1.1597 | 1.1785 |
| 圓柱 | 1.5994 | 2.5582 | 1.3384 |
| 球體 | 2.0288 | 4.1159 | 1.4793 |
λ₁² 是指數裡的衰減速率，形狀把它從 1.16 拉到 4.12。而 A₁ 只從 1.18 走到 1.48。
這兩個變化的性質完全不同。A₁ 出現在對數裡，λ₁² 出現在分母上：
實算一次。取 θ = 0.1（走完九成路程）：
- 平板：Fo = ln(1.1785/0.1) / 1.1597 = ln(11.785) / 1.1597 = 2.4669 / 1.1597 = 2.127
- 圓柱：Fo = ln(1.3384/0.1) / 2.5582 = 2.5941 / 2.5582 = 1.014
- 球體：Fo = ln(1.4793/0.1) / 4.1159 = 2.6941 / 4.1159 = 0.655
平板／球體 = 3.25 倍。
再取 θ = 0.00909（昨天用的「距浴溫 0.5°C」標準）：
- 平板：Fo = ln(1.1785/0.00909) / 1.1597 = 4.8646 / 1.1597 = 4.195
- 圓柱：Fo = ln(1.3384/0.00909) / 2.5582 = ln(147.24) / 2.5582 = 4.9921 / 2.5582 = 1.952
- 球體：ln(1.4793/0.00909) / 4.1159 = ln(162.74) / 4.1159 = 5.0921 / 4.1159 = 1.237
平板／球體 = 3.39 倍。
注意這個比值隨著 θ 變小而變大：3.25 → 3.39 → 極限是 λ₁²(球)/λ₁²(板) = 4.1159 / 1.1597 = 3.549。
這就是答案。A₁ 的差異在對數裡，它只造成一個固定的 Fo 偏移；而且更諷刺的是，A₁ 的方向是反向的——球的 A₁ 比較大（1.48 > 1.18），意思是球心在一項近似下「起點看起來更落後」，這會讓球稍微變慢。但 λ₁² 在分母，它是乘性的，而且隨著時間拉長，對數項的相對重要性遞減，係數的影響被稀釋到零。
形狀改變的是衰減的速率常數，不是起跑點的位置。速率常數在指數裡，指數永遠贏。
那 λ₁ 的差異又是從哪來的？從特徵方程的上界。看回上面那張表：λ₁ 的上界分別是 π/2、2.4048、π。這三個數字有物理意義——在 Bi → ∞（表面瞬間到達環境溫度）的極限下，λ₁² 分別是：
比值 1 : 2.344 : 4.000。球體剛好是平板的整整四倍（π² 對 (π/2)²）。
另一頭，Bi → 0（集總熱容，內部溫度均勻，全部熱阻在表面）的極限也很有意思。這時 λ₁² ≈ (m+1)·Bi，其中 m = 0、1、2 對應板／柱／球。驗證一下 Bi = 0.25 的數值解：λ₁² = 0.2305 / 0.4700 / 0.7136，比值 1 : 2.04 : 3.10——確實逼近 1 : 2 : 3。
而 1 : 2 : 3 是什麼？是表面積對體積比。
- 半厚 R 的平板（中心絕熱）：A/V = 1/R
- 半徑 R 的長圓柱：A/V = 2πRL / (πR²L) = 2/R
- 半徑 R 的球：A/V = 4πR² / (4πR³/3) = 3/R
代進去確認：λ₁²·Fo → (m+1)·Bi·Fo = (m+1)·(hR/k)·(αt/R²) = (m+1)·h·t/(ρC_p R) = (A/V)·h·t/(ρC_p)。這正是集總熱容的標準指數。低 Bi 時形狀效應純粹是幾何的（面積比），高 Bi 時是波動的（特徵值比）。中間的 Bi = 2 落在 1:2.21:3.55，就在兩者之間。
這也順帶解釋了一件廚房裡的事：為什麼肉丸比同樣「粗細」的肉條熟得快這麼多，而不是快一點點。同樣半徑，球是柱的 1.6 倍快（4.1159/2.5582），而且這個倍數作用在指數上。把東西揉成球，是所有形狀操作裡回報最高的一個。
三、立方體怎麼辦：Newman 乘積解，以及 β 內插的代價
實際的肉塊不是無窮大的平板也不是無窮長的圓柱。Baldwin 給了一組實用的分數 β：立方體取 β = 1.25、方柱取 β = 0.70、2:3:5 的長方磚取 β = 0.28。
這是內插法，不是解出來的。真正的立方體有精確解，叫 Newman 乘積法則：三維的解等於三個一維平板解的乘積。
也就是說，等效的 λ₁² 是平板的三倍，等效的 A₁ 是平板的三次方。Bi = 2 時：λ₁²_eff = 3 × 1.1597 = 3.479，A₁_eff = 1.1785³ = 1.637。
拿它跟 β = 1.25 的線性內插比一比（內插值：λ₁² = 2.5582 + 0.25 × (4.1159 − 2.5582) = 2.948，A₁ = 1.3384 + 0.25 × (1.4793 − 1.3384) = 1.374）：
| 目標 θ | Newman 精確解 Fo | β = 1.25 內插 Fo | 差異 |
|---|---|---|---|
| 0.1 | 0.804 | 0.889 | +10.6% |
| 0.00909 | 1.493 | 1.702 | +14.0% |
**內插法算出來的時間比精確解長 10–14%。**方向是保守的（會叫你煮更久），所以拿來當安全邊際沒問題；但如果你在做的是逆算——例如想從表格反推「我這塊肉到底達成了幾個 log」——那 14% 的時間差就會變成致死量的高估。
（已標存疑：這組比較是我自己用 Newman 乘積與線性內插算的，Baldwin 本人沒有把 β = 1.25 的推導寫出來，我不確定他是用線性內插還是某種擬合。如果他用的是對 λ₁² 而非對 β 的擬合，誤差可能小得多。這個 10–14% 應該當成「內插法在這個 Bi 附近的量級」，不是他表格的實際誤差。）
四、h 到底是多少：從 8 到 10000 的六個數量級
上面所有推導裡，唯一需要「查外面」的數字就是 h。它不是食材的性質，是環境的性質——流體是什麼、流得多快、有沒有相變。工程手冊給的典型範圍大致是這樣（自然對流的空氣 5–25 W/m²·K，水與液體的自然對流 50–3000，冷凝蒸氣可以到數千甚至上萬）：
| 環境 | h（W/m²·K） | 說明 |
|---|---|---|
| 靜止空氣（冰箱內、流理台上） | 5–25，Baldwin 取 8 | 完全靠自然對流，最差的情況 |
| 家用旋風烤箱 | 14–30 | 強制對流但流體是空氣，密度太低 |
| 水浴（有循環） | 95–155（Baldwin） | 這是 sous vide 的工作區間 |
| 冰浴（有水、有攪動） | ≈ 100 | 昨天算冷卻用的值 |
| 冷凝蒸氣 | 1000–10000+ | 相變放出潛熱，數量級的跳躍 |
| 沸水 | 2500–25000 | 同樣是相變，反向 |
（已標存疑：這張表是工程手冊的範圍，不是你家設備的實測值。水浴的 h 對「循環泵流量／鍋子大小／袋子有沒有貼到別的袋子／有沒有用架子隔開」極度敏感；袋子疊在一起時局部 h 可能掉到 20 以下。烤箱的 h 隨風扇檔位、爐內位置、烤盤材質變化可以差三倍。這裡的數字只用來看數量級與趨勢。）
這張表的關鍵訊息是那個相變的斷層。空氣、水、蒸氣之間不是線性的差別，是數量級的差別。冷凝蒸氣之所以那麼猛，是因為它不是靠溫差傳熱，是靠相變潛熱——水的汽化潛熱在 100°C 約 2257 kJ/kg，而同樣一公斤水降 1°C 只放 4.18 kJ。所以蒸氣一碰到冷表面就凝結，把潛熱一次倒出來。
但這裡有一個很少被講的陷阱：冷凝的 h 對不凝結氣體極度敏感。蒸氣裡混進來的空氣會在液膜表面堆積成一層氣體邊界層，蒸氣要先擴散穿過它才能凝結。文獻上常見的說法是，體積分率百分之幾的空氣就能讓冷凝 h 掉一半以上。這就是為什麼萬能蒸烤箱在「純蒸氣」模式與「熱風＋蒸氣」模式下的表現差這麼多——後者本質上是在冷凝介面塞了一堆空氣。
還有一層更根本的限制，跟蒸氣多猛沒關係：**sous vide 的水浴溫度是 55–65°C，這個溫度的蒸氣壓遠低於大氣壓，你根本得不到純蒸氣環境。**低溫蒸烤箱裡跑的是「濕熱空氣」，h 頂多是幾十的量級，不是幾千。真正能拿到冷凝 h 的只有 100°C 的蒸鍋。所以「用蒸的比水浴快」在低溫域是完全錯的。
五、加大馬力的邊際報酬：一個很掃興的計算
現在把上面那張 h 表整個代進一項級數解，看看差多少。取一塊 40 mm 厚的平板肉（半厚 R = 0.02 m），k = 0.48 W/m·K，α = 1.4 × 10⁻⁷ m²/s，目標是中心到達「距環境溫度 0.5°C 以內」（θ = 0.00909，跟昨天同一個標準）。R²/α = 4 × 10⁻⁴ / 1.4 × 10⁻⁷ = 2857 秒。
| 環境 | h | Bi | λ₁ | λ₁² | A₁ | Fo | 時間 |
|---|---|---|---|---|---|---|---|
| 靜止空氣 | 8 | 0.33 | 0.5472 | 0.2994 | 1.0495 | 15.86 | 12 小時 35 分 |
| 旋風烤箱 | 25 | 1.04 | 0.8732 | 0.7624 | 1.1225 | 6.32 | 5 小時 1 分 |
| 水浴（弱） | 95 | 3.96 | 1.2621 | 1.5930 | 1.2281 | 3.080 | 2 小時 27 分 |
| 水浴（強） | 155 | 6.46 | 1.3628 | 1.8573 | 1.2505 | 2.651 | 2 小時 6 分 |
| 冷凝蒸氣 | 2000 | 83.3 | 1.5522 | 2.4092 | 1.2730 | 2.051 | 1 小時 38 分 |
| 完美傳熱 | ∞ | ∞ | 1.5708 | 2.4674 | 1.2732 | 2.003 | 1 小時 35 分 |
從這張表可以讀出三件事。
**第一，h 從 95 拉到 155（增加 63%），時間只從 147 分鐘降到 126 分鐘，快了 14%。**買一台流量大一倍的循環機，換來的是這個。
第二，h 從 95 直接拉到無窮大（花無限多錢，買到物理上不可能的設備），時間從 147 分鐘降到 95 分鐘，只快了 35%。
這是 h 這條路的理論天花板。因為 Bi = 3.96 的時候，內部熱阻已經是主導項；再怎麼改善表面，也只能拿走那 35%。
**第三，反過來，從靜止空氣（h = 8）升級到水浴（h = 95），時間從 755 分鐘掉到 147 分鐘——快了 5.1 倍。**同樣是「改善 h」，在 Bi 小的那一端回報是巨大的，在 Bi 大的那一端幾乎沒有。
這就是 Bi 這個數的全部意義：它告訴你錢該花在哪一邊。Bi ≪ 1，瓶頸在表面，改善 h 立刻見效；Bi ≫ 1，瓶頸在內部，唯一的出路是改幾何——切薄、切小、揉成球。
再看薄件。同樣的計算對 20 mm 平板（R = 0.01，Bi = 1.98，R²/α = 714 秒）：水浴 h = 95 是 Fo = 4.218 → 50.2 分鐘；完美傳熱是 Fo = 2.003 → 23.8 分鐘。省下 53%。薄件的 Bi 小，所以表面熱阻占比高，改善 h 的回報就大。
六、順便，昨天那個 2.94 的完整解釋
昨天算過：Baldwin 的升溫表上 20 mm 平板是 50 分鐘、40 mm 是 150 分鐘，比值 2.94，不是「厚度加倍時間四倍」的 4。當時的解釋是「Bi 變大之後所需 Fourier 數從 4.195 掉到 3.083」。現在可以說得更完整。
在 Bi → ∞ 的極限下，λ₁ 與 A₁ 不再隨 R 變化（都固定在 π/2 與 4/π = 1.2732），所以所需的 Fo 對兩個厚度是同一個數 2.003。這時：
實算驗證：95.4 / 23.8 = 4.01（差在四捨五入）。
**「厚度加倍時間四倍」是一條在 h = ∞ 才成立的定律。**它在有限 h 下必然低於 4，因為厚件的 Bi 比較大，所需 Fo 比較小，部分抵銷了 R² 的成長。水浴（h = 95）下是 2.94，靜止空氣（h = 8）下會更低——20 mm 在 h = 8 的 Bi 只有 0.167，λ₁² = 0.1622，A₁ = 1.0273，Fo = ln(1.0273/0.00909)/0.1622 = 4.7284/0.1622 = 29.15，t = 29.15 × 714 = 20813 秒 = 347 分鐘；40 mm 是 755 分鐘。比值 2.18。
整理成一句話：厚度加倍的時間倍數，隨著 h 從 0 到 ∞，從 2 單調爬到 4。（h → 0 的極限是集總熱容，t ∝ V/A ∝ R，所以是 2 倍。）水浴的 2.9 就落在中間偏上。
七、Heisler 圖，與 Fo > 0.2 那條護欄
在有電腦之前，上面這些超越方程沒辦法隨手解。1947 年 M. P. Heisler 把平板、圓柱、球的解畫成一組對數座標的圖表（Trans. ASME 69 期），後來配上 Gröber 的累積熱量圖，成為傳熱學教科書裡的標準工具。用法是：先算 Bi，在圖上找那條 1/Bi 的曲線，再用 Fo 讀出中心的 θ。
Heisler 圖畫的其實就是一項級數解。所以它繼承了同一個限制：**Fo > 0.2 才準。**這條線的來由是 λ₂ 那一項的殘量——Fo = 0.2 時第二項的貢獻已經降到 2% 以下，可以忽略。Fo < 0.2 的時候，第二、三、四項都還醒著，一項近似會失真。
這條護欄在廚房裡對應到什麼？把 Fo = 0.2 換算成實際時間（α = 1.4 × 10⁻⁷）：
| 半厚 R | Fo = 0.2 對應的時間 | R²/α |
|---|---|---|
| 5 mm | 36 秒 | 179 秒 |
| 10 mm | 2.4 分鐘 | 714 秒 |
| 15 mm | 5.4 分鐘 | 1607 秒 |
| 20 mm | 9.5 分鐘 | 2857 秒 |
| 30 mm | 21.4 分鐘 | 6429 秒 |
物理意義比數學意義更有用：在 Fo < 0.2 的那段時間裡，中心根本還沒「感覺到」表面。熱擾動從表面傳到中心需要時間，在那之前中心的溫度基本上一動不動，材料表現得像半無窮大固體（用 erfc 解，跟 R 完全無關）。
所以一塊 60 mm 厚的肉下水之後，前面 21 分鐘中心是不動的。這解釋了很多人插針量溫時的困惑：下水二十分鐘，探針讀數幾乎沒變，於是懷疑機器壞了。沒壞，這是 Fourier 數還沒到 0.2。
昨天那五個計算的 Fo 分別是 4.195、3.083、0.624、3.624、2.028，全部大於 0.2，所以一項解都是有效的。但今天要算的 come-up time 會撞到這條線，下一節就是。
八、CUT 實算：一塊 60 mm 圓柱吃掉多少「六小時預算」
昨天提到 FSIS Appendix A（FSIS Cooking Guideline for Meat and Poultry Products, Revised Appendix A，2021 年 12 月）在致死時間—溫度表的註腳裡建議：產品溫度停留在 50–130°F（10–54.4°C）之間的總時間限制在 6 小時以內，目的是控制金黃色葡萄球菌的耐熱腸毒素。當時只說了「很容易吃光」，沒算。今天算。
（已標存疑：FSIS 官網上這份指引的文件編號，我查到的字串在 2021-0014 與 FSIS-GD-2021-14 之間不一致；昨天引的 Appendix B 是 FSIS-GD-2021-13。兩份是同一批於 2021 年 12 月 14 日在 Federal Register 公告的（86 FR 71007），實施驗證從 2022 年 12 月 14 日起。編號請以官網頁面為準。）
情境：一塊 60 mm 直徑的圓柱狀肉（想像一條綁好的豬里肌或整條雞胸捲），從冰箱 4°C 直接下 55°C 水浴，h = 95。
先算參數：R = 0.03 m，Bi = 95 × 0.03 / 0.48 = 5.94。解 λ·J₁(λ)/J₀(λ) = 5.94 得 λ₁ = 2.0458，λ₁² = 4.1853，A₁ = 1.5242。R²/α = 9 × 10⁻⁴ / 1.4 × 10⁻⁷ = 6429 秒。
中心進入管制區間（到 10°C）：θ = (55 − 10)/(55 − 4) = 45/51 = 0.8824。 Fo = ln(1.5242 / 0.8824) / 4.1853 = 0.5468 / 4.1853 = 0.1306。
**Fo = 0.131 < 0.2，一項解在這裡無效。**所以我另外跑了一次一維有限差分數值解（80 個網格點，顯式格式），結果是 12.3 分鐘；一項解算出來是 14.0 分鐘。一項解高估了 14%。方向是「說它比較慢進入危險區」，對計算區間停留時間來說是反向的樂觀，要注意。
中心離開管制區間（到 54.4°C）：θ = (55 − 54.4)/51 = 0.6/51 = 0.01176。 Fo = ln(1.5242 / 0.01176) / 4.1853 = ln(129.6) / 4.1853 = 4.8646 / 4.1853 = 1.1623。 t = 1.1623 × 6429 = 7472 秒 = 124.5 分鐘。數值解給 124.5 分鐘，完全吻合（這裡 Fo 遠大於 0.2）。
途中的檢查點（數值解）：到 20°C 是 20.0 分鐘、到 26.7°C 是 25.7 分鐘（一項解給 25.9，Fo = 0.241，剛過門檻就對上了）、到 40°C 是 42.1 分鐘、到 50°C 是 70.2 分鐘。最後那 4.4°C 用掉了 54 分鐘——因為驅動溫差已經只剩幾度，指數尾巴拖得很長。
升溫段的管制區間停留時間 = 124.5 − 12.3 = 112.2 分鐘 = 1.87 小時。
六小時預算已經用掉 31%，而肉才剛煮好。
接著算冷卻段。用昨天冰浴那組參數（α = 1.1 × 10⁻⁷，因為熟肉的含水率與熱性質已經改變），從中心 54.4°C 開始降：
冰浴（h = 100，浴溫 0°C）（數值解）：
- 到 26.7°C（80°F）：36.4 分鐘
- 到 10°C（50°F）：68.0 分鐘
- 到 7.2°C（45°F）：78.5 分鐘
- 到 4.4°C（40°F）：94.3 分鐘
冰箱靜止空氣（h = 8，箱溫 4°C）（數值解）：
- 到 26.7°C：139.6 分鐘（2.33 小時）
- 到 10°C：344.6 分鐘（5.74 小時）
- 到 7.2°C：441.4 分鐘
- 到 4.4°C：永遠到不了（箱溫就是 4°C，θ 永遠 > 0）
現在對照法規。
FSIS Appendix B Option 1.1（130→80°F ≤ 1.5 小時、80→40°F ≤ 5 小時）：
- 冰浴：第一段 36.4 分鐘 ✅、第二段 94.3 − 36.4 = 57.9 分鐘 ✅ → 通過
- 冰箱：第一段 2.33 小時 > 1.5 小時 ❌ → 不通過。而且也過不了 Option 1.5（2 小時）；1.6–1.8 是從 126°F 起算的、1.3–1.4 需要亞硝酸鹽或 a_w 柵欄，這塊肉都沒有。八個選項全滅。
Appendix A 的 50–130°F 六小時：
- 冰浴路線：升溫 1.87 小時 + 冷卻 68.0 分鐘 = 1.87 + 1.13 = 3.00 小時 ✅（用掉預算的一半）
- 冰箱路線：升溫 1.87 小時 + 冷卻 5.74 小時 = 7.61 小時 ❌（超標 27%）
**兩本帳同時掛掉。**而這只是一塊六公分粗的肉、一個完全常見的家庭做法。
這裡有個值得停一下的觀察：**升溫那 1.87 小時是你付不掉的固定成本。**它由 R²、α、Bi 決定，除了切薄之外沒有辦法縮短——你把水浴溫度調高會加快，但那就改變了你要的熟度。所以六小時的預算裡，光是「把肉煮熟」就先扣掉三分之一，剩下四小時全部得留給冷卻。這就是為什麼冰浴不是可選項，是必需品。
順帶修一個昨天講得不夠精確的地方：昨天說「55°C 水浴的整段升溫都在管制區間裡」。嚴格來說不對——區間下緣是 10°C，而肉從冰箱出來是 4°C，所以最前面那 12 分鐘（中心 4 → 10°C）不算在內。這 12 分鐘是白送的。對六小時預算來說是 3.4%，不多，但算帳就要算準。
九、再加熱那本帳：74°C、兩小時、與物理上的厚度上限
冷卻完之後還有一段：吃之前的再加熱。這一段的法規在 FDA Food Code 3-403.11（華盛頓州對應條文 WAC 246-215-03440）：經過烹調、冷卻的 TCS 食品，若要重新熱持存，必須讓所有部位在 2 小時內達到 165°F（74°C）持續 15 秒。微波再加熱另有規定：達 165°F、並且要翻攪或旋轉、加蓋、加熱後加蓋靜置 2 分鐘。商業預製的即食品門檻較低，135°F（57°C）即可。
三個數字，三個不同的理由。
74°C／15 秒是把致死量拉滿到瞬時區。對照 Appendix A：70°C 以上所需的 log 削減是瞬間達成的，74°C 留了 4°C 的邊際，涵蓋了探針誤差與產品內的溫度分布。
2 小時管的是再加熱過程本身穿過危險帶的時間。這一段跟前面 CUT 是同一本帳——再加熱時產品又要從 5°C 爬到 74°C，中間整段都在 10–54.4°C 裡面。
微波要靜置 2 分鐘，是因為微波加熱的能量沉積極度不均勻（駐波、穿透深度、介電常數隨溫度變化），關機瞬間的溫度分布可以差二三十度。靜置讓熱傳導把峰谷抹平——這是一段純粹的、寫進法規裡的傅立葉方程。
現在用今天的工具算一件實用的事：兩小時的上限，物理上卡住了多厚的東西？
40 mm 平板（R = 0.02，Bi = 3.96，λ₁² = 1.5930，A₁ = 1.2281，R²/α = 2857 秒），從 4°C 加熱到中心 74°C：
- 80°C 水浴：θ = (80 − 74)/(80 − 4) = 6/76 = 0.07895。Fo = ln(1.2281/0.07895)/1.5930 = ln(15.556)/1.5930 = 2.7444/1.5930 = 1.723。t = 1.723 × 2857 = 4922 秒 = 82.0 分鐘
- 85°C 水浴：θ = 11/81 = 0.1358。Fo = ln(9.043)/1.5930 = 2.2020/1.5930 = 1.382。t = 65.8 分鐘
- 90°C 水浴：θ = 16/86 = 0.1860。Fo = ln(6.602)/1.5930 = 1.8873/1.5930 = 1.185。t = 56.4 分鐘
三個都在兩小時內，但注意浴溫從 90 掉到 80，時間就從 56 分鐘漲到 82 分鐘——漲了 45%，只因為驅動溫差從 86°C 縮到 76°C。θ 在對數裡，而它對浴溫的敏感度在接近目標溫度時會爆炸。這就是為什麼「用 75°C 的水浴回溫到 74°C」在紙上看起來只差 1°C，實際上是永遠到不了。
解出兩小時剛好用完的厚度（用 brentq 反解 t(R) = 7200）：
| 加熱環境 | 最大半厚 | 最大全厚 |
|---|---|---|
| 80°C 水浴（h = 95） | 25.1 mm | 50 mm |
| 85°C 水浴（h = 95） | 28.5 mm | 57 mm |
| 175°C 旋風烤箱（h = 25） | 42.0 mm | 84 mm |
（已標存疑：175°C 烤箱那一列的物理假設很粗糙——它把爐溫當成常數表面驅動溫度，忽略了表面蒸發冷卻與結殼。實際上肉表面在水分蒸發完之前會被鎖在接近 100°C，這會大幅降低有效驅動溫差；結殼之後 k 與 α 又都改變。這個 84 mm 應該當上界看，不是可用的數字。水浴那兩列因為沒有蒸發、驅動溫差明確，可信度高得多。）
換算成廚房語言：**一塊超過 5–6 公分厚的冷藏熟肉，用 80–85°C 的水浴回溫，會撞到兩小時的紅線。**要嘛切開，要嘛提高浴溫，要嘛不要走再加熱這條路（改成冷食，或改成小份量分裝）。
十、log 削減可以相加嗎
最後一個問題，也是最容易搞錯的一個。既然第一次烹調做了 7-log，再加熱到 74°C 又是一次「瞬時致死」，那總共是不是 14-log？或者說，第一次沒煮夠、再加熱補一次，是不是就補回來了？
**不能這樣加。**至少有三個獨立的理由。
**理由一：存活下來的不是原來那群菌。**對數線性的 D 值模型假設族群同質——每個細胞的死亡機率一樣。真實族群不是。第一次熱處理殺掉的是最脆弱的那 99.99999%，活下來的是最耐熱的那條尾巴。對這個新族群再做一次熱處理，它的 D 值比原始族群大。這就是存活曲線上「拖尾」（tailing）的來源，也是為什麼 Weibull 模型（log N = −b·t^n）在很多資料集上比對數線性擬合得更好。把兩次處理的 log 直接相加，等於假設 n = 1 且族群同質，這兩個假設都不成立。
理由二：熱休克誘導耐熱性。細菌在亞致死溫度（42–54°C）停留一段時間會啟動熱休克反應——GroEL、DnaK 這些分子伴護蛋白大量表現，rpoE、rpoH、htrA 這些基因被誘導，膜流動性下降。文獻報導經過熱休克的沙門氏菌，在雞肉、火雞肉與實驗培養基中的 D 值都顯著高於未休克的對照組。而慢速的升溫與慢速的冷卻，本身就是一場長時間的熱休克。這是低溫慢煮相對於快速烹調的一個結構性劣勢：你在殺菌的同時，也在訓練倖存者。
（已標存疑：我沒有查到一組能直接引用的「D₆₀ 上升幾倍」的數字，文獻報導的幅度隨菌株、基質、休克溫度與時間差異很大。這裡只主張方向——熱休克提高耐熱性、且效果顯著——不主張倍數。要做製程驗證的話必須用你自己基質的實測資料，不能套文獻值。）
**理由三，也是最重要的：再加熱對毒素與芽孢幾乎完全無效。**這一點必須分開講，因為三種毒素的行為完全不同：
- 肉毒桿菌毒素是耐熱性差的蛋白質。Woodburn 等人 1979 年在《Journal of Food Science》上的研究測了 A、B、E、F 四型在數種食品與緩衝液中的失活速率，結論是 85°C 5 分鐘或 79°C 20 分鐘可以破壞每公克 10³ LD₅₀ 的毒素；CDC 的公衛指引也採用 85°C 至少 5 分鐘這個數字。所以再加熱確實能處理肉毒毒素——但要 85°C 而不是 74°C，而且要所有部位都到。同一份研究還指出毒素在酸性食品（pH 4.2 的番茄湯）中比在低酸食品（pH 6.2 的玉米罐頭）中更穩定，這跟直覺相反。
- 金黃色葡萄球菌腸毒素是耐熱的。它能撐過一般的烹調、再加熱，甚至商業殺菌。這就是 Appendix A 那條 50–130°F 六小時規則的全部理由——產毒之後就沒有補救手段了，只能事前不讓它產。
- 仙人掌桿菌的嘔吐毒素 cereulide 是一個環狀十二胜肽，耐熱耐酸耐蛋白酶到近乎荒謬的程度，一般烹調與再加熱對它完全沒有作用。這是為什麼隔夜炒飯的問題不能靠「熱透一點」解決。
而芽孢在 74°C 面前毫髮無傷——昨天算過，Group II 肉毒桿菌芽孢的 6D 在 75°C 需要 520 分鐘。
所以正確的心智模型是：再加熱是熱持存的入場券，不是第二次殺菌。它處理的是冷卻後可能重新長起來的營養細胞，以及（在 85°C 以上時）已經產生的肉毒毒素。它不處理芽孢，不處理耐熱腸毒素，也不能拿來補救一次不足的烹調。
更正 1：「厚度加倍時間四倍」——四倍是 h = ∞ 的極限，實際落在 2 到 4 之間
昨天已經指出實際比值是 2.94 不是 4，但只給了「Fo 從 4.195 掉到 3.083」這個現象描述。完整的說法是：
在 Bi → ∞ 的極限下，λ₁ 與 A₁ 不再隨 R 變化（固定在 π/2 與 4/π），所需 Fo 對所有厚度都是同一個數，於是 t ∝ R² 精確成立，比值恰好 4.00。實算驗證 95.4/23.8 = 4.01。
在 Bi → 0 的極限下，退化成集總熱容，t ∝ ρVC_p/(hA) ∝ R，比值是 2.00。
所以這個倍數是單調地從 2 爬到 4 的。水浴（h = 95）：20→40 mm 是 2.94。靜止空氣（h = 8）：347 分鐘 → 755 分鐘，是 2.18。烤箱大約在 2.5 附近。
實務推論：在冰箱裡冷卻時，切薄的回報比在水浴裡加熱時小。因為冰箱的 Bi 小，時間對厚度只有接近一次方的敏感度。所以冰箱降溫真正該做的第一件事不是切薄（那只買到 2 倍），而是把 h 從 8 拉到 100（那買到 5 倍以上）——也就是換成冰浴。這修正了一個很常見的優先順序錯誤。
更正 2：「買更強的循環機能明顯加快」——上限是 35%，而且你摸不到
這條在網路討論區裡經常出現：換更大流量的棒子、加裝循環泵、用更大的鍋。實算的結果很掃興。
40 mm 平板，h 從 95 拉到 155（流量提升六成以上），時間 147 → 126 分鐘，只快 14%。h 拉到無窮大（物理極限），147 → 95 分鐘，只快 35%。
原因是 Bi = 3.96 已經進入內部熱阻主導區。表面熱阻占總熱阻的比例大約是 1/(1+Bi) 這個量級，Bi ≈ 4 時就只剩兩成上下，你再怎麼優化都是在優化那兩成。
但這條規則有一個重要的例外，而且例外的方向是往下的。上面所有計算都假設袋子完全泡在流動的水裡。實際情況常常不是——袋子浮起來（部分表面在空氣中，那一面的 h 掉到 8）、袋子疊在一起（接觸面幾乎沒有對流）、鍋子太滿導致循環死角。這些情況下有效 h 可能掉到 20 甚至更低，而 Bi = 0.83 那一區，時間對 h 是敏感的。
所以正確的操作重點不是「加大馬力」，而是確保每一個袋子的每一面都真的接觸到流動的水：用架子隔開、壓重物防止浮起、不要一次塞太多、留出循環通道。這些是免費的，而且比換設備有效得多。
更正 3：「形狀只是個小修正」——球比板快 3.5 倍，而且差在指數
這是本文的核心結論，值得單獨列出來，因為它跟大多數人的直覺差距最大。
同樣的特徵長度 R、同樣的 Bi = 2，中心到達 θ = 0.00909 所需的 Fourier 數：平板 4.195、圓柱 1.952、球體 1.237。板／球 = 3.39，而且這個比值隨著要求越嚴格（θ 越小）越接近極限值 λ₁²(球)/λ₁²(板) = 4.1159/1.1597 = 3.55。
關鍵在於這個差異的性質：λ₁² 在指數裡，A₁ 在對數裡。A₁ 的差異（1.18 → 1.48）只造成一個固定的 Fo 偏移，而且方向還是反的；λ₁² 的差異是乘性的、永久的。時間拉得越長，係數越不重要，指數越支配一切。
λ₁² 的來源本身也分兩段：低 Bi 時 λ₁² → (m+1)Bi，比值 1:2:3，就是表面積體積比；高 Bi 時比值 1:2.344:4.000，是特徵方程根的比值（(π/2)²、2.4048²、π²）。Bi = 2 落在中間，是 1:2.21:3.55。
順帶修正一個技術細節：Baldwin 給的分數 β（立方體 1.25、方柱 0.70、長方磚 0.28）是內插近似，不是解出來的。立方體有精確解——Newman 乘積法則，三個一維平板解相乘，等效 λ₁² 是平板的三倍、等效 A₁ 是平板的三次方。Bi = 2 時精確解給 Fo = 1.493，β = 1.25 的線性內插給 1.702，內插高估 14%（方向保守）。
更正 4：一項級數解不能隨便用——Fo < 0.2 時它會說謊
昨天連續用了五次 θ = A₁·exp(−λ₁²Fo)，沒有檢查適用範圍。今天補上：這個近似只在 Fo > 0.2 成立，因為它把第二項以後全部丟掉了，而 Fo < 0.2 時那些項還醒著。
實際撞到的例子就在本文第八節：60 mm 圓柱中心從 4°C 升到 10°C，一項解算出 Fo = 0.131，落在門檻以下。用有限差分數值解對照，真實答案是 12.3 分鐘，一項解給 14.0 分鐘，高估 14%。
這個誤差的方向要特別注意。在升溫問題裡，一項解說「比較慢」聽起來是保守的；但如果你在算的是「產品在管制區間停留多久」，一項解把進入時刻往後推，會讓你低估停留時間——那就是反向的樂觀。本例中一項解算出的區間停留是 110.5 分鐘，數值解是 112.2 分鐘。差 1.7 分鐘不致命，但如果換成薄件或大溫差的情境，這個偏差會放大。
好消息是昨天那五個計算的 Fo（4.195、3.083、0.624、3.624、2.028）全部安全過關。壞消息是這件事必須每次都檢查，而檢查方法很簡單：算完 Fo 之後看一眼是不是大於 0.2；不是的話，那段時間裡物體行為接近半無窮大固體，該用 erfc 解或直接跑數值。
至於 Fo < 0.2 的物理意義，比數學意義更值得記：**熱擾動還沒傳到中心。**60 mm 的肉塊下水後前 21 分鐘、40 mm 前 9.5 分鐘、20 mm 前 2.4 分鐘，中心溫度幾乎不動。這不是設備問題。
更正 5：「再加熱到 74°C 等於重新殺一次菌」——三重錯誤
這是最危險的一條誤解，因為它會被拿來合理化「第一次沒煮夠沒關係」。
**錯誤一，log 不能相加。**存活者是原族群裡最耐熱的尾巴，對它們再做一次處理，適用的 D 值更大。這是拖尾現象，也是 Weibull 模型比對數線性模型更貼合真實資料的原因。
**錯誤二，緩慢的升降溫會誘導耐熱性。**亞致死溫度停留會啟動熱休克反應（GroEL、DnaK、rpoE/rpoH/htrA），文獻一致報導沙門氏菌、李斯特菌、大腸桿菌 O157:H7 在熱休克之後的 D 值上升。sous vide 的慢升溫與冰箱的慢降溫，都是熱休克。
錯誤三，也是最關鍵的：毒素與芽孢的行為完全不同。
| 對象 | 74°C 再加熱 | 85°C 5 分鐘 | 121°C |
|---|---|---|---|
| 營養細胞（沙門氏菌等） | 有效 | 有效 | 有效 |
| 肉毒桿菌毒素 | 不足 | 有效（Woodburn 1979；CDC 指引） | 有效 |
| 金黃色葡萄球菌腸毒素 | 無效 | 無效 | 大致無效 |
| 仙人掌桿菌 cereulide | 無效 | 無效 | 無效 |
| 各種芽孢 | 無效 | 無效（Group II 6D 需 90°C×10 min） | 有效（12D botulinum cook） |
這張表的實務讀法：再加熱是熱持存的入場券，不是補考機會。Food Code 3-403.11 那條 74°C／15 秒／2 小時，管的是「冷卻後可能重新長起來的營養細胞」，並且用 2 小時的上限限制再加熱過程本身通過危險帶的時間。它預設的前提是第一次烹調與冷卻都已經合格。
也因此，**產毒是不可逆的。**這句話應該和昨天那句「要嘛留在 55°C 以上，要嘛快速衝到 5°C 以下」放在一起記——後者講的是不要給它機會，前者講的是給了就沒得反悔。
🧠 記
- 無因次化之後，整個熱傳導問題只剩兩個參數：Fo = αt/R²（無因次時間）與 Bi = hR/k（內外熱阻比）。Bi 相同的兩個系統，θ–Fo 曲線完全一樣。
- 分離變數必然給出 exp(−λ²Fo)：時間那一半只能是指數衰減。λ 由空間那一半的特徵方程挑選——平板 λ tan λ = Bi、圓柱 λJ₁(λ)/J₀(λ) = Bi、球體 1 − λ cot λ = Bi。
- **形狀改變的是指數（λ₁²），不是係數（A₁）。**Bi = 2 時 λ₁² 是 1.16 / 2.56 / 4.12（板/柱/球），A₁ 只有 1.18 / 1.34 / 1.48。板比球慢 3.39 倍，極限 3.55 倍。
- λ₁² 的形狀比在 Bi → 0 時是 1:2:3（等於表面積體積比 (m+1)/R），在 Bi → ∞ 時是 1:2.344:4.000（等於 (π/2)²、2.4048²、π²）。
- 立方體有精確解：Newman 乘積，λ₁²_eff = 3λ₁²(板)、A₁_eff = A₁³(板)。Baldwin 的 β = 1.25 內插高估時間約 10–14%（保守方向）。
- h 的量級：靜止空氣 5–25、旋風烤箱 14–30、水浴 95–155、冷凝蒸氣 1000–10000+。相變造成的斷層是數量級的，但 sous vide 溫域拿不到冷凝蒸氣。
- 加大馬力的天花板是 35%。40 mm 平板：h 95→155 只快 14%，h 95→∞ 只快 35%。Bi ≫ 1 時瓶頸在內部，錢要花在幾何上；Bi ≪ 1 時才該花在 h 上（8→95 是 5.1 倍）。
- 「厚度加倍時間四倍」的倍數隨 h 從 2 單調爬到 4：靜止空氣 2.18、水浴 2.94、完美傳熱 4.01。
- 一項級數解只在 Fo > 0.2 有效（Heisler 圖同此限制）。實例：60 mm 圓柱中心到 10°C，一項解 14.0 分鐘、數值解 12.3 分鐘，高估 14%。Fo < 0.2 的物理意義是「熱還沒傳到中心」——60 mm 前 21 分鐘中心不動。
- CUT 實算：60 mm 圓柱、4°C 下 55°C 水浴，中心 12.3 分鐘進 50°F、124.5 分鐘出 130°F，管制區間停留 1.87 小時，先吃掉 Appendix A 六小時預算的 31%。
- 同一塊肉的冷卻：冰浴第一段 36.4 分鐘（過 Option 1.1）、區間停留 1.13 小時，總計 3.00 小時 ✅；冰箱第一段 2.33 小時（八個選項全滅）、區間停留 5.74 小時，總計 7.61 小時 ❌。
- Food Code 3-403.11：冷卻後再熱持存需 2 小時內、所有部位 74°C／15 秒；微波另需翻攪、加蓋、靜置 2 分鐘（抹平微波的不均勻沉積）。實算兩小時上限：85°C 水浴約 57 mm、80°C 水浴約 50 mm。
- log 削減不能相加：存活者是耐熱尾巴（拖尾／Weibull）、慢升降溫會誘導熱休克耐熱性、而且 74°C 對芽孢與耐熱毒素完全無效。肉毒毒素要 85°C 5 分鐘（Woodburn 1979）；金黃色葡萄球菌腸毒素與 cereulide 煮不掉。
✍️ 實踐
一、算出你家最常做的那塊肉的 Bi，然後決定錢往哪裡花 量最厚處的一半當 R（公尺），乘以 95，除以 0.48。得到的數字如果大於 3，別再研究設備了——你的瓶頸在肉裡面，唯一有效的槓桿是切薄或切小。如果小於 1（薄片、小塊、或你是在冰箱裡冷卻），那改善 h 的回報很大，優先處理循環與接觸。這一個除法決定了你接下來所有的努力方向。
二、用有沒有「動」來檢驗你的鍋，而不是用機器規格 下次水浴運轉時，丟一小撮乾燥香草或一小片蔬菜屑進去，看它有沒有沿著鍋壁繞一圈回來。如果它卡在某個角落打轉，那個角落就是 h ≈ 20 的死區。把袋子架起來、彼此隔開、壓上重物防浮，然後再丟一次。這比買新機器有效，而且免費。
三、把「起鍋後 21 分鐘」寫進你的計時器 拿你最常做的厚度算一次 Fo = 0.2 的時間（R² / α × 0.2，α 取 1.4 × 10⁻⁷）。5 mm 半厚是 36 秒、10 mm 是 2.4 分鐘、15 mm 是 5.4 分鐘、20 mm 是 9.5 分鐘、30 mm 是 21.4 分鐘。在這段時間之內插針讀到的中心溫度沒有資訊量，不要據此判斷任何事。設一個鬧鐘，時間到了再開始認真讀數。
四、把再加熱的厚度上限量出來，寫在冰箱上 如果你有分裝冷藏熟食的習慣，今天做一件事：規定所有要再加熱的分裝品厚度不超過 4 公分。理由是上面那個 50–57 mm 的計算——4 公分留了三成邊際，涵蓋你家水浴溫度不穩、袋子貼在一起、冰箱溫度比想像中低這些現實。超過 4 公分的東西，改成冷食或當天吃完。
五、實測一次你的圓柱與球 拿同一份絞肉，一半捏成直徑 5 公分的球，一半塑成直徑 5 公分的柱（長度至少 15 公分），一起下 65°C 水浴，各插一支探針。記下兩者到達 60°C 的時間。理論預測球會快約 1.6 倍（λ₁² 之比 4.19/2.56，Bi ≈ 5.9）。這個實驗會把「形狀在指數裡」這件事變成你手上的一條曲線，比任何公式都好記。
六、給自己一條再加熱的鐵律 再加熱只用在「第一次已經合格烹調並合格冷卻」的東西上。任何「上次好像沒煮透」「放了幾天不確定」的東西，不要靠再加熱補救——芽孢與耐熱毒素不吃這一套。這條規則沒有例外，也不需要判斷力，這正是它的價值。
🔗 延伸學習
- Douglas Baldwin, A Practical Guide to Sous Vide Cooking：附錄 A 完整公開了帶 β 幾何因子的熱傳導模型、h 與 α 的取值、以及所有升溫表與巴氏殺菌表的計算方法。本文所有反算與比對都是用他公開的參數做的；本文對 β = 1.25 立方體近似的檢驗也是針對這份文件。
- FSIS Cooking Guideline for Meat and Poultry Products (Revised Appendix A)：50–130°F 六小時的建議、按脂肪含量與禽種分欄的致死時間—溫度表、以及相對濕度作為關鍵操作參數的規定，都在這份文件裡。與昨天引用的 Appendix B（FSIS-GD-2021-13）是同一批於 2021 年 12 月 14 日公告（86 FR 71007）、2022 年 12 月 14 日起適用的指引。
- WAC 246-215-03440（華盛頓州對應 FDA Food Code 3-403.11 再加熱條文）：165°F／15 秒、2 小時上限、微波需翻攪加蓋並靜置 2 分鐘、商業預製即食品 135°F 等條文的完整文字。
- Engineering ToolBox: Convective Heat Transfer：本文 h 值範圍表的來源（空氣自然對流 5–25、液體自然對流 50–3000、冷凝蒸氣的量級）。（已標存疑：這是工程手冊的通用範圍，涵蓋的是一般工業表面，不是食品幾何的實測值；用來看數量級可以，代進計算之前應該用自己設備的降溫曲線反推真實 h。另外，FDA 已把新版 Food Code 列入 2026 年的優先交付項目，但截至撰稿時我沒能查證到它是否已正式發布；目前可確認的完整版本仍是 2022 Food Code 加 2024 年 11 月的 Supplement。）
💬 問 AI
我要把一個實際的烹調流程從頭到尾用熱傳導與微生物動力學算一遍，
請每一步都寫出公式、代入的數值、以及數值的出處。
我的情境（請用我給的值，不要替換成典型值）：
- 食材與幾何：____（形狀：平板／圓柱／球／立方體，最厚處 ____ mm）
- 起始溫度 ____ °C，水浴／烤箱溫度 ____ °C
- 冷卻方式 ____，冷卻介質溫度 ____ °C
- 之後要冷藏 ____ 天，食用前是否再加熱：____
請依序處理：
1. 參數確認
- 列出你要用的 k、ρ、C_p、α，標明文獻出處與適用溫度範圍。
- 列出升溫段與冷卻段各自的 h 估計值與依據，並說明不確定度。
- 算出 Bi = hR/k，明確說出瓶頸在表面還是在內部。
2. 特徵值
- 依幾何選擇正確的特徵方程（板 λtanλ=Bi；柱 λJ₁/J₀=Bi；球 1−λcotλ=Bi），
解出 λ₁ 與 A₁，並列出你用的求根過程。
- 若是立方體或長方體，請用 Newman 乘積解，不要用分數 β 內插；
若你選擇用內插，請同時給出乘積解並說明兩者差多少。
3. 時間計算（每一步都要檢查 Fo > 0.2）
- 中心到達目標溫度的時間。
- 中心進入與離開 10–54.4°C（50–130°F）的時刻，算出停留總時數。
- 冷卻段：中心走完 54.4→26.7°C 與 26.7→4.4°C 各需多久。
- 任何 Fo < 0.2 的步驟，請明確標示一項級數解不適用，
並改用半無窮大解或數值解，說明你換了方法。
4. 法規對照
- 對照 FSIS Stabilization Guideline (Revised Appendix B) Table 1 的八個選項，
逐一判定通過或不通過，指出是哪一段超時。
- 對照 Revised Appendix A 的 50–130°F 六小時建議，算出我用掉幾成預算。
- 若要再加熱，對照 FDA Food Code 3-403.11，
算出在我的加熱環境下中心到 74°C 需要多久，是否在 2 小時內完成。
5. 敏感度與槓桿
- 分別算出「h 加倍」「厚度減半」「換成球形」三種改動各能省下多少時間，
排出優先順序，並說明為什麼是這個順序（要提到 Bi 落在哪一區）。
- 給出在我這個 Bi 下，h→∞ 的理論最大節省百分比。
要求：
- 不要為了讓結論好看而縮短任何官方時程；不通過就直接說不通過。
- 明確區分「法規要求」「官方指引（safe harbor）」「文獻建議」「你的估算」四種強度。
- 任何條號、版本、文件編號不確定的地方，標註「不確定」，不要編造。
- 不要用「再加熱可以補救」來解決任何一個不通過的環節，
並在最後說明為什麼 log 削減不能相加。