关于“和细长结构有关的重叠网格技术”的问题

apfel

anyone介绍过的“和细长结构有关的重叠网格技术”：处于多维空间中的一维细长结构，多维空间和一维细长结构分别用两套独立的网格；两个区域解的耦合是通过“交换源项”来实现，即先找到两套网格交叉的部分，用q＝h×A×dT求出他们之间的热交换（或其他物理量交换），这个q即为两套网格中交叉网格部分的源项，这样分别算出两个区域内的温度分布，温度更新之后，又可计算出q的新值。。。如此反复，迭代求出两个区域内的温度分布。
嗯，希望我说清楚了。现在的问题是结果不收敛，一直变一直变！
各位有处理过此类问题的吗？帮我看看思路有没有问题？
anyone路过的话请一定帮忙看看。
多谢了！

anyone

不收敛可能原因可以有是很多的. 如果只从区域解的耦合来讲:
1. 首先保证守恒性, 也就是加在2边的Q数值相等,符号相反,总能量不变.
2. 你现在是显式处理, 收敛会变慢, 隐式对代数求解器有一定要求不是很好做. 显式可以把源项线化, 可以改进收敛性.
3. 如果 H 的值比较小收敛会容易一点. 你可以人为把H变小, 看能不能收敛, 验证一下其他部分是不是做对了.

apfel

我现在怀疑它是不是能收敛：只要有温差，加在两边的Q就不为0，那两边的温度就要一直变下去啊？

anyone

如果高温一侧有热源, 低温一侧有散热就会达到平衡.
如果低温一侧是绝热, 高温是给温度值, 最后温差会变0.
难到你条件给错了?

apfel

做的是组织与血管的传热，血管先给定一个温度分布，求出组织温度作为初始温度；血管比组织温度高；每经过一次迭代血管和组织的温度都下降，血管是由于向组织散热而温度降低，而组织由于这一次迭代的Q比上一次要小所以温度也下降，所以他俩的温度以几乎相等的趋势一直降，看不出有收敛的趋势。

anyone

看起来你传热的概念有一些问题.
你把问题说完整了, 边界条件/初始条件, 物性, 还有血管/组织间换热系数(你怎么得到的), 我可以帮你看看.

apfel

[这个贴子最后由apfel在 2005/06/29 12:16pm 第 1 次编辑]

一写就多了，麻烦你帮我看看吧，多谢多谢^_^

anyone

看了你写的东西, 有以下看法:
1. 根据你的边界条件, 组织和血管内的温度应该在37℃-37.5℃之间, 如果出界了, 肯定算错了, 如果没出界, 计算也许是对的, 如果温度分布趋势也对, 具体温度值和你给的参数有关.
2. 影响组织和血管内的温度值的主要参数是
  A. 血管的流速,密度,和热容量
  B. 组织的传热系数.
3. 既然用已知的努谢尔数来做, 血管用一维模型就足够了, 解起来也方便(用TDMA), 而且流速恒定不用解流动问题.
4. 血管进出口的边界条件是容易出错的地方, 建议用简化模型验证一下:
  A. 把组织拿掉,只解血管, 把热源也拿掉, 结果(不管什么初始条件)温度应该都是37.5℃
  B. 只解血管, 加上热源, 假设组织温度不变37℃,看血管温度分布. 这个问题用一维可以比较容易得到解析解, 比较一下就知道对错了.