旋轉、飛躍、球進了耶！－香蕉球物理學

《論語》〈鄉黨〉裡，關於吃飯的規定，超—級—多—！！「食不厭精，膾不厭細。食饐而餲，魚餒而肉敗，不食。色惡不食，臭惡不食。失飪不食，不時不食。割不正不食，不得其醬不食。肉雖多，不使勝食氣。唯酒無量，不及亂。沽酒市脯不食。不撤薑食，不多食。祭於公，不宿肉。祭肉不出三日，出三日，不食之矣。食不語，寢不言。雖疏食菜羹瓜祭，必齊如也。」^[1]吼～這麼囉嗦，有本事自己來啊！

有些男人激不得。

物理學家 Mathieu Sellier 與 Edouard Boujo 就因為前者的妻子提出挑戰，而用電腦運算出最佳烹調模型，還在 2019 年的《物理評論流體》（Physical Review Fluids）期刊上，分享成果，造福饕客。^{[2, 3, 4]}全文第一句話，是這麼說的：「本論文研究固化流體薄膜，受制於複雜的運動學，在固體表面的流動…」，^{[5, 6]}意思是「我們要教大家做可麗餅。」

完美可麗餅的定義

撇開二位科學家基於品味差異，而無法達成共識的餡料不談，^[2]這個研究中，可麗餅的最高境界，被定義為「厚度均勻，無孔洞，且呈現完美圓形」。^{[6, 7]}要在自家廚房，達成此終極目標，通常會遇上難題：當麵糊在鍋底鋪開，同時也會逐漸被煮熟。如果水平放置鍋子，麵糊便在平均地觸及鍋緣之前固化。為避免這個問題，一般有兩種常見的解決辦法：第一種是用刀具迫使麵糊在鍋中散開，類似刮刀塗層的動作；另個做法則是將鍋子傾斜旋轉，令麵糊往低處流動。^[6]

運算製作可麗餅的模型

二位科學家採用「伴隨優化」（adjoint optimisation）的數學方法，描述流體在活動容器中的運動，模擬以最小施力，鋪出最平坦的可麗餅。^{[7, 8]}其中考慮的因素，包含：以通過鍋子圓心的縱軸為中心運動；藉重力鋪開麵糊；以及隨溫度改變濃度的麵糊與旋轉中的鍋子的互動。^[3]經過一番（讓人反胃的）計算過程，他們找到最佳的做法：先快速翻動鍋子，然後在煮的過程中，緩慢旋轉。^[7]

詳細的技巧，如下：一倒入麵糊，就馬上以陡峭的角度傾斜鍋子，把液體逼到邊緣。接著，手腕輕扭，轉鍋子一圈，確保麵糊完整鋪滿底部。在轉動的時候，傾斜的角度得逐漸縮小，轉速也隨著麵糊固化而趨緩。當覆蓋動作完成，鍋子也恢復水平狀態。^{[4, 7, 8]}

圖中，深紅處麵糊最厚，深藍則最薄。可麗餅的製作流程，由左上開始，先向下，再依序往中、右二欄進行。^[7]

起初濃厚的（紅色）麵糊被推向鍋子的右上緣，把稀薄的（淺藍）剩料拋在後頭。然而隨著順時鐘的轉動，麵糊逐漸勻稱地分佈於整個鍋底。^[7]

_{圖／參考資料 6，figure 6}

做鬆餅救眼疾

科學家們之所以對餅皮類食物的製作如此著迷，是因為類似的手法不僅能生產巧克力，幫智慧型手機螢幕鍍膜，^[4]還可以懸壺濟世。2016 年倫敦大學學院（University College London）在 YouTube 上，也發佈了一個看似玩物喪志的作品。全長約 5 五分鐘的影片裡，前 4 分鐘幾乎都在以嚴謹的態度，講述鬆餅（此指 pancake，而非 waffle）的製作。到了最後卻話鋒一轉，道出製餅技術與眼疾治療的關係。原來手術中控制眼睛內部液體外流的皮瓣（surgical flaps），就要倚靠類似的原理來研發。^[9]

救世的精神，於是賦予了科學家一個精進廚藝的學術使命。

參考資料

1. 論語/鄉黨第十（維基文庫）
2. Making the Perfect Crêpe (APS Physics, 2019)
3. The hard-hitting science behind crepes and beauty pageants (Chemical & Engineering News, 2019)
4. Using fluid dynamics to perfect crêpe cooking techniques (Phys.org, 2019)
5. Boujo E and Sellier M. (2019) ‘Pancake making and surface coating: Optimal control of a gravity-driven liquid film’. Physical Review Fluids, 4, 064802.
6. Boujo E and Sellier M. (2019) ‘Pancake making and surface coating: Optimal control of a gravity-driven liquid film’. arXiv
7. Physicists Think They’ve Finally Found the Trick to Making Perfect Crepes (Science Alert, 2019)
8. A computer model explains how to make perfectly smooth crepes (Science News, 2019)
9. Understanding the physics of pancakes to save sight (University College Lodon on YouTube, 2016)

疫情期間學烹飪，再拍照上傳社群網站，是凡夫俗子的成果發表；將煮義大利麵的心得筆記，發表在《流體物理學》（Physics of Fluids）期刊上，^[1]則是科學家的華麗炫技。

煮麵的動機

美國伊利諾大學Sameh Tawfick副教授的實驗室，專攻靈活可變形的纖維和有彈性的結構，所產生的「流體結構交互作用」（fluid structure interaction）。「過去幾年老是開玩笑，說義大利麵的黏著力與我們的研究息息相關」，他說團隊發覺分析麵條力學質地的轉變，「可以體現黏著力、力學質地和烹煮熟度的關聯。」^[2]

以上研究動機有聽沒懂，無所謂。煮麵要緊。

煮麵的方法

當 COVID-19 疫情襲來，學生與博士後研究人員，開始在家中和實驗室裡煮義大利麵。^[2]

他們在實驗室裡，用的不是鍋子和瓦斯爐等通俗的烹飪設備，而是燒杯與加熱板。^[3]加熱板是一種有發熱平面的機器，比開放式的火源安全。^[4]

依照一般食譜的建議，煮義大利麵的步驟，是先把麵條扔進大滾的鹽水中，等熟了再撈出來瀝乾。^[5]不過，實驗的這個部份沒加鹽，大概是想把變數降到最低，後面再討論食鹽的作用。以滾水煮熟的麵條，在脫離水的時候，殘留的水份會在麵條之間形成「彎液面」，以「表面張力」把一根根的麵條黏在一起。^{[2, 3]}

_{左圖／參考資料3}

整個過程就像圖中所示，左側為二條（橙色）義大利麵條，被置在（粉色）加熱板上的燒杯內烹煮；中央是麵條離水；右邊則為局部放大圖，呈現義大利麵條之間的彎液面。

表面張力之所以會產生，是因為水分子與水分子之間，每個方向的引力本來都一樣，可以相互抵銷，然而到了水面就失衡了，因此有儲存能量的張力。 ^{[6, 7]}

_{右圖／參考資料3}

傳統義大利麵講求的口感，叫做「al dente」，意思是「煮到內硬外軟」，恰到好處。從物理的角度來看，水份由麵條表面，逐漸擴散進入內部，所以首先軟化的當然是最外層。吸水的過程中，麵條體積會隨之膨脹。煮愈久，效果愈明顯。下圖是研究團隊在觀察義大利麵「吸濕膨脹」（hygroscopic swelling）時，進行的量化紀錄。^[3]

圖片上方中間的麵條剖面，從圓心向外，由深至淺，有輕微的色彩變化。然而，水份達到飽和後的右上剖面，便不再有任何的顏色漸層，大概就是所謂煮過頭的狀態。^[3]

至於縱向的吸濕膨脹，則展示在下方。以最左邊的生麵條為比例尺，對照用 100 度 C 的滾水，分別烹煮 12、18、24、30 分鐘後的膨脹變化。整體而言，麵條剖面周長增加的比例，大於長度的成長比例。這是基於內部沒有與水接觸的核心，限制了麵條縱向的延展。^[3]

_{左圖／參考資料3}

所以，到底要怎麼做，才能擁有al dente的口感？

研究團隊發現義大利麵條達到 al dente 前，其周長與長度分別的膨脹率相比，所得的比率是3.5比1。一旦超過了，就會軟爛。^[2]此外，由於麵條煮愈久，離水時彼此相黏的部份就愈長。研究團隊認為，專業廚師也可以測量相黏長度，來推論起鍋時間，以後再將成熟的技術，推廣至普通民家…^[3]（原文口氣意外地認真。）總之，要知道煮好沒，不是用嘴試吃，也別拿錶計時，科學家的建議竟是用尺測量！

煮麵的鹽和光

明明義大利麵條的包裝上，都有建議的烹煮時間。為什麼科學家不直接告訴大家，煮多久能起鍋？原來如果照正常的煮法，在水中加食鹽，麵條的化學和力學特性都會起變化。比起用蒸餾水，鹽水不僅有助麵條膨脹，而且會增添嚼勁。^[3]此外，Sameh Tawfick 副教授解釋，滾水中的鹽量，會改變達到 al dente 所需的時間。有鑑於此，他未來要探討食鹽，在義大利麵膨脹時所扮演的角色。^[2]同時，這個研究正如一道照亮前程的光，或許會引領其他人，也來嘗試用簡單的方式，研究軟物質的特性。

參考資料：

1. Hwang J., Ha J., Siu R., Kim Y. S., and Tawfick S. (2022) ‘Swelling, Softening, and Elastocapillary Adhesion of Cooked Pasta’, Physics of Fluids, 34 (042105)
2. Physics Models Better Define What Makes Pasta Al Dente (Physics of Fluids, 2022)
3. Hwang J., Ha J., Siu R., Kim Y. S., and Tawfick S. (2022) ‘Swelling, Softening and Elastocapillary Adhesion of Cooked Pasta’, arXiv
4. Hot Plate Use and Safety in Laboratory (University of Wisconsin-Madison, 2013)
5. Spaghetti and Meatballs (Gourmet Traveller, 2018)
6. 第11章  有趣的界面現象（國立成功大學化學工程學系）
7. Chemical Science Lesson Plan: Hydrogen Bonding and Surface Tension (University of Illinois, 2010)
8. Enjoy better-cooked pasta with…physics and a ruler? (University of Illinois, 2022)

下雨天的時候走在路上，天氣涼涼的，聽著雨聲的感覺非常好。但是你有沒有想過，為什麼雨滴會從天上掉下來？

「啊！就像蘋果會掉到地面一樣，會受到重力的作用嗎？」你可能會這麼說。

好，那我們這邊就來帶大家算一下，一滴雨從高空落到地面，純粹只有受到重力時，應該是什麼樣子的感覺吧！

只有受到重力作用雨滴的運動分析

當不考慮空氣阻力時，由高空落下的物體全程會受到重力加速度值 g 的作用，而因為地表的重力加速度約為定值，以海平面且緯度 45º 為標準，其數值為 9.8m/s^{2 [1]}。因此雨滴從高空落下時，可以視為一個單純的等加速度運動，而這個運動我們又稱之為自由落體。

假設雨滴是靜止落下且受到重力加速度值 g 作用，即可根據等加速度運動公式，求得雨滴從高度 h 自由落下時的末速度值：

然而，在探討雨滴落下的末速度之前，我們必須對於雲的分類以及大致上的高度有一個基本的了解，才能比較明確地知道我們要探討的雨滴大概是從什麼樣的高度落下來的。

氣象學家 Luke Howard 於 1803 年中的著作《論雲的變形》（The Essay on the Modification of Clouds）中，按照不同雲的形狀、組成、形成原因，將雲分為 10 大雲屬，並且將這 10 大雲屬劃為三個雲族，分別為：位於距地表 6,000 至 7,000 公尺的高雲族，位於距地表 2,000 至 6,000 公尺的中雲族，以及位於距地表 0 至 2,000 公尺的低雲族^[2]。另外，則還有橫跨了三個不同雲族高度的直展雲族，常常造成短暫但是相當豐沛的降雨量^[3]。

按照國際氣象組織所提供的分類，以及 Luke Howard 的定義，天空中主要的降雨來源為積雨雲（cumulonimbus）以及雨層雲（nimbostratus），降雨來源以雨層雲較為常見，且其雲底多為 1,200 公尺以下。故我們這邊計算雨滴的高度時，便以 1,200 公尺作為高度的參考依據。

因此，當一滴雨從高空落下，代入前述自由落體公式，即可計算出雨滴理論上應該要有的末速度：

根據上述的計算式子可以知道，當雨滴從高處落下時，如果沒有任何的空氣阻力，雨滴落到地面的速度大約會是 153 m/s。

對於這個數字沒有感覺嗎？那這邊簡單地計算給你看一下，讓你有點 fu。但是在這個計算之前，首先我們要先對於雨滴的大小有個概念。

依照 2009 年的相關研究^[4]顯示，小雨滴在降落時幾乎是圓形，可是隨著體積越大，就會變得越扁平，受到空氣的影響也會越明顯。當雨滴達到特定的大小時，就會被切割為較小的雨滴，也因此最大的雨滴直徑會被限制在 6 mm 左右。

而按照另一個研究^[5]對於雨滴粒徑的分布探討，發現雨滴的直徑多數是落在 0.5 mm 至 4 mm 之間，也就是半徑 0.25 mm 至 2 mm 之間。

這邊先姑且不論雨滴本身的化學成分所帶來的密度差異，以及落下過程中的密度和質量變化。因此我們可以簡單的利用密度、質量和體積的關係式，假設有一顆雨滴的成分皆為水，密度為 1 g/cm³，半徑 2 mm，且為均勻球體的情況下，計算這顆雨滴的質量如下：

接著，我們利用牛頓第二運動定律和動量衝量的概念，來計算平均一顆雨滴所造成的衝擊力大小。這邊，我們假設你是淋雨的狀態，雨滴跟你的腦袋接觸的時間大約為 0.001 秒，且雨滴最後會完全靜止在你的腦袋上，也就是末速度為 0。

此時，造成雨滴會有速度變化的作用力有二，一為雨滴所受到的重力、二為腦袋給雨滴的正向力。根據牛頓第三運動定律，腦袋給雨滴的作用力，與雨滴給腦袋的作用力，為「作用力與反作用力」之間的關係。

那我們要怎麼知道雨滴對於腦袋的衝擊力有多少呢？

根據前面的假設，我們假設腦袋給雨滴的作用力使用變項為 N，可以列式如下：

雖然我們前面說，在計算正向力 N 時，應該要將重力納入考量，不過實際計算後會發現雨滴本身重量也不算大，相較之下，後面的重力項是可以忽略的，因此計算出來的衝擊力約為 0.52 kgw。

嗯？你說你還是沒有感覺嗎？再說白話一點好了，這個重量就差不多是一瓶 500 ml 的礦泉水壓在你身上的感覺。這只是單一顆雨滴，平常在下雨的時候絕對不可能只有一顆雨滴。一瓶礦泉水壓在身上其實是有感覺的，那很多雨滴下在身上，等同於很多很多瓶礦泉水壓在身上，那肯定也是非常有感。

修但幾勒，這個結論跟我們平常淋雨的感覺完全不同吧！那到底問題出在哪裡？

其實雨滴不只受到重力的作用

雨雲本身存在於大氣層的對流層內，而對流層內充滿很多空氣分子。當雨滴在這些空氣分子所形成的「流體」裡面移動的時候，會使得雨滴本身除了受到重力以外，還會額外受到空氣阻力（drag force）的作用。

在流體動力學中，在流體中移動的物體會受到一個和運動方向相反的阻力。這個阻力來自流體，會存在於兩個流體層之間，或者是流體與固體之間。可是，這和以往我們所學的固體和固體之間的摩擦力不同，因為物體在流體中受到的阻力其實是和物體移動的速度有關^[7][8]。

物體在流體中所受到的阻力，會受到物體大小、形狀、特性，以及流體性質的影響。阻力方程式（drag equation）概括了這些因素，描述如下^[7]：

其中，ρ 為流體的密度（如果是在空氣中，則是空氣的平均密度）、A 為物體在流體中的有效面積、v 為物體在流體中之速度；C_D 則是阻尼係數，是一個沒有因次的數字，一般來說會跟物體的形狀以及雷諾數（Reynolds number）有關。

而雷諾數則是在流體動力學之中，流體慣性力（inertial force）和黏性力（viscous force）的比值，用來預測流體狀態的無因次物理量。對於不同的流體來說，雷諾數會有很多不同的表達方式，但一般來說都會包含流體的密度（density）、黏滯性（viscosity）、流體的流速，以及特徵長度或尺寸。

最基本的雷諾數可以表示如下^[9]：

其中，ρ 為流體的密度，v 為流體的平均流速、D 為特徵長度，而 μ 則為流體的黏滯性。

雷諾數低的時候，流體會呈現層流（laminar flow）的狀態。流體分子會在每一層中平順流動，相鄰層之間就像堆疊的紙牌，鮮少或甚至幾乎沒有混合，當然也不會產生漩渦^[10]。

相反地，在雷諾數高的時候，流體則是會呈現紊流（turbulent flow）的狀態，流體的流速跟壓力沒有一定的變化規律，流體分子也沒有明顯的平行層，很常會互相混合在一起^[11]。

扯遠了扯遠了，我們還是繼續回到原本的阻力方程式。

根據實驗觀察，在雷諾數較高，也就是流體的密度較大、流速較快，而且黏滯性較小時，阻力係數可以幾乎視為定值。此時，阻力就會跟流體流速的平方成正比，公式如下：

而在雷諾數低，也就是流體密度較小、流速較慢且黏滯性較大時，阻力係數會和雷諾數的倒數成正比，因此我們結合雷諾數本身的定義以及阻力方程式，就可以知道「在雷諾數較低時，阻力與流速之間的關係為線性關係」，公式如下：

依照前面講過的阻力方程式和流速之間關係的背景知識，讓我們回到最一開始遇到的雨滴問題。

之前在分析雨滴的受力時，只有考慮到重力的作用，計算出雨滴自 1200 m 高的雨雲雲底落到頭上時，速度約為 153 m/s。在考慮到空氣阻力時，由於阻力與雨滴的運動方向恆相反，因此我們可以將雨滴的質量先以 m 作為變項，假設雨滴為正球形且半徑為 R，繪製雨滴所受到的力圖如下：

因為空氣阻力恆與物體運動的速度反向，而雨滴在落下的時候，速度一定是向下的，加速度也向下，故空氣阻力會向上。

阻力方程式中的 A 是投影的等效面積，在球形的雨滴中，即為上圖斜線部分，可以用半徑 R 和圓面積的公式來計算。此時，我們利用牛頓第二運動定律計算雨滴運動過程中所受到的加速度量值，來觀察雨滴運動的情形：

如果今天的流體狀況是屬於高雷諾數的情況（流體的密度較大、流速較快且黏滯性較小）時，則前述的式子可以下表示，並計算出加速度的關係式：

反之，如果是低雷諾數的情形（流體的密度較小、流速較慢且黏滯性較大），則前述的式子可以下表示，也順手計算出加速度的關係式：

從前面的兩條化簡式子，可以看出雨滴掉落時，不論雷諾數如何，速度漸大都將造成阻力漸大，並使得加速度漸小。當達到一定的速度時，雨滴就不再會有加速度，而是改以等速度的方式落下。此時，雨滴所具有的速度即終端速度（terminal velocity, vt）。在終端速度時，我們可以知道雨滴所受到的重力與拖曳力達到力平衡，因此可以根據不同的雷諾數而列式。高雷諾數的情況下所計算出的終端速度如下：

低雷諾數的情況下所計算出的終端速度如下：

我們這邊以高雷諾數的流體情形來考量大氣中的情況，與前面的條件相同假設，也就是雨滴為半徑是 2 mm 的正球體，雨滴密度主要成分為水，因此密度為 1000 kg/m³，而阻尼係數這邊我們根據雨滴的形狀和經驗公式簡單取 0.6 來概略估算^[14]。

利用高雷諾數的情況計算終端速度實際值時，會需要流體的密度。在這裡，我們討論的對象是空氣中的雨滴，故理想上（當然，這是很理想的情況下）可以使用理想氣體方程式來求出於 1 大氣壓、20ºC 時候的空氣密度，來代入終端速度的公式。

代入我們目前空氣的條件，也就是 1 大氣壓、20ºC 的情形，而這邊務必將所有單位都轉為 SI 制，加上理想氣體常數，此時使用的是 8.314。其中，M 為空氣的分子量，我們這邊使用 28.97 g 配合以上的條件代入計算^[15]。

將前述所得到的空氣密度數值，結合前面的其他條件，代入高雷諾數情況的終端速度公式，即可計算終端速度：

由計算結果可以知道，當考慮到空氣阻力時，雨滴會以 8.52 m/s 的終端速度落下，比起之前純粹考慮重力時，求出的 153 m/s 來說小了非常多，是原本的二十分之一。按照牛頓第二運動定律，這樣的雨滴打到腦袋時，對於腦袋瓜的正向力也會減為原本的二十分之一。如此一來，就比較像我們平常淋雨的情況了。

由前面的計算過程，我們可以明白從高空落下的雨滴不只有受到重力。能夠讓我們下雨天走在路上不被雨滴狠狠槌死的最重要因素，其實就是空氣阻力的功勞。同時，我們可以知道，造成雨滴落下的運動過程並非等加速度，而是變加速度運動。利用牛頓第二運動定律得出加速度的關係式後，也知道速度越來越大，加速度就會越來越小。在加速度為 0 時，則會以終端速度等速落下。

最後，讓我們來感謝空氣阻力，讓每一個人在下雨天的時候都能安心走在路上。

註解

• 註 1：Halliday, D., Resnick, R. and Walker, J. (2011) Fundamentals of physics. 9th Edition, John Wiley & Sons, Inc., New York.
• 註 2：Classifying clouds | World Meteorological Organization
• 註 3：中央氣象局，天空的魔術師—千變萬化的雲朵
• 註 4：Villermaux, E., Bossa, B. Single-drop fragmentation determines size distribution of raindrops. Nature Phys 5, 697–702 (2009).
• 註 5：McFarquhar, G. M. (2010). Raindrop size distribution and evolution. In Rainfall: State of the Science (pp. 49-60). (Geophysical Monograph Series; Vol. 191). American Geophysical Union.
• 註 6：Drop (liquid) – Wikipedia.
• 註 7：Drag (physics) – Wikipedia.
• 註 8：NASA (2010). What is drag?
• 註 9：Reynolds number – Wikipedia.
• 註 10：Laminar flow – Wikipedia.
• 註 11：Turbulence – Wikipedia.
• 註 12：What is Laminar Flow?
• 註 13：Viscosity – Wikipedia.
• 註 14：Drag coefficient – Wikipedia.
• 註 15：科學Online，大氣組成的檢測(Air)-上