昨天留下的限制清單裡有一條最容易被忽略:算供給時我把驅動壓力 ΔP 固定在 09-30 的值,而 Larson & Zok 2018 全文裡剛好有一條經過實測驗證的軸向毛細壓力公式,可以拿來對帳。對完之後的結果是,09-30 的 ΔP 比文獻公式高 2.9 到 4.5 倍,把它改掉,交叉速度 V* 會縮到原來的 0.41 至 0.53 倍(φ = 0.5 時由 26.5 降到 14.0 公分/秒);再把昨天的叢集效應(K_L 升至 1.67 至 3.87 倍)疊上去,V* 回到 19 至 32 公分/秒。兩個修正方向相反、大小相近,淨結果是「V* 在十幾到三十幾公分/秒」這個量級沒有動,但 09-30 的 26.5 是靠兩個沒對帳的假設湊出來的。今天另外確認了昨天的方向疑點(Chen & Papathanasiou 的摘要明寫軸向滲透率隨不均勻而上升),並把毛筆內部的毛細數對到 Larson & Zok 的實驗範圍:一般運筆速度落在他們資料的邊緣或以外,不飽和折減只出現在很慢的筆速,而那時供給遠大於需求,折減無傷。
📖 學
一、昨天兩件懸案:方向定案,Seuffert 的公式仍沒拿到
第一件是 Chen & Papathanasiou 2007 的方向。昨天只能從 Larson & Zok 的引言轉引,今天另外查到該篇在 Composites Science and Technology 第 67 卷第 7–8 期的摘要轉述(搜尋工具回傳的摘要內容;ScienceDirect 原頁與 ResearchGate 頁我這邊都回 403,讀不到全文):研究用約 600 根隨機放置的纖維截面、以邊界元法解軸向 Stokes 流,結果是軸向滲透率隨微結構由均勻變成不均勻而上升,並提出「平均最近纖維間距」當不均勻度的參數。這與 Larson & Zok 的引言方向相同,昨天二手摘要裡「下降」的說法,現在我判為誤轉述,仍屬推論,因為我沒讀到原文。Shou, Ye & Fan 2015 的摘要(SAGE 頁)用 Voronoi 鑲嵌量化纖維分布的隨機度、以單一參數表示,但摘要沒有明寫方向,方向仍然只是 Larson & Zok 的轉述。
第二件是 Seuffert 2021 圖 15 的定義。今天重讀全文網頁,確認原文是「normalized difference between the averages of the SVE and MVE simulations … calculated with a normalization using the geometric mean of the values」,但公式本體在我用的網頁抽取工具裡沒有出現,所以昨天那個「(MVE − SVE)/√(MVE × SVE) 的話最大值 2 對應 5.8 倍」仍是我猜公式後的換算,不能當論文數字,維持昨天的標示。倒是有兩句新的原文值得記:軸向流動「depends on the flow inside large gaps, whereas the flow between very close fibers can be neglected」,這與上一篇的凸函數論證是同一個方向(密處貢獻小、疏處貢獻大);而橫向的結果正好相反,「the transverse permeability of the SVEs is higher than the one predicted using MVEs」,真實切片的橫向滲透率比人工隨機低。對毛筆的意義是推論:墨沿軸向送到筆尖受惠於聚束,但墨從筆心橫向擴散到外圈毛,在真實的叢集筆毫裡反而更慢。
二、對帳毛細壓力:09-30 的 ΔP 高了 2.9 到 4.5 倍
我今天把 Larson & Zok 的 PDF 文字層抽出來讀(公式符號在抽取後有錯位,下列公式我都用文中給的數值重算驗證過)。他們量的是 SiC 纖維(平均半徑 6.3 微米)束,孔隙率 0.32–0.38,流體黏度 68–82 mPa·s、表面張力 30 mN/m,接觸角實測 26.3 ± 8.2°。軸向毛細壓力用這個式子:
P_c = F γ cosθ (1−ε) / (2 ε r̄_f),軸向流動 F = 4(文中引自前人)
代入 ε = 0.35、cosθ = cos 26.3° = 0.897,算出 15.9 千帕,與他們另行報告的平均毛細壓力 16 ± 1 千帕吻合,所以這條式子的數字沒有抄錯,也確實是實測量級。我自己再推一次它的來歷:床的水力半徑 r_h = ε r / (2(1−ε)),毛細壓力 = γ cosθ / r_h,代進去正好就是 F = 4。換句話說,文獻用的是整個床的平均水力半徑,而 09-30 用的是 gap = r(√(φ_max/φ) − 1) 這個最窄間隙當半縫(ΔP = 2γ/gap,等於把水力半徑取成 gap/2)。兩者差一個「最窄通道」與「平均通道」。
毛筆參數沿用 09-30(r = 42.5 微米、γ = 0.04 牛頓/公尺、接觸角取 0 度即 cosθ = 1,這是上界):
| 填充率 φ | 09-30 的 ΔP(最窄間隙) | 文獻式 P_c(平均水力半徑) | 倍數 | V*(09-30) | V*(改用文獻式) |
|---|---|---|---|---|---|
| 0.1 | 936 Pa | 209 Pa | 4.47 | 183 cm/s | 74.6 cm/s |
| 0.3 | 2548 Pa | 807 Pa | 3.16 | 56.9 cm/s | 28.5 cm/s |
| 0.5 | 5428 Pa | 1882 Pa | 2.88 | 26.5 cm/s | 14.0 cm/s |
V* 欄的算法沿用 09-30 的關係:Q_supply = K_L · ΔP · A / (μ ℓ),Q_demand = 4.00×10⁻⁷ V^(5/3),所以 V* ∝ (K_L ΔP)^(3/5)。φ = 0.5 時:ΔP 比值 2.88,2.88^0.6 = 1.89,26.5 ÷ 1.89 = 14.0。同一支筆的 K_L 一個都沒動(1.36×10⁻¹⁰ m²),變的只有驅動壓力。
這不等於 09-30 的數字「錯了」。Larson & Zok 同時報告:在整個局部孔隙率範圍內,局部毛細壓力落在 6 至 40 千帕,而且液體最先出現在最窄的通道(局部毛細壓力最高處)。所以最窄間隙的值是真實存在的局部吸力,文獻式是床平均吸力;供給整體該用哪一個,取決於墨是走窄通道還是走全部通道,我沒有查到對毛筆的直接答案。保守的讀法是:09-30 的 V* 屬上界,文獻式的 V* 屬較合理的中心值。
三、三個修正疊起來:誤差預算
昨天叢集效應的 V* 欄(36.1、46.2、59.7)是「K_L 變大、ΔP 不變」算的。把基準換成文獻式 P_c 後重算,φ = 0.5 的基準 V* 為 14.0,叢集因子 V* ∝ K^0.6:
| 情境(平均 φ = 0.5) | K_L 倍數 | V*(cm/s) |
|---|---|---|
| 09-30 原模型 | 1.00(ΔP 取最窄間隙) | 26.5 |
| 改用文獻 P_c,均勻分布 | 1.00 | 14.0 |
| 30% 截面叢集成 0.7 | 1.67 | 14.0 × 1.67^0.6 = 19.1 |
| 30% 截面叢集成 0.8 | 2.53 | 14.0 × 2.53^0.6 = 24.5 |
| 50% 截面叢集成 0.7 | 3.87 | 14.0 × 3.87^0.6 = 31.6 |
結論有兩層。數量級穩:不論怎麼修,V* 在十幾到三十幾公分/秒,也就是仍然在正常運筆速度的範圍邊緣,「正常書寫是否饑餓要看筆有多濕」這句話不需要收回。精確數字不穩:任何一個單一修正都能讓 V* 變動將近一倍,三個修正互相抵消只是巧合,不是證明模型準。兩個尚未處理的方向在表外:叢集區的毛細壓力也會變(局部更密,吸力更大,與整體攤平的假設不相容),以及叢集讓出口截面 A 縮小(昨天已標),我都沒算進去。
四、內部毛細數對上 Larson & Zok 的實驗範圍
Larson & Zok 的主要結果之一是不飽和滲透率隨毛細數 Ca 上升:對數對數迴歸斜率 B = 0.37 ± 0.04(p = 1.4×10⁻⁹),樣本為 22 支試片的 31 個位置,孔隙率 0.32–0.38,Ca_ff 範圍約 10⁻⁵ 到 10⁻³。他們文字說這個範圍內滲透率上升「nearly an order of magnitude」,但迴歸斜率換算出的結果是 100^0.37 = 5.5 倍,兩者不完全吻合,我沒找到文中的解釋,猜是端點與散布的差別,這是我的推論。實驗分三種狀況:無外加壓力、純毛細吸入(Ca_ff ≈ 4×10⁻⁵),穩態飽和度只有 0.81,按他們的 GPE 方法算出的 κ_u/S 比飽和值低約 30%;中等壓力 138 千帕(Ca_ff ≈ 3×10⁻⁴),飽和度 0.93,滲透率只低約 8%(0.0214 對 0.0196,無因次值);高壓(超過 276 千帕,Ca_ff ≥ 7×10⁻⁴),飽和度超過 0.99,兩者幾乎相同。機制是慢流時液體偏走最窄通道,窄通道的流量低於平均,所以整體滲透率偏低。
把毛筆內部的流動換成同一個毛細數:用通過 A = 1.18×10⁻⁶ 平方公尺的達西流速 q = Q_demand / A,示蹤流速 ν = q/ε,Ca = μν/γ(μ = 0.01 Pa·s、γ = 0.04 N/m,與 09-29 一致),Q_demand 取 09-30 的需求曲線。供給/需求比用文獻式 P_c 算,以 φ = 0.5 為例:
| 運筆速度 | Q_demand | 內部 Ca(φ = 0.3 / 0.5) | 對應 Larson & Zok 的區域 | 供給/需求(φ = 0.5) |
|---|---|---|---|---|
| 1 cm/s | 1.86×10⁻¹⁰ m³/s | 5.6×10⁻⁵ / 7.9×10⁻⁵ | 近純毛細區,κ_u/S 約 −30% | 82 |
| 5 cm/s | 2.72×10⁻⁹ m³/s | 8.2×10⁻⁴ / 1.15×10⁻³ | 近飽和區 | — |
| 10 cm/s | 8.62×10⁻⁹ m³/s | 2.6×10⁻³ / 3.6×10⁻³ | 超出他們資料上緣(2×10⁻³) | 1.76 |
| 20 cm/s | 2.73×10⁻⁸ m³/s | 8.3×10⁻³ / 1.16×10⁻² | 超出資料約 4 至 6 倍 | — |
(10 公分/秒那一列:φ = 0.5 時供給 1.51×10⁻⁸ ÷ 需求 8.62×10⁻⁹ = 1.76。)
兩個含意。第一,慢筆(約 4 公分/秒以下,也就是 Ca 低於 7×10⁻⁴;小楷的起收筆、行筆很慢處)的內部流速落在他們發現「偏走窄通道、滲透率降低」的範圍,最明顯的折減約三成出現在 1 公分/秒以下。這個折減若套到 V* 上只讓它變成 0.7^0.6 = 0.81 倍,而且慢筆時供給遠大於需求(1 公分/秒時 82 倍,4 公分/秒時約 8 倍),折減不傷。第二,V* 本身(Ca 在 0.006 到 0.06)離他們的實驗範圍至少三倍以上,在 V 附近滲透率對 Ca 的行為沒有任何直接資料*,我只能說高壓端飽和度已趨近 1,往更高 Ca 外推沒有看到折減的理由,但那是外推。
五、類比的邊界
這個對照不是等價的,先列清楚。流體:他們的前驅物聚合物黏度 68–82 mPa·s、γ = 30 mN/m;墨汁假設 10 mPa·s、40 mN/m,墨還是含碳黑顆粒的懸浮液。纖維:半徑 6.3 微米、孔隙率 0.32–0.38(填充率 0.62–0.68),毛筆半徑 42.5 微米、填充率 0.1–0.5,兩者的填充率範圍沒有重疊;GPE 的校準也只對六角排列 ε = 0.2–0.7,毛筆的疏鬆端(ε 到 0.9)在校準範圍外。過程:他們量的是液體侵入乾床的前緣;毛筆蓄墨區是預先濕透的,不飽和只會出現在乾掉的筆尖,與寫到墨乾時的情況相近但不是同一件事。此外他們觀察到毛細力把密處的纖維拉得更近(彈毛細效應),使床變得更不均勻,並預期這個效應在「較低纖維體積分率」更大(原文為預期,不是量測結果),這與 09-22 的毛細聚束同向,但在毛筆上沒有直接測量。
🧠 記
**一、Chen & Papathanasiou、Larson & Zok、Seuffert 三條路都指向同一方向:軸向滲透率隨纖維分布不均勻而上升;橫向則相反。**昨天二手摘要的「下降」判為誤轉述,但我仍沒讀到 Chen 與 Shou 的全文,所以是轉述層級的結論。對毛筆,軸向供墨受惠於聚束,橫向把墨擴散到外圈毛則可能受損(推論)。
**二、09-30 的 ΔP 偏高 2.9 到 4.5 倍。文獻的軸向毛細壓力式 P_c = Fγcosθ(1−ε)/(2εr)(F = 4)在 Larson & Zok 的材料上算出 15.9 千帕,與實測 16 ± 1 千帕吻合;φ = 0.5 的毛筆改用它之後 V 由 26.5 降到 14.0 公分/秒。疊上叢集效應(K_L 升 1.67 至 3.87 倍)後回到 19 至 32,所以「V 在十幾到三十幾公分/秒」的量級穩,精確值不穩。
*三、不飽和折減(約 30%)只出現在慢筆(Ca 約 10⁻⁵ 至 10⁻⁴),那時供給至少是需求的八倍,對書寫沒有影響;V 附近比他們的資料上限高三倍以上,沒有直接實驗支持,只能外推。**邊界:墨汁與毛筆的填充率、黏度、纖維直徑都和他們的材料不同,對照只能當量級參考。
✍️ 實踐
今天的模型有一個可以直接量的預測:毛筆從乾的狀態吸墨時,前緣位置的平方正比於時間(達西吸入),斜率 D = 2κP_c / (με)。Larson & Zok 的式 5 就是這個形式(無外加壓力,Pm = 0)。二十分鐘內做完。
預測值(cosθ = 1、K 與 P_c 用上面的文獻式;因此是上界):
| 填充率 φ | D(m²/s) | 前緣到 5 mm | 前緣到 16 mm |
|---|---|---|---|
| 0.3 | 2.4×10⁻⁴ | 0.10 秒 | 1.1 秒 |
| 0.5 | 1.0×10⁻⁴ | 0.24 秒 | 2.5 秒 |
(φ = 0.1 時重力與毛細壓力同量級,不列。16 毫米處的重力在 φ = 0.5 約佔 P_c 的 8%,筆桿須水平,避免重力干擾。)
器材:洗淨並完全晾乾的小楷羊毫(已用過、膠已退掉的筆,新筆膠硬不適用)、玻璃片或光滑塑膠片、墨汁約 0.1 毫升(兩滴)、手機慢動作(120 或 240 幀/秒)、尺貼在筆旁。
**第一步(5 分鐘):架設。**筆桿水平固定在桌邊(膠帶黏住即可),筆尖輕輕貼在玻璃上的一滴墨,不要壓。尺沿筆毫貼好,從上方拍,畫面包含尺與筆毫。羊毫偏白,墨前緣在毛上看得出來。
**第二步(10 分鐘):拍三次。**每次換一個乾的位置(或把筆洗淨晾乾),只拍墨從筆尖爬上去的前十秒。逐幀看,記五個時間點的前緣位置 x。
第三步(5 分鐘):算斜率。把 x² 對 t 畫成圖,看是否近似一條直線(達西行為),斜率就是 D。與上表比較:實測 D 比預測小多少倍。實測偏慢可以有三個來源,這個實驗分不開:接觸角大於零(cosθ < 1)、有效填充率比假設的高(筆毫黏緊)、墨裡碳黑顆粒落後液體前緣(看到的是顏色前緣,低估液體前緣)。如果實測 D 比預測小超過二十倍,代表我的 K_L 或 P_c 至少有一個錯得離譜,不能再當作量級使用。
**明天接:**拿今天量到的 D 反推「有效 cosθ × K」的乘積,與 09-30 的最窄間隙版本比較,看實測落在文獻式與最窄間隙式之間的哪一側;並去查 Um & Lee 1997 與 Yang & Jin 2014 兩篇的原文,今天只在 Larson & Zok 的引用裡看到名字。
🔗 延伸學習
- Larson, N.M. & Zok, F.W., “Insights from in-situ X-ray computed tomography during axial impregnation of unidirectional fiber beds”, Composites Part A 107 (2018)(OSTI 開放全文) — 今天第二至五節的主要依據,我讀的是 PDF 文字層(公式符號有錯位,已用數值重算);毛細壓力式、Ca 與不飽和滲透率關係、GPE 校準都在這篇。
- Seuffert, J. et al., “Micro-Scale Permeability Characterization of Carbon Fiber Composites Using Micrograph Volume Elements”, Frontiers in Materials 8 (2021) — 第一節,軸向與橫向的方向相反結果與「大間隙主導軸向流動」這句,圖 15 的公式本體今天仍沒讀到。
- Chen, X. & Papathanasiou, T.D., “Micro-scale modeling of axial flow through unidirectional disordered fiber arrays”, Composites Science and Technology 67 (2007)(ScienceDirect 摘要頁) — 第一節方向定案所依的摘要轉述,我這邊該頁回 403,全文未讀。
- Shou, D., Ye, L. & Fan, J., “On the longitudinal permeability of aligned fiber arrays”, Journal of Composite Materials 49 (2015)(SAGE 摘要頁) — 用 Voronoi 量化隨機度的單元模型,摘要未寫方向;是 Larson & Zok 評比的兩個局部公式之一(昨天說的 Voronoi 配 Shou)。
💬 問 AI
我在檢查一個毛筆供墨模型的兩處修正,請逐項回答,並嚴格區分三類:(一)文獻明確寫了的、(二)你能查證或憑訓練資料確定的、(三)純屬推論的,每項標明把握程度(高/中/低);不確定或查不到就直說,不要為了流暢而編造文獻、公式或數字。
【A 組:毛細壓力公式】
1. Larson & Zok 2018(Composites Part A 107)用軸向毛細壓力 P_c = F γ cosθ (1−ε) / (2 ε r),F = 4。我自己推導的來歷是:床的水力半徑 r_h = ε r / (2(1−ε)),P_c = γ cosθ / r_h,正好對上 F = 4。
(a) 這個推導與 F = 4 的來歷對嗎?文獻的 F 在什麼條件下會不是 4?
(b) 我之前用「最窄間隙 gap = r(√(φ_max/φ)−1)」的半縫當水力半徑,得到的 ΔP 比平均水力半徑高 2.9 到 4.5 倍。在「墨從蓄墨區流向紙上接觸區」的供給問題裡,該用最窄通道的吸力、床平均吸力,還是別的?為什麼?
(c) 我取 cosθ = 1(上界)。羊毫(角蛋白纖維)與墨汁的接觸角大概是多少?你有沒有可靠數據?
【B 組:誤差預算】
2. 我把三件事疊在一起:毛細壓力改用文獻式(V* 變 0.4 至 0.53 倍)、叢集使 K_L 升 1.67 至 3.87 倍(V* ∝ K^0.6,升 1.36 至 2.25 倍)、不飽和折減約 30%(V* 乘 0.81)。請檢查:
(a) V* ∝ (K_L · ΔP)^(3/5) 這個關係成立的前提是什麼?需求端 Q_demand ∝ V^(5/3) 在 Ca 遠大於 0.003 時還適用嗎?
(b) 叢集區域的局部毛細壓力也會改變,我沒有算進去。方向和量級大概是什麼?
(c) 這三項修正互相抵消,是巧合還是有結構性原因?
【C 組:不飽和滲透率的適用性】
3. Larson & Zok 量到的 κ_u/S 隨 Ca 上升,對數迴歸斜率 0.37 ± 0.04,Ca 約 10⁻⁵ 到 10⁻³,孔隙率 0.32–0.38。文字說「nearly an order of magnitude」,我換算斜率只有 5.5 倍。這個差距你怎麼解釋?
4. 我把毛筆內部的流動換成同一個毛細數(Ca = μ q / (ε γ),q 取需求曲線除以出口截面),結論是一般運筆速度(5–20 公分/秒)落在他們實驗資料的邊緣或以外。這個換算合理嗎?他們的 Ca 是「流動前緣速度」定義,和我的「出口達西流速」定義是同一種量嗎?
5. 他們研究的是液體侵入乾床,毛筆蓄墨區是預先濕透的。不飽和滲透率的結論,在毛筆「寫到墨乾」時適用到什麼程度?
【D 組:反駁與驗證】
6. 請站到反面,攻擊我最核心的論點:「V* 在十幾到三十幾公分/秒這個量級是穩的」。證據是一篇讀到全文的論文、一個自己重算驗證過的公式、以及三個互相抵消的修正。
7. 我設計了一個吸墨實驗:乾筆水平,筆尖碰一滴墨,慢動作拍前緣位置 x 對時間 t,看 x² 對 t 是否為直線,斜率 D = 2κP_c/(με)。請指出這個實驗最容易出錯的地方,以及碳黑顆粒落後液體前緣會造成多大誤差(若有文獻請標明)。
8. 最後請把所有存疑項目分成三堆:(一)你憑訓練資料就能判斷的;(二)必須查一手文獻的(請列出你認為該先查哪幾篇);(三)沒人做過、需要用毛筆實驗才能回答的。