Re: [請益] 暫態移動熱源
如果我是你,我會先去研究解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
討論串 (同標題文章)
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章