作者aicwebb (peter)
看板Fortran
标题Re: [请益] 暂态移动热源
时间Sun Apr 15 13:27:12 2012
如果我是你，我会先去研究解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
1F：推 adiadiadi:谢谢你！ 其实我现在也是这样写...只是想知道有没有其 04/15 13:40
2F：→ adiadiadi:更简单的表示法。 引为这样写感觉很像暴力展开.. 04/15 13:41
3F：推 philodi:威斯康辛一哥 06/10 17:23