[請益] 誠心請教-目前遇到的問題

看板Mechanical (機械)作者gmichen (G_Mi)時間13年前 (2013/01/17 13:15)推噓1(1推 0噓 1→)

留言2則, 2人參與討論串1/1

各位大大們好!! 感謝上次有許多大大熱心的幫助，尤其是兩位M大!! 不過這次又遇到了點困難想請教各位大大.. 先感謝各位大大抽空看完並回復了!! 我目前把網格建立起來，如圖 (0,1) ┌───┬───┐ │ │ │ │ │(0,0) │ (-1,0)├───┼───┤(1,0) │ │ │ │ │ │ └───┴───┘ (0,-1) 若是設定Delta x和Delta y都為1um，所以Fo可以用公式Fo=[(a*Delta t)/(Delta x)^2]以及2-D的邊界(Fo小於等於1/4)代入，求得最大Delta t值為3.125 nano-sec. 接著由初始值(P=0)溫度求出P=1的溫度值，但是那求出的是在(0,0)下一刻時間的溫度，那麼依照我想要的情況，我在上面的(0,1)設一個較高溫度點，假設它溫度為26度，其他點都設為25度的話，我就得到(0,0)他下一刻時間的溫度了。但是，發現Fo因為我先前將它定為0.25(就是1/4)，因此explicit方式的(1-4Fo)T(0,0)這項就會一直為0，所以不管計算幾次它沒辦法影響到溫度改變... 所以第一個想請問該怎麼去定義Fo的值? 第二個問題是：若是我想要求其他點(0,1)、(0,-1)、(-1,0)的下一刻時間的溫度，是不是要把網格畫的多格一點，再由周遭溫度利用同樣方式去計算它的下一刻的溫度呢? 如下圖 ┌───┬───┐ │ │ │ │(0,1) │ │ ┌───┼───┼───┼───┐ │ │ │ │ │ │ │ │(0,0) │ │ ├───┼───┼───┼───┤ │(-1,0)│ │ │(1,0) │ │ │ │ │ │ └───┼───┼───┼───┘ │ │(0,-1)│ │ │ │ └───┴───┘ 以上是我目前的解法跟遇到的問題，再麻煩大大們幫我解答!! 真的感激不盡... -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.125.35.181