一维热传导：自行编程计算与comsol计算结果对比

talkingbird · 发表于 2013-5-27 21:09:56

一维热传导方程，无内热源，限定边界条件，两端温度都是定值。
使用差分法，将将一维热传导方程进行离散化，编写matlab程序，根据初值和边界条件，进行迭代，计算了1秒内的温度变化。
comsol中材料使用的参数与matlab中计算用的参数一致。
对比结果见下图

talkingbird · 发表于 2013-5-27 21:13:47

计算结果相差很大，是什么原因呢？

talkingbird · 发表于 2013-5-27 21:27:28

而且，comsol软件的计算结果，在刚开始的时候，曲线还出现了一个尖角，我缩短了步长，发现仍然，存在这样的尖角，这是因为comsol软件的计算能力不行吗？

talkingbird · 发表于 2013-5-27 21:30:03

如图~~~~~

青梅煮酒28 · 发表于 2013-5-27 21:39:26

多做几个例子再试试，比如用二维热传导，弹性波，很多书有现成的偏微分方程，matlab程序

TBE_Legend · 发表于 2013-5-28 10:22:49

这个问题和comsol计算浓度时出现负值很类似。

若方便，把comsol文件传上来，大家帮你看看。

105032572 · 发表于 2013-5-29 22:46:48

众人拾柴火焰高

jinxiaoli_1981 · 发表于 2013-5-30 08:41:08

这跟comsol的网格节点分布有关，你把端部网格加密，这种数值突变就会改善

talkingbird · 发表于 2013-5-30 17:17:37

TBE_Legend 发表于 2013-5-28 10:22
这个问题和comsol计算浓度时出现负值很类似。

若方便，把comsol文件传上来，大家帮你看看。 ...

只能传几百K的文件呢
传不了呢

hanlianhuan · 发表于 2013-6-5 15:32:17

大概试了一下可能跟网格划分有关，网格密一点就解决了。

TBE_Legend · 发表于 2013-6-12 01:03:34

talkingbird 发表于 2013-5-30 17:17
只能传几百K的文件呢
传不了呢

用网盘，金山快盘等。

maxiao0 · 发表于 2013-6-27 19:59:21

这个存在一个阶跃，所以高阶形函数在网格点间会有一个波动。如果加密下网格，同时选用一阶线性形函数，便可以消除这个波动
但是从计算本事来看，选用高阶形函数在时间步后期，相同网格密度下，误差会更低些。

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

[05.传热] 一维热传导：自行编程计算与comsol计算结果对比