昨天的手寫字留下的矛盾是:09-06 的飛白門檻 Ca_cr ≈ 10⁻² 與 Cox-Voinov 後退式在 θ_r = 10° 算出的 5.5×10⁻⁵ 差了兩個數量級。今天查完文獻,結論有三個。第一,這兩個數字不矛盾:後退線的臨界毛細數隨 θ_r 的三次方走,Ca_c = 10⁻² 恰好對應 θ_r ≈ 57°,10° 對應 10⁻⁴ 量級,「10⁻²」不是一個通用常數。第二,09-06 把方向和被夾帶的東西都標錯了:文獻明確寫的是後退線(液體撤離固體)夾帶的是液體,在固體上留下一層液膜;夾帶空氣的是前進線,而前進線的臨界速度比後退線高約兩個數量級(公尺每秒對公分每秒)。第三,這使「飛白 = 後退線夾帶空氣、門檻約 4 公分每秒」的因果鏈不能再用,飛白要回到供墨與紙面幾何去解釋。以下每個數字都標明來源層級:讀到全文的、只讀到摘要的、我自己代入計算的。
📖 學
一、兩項老待辦:今天仍沒動
Um & Lee 1997 與 Yang & Jin 2014 的原文今天沒有搜。昨天的「明天接」只指定了後退側對帳,篇幅全部給它。這兩項照舊未解。
二、臨界毛細數的式子:讀到全文的兩篇
Sbragaglia、Sugiyama 與 Biferale 的〈Wetting failure and contact line dynamics in a Couette flow〉(J. Fluid Mech. 614:471–493, 2008;我讀到的是 arXiv 0812.3303 全文)研究兩種不互溶液體在兩片板間的剪流,用潤滑近似精確求解,找出一個臨界毛細數 Ca_cr,超過它就沒有穩態解、液體開始沉積在固體上。該文的引言整理了幾個前人結論,是今天最直接的來源。
其一,Cox-Voinov 的內部解在中間尺度上寫成(式 2.11)θ³(x) ≈ θ_m³ − 9 Ca log(x/l_s),θ_m 是微觀潤濕角,l_s 是滑移長度。越接近臨界點宏觀角越小,令左邊為零就得到尺度關係:Ca_cr 正比於 θ_m³,反比於 log(1/λ),λ 是滑移長度與巨觀尺度的比(式 2.12)。其二,Eggers(2004、2005)對「平板從液槽往上拉」的完整流體力學計算給出(式 2.13)Ca_cr = (θ_m³/9) × [log(…)]⁻¹,括號內含平板傾角、Airy 函數的極值與 λ,論文說它「在微觀角 θ_m 上有普適性(Ca_cr ∼ θ_m³),但幾何前因子不普適」。其三,Hocking(2001)窄縫重力流的結果(式 2.14)Ca_cr ≈ 0.6 / (3 log(λ⁻¹)) × θ_m³,該文的圖 8 顯示它比 Eggers 與他們自己的數值都偏高。其四,Kistler(1993)的假設是動態角降到 0 時潤濕失效,而實驗上接觸線在失效附近會變成 V 形(Blake & Ruschak 1979);該文也提到 Snoeijer 等人 2006 的提拉實驗沒有在 Eggers 預測的 Ca_cr 發生轉變,原因是形成了毛細脊,並推測與接觸角遲滯有關,連續模型沒處理遲滯。
這個式子的重點只有一個:臨界毛細數不是常數,是 θ³ 乘上一個對數。對數裡的 λ 動 100 倍(10⁻³ 到 10⁻⁵),Ca_cr 只降約四成,而 θ 動兩倍,Ca_cr 動八倍。
三、大角度也成立嗎:Kansal 等人 2024
Kansal、Datt、Bertin 與 Snoeijer 的〈Viscoelastic wetting transition: beyond lubrication theory〉(arXiv 2407.07833,我讀到全文)處理平板提拉,把潤滑理論推廣到任意接觸角。本文的牛頓流體極限(式 47,黏彈參數 N0 = 0)給出
Ca_c / θ_e³ = 1 ÷ [ 9 ln( c₃ · (Ca_c/θ_e³)^(−1/3) ) ],c₃ = 0.423 θ_e² ÷ (λ̄_s θ_p)
(我從 PDF 抽出的文字排版有些亂,式子是我依其定義重讀的;小板傾角取 θ_p,垂直板取 θ_p → √2。)λ̄_s 是滑移長度除以毛細長度,該文固定取 10⁻⁴。文中的圖給了三個可核對的量級:θ_e = 10° 的小傾角案例,Ca_c 的縱軸是 10⁻⁴ 量級;θ_e = 90° 垂直板的縱軸是 10⁻² 量級;Ca_c/θ_e³ 在 θ_e 從 30° 到 150° 的曲線在縱軸 10⁻² 量級上收斂,論文稱潤滑結果「在大到 150° 的接觸角依然異常準確」,理由是 Cox 的大角度展開含的積分和 θ³ 在 150° 內只差幾個百分點。這給了 09-06 與昨天的 θ³ 外推一個正面的支持:在 57° 用立方律,不算外推過界。
四、對帳:同一條線,兩個不同的 θ_r
墨 γ = 0.04 N/m、μ = 0.01 Pa·s(沿用前幾天的值,黏度是從 D 反推的假設),γ/μ = 4 m/s,所以 U_c = 4 × Ca_c。兩種算法並列:簡單版取昨天的式子令動態角為零,Ca_c = θ_r³ ÷ (9 × 10.7),即 ln(L/L_S) = 10.7;Kansal 版是解上面的式 47(λ̄_s = 10⁻⁴)。
| θ_r | 簡單版 Ca_c | 簡單版 U_c | Kansal 式 Ca_c | Kansal 式 U_c |
|---|---|---|---|---|
| 10° | 5.5×10⁻⁵ | 0.022 cm/s | 7.4×10⁻⁵(小傾角板) | 0.030 cm/s |
| 30° | 1.5×10⁻³ | 0.60 cm/s | 2.0×10⁻³ | 0.78 cm/s |
| 57° | 1.0×10⁻² | 4.1 cm/s | 1.15×10⁻² | 4.6 cm/s |
| 70° | 1.9×10⁻² | 7.6 cm/s | 2.05×10⁻² | 8.2 cm/s |
| 90° | 4.0×10⁻² | 16 cm/s | 4.1×10⁻² | 16.5 cm/s |
(驗算:57° = 0.9948 rad,立方 0.9846,÷ 96.3 = 1.02×10⁻²,× 4 = 4.1 cm/s。)兩種算法差不到 30%,但這是兩個只差在 ln 取值的同類式子,不是獨立驗證。敏感度:57° 時 ln 取 5 與 15,U_c 分別 8.8 與 2.9 cm/s;λ̄_s 取 10⁻³ 至 10⁻⁵,Kansal 式的 U_c 為 6.2 至 3.7 cm/s。
所以昨天的矛盾這樣解:09-06 的 10⁻² 與 4 公分每秒,等價於假設 θ_r ≈ 57°;昨天引用的「後退角近 0」若成立,臨界速度就是每秒零點幾毫米到零點幾公分,不是 4 公分。兩者的差不是某個機制算錯了,是輸入的 θ_r 不同。反過來,若書寫速度取 10-03 的 V* 約 14 公分每秒,Ca = 0.035,在簡單版下需要 θ_r ≈ 86° 才不被夾帶。09-06 那條「Ca_cr ≈ 10⁻² 來自平板塗佈實驗」我今天回頭找,找不到對應的原始出處,搜尋結果裡有一段摘要寫「前進線可到 Ca 約 50、後退線約 10⁻²」,但它沒指向可打開的論文,所以那兩個數字我不再引用。
五、更正:後退線夾帶的是液體,前進線夾帶的才是空氣
09-06 寫「後退接觸線…開始夾帶空氣」。這個方向是錯的。Ledesma-Aguilar、Hernández-Machado 與 Pagonabarraga 的〈Theory of wetting-induced fluid entrainment by advancing contact lines on dry surfaces〉(arXiv 1306.2569,我讀到全文)引言寫得很直接:液體前沿被強迫從固體撤退時,在親水基材上能完全脫濕的速度上限是 U_rec ∼ Aγθ_e³/η_L(A 是幾何前因子),超過它液體被固體帶走,表面留下一層薄膜;前進線的情況是氣體被排開,臨界速度可以大得多,量級約公尺每秒,後退約公分每秒,他們的解釋是前進線排開的是氣體、後退線撤走的是液體,氣體的黏滯摩擦小,所以前進線的失穩門檻較高。他們也寫到前進臨界速度與接觸角的關係複雜,簡單的 θ³ 外推不能涵蓋。
Marchand、Chan、Snoeijer 與 Andreotti 的〈Air Entrainment by Viscous Contact Lines〉(Phys. Rev. Lett. 108, 204501, 2012;arXiv 1110.5274,我讀到全文)把平板插進黏性矽油,在超過臨界速度時看到厚約 10 微米的空氣膜,然後碎成氣泡。對 09-06 的推論有兩處傷害。其一,黏度從 0.02 變到 5 Pa·s(約 250 倍),夾帶速度只降約 10 倍,遠弱於預期的 1/η,所以前進線的空氣夾帶不能套 U_c ∝ 1/η;09-06 把「膠加倍、飛白門檻砍半」寫成通則,對後退線的液膜夾帶依 U_rec ∼ γθ³/η 才成立。其二,他們量到同一種液體「重新潤濕」的洞擴張速度(每秒 1 到 10 公分)比反向的「脫濕」洞(每秒 40 微米到 1 公分)快幾個數量級,最高約一千倍,正是前進與後退的不對稱。
把這三條線對照,針對墨與毛筆:
| 接觸線 | 誰被夾帶 | 門檻量級 | 對書寫的意義 |
|---|---|---|---|
| 墨在紙上往前推(前進,排開空氣) | 空氣 | 公尺每秒(Ledesma 引言) | 書寫 1 至 30 公分每秒低兩個數量級,一般不會因此留白 |
| 墨從毛或從紙面撤退(後退) | 墨(留下液膜) | γθ_r³/η,見上表,零點幾至十幾公分每秒 | 不是「斷」,是「帶膜」;膜厚歸 Landau–Levich(09-06 的 h = 1.34R·Ca^(2/3)) |
| 毛束內部的流動 | 不適用 | 由 V*(10-02 至 10-03)的供給對需求決定 | 飛白主要在這裡 |
這張表的前兩列是文獻明說的;哪一條線對應筆的哪個位置,是我的推論,沒有一篇文獻是針對毛筆的。若推論成立,飛白就不是「接觸線被扯斷」,而是供墨不足或紙面粗糙造成覆蓋不全。這與 09-29 到 10-03 的供給模型一致,那條線從頭到尾沒有用到 Ca_cr。Ledesma 的引言還提到,表面的非均質(他們引用超疏水表面的例子)可以降低前進線的夾帶門檻,我只讀到引言層級的一句話,紙面粗糙是否同樣有效,今天沒有答案。
六、θ_r 趨近 0 時會怎樣
Ca_c ∝ θ_r³ 的直接推論是 θ_r → 0 時門檻趨近 0,意思是對完全潤濕的表面(毛髮被墨完全潤濕),只要有速度就會帶出液膜,沒有「安全速度」。Sbragaglia 的引言轉述 Eggers 的結論是「不穩定發生在外部靜態彎月面接近完全潤濕流體的形狀(表觀角趨近 0)時」,與這個圖像一致。但 Devauchelle、Josserand 與 Zaleski 的〈Forced dewetting on porous media〉(arXiv cond-mat/0509349)研究零微觀接觸角的多孔板提拉,摘要說仍存在一個臨界毛細數,且其對板傾角的依賴在大控制參數下是 3/2 次方;我只讀到摘要,沒有全文,不知道他們的式子與量級,也不知道孔隙吸液如何改變了 θ→0 的極限。毛筆毛束與宣紙都是多孔體,這一點在今天的模型裡完全沒有處理。
🧠 記
**一、後退線的臨界毛細數 Ca_c ≈ θ_r³ ÷ [9 ln(…)],不是常數:**θ_r = 10°、30°、57°、90° 依序約 5.5×10⁻⁵、1.5×10⁻³、1.0×10⁻²、4.0×10⁻²;換成墨(γ/μ = 4 m/s)是 0.022、0.60、4.1、16 公分每秒。09-06 的 10⁻² 與 4 公分每秒等價於 θ_r ≈ 57°,和昨天的「θ_r 近 0」不矛盾,只是輸入不同。
**二、09-06 標錯了:後退線夾帶液體(留下液膜),前進線才夾帶空氣,且前進門檻約公尺每秒、後退約公分每秒;前進線的夾帶速度對黏度的依賴遠弱於 1/η。**飛白不能再寫成「後退線被空氣扯斷、門檻 4 公分每秒」。
**三、毛筆與宣紙是多孔體,今天的式子全來自平板或平行板的非多孔幾何;多孔板的摘要顯示零接觸角仍有臨界值,但式子沒讀到。**所有 θ_r 的數值是假設,沒有墨在羊毛或宣紙上的實測。
✍️ 實踐
做一個可以在廚房完成的後退線門檻測試,檢驗 Ca_c ∝ θ³ 與 U_c ∝ 1/η 兩件事。器材:載玻片或塑膠片(PET、投影片)各一、燒杯、清水、甘油或糖漿、尺、手機慢動作、一個能等速提拉的辦法(把片子夾在釣魚線上用手轉動捲筒,或用步進馬達;手拉也可以,用尺在影片裡逐幀讀位置)。
第一步,量接觸角。把片子浸入液體再慢慢拉出,從側面拍照,量液面貼著片子處的後退角 θ_r,寫下角度與影像。
第二步,逐步加速提拉,看哪個速度開始在片子上留下完整的液膜(膜面平整反光),而不是乾乾淨淨離開。記下門檻速度 U_c。水的 γ/μ 約 72 m/s,按簡單版,θ_r = 20°、30°、45° 的預測是 3.2、11、36 公分每秒。若量到的 U_c 與預測差一個數量級以上,先檢查 θ_r 是否量錯,再檢查片子表面是否乾淨。
第三步,把水換成含一半甘油的水(黏度升高,預期約 6 倍,需用黏度計或查表確認),θ_r 近似不變,預測 U_c 降為約 1/6。這測的是 1/η,但那是後退線的標度,不是前進線。
限制:毛筆不是平板;接觸角遲滯會讓門檻拖到毛細脊出現之後(Sbragaglia 引言提到的現象);玻璃片若有油污,θ_r 會偏大;這個實驗只驗證平板上的標度,不驗證毛束。
明天接:(一)去讀 Devauchelle 等人 2005 的全文,取出多孔板 θ→0 的臨界毛細數式子,檢查是否能代入毛束的孔隙率與滲透率;(二)重寫 09-06 的飛白段落,把「接觸線斷裂」改成「供墨與覆蓋」的說法,並與 V* 的供給模型對上;(三)Um & Lee 1997 與 Yang & Jin 2014 原文。
🔗 延伸學習
- arXiv 0812.3303:Sbragaglia, Sugiyama & Biferale, Wetting failure and contact line dynamics in a Couette flow(J. Fluid Mech. 2008) — 今天式 2.11 至 2.14 與 Ca_cr ∼ θ_m³ 的來源,我讀到全文。
- arXiv 2407.07833:Kansal, Datt, Bertin & Snoeijer, Viscoelastic wetting transition: beyond lubrication theory — 牛頓極限的式 47 與大角度仍準確的論述,讀到全文;主題是黏彈流體,我只用其牛頓極限。
- arXiv 1306.2569:Ledesma-Aguilar et al., Theory of wetting-induced fluid entrainment by advancing contact lines on dry surfaces — 引言明寫後退線夾帶液體、前進線夾帶氣體,以及兩者臨界速度的量級差,讀到全文。
- arXiv 1110.5274:Marchand, Chan, Snoeijer & Andreotti, Air Entrainment by Viscous Contact Lines(Phys. Rev. Lett. 108, 204501, 2012) — 前進線空氣夾帶對黏度的弱依賴,以及潤濕與脫濕速度的不對稱,讀到全文。
💬 問 AI
我在修正一個毛筆吸墨與飛白的模型,請逐項回答,並嚴格區分三類:(一)文獻明確寫了的、(二)你能查證或憑訓練資料確定的、(三)純屬推論的,每項標明把握程度(高/中/低);不確定或查不到就直說,不要為了流暢而編造文獻、公式或數字。
【A 組:後退線臨界毛細數】
1. 我引 Cox-Voinov 令動態角為零得 Ca_c = θ_r³ / (9 ln(L/L_S)),並用 Kansal 等人 2024(arXiv 2407.07833)的式 47 牛頓極限 Ca_c/θ_e³ = 1/[9 ln(c3 (Ca_c/θ_e³)^(-1/3))]、c3 = 0.423 θ_e²/(λ_s θ_p) 交叉檢查。請確認這兩式與我的代入是否正確,並指出 Eggers 2004/2005 的原式裡 log 的內容。
2. 我取墨 γ = 0.04 N/m、μ = 0.01 Pa·s,算得 θ_r = 57° 時 Ca_c ≈ 10⁻²、U_c ≈ 4 cm/s;θ_r = 10° 時 Ca_c ≈ 5.5×10⁻⁵、U_c ≈ 0.022 cm/s。請核算,並說明這個 θ³ 標度在 θ_r 大於 60° 時還可信到什麼程度。
3. 接觸角遲滯(前進角與後退角不同,表面粗糙或化學不均)如何改變這個門檻?Snoeijer 等人 2006 提拉實驗沒有在 Eggers 預測的 Ca_cr 轉變,原因是什麼?
【B 組:前進與後退的夾帶物】
4. 請確認:後退線夾帶的是液體(留下液膜),前進線夾帶的是氣體。我原先把後退線寫成「夾帶空氣」,請指出我這個錯誤在文獻裡對應的正確表述。
5. Ledesma-Aguilar 等人 2013 說前進臨界速度約公尺每秒、後退約公分每秒;Marchand 等人 2012 說前進線夾帶速度對黏度的依賴遠弱於 1/η。對黏度 10 mPa·s、表面張力 40 mN/m 的墨,前進線的空氣夾帶 Ca 大約在什麼範圍?文獻裡有沒有毛細數的實測表?
6. 我曾引用「前進線可到 Ca 約 50、後退線約 10⁻²」這組數字,但找不到原始論文。你知道這組數字出自哪裡嗎?如果不確定,請直說。
【C 組:毛筆與多孔體】
7. Devauchelle、Josserand 與 Zaleski 2005(arXiv cond-mat/0509349)研究零接觸角的多孔板提拉,摘要稱仍有臨界毛細數,對板傾角的依賴是 3/2 次方。請說明這個結果的物理圖像,以及能否搬到纖維束(毛筆)的情況。
8. 毛筆「飛白」的成因,在你所知的流體力學文獻中,有哪些候選機制(供墨不足、毛束間隙、紙面粗糙、接觸線失穩)?哪些有定量證據、哪些只是推論?
9. 請把所有存疑項目分成三堆:(一)憑訓練資料就能判斷的;(二)必須查一手文獻的;(三)需要用毛筆與宣紙實驗才能回答的。