上传我写的DG源码

ch06 · 发表于 2012-6-15 14:47:05

没用过DG，帮不上忙了

shirazbj · 发表于 2012-6-15 15:07:45

已经算到1300了，看来要算到3000.
一个结果还没有，要是断电就玩完了。
要把我的老爷机算死了，已经一下午了呀。

shirazbj · 发表于 2012-6-15 18:24:00

果然是3000. 程序的输出做的好简单。
原来在main的开头可以改。
结果文件是存的三角形每个顶点上的值？

print *,t/(0.018/3.0/(um+sqrt(gm*pm/rm))),tn

是不是要给的tn要大于这个算出来的值？

[ 本帖最后由 shirazbj 于 2012-6-16 10:18 编辑 ]

zqb138 · 发表于 2012-6-15 23:34:57

一般是要大于这个算出来的值，这个值是CFL条件的一个估计，但是有时估计不是很准确，按照DG稳定性理论要求|Lamda|max*dt/dx<=1/(2*k+1)来估计dt

shirazbj · 发表于 2012-6-16 08:00:34

终于算出数来了。

[ 本帖最后由 shirazbj 于 2012-6-16 12:51 编辑 ]

zqb138 · 发表于 2012-6-16 15:52:42

不用算那么夸张的长时间吧，你的电脑什么配置呀？
还有就是我觉得不看三维立体图是看不出实际的计算结果的，毕竟云图只是一个视觉上的大概结果，真正的细节问题必须看三维立体效果图。

shirazbj · 发表于 2012-6-16 22:30:12

Gfortran比较奇怪，有时编译完运行极慢。改改代码再编译，运行又快了。

光dg，不加限制器是不是不行？那样好像很容易溢出。

这个是不是指能算可压的？

[ 本帖最后由 shirazbj 于 2012-6-17 09:46 编辑 ]

zqb138 · 发表于 2012-6-17 16:59:30

这个就是算可压缩的呀，对具有强间断的非线性问题，涉及激波，没限制器常常会不稳定。不可压缩的NS方程也可用DG，不过没有尝试过。

fly100fly100 · 发表于 2012-6-18 19:05:27

怎么一直是Euler方程的？求NS

shirazbj · 发表于 2012-6-18 21:11:37

原帖由 zqb138 于 2012-6-13 22:53 发表

 登录/注册后可看大图

从结果看在激波处有较多的刺，应该是限制器的问题，其他几个算例也有这样的问题，实在不明白啊。

是不是网格太粗了，显得有刺？

zqb138 · 发表于 2012-7-1 08:29:15

这个网格数应该不少了，这个程序中可能还是有些问题，只是我已经形成思维定式了，怎么也看不出来了。

iceflying · 发表于 2012-7-13 16:10:11

居然被我发现了

金蛇狂舞 · 发表于 2012-7-18 15:46:55

好人啊。赞一个。

金蛇狂舞 · 发表于 2012-7-19 15:43:46

请问你的DG程序中，单元界面上流通量是怎么取的？

zqb138 · 发表于 2012-7-19 16:39:44

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

上传我写的DG源码

回复 23# zqb138 的帖子

回复 1# zqb138 的帖子

站点信息