推 adiadiadi:謝謝你！ 其實我現在也是這樣寫...只是想知道有沒有其 04/15 13:40
→ adiadiadi:更簡單的表示法。 引為這樣寫感覺很像暴力展開.. 04/15 13:41
推 philodi:威斯康辛一哥 06/10 17:23
如果我是你，我會先去研究解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