Re: [請益] 暫態移動熱源

看板Fortran作者 (peter)時間12年前 (2012/04/15 13:27), 編輯推噓2(201)
留言3則, 2人參與, 最新討論串1/2 (看更多)
如果我是你,我會先去研究解2-D heat equation的演算方法為何, 這應該google或是查一些CFD的書都可以查到,知道怎麼做了以後,先假設一個最 簡單的邊界條件,例如在y = 0 處為等熱通量來coding,程式跑了沒問題之後 再去修改成你要的邊界條件。以你所說的邊界條件,可能可以這樣寫: do it = 1, nt !time loop q_dot = 0 if (t < tau) then do i = 0, x1 q_dot(i) = constant end do else if (tau < t < 2*tau) do i = x1,x2 q_dot(i) = constant end do else if (.....) . . . end if end do 這寫法應該是笨笨的寫法,但目前想到就這樣了。 ※ 引述《adiadiadi (笑一笑沒什麼)》之銘言: : 各位好 : 最近在寫一個暫態二維的熱擴散方程式 : 小弟才疏學淺,在移動熱源的地方一直卡住 : 這熱源由雷射加熱(加熱寬度s) : 每隔時間 tou 會在下一個間隔landa的位置加熱一次 : 加熱位置皆為j=0的地方(即x軸) : 下面為簡單圖示:XD : tou=1 : ^ y : | : | : | : | : | : |___________> x : ^ : | : power : tou=2 : ^ y : | : | : | : | : | : |___________> x : ^ : | : power : tou=3 : ^ y : | : | : | : | : | : |___________> x : ^ : | : power : 現在腦殘中...希望可以得到一個思考方向 : 或是有code可以參考 : 感激不盡 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 99.126.178.148

04/15 13:40, , 1F
謝謝你! 其實我現在也是這樣寫...只是想知道有沒有其
04/15 13:40, 1F

04/15 13:41, , 2F
更簡單的表示法。 引為這樣寫感覺很像暴力展開..
04/15 13:41, 2F

06/10 17:23, , 3F
威斯康辛一哥
06/10 17:23, 3F
文章代碼(AID): #1FYbmoa7 (Fortran)
文章代碼(AID): #1FYbmoa7 (Fortran)