昨天把供給端的模型釘在兩個隨手假設上：縫隙寬 g = 5 微米、傳輸距離 ℓ = 2 毫米，兩個都標了「自行假設,已標存疑」。今天去查文獻,想把這兩個數字換成真正有出處的東西,結果查到的不是修正,是整個模型骨架都要換——昨天那個「兩片平行板間的潤滑流」,原來在複合材料製造工程裡有一個三十多年歷史的正式解,而且套上真正的公式之後,供給量暴增了三到五個數量級,交叉速度 V* 從每秒 0.18 公分,一路推到每秒 26 到 183 公分——不再是遠低於正常運筆速度,而是卡在正常運筆速度的範圍之內或之上。這不是把昨天的結論修得更精確,是把「正常書寫是否全程饑餓」這句話,從「幾乎確定」拉回「要看筆有多濕」。
📖 學
一、「饑餓區」不是這個系列發明的詞,塗佈工程裡它有正式名字和參數
先把好消息說完。查滾筒塗佈(roll coating)文獻,饑餓入口與供墨過剩入口的區分是這個領域的基本功課,正式名稱是 flooded inlet 與 starved（或 meniscus）inlet,用一個無量綱的流量比 λ 描述——λ 大於 1 是過剩供給,λ 小於 1 是饑餓供給,這正是昨天 Q_demand 對 Q_supply 那個比值想表達的東西,只是換了記號。Leeds 大學 Gaskell 等人那篇 meniscus roll coating 的實驗研究把 λ 定義成入口流量與「滾筒速度 × 最小半縫隙」的比值,並指出饑餓入口的上游彎液面會產生一個低於環境壓力的壓力場——這和昨天用毛細吸力 ΔP = 2γ/g 驅動供給的圖像,方向完全一致。
也就是說,昨天用能量與流量平衡重新推出來的東西,不是重新發明輪子,是撞上了一個已經被命名、被參數化過的現象。這件事本身值得記一句:「饑餓區」這個框架是對的方向,錯的是供給那條曲線的具體算法。
二、供給的真正模型：Gebart 1992,同一個潤滑理論,不同的幾何
昨天那句「兩片平行板間的潤滑流」,查了才知道它在纖維複合材料的樹脂灌注成型(RTM, resin transfer molding)工程裡有一個標準解,叫 Gebart 模型(Gebart, B.R., “Permeability of Unidirectional Reinforcements for RTM”, Journal of Composite Materials 26 (1992), 1100–1133)。這個模型處理的問題和毛筆的供墨機制在幾何上幾乎一模一樣:一束平行排列的圓柱纖維,樹脂(對應到這裡就是墨)要滲過纖維間的窄縫,而窄縫的潤滑流正是主導阻力的機制——這句話直接對上昨天的假設,只是昨天是憑空猜的,這裡是三十多年來被大量實驗校驗過的。
Gebart 給出兩個方向的滲透率,關鍵是要分清楚哪一個對應毛筆:
橫向滲透率(K_T,垂直於纖維軸,流體要繞過一根根纖維前進)：K_T = C₁ r² (√(φ_max/φ) − 1)^2.5
縱向滲透率(K_L,沿著纖維軸方向,流體走的是纖維之間的長條縫隙)：K_L = (8r²/c) × (1−φ)³/φ²
r 是纖維半徑,φ 是纖維體積分率(填充率),φ_max 是該排列方式的最密堆積填充率。方形排列 c = 57、φ_max = π/4;六角排列 c = 53、φ_max = π/(2√3) ≈ 0.907,C₁ = 16/(9π√6) ≈ 0.231(方形是 16/(9π√2))。
毛筆裡墨要走的路,是從蓄墨區沿著筆毛的方向送到接觸區——這是沿纖維軸方向的流動,對應的是 K_L,不是 K_T(K_T 描述的是墨要「橫著」穿過一整排纖維,那是紙吸墨的洇散模型,不是筆供墨的模型)。這一點昨天完全沒有分清楚——昨天的兩片平行板模型其實混著算,沒有明確是哪個方向。
三、重算供給:三個填充率,三個天差地遠的答案
取六角排列(c = 53、φ_max = 0.907)、r = 42.5 微米(09-22 釘死的羊毛半徑)。填充率沿用 09-27 第五節 Kuwabara 模型用過的同一組值(0.1、0.3、0.5),這樣兩天的推導可以互相對照。
Gebart 公式裡有一個副產品:(√(φ_max/φ) − 1) 這一項本身就是纖維間隙相對半徑的無量綱寬度,不必再像昨天那樣憑空猜一個 g。用它反推真正的間隙:
gap = r × (√(φ_max/φ) − 1)
| 填充率 φ | 間隙 gap | 驅動壓力 ΔP = 2γ/gap | K_L | Q_supply(ℓ = 2mm、A = πa² 不變) |
|---|---|---|---|---|
| 0.1 | 85.5 μm | 936 Pa | 1.99×10⁻⁸ m² | 1.10×10⁻⁶ m³/s |
| 0.3 | 31.4 μm | 2548 Pa | 1.04×10⁻⁹ m² | 1.56×10⁻⁷ m³/s |
| 0.5 | 14.7 μm | 5428 Pa | 1.36×10⁻¹⁰ m² | 4.37×10⁻⁸ m³/s |
（A 沿用 09-22 δ = 2mm 那一列的接觸半寬 a = 0.613mm,取 A = πa² ≈ 1.18×10⁻⁶ m²,作為墨離開蓄墨區、進入接觸區那個「出口截面」的估計；ℓ = 2mm 沿用昨天的假設,今天只換供給模型,不換傳輸距離,好讓兩天的差異只來自一個變數。）
先看間隙。真正由 Gebart 幾何反推出來的間隙是 14.7 到 85.5 微米,比昨天憑空猜的 5 微米大了 3 到 17 倍——猜錯的方向是把間隙猜得太窄。但即使 ΔP 因此比昨天小(5428 對 16000 帕,以及更小的填充率下更小),Q_supply 依然比昨天暴增,原因不在 ΔP,在於流動的截面：昨天把供給算成一條寬 w、縫隙 g 的窄槽,截面積只有 w×g 這麼小(數量級 6×10⁻⁹ m²);今天用整個蓄墨柱的截面 πa² ≈ 1.18×10⁻⁶ m²,大了兩百倍。這才是暴增的真正來源——昨天的模型讓墨只能從一條縫擠過去,今天的模型讓墨從整根蓄墨柱的橫截面一起滲透,而後者才是纖維束滲流的真實幾何。
對照昨天算出的 Q_supply ≈ 1.03×10⁻¹¹ m³/s(那個數字與速度無關,也與填充率無關,因為根本沒有把填充率放進模型):今天在 φ = 0.5(密實)時是昨天的 4243 倍,φ = 0.3 時是 15,150 倍,φ = 0.1(疏鬆)時是 106,800 倍。
四、交叉點暴增:V* 從 0.18 公分/秒變成 26 到 183 公分/秒
沿用昨天的需求公式 Q_demand ≈ 4.00×10⁻⁷ × V^(5/3)(V 以 m/s 計),解交叉點:
| 填充率 φ | V*(m/s) | V*(cm/s) | 對應 Ca = 0.25V* |
|---|---|---|---|
| 0.1(疏鬆,新開鋒或墨極飽) | 1.83 | 183 | 0.458 |
| 0.3 | 0.569 | 56.9 | 0.142 |
| 0.5(密實,半乾或陳墨黏稠殼層) | 0.265 | 26.5 | 0.0662 |
昨天的 V* ≈ 0.18 公分/秒,今天最保守(填充率最密)的版本也有 26.5 公分/秒,是昨天的 147 倍;最疏鬆的版本是 1017 倍。
這不是誤差範圍內的修正。昨天的結論是「正常運筆速度(幾到幾十公分/秒)全部落在饑餓區,而且遠遠落在裡面」;今天的結論是「饑餓區的交叉點,就卡在正常運筆速度的範圍邊緣或裡面」——**一支墨飽而鬆散的筆,在中等速度下可能根本還沒進入饑餓區;一支半乾而毛已經開始局部黏合(填充率被有效拉高)的筆,饑餓區的門檻低得多,寫兩三下就會進去。*填充率不是一個抽象參數,它就是「這支筆現在有多濕」的物理讀數——09-22 已經算過,墨越少,毛細力把毛黏得越緊,有效填充率上升。也就是說,同一支筆從飽墨到半乾的過程中,自己的 V 會一路下降——這剛好是王鐸涨墨轉渴筆的另一個機制,而且是今天新找到的,和昨天講的「需求隨速度上升」是兩個獨立、方向相同、會疊加的效應。
五、誠實的邊界:新交叉點早就跑出模型自己的有效範圍
上面這張表算得漂亮,但要立刻潑自己冷水。09-29 那篇查到 Bretherton/Landau-Levich 膜厚公式嚴格成立的範圍是 Ca < 0.003(對應 V < 1.2 公分/秒)。今天算出的三個 Ca* 是 0.066、0.142、0.458——分別是那個門檻的 22 倍、47 倍、153 倍。
換句話說:今天算出的「交叉點」,落在理想膜厚公式(h ∝ Ca^(2/3))早已失效的區域。 山峰模型的兩條曲線(先升的理想曲線、後降的饑餓曲線)要嘛在山頂交會才有意義,但如果山頂本身落在理想曲線失效的地方,那個交會點的精確位置就沒有意義,只能當「量級」看。這條存疑必須寫在最前面,而不是塞在附錄——今天所有的 V* 數字都只能說「大約在每秒十幾到接近兩公尺這個量級」,不能再精確到個位數的公分。
好消息是:結構性的結論不受這條存疑影響。不管理想曲線在哪裡失效,需求隨 V^(5/3) 成長、供給大致不隨 V 變化(只隨 φ、ℓ、A 這些幾何量變),這個「兩條曲線形狀不同、必然交叉」的骨架是穩的。變得不確定的只是交叉點的精確座標。
六、填充率是一個真實存在的工程旋鈕:一支馬克筆的專利當證據
如果「筆有多濕決定饑餓門檻在哪裡」這件事聽起來像是為了讓模型自洽而編出來的,查一下現代筆類工程就知道這不是憑空想像。氈頭筆(felt-tip pen/marker)的墨水儲存結構,幾十年來一直是壓實的纖維束,靠纖維間的毛細通道儲墨釋墨,美國專利 US7470079B2描述的正是這種設計:主蓄墨區與輔助蓄墨區刻意做成不同的纖維密度(範例給出主蓄墨區約 0.234 g/cm³、輔助區約 0.20 到 0.23 g/cm³),密度差造成的毛細力差,用來主動驅動墨水從一區流向另一區、維持穩定供墨。
這給了今天的模型一個獨立於書法史料之外的驗證:用纖維填充密度來控制毛細供給速率,是一個被工業界實際採用、而且刻意調校的設計手段,不是這個系列自己想像出來的類比。毛筆製作裡「兼毫」(狼毫與羊毫混紮,兩種毛的粗細與硬度不同,等效填充率也不同)其實做的是同一件事,只是沒有人用「填充率控制饑餓門檻」這樣的物理語言描述過。
七、還沒算進去、但方向明確的一條:傳輸距離 ℓ 會不會隨著蘸墨後的行程增加
今天固定 ℓ = 2mm 沒有動,只換了供給的幾何模型。但有一個結構性的問題今天沒有時間展開:一次蘸墨之後,墨最初濕潤整根筆頭,隨著書寫消耗,濕潤區的邊界應該會往筆根方向退縮,墨要從蓄墨區走到接觸區的距離 ℓ 因此會隨著已經寫掉的字數增加而變長。Q_supply ∝ 1/ℓ(從 Darcy’s law 直接讀出來),ℓ 變長,供給進一步下降——這是一個和「速度上升」完全獨立、但方向相同的機制:**同一次蘸墨,寫得越久,即使速度不變,也會越寫越乾。**這比昨天講的「加速轉乾」涵蓋的範圍更廣,可能才是王鐸「濃淡燥潤相映」裡「不需要換速度也會轉乾」那部分的真正機制。今天沒有把 ℓ(字數)這個變數量化,留著。
八、依然查不到的三件事
誠實列出來,不要假裝解決了。第一,飛白斷裂門檻 Ca ≈ 3 這個數字的原始出處,這幾天查了幾次,沒有找到蔡邕堊帚飛白與毛細數 3 這個具體對應的一手文獻,只能繼續當作系列內部的推論,不是外部驗證過的數字。第二,渴筆與飛白是否為兩種不同物理機制,書畫史或紙墨科學文獻裡沒有查到明確區分的研究。第三,直接用高速攝影或流變儀測量真實毛筆書寫時供墨過程的實驗研究——查到的都是機器人書法的軌跡視覺辨識、可微分書寫模型這類計算機圖形/視覺研究,沒有查到針對墨水在筆毫間實際流動速率的物理量測。這三條留著,不編。
🧠 記
**一、「饑餓區」在塗佈工程裡有正式名字:starved / meniscus inlet,用流量比 λ 參數化。**這驗證了昨天的框架方向對,不是自己發明的概念。
**二、供給模型從「兩片平行板」換成 Gebart 1992 的纖維束滲透率公式,沿纖維軸方向要用縱向滲透率 K_L = (8r²/c)(1−φ)³/φ²,不是橫向的 K_T。**這是複合材料樹脂灌注成型工程裡驗證過三十多年的標準解,幾何上與毛筆供墨幾乎一致。
**三、換模型之後,Q_supply 暴增三到五個數量級。**填充率 φ = 0.5 時是昨天估計的 4243 倍,φ = 0.3 時 15,150 倍,φ = 0.1 時 106,800 倍——主因不是驅動壓力變大(反而變小了),是流動截面從「一條窄縫」換成「整根蓄墨柱的橫截面」,大了約兩百倍。
四、交叉速度 V 從每秒 0.18 公分,變成每秒 26.5 到 183 公分,取決於筆的填充率。這個範圍卡在正常運筆速度之內或之上,不再是遠低於它。「正常書寫是否全程饑餓」不再是幾乎確定的結論,要看筆有多濕。*
**五、新交叉點的 Ca 值(0.066 到 0.458)比昨天算出的模型有效邊界(Ca < 0.003)高了 22 到 153 倍。**山峰模型的骨架(需求冪次上升、供給大致平坦、必有交叉)仍然穩,但交叉點的精確位置只能當量級看,不能當精確值用。
**六、填充率是一個真實的工程旋鈕,有專利為證。**美國專利 US7470079B2 描述的氈頭筆蓄墨結構,刻意用不同纖維密度(0.234 對 0.20–0.23 g/cm³)控制毛細供墨速率——這不是這個系列發明的類比,是被工業界實際採用的設計手段。
**七、還沒算的一條:傳輸距離 ℓ 可能隨一次蘸墨後的行程增加而變長,Q_supply ∝ 1/ℓ,這會讓「越寫越乾」在速度不變時也發生。**這條方向明確但今天沒有量化。
**八、三件事依然查不到:飛白 Ca≈3 的原始出處、渴筆與飛白是否為不同機制的文獻區分、任何直接量測毛筆書寫供墨速率的高速攝影或流變儀研究。**繼續標記存疑,不編造。
✍️ 實踐
今天的推導改變很大,但可驗證的落點很具體:同一支筆,飽墨與半乾兩種狀態下,饑餓門檻應該明顯不同。
器材:一支中楷或小楷、生宣、墨、手機(能拍慢動作更好)、尺。
**第一步(6 分鐘):飽墨狀態下找饑餓門檻。*飽蘸墨,固定一個能穩定重複的壓深,依序用大約 5、15、30、60、100 公分/秒的速度(用「一秒內畫完多長」來估算,不必精準)各畫一條 10 公分的橫。預測:如果今天的模型方向對,飽墨(填充率低,供給旺盛,V 較高)時,前面幾個速度墨色應該相對穩定,只有最快的一兩條才明顯轉乾。
**第二步(6 分鐘):同一支筆,先寫到接近半乾,再重複同一組速度。*故意先把筆寫到明顯感覺變澀(這代表墨已經不多,填充率被有效拉高、V 下降),再用同一組速度重複第一步的橫線實驗。預測:半乾狀態下,轉乾出現得更早——中等速度就已經開始出現渴筆絲狀紋理,不需要等到最快那一條。
第三步(8 分鐘):單次蘸墨、固定速度,畫到底,只看行程長度。飽蘸墨,選一個中等偏慢、你確定不會單靠速度就轉乾的速度,持續畫一條盡量長的線(可以繞著寫、或分成好幾段接力,重點是不重新蘸墨、速度盡量恆定)。預測:即使速度全程不變,墨色也會隨著行程增加逐漸變淡,最後出現渴筆——這是在驗證第七節「ℓ 隨行程增加」這個假設。如果速度恆定、墨色卻不隨行程惡化,那條假設就不成立。
**第四步(5 分鐘):比較記錄。**把六到八段線索拍照,不當場判斷濃淡,晚一點或找人排序比較,避免期望造成的偏誤(09-29 已經強調過這一點,今天延用)。
🔗 延伸學習
- Gaskell, P.H. et al., “An experimental investigation of meniscus roll coating”（Leeds 大學 eprints,開放取用 PDF） — 第一節「饑餓入口 / 過剩入口」λ 參數的出處,塗佈工程對這個系列自己重新發明的「饑餓區」給出了正式名稱與量化框架。
- Gebart, B.R., “Permeability of Unidirectional Reinforcements for RTM”, Journal of Composite Materials 26 (1992) — 第二、三節整段推導的公式來源,樹脂灌注成型工程裡纖維束滲透率的標準解,今天用它取代昨天憑空假設的「兩片平行板」模型。摘要頁可讀,全文需訂閱。
- US Patent 7470079B2, “Writing article having a capillary reservoir with…”（Google Patents,開放取用全文） — 第六節填充率作為真實工程旋鈕的證據,氈頭筆蓄墨結構用不同纖維密度(0.234 對 0.20–0.23 g/cm³)主動控制毛細供墨速率,給出具體數字。
- Jackson, G.W. & James, D.F., “The permeability of fibrous porous media”, Canadian Journal of Chemical Engineering 64 (1986) — 纖維多孔介質滲透率這個領域最常被引用的回顧文獻之一,涵蓋從高分子鏈到玻璃纖維的大範圍實驗數據,可以用來檢查 Gebart 模型在毛筆這種相對疏鬆填充率(φ 0.1–0.5)下是否仍然可靠——今天沒有時間查證這一點,留給下一步。
💬 問 AI
我把毛筆的供墨模型從「兩片平行板間的潤滑流」(自己假設的)換成 Gebart 1992 的纖維束滲透率公式(複合材料 RTM 工程的標準解),想請你逐項檢查這個替換是否合法。請嚴格區分「有文獻或實測支持」與「你自己的推論」,每一項標明把握程度(高／中／低);不確定就明說,不要為了答得流暢而編造數據、公式或文獻。
【A 組:模型替換是否合法】
1. 我認為墨從蓄墨區送到接觸區是「沿纖維軸方向」的流動,所以該用 Gebart 的縱向滲透率 K_L = (8r²/c)(1−φ)³/φ²,而不是橫向滲透率 K_T = C₁r²(√(φ_max/φ)−1)^2.5。
（a）這個方向判斷對嗎?毛筆筆毫的排列、以及墨實際的流動路徑,真的是「平行於纖維軸」這麼乾淨嗎?會不會其實是橫向與縱向混合的三維滲流,兩個滲透率都低估或高估了真實值?
（b）Gebart 模型是為連續纖維複合材料(纖維幾乎完美平行、長度遠超過流動距離)設計的。毛筆筆毫是漸細的錐形排列、長度只有 1.6 公分,而且靠近筆尖的地方纖維密度和方向都會變化。這個模型套在這種幾何上,誤差可能有多大?
【B 組:數字本身】
2. 我用 φ = 0.1、0.3、0.5 三個填充率(沿用 09-27 Kuwabara 模型用過的同一組值),算出間隙 gap = r(√(φ_max/φ)−1) 分別是 85.5、31.4、14.7 微米,再用 ΔP = 2γ/gap 當驅動壓力。
（a）用 Gebart 公式內部那個幾何項反推「間隙」再代回拉普拉斯壓力公式,這個做法合理嗎?還是這是把兩個不同尺度的東西硬湊在一起?
（b）我把出口截面積取成 A = πa²(a 是 09-22 算出的接觸半寬,δ=2mm 時 a=0.613mm)。這個截面選得對嗎?真實的「供墨出口」應該多大?
3. 換了模型之後,Q_supply 從昨天的 1.03×10⁻¹¹ m³/s 暴增到今天的 4.37×10⁻⁸ 到 1.10×10⁻⁶ m³/s(三到五個數量級),交叉速度 V* 從每秒 0.18 公分變成每秒 26.5 到 183 公分。
（a）這個暴增幅度合理嗎?會不會是我在套用 Gebart 公式時單位或幾何因子出錯,才會差這麼多個數量級?
（b）如果暴增是真的(不是算錯),它對「正常書寫是否處於饑餓區」這個結論的衝擊有多大?我的說法是「不再是幾乎確定,要看筆有多濕」,這個修正合理嗎,還是我又走過頭了?
【C 組:自己標出的存疑】
4. 我發現新算出的交叉點 Ca 值(0.066 到 0.458)比昨天查到的 Bretherton 公式有效邊界(Ca < 0.003)高了 22 到 153 倍,並主張「山峰模型的骨架還穩,但交叉點精確位置只能當量級看」。
（a）這個補救說法站得住嗎?還是說,一旦交叉點落在理想曲線失效的地方,連「存在一個先升後降的山峰」這個定性結論都動搖了?
（b）在 Ca 這麼大的範圍,有沒有別的、真正適用的膜厚理論(不是 Bretherton 的小 Ca 近似)可以用來檢查我的結論方向?
【D 組:新開的線索】
5. 我提出一個沒有量化的假設:同一次蘸墨之後,濕潤區邊界隨書寫消耗往筆根退縮,傳輸距離 ℓ 因此隨行程增加而變長,Q_supply ∝ 1/ℓ,所以「速度不變也會越寫越乾」。
（a）這個物理圖像合理嗎?濕潤區退縮的機制,和多孔介質乾燥(drying front)的標準理論是同一件事嗎?
（b）如果要把這個效應量化,最少需要哪些額外的假設或量測?
6. 我用一個氈頭筆蓄墨的美國專利(US7470079B2,兩個蓄墨區纖維密度 0.234 對 0.20–0.23 g/cm³)當作「填充率是真實工程旋鈕」的獨立證據。
（a）這個類比有多牽強?氈頭筆的蓄墨結構和毛筆筆毫的物理有沒有根本性的差異(例如材質是合成纖維不是角蛋白、通道更規則)?
（b）現代書畫用品或製筆工藝文獻裡,有沒有更直接討論「兼毫」(不同獸毛混紮)如何影響蓄墨與供墨速率的資料?
【E 組:反駁與取捨】
7. 請站到反面,給出你能想到的、對「換成 Gebart 模型後供給暴增,饑餓門檻大幅上修」這個結論的最強反駁。特別攻擊——
（a）Gebart 模型假設纖維是規則排列(方形或六角),真實毛筆筆毫的排列是隨機的、還帶著毛細力造成的局部聚束(09-22 算過的協同硬化)。隨機排列會不會讓真實滲透率遠低於規則排列的預測,把結論打回昨天那個量級?
（b）我固定截面積 A 不隨速度變化,但真實書寫時筆頭形變、接觸區大小都在動態變化,這個簡化會不會系統性地高估或低估供給?
8. 最後請把上面所有條目分成三堆:(一)我可以用手機拍照和簡單器材自己驗證的;(二)必須查一手文獻的;(三)目前沒人知道、需要真的做實驗(顯微鏡、流變儀、或纖維滲透率量測儀)才能回答的。