参多角形台风眼与旋转弹簧的关系 - 动力与模拟研究室.PPT

* 　　本報告利用平面上的旋轉彈簧了解多角形颱風眼中心空氣塊的受力情形。旋轉彈簧是一個被輕彈簧繫著的質點在水平面上旋轉的動力系統，它將表現一個多重時間尺度的運動、呈現多角星形的軌跡，我利用空間和時間的無因次量總結數值實驗的結果。旋轉彈簧質點的受力情形解釋了空氣塊在颱風眼中受到小型低壓中心的影響，以及平均氣壓梯度力隨與颱風中心距離增加的趨 勢，也凸顯颱風中心小渦旋的重要性。 　　衛星觀測資料顯示，颱風眼並不一定是圓形卻可能是多角形的結構，例如2001年侵台的納莉颱風、2003年襲美的Isabel颶風，都被觀測到明顯的三角形、五角星形的颱風眼結構（圖一）。這個現象在被觀測到之前，就已經被科學家以正壓模式預測出(Kossin和Schubert 2001)。 　　描述一種質點軌跡形狀，需要兩個無因次的量： 壹、前言 圖一、颱風眼的可見光衛星影像，上圖為2001年納莉颱風，下圖為2003年Isabel颶風 角形的軌跡。它與多角形颱風眼的關係引起我的好奇心，希望找出颱風眼中類似的作用力，了解氣壓梯度力的分布情形。 貳、旋轉彈簧的行為 大開眼界──利用彈簧探討颱風眼形狀的秘密 一、旋轉彈簧系統簡介 　　質點對原點的總能量、角動量L守恆，會表現出二維的運動，並帶有旋轉、伸縮週期兩種時間尺度。當初始時彈簧所給予的拉力恰提供旋轉所需的向心力，且徑向速度為零時，質點處於平衡狀態、其軌跡會是圓形的，設平衡 半徑為r*，速度為v*，滿足 　　旋轉彈簧系統由一個力常數為k 、原長為r0且沒有質量的理想彈簧，一端固定在座標原點，另一端繫著質量為m的質點在光滑水平面上運動（圖 二）。 　　 L = mr*v*　、 　　若不處於平衡狀態時，質點會呈現週期性的運動，產生多角形的軌跡。 　　設質點座標為(x,y)，與原點距離為r，運動時間為t，由虎克定律可以列出以下運動方程 式： 　　給定初始x, y位置及速度dx/dt, dy/dt即可以4階Runge Kutta數值積分方法透過MATLAB軟體 模擬質點的軌跡。 二、系統變數分析 　　本報告中將旋轉後相同的軌跡形狀視為同一種，因此可將初始徑向速度、y座標設為零不失一般性。 　　總共只有五個獨立的變數決定旋轉彈簧的行為：L、xi、m、r0、k。根據Buckingham Pi定理，系統可被含有兩個無因次變數的方程式所 描述。由平衡狀態時的受力方程式可推得 為無因次的時間 另一個無因次變數則是無因次的空間 兩者互相獨立。 1. 角數W：質點旋轉一圈時經歷的振盪次數，可能不是整數，反應在軌跡圖形上就是角數。 2. 振盪幅度A：　　　　　　　　，代表質點運動時遠離平衡位置的程度。 圖三、數值模擬旋轉彈簧質點軌跡與各種變數的關係 　　一個被彈簧繫著的旋轉質點在運動過程中，會受到拉近、推出的力量，產生多 Typhoon asymmetries and the rotating spring with multiple time scales 圖二、旋轉彈簧系統 　　圖四顯示角數W與K、N的關係，當K愈大、N愈小時，W與K的值會很接近，這是由於離開平衡半徑很少、軌跡仍接近圓形，同時旋轉週期遠大於振盪週期。因此旋轉一圈的週期 可以用平衡半徑的旋轉週期近似。 三、實驗結果 圖四、角數W與N、K的關係　(a)W(N,K)　(b)W(N,K)/K，其中綠色部分的值很接近1 (a) (b) 　　圖五顯示振盪幅度A與N、K的關係，同樣當K愈大、N愈小時，A會與N的值很接近。若N太大，質點能向外拉伸彈簧的距離比向內壓 縮的距離大的多，造成A與N的值不接近。 圖五、振盪幅度A與N、K的關係　(a)A(N,K)　(b)A(N,K)/N，其中橘色部分的值很接近1 (a) (b) 参、多角形颱風眼與旋轉彈簧的關係 圖六、正壓模式模擬的多角形颱風眼，填色部分代表氣壓，數值是與外圍的壓力差(hPa)，如右方色標所示，黑線代表流 線，白色箭號代表氣壓梯度力 　　圖七是質點在旋轉彈簧中和空氣塊在颱風眼中的受力比較圖，在多角形颱風眼中空氣塊受到的氣壓梯度力受到不同尺度擾動的影響，表現得相當複雜。 (a) (b) (c) (d) 圖七、受力比較圖　(a)旋轉彈簧的質點軌跡，藍圈是質點每轉過固定角度時的位置　(b) 藍圈處的加速度，由左而右對應到逆時針方向的藍圈　(c)填色部份為氣壓，紅線為一條流線，紅圈是對應到(a)圖藍圈的位置　(d)紅圈處的氣壓梯度 力，由左而右對應到逆時針方向的紅圈 一、正壓模式中的多角形颱風眼 二、氣壓梯度力分析 　　以距中心最遠處為中點，沿著流線平均兩側的氣壓梯度力，對與中心的夾角做圖。並以同樣的方式計算旋轉彈簧質點的平均加速度，與多角形颱風眼做比較： 圖八、平均方法示意圖　(a)旋轉彈簧平均加速度　(b)白色箭號代表氣壓梯度力，