| 查看: 1580 | 回复: 6 | |||
[交流]
simple算法-数值传热学-方形空腔内自然对流换热计算发散-解决方法求助-MATLAB编写已有2人参与
|
|
暑期自学数值传热学，与研究课题相关，用simple算法自编计算一个方型空腔内自热对流的换热，程序matlab语法上应该没问题，但是始终连续性方程发散，具体见下： 方腔尺寸：4cm*4cm 流体：水 边界温度条件：W（west）:100度 S:25度 N:绝热 S:绝热 水初始温度：25度 网格数量：10*10、15*15、20*20。40*40 迭代收敛判据（各方程余量的范数为err） 动量方程u:err<0.001 动量方程v:err<0.001 能量方程t:err<0.001 质量方程m:err<0.001 各方程迭代方法：高斯赛德尔迭代 对流项采用一阶迎风差分 时间步长tt：tt=0.01s 结果：网格采用20*20时，30步时间步长（0.3s）后发散，连续性方程发散； 采用网格40*40时，3步（0.03s）即发生发散，同样连续性方程发散； 索性减小网格尺寸10*10/15*15，反而可以完成一分钟的计算不发散，但是结果流场不太对； 将对流项差分改为混合格式、将err提高为0.0001也没有效果，代码写完自己调试了两周，毫无进展感觉还是需要指点，请大伙帮帮忙，指指路。 迭代的最后一次的计算结果（流线图和温度场），图片形式插入附件了，名称包含相关设置，例如20-20-30-l表示网格20*20-第30步-流线图  15-15-6000-l.jpg  15-15-6000-t.jpg  20-20-30-l.jpg  20-20-30-t.jpg  40-40-3-l.jpg  40-40-3-t.jpg |
回复此楼» 猜你喜欢
香港科技大学全奖招收博士/博士后，光子、纳米材料、液晶显示方向
已经有15人回复
-
法国博士招聘
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
法国博士招聘
已经有0人回复
-
超导三元凸包图
已经有0人回复
-
最近遇到这样一个问题
已经有10人回复
-
论文接收，散金币了！
已经有167人回复
-
vasp算居里温度
已经有0人回复
2楼2019-08-04 00:49:02
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
3楼2019-08-04 09:00:12
4楼2019-08-06 23:44:23
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
5楼2019-08-07 08:45:23
kingspin
木虫 (著名写手)
- 应助: 254 (大学生)
- 金币: 4472.5
- 红花: 74
- 帖子: 1635
- 在线: 315.5小时
- 虫号: 2585875
- 注册: 2013-08-07
- 性别: GG
- 专业: 聚合物共混与复合材料
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
涡的尺寸和形态一般来说与能量有关，外载荷很小时，一般形成一系列小尺寸的涡，当体系能量提高后，小涡会融合成大涡。但如果再提高，就跨入了湍流的范畴，涡结构就被破坏了，流场不再具有稳态形式。如果我的记忆没错的话 发自小木虫IOS客户端 |
6楼2019-08-07 08:50:35
|
666 发自小木虫Android客户端 |
7楼2019-10-16 15:04:16
15