找回密码
 注册
流体中文网»论坛 › 教学科研 › CFD基础理论 › DG方法求解二维Euler方程的计算效果图 速度 u
返回列表 发新帖
查看: 3064|回复: 11

DG方法求解二维Euler方程的计算效果图 速度 u

[复制链接]
发表于 2011-11-8 10:58:23 | 显示全部楼层 |阅读模式

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？注册

x
这是速度u的contour图
u.jpg
回复

使用道具 举报

lwd1981
发表于 2011-11-23 12:39:04 | 显示全部楼层

回复 1# glandetian 的帖子

计算的应该是很低马赫数情况下的流动，算的不太好啊，出现尾迹了，数值粘性过大，关键是还特意把涡弄出来，呵呵。看来楼主需要补补流体力学的课啊，哈哈。按照达朗贝尔杨缪，压力分部应该前后左右，上下都对称！既没有阻力也没有升力的，呵呵。还有就是如果你用的是一阶格式，应该是可以捕捉激波的。二阶格式也不是任何时候捕捉激波都要限制器的，至少我的程序是这样，呵呵。
    要算好极低马赫数流动，构造合理的格式数值耗散是十分重要的，还有就是为了克服数值刚性，最好是采用预处理技术。

[ 本帖最后由 lwd1981 于 2011-11-23 04:55 编辑 ]
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:36:11 | 显示全部楼层

回复 2# lwd1981 的帖子

谢谢回复啊 ^_^
我在帖子“发一个自己写的间断有限元程序---求解二维Euler方程”中已经说明了出现涡的原因了----没有做边界处理。
现在，我已经将边界由直的改为二次插值的了 计算结果还算过得去
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:37:20 | 显示全部楼层
请看附图 再没有涡出现了

[ 本帖最后由 glandetian 于 2011-11-24 10:24 编辑 ]
u.jpeg
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:41:13 | 显示全部楼层
速度v
v.jpeg
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:41:45 | 显示全部楼层
压力p
p.jpeg
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:42:30 | 显示全部楼层
马赫数Ma
Ma.jpeg
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-24 09:50:45 | 显示全部楼层

回复 2# lwd1981 的帖子

非常感谢你的指点！ ^_^
上面的计算，来流马赫数都是0.3    形函数都是二次多项式
我的形函数是一次多项式时，如果不加限制器，也不能捕捉激波，因为总会出现负压力、负密度。当出现这种情况时，我总是让程序自动停止。用二次多项式时就更是如此了。
请问你是怎么处理负压力负密度问题的？

[ 本帖最后由 glandetian 于 2011-11-24 09:54 编辑 ]
回复 支持 反对

使用道具 举报

zqb138
发表于 2011-11-26 15:27:52 | 显示全部楼层
分片常数的一阶格式耗散那么强，几乎磨平激波。分片线性的二阶格式什么时候能在不用限制器的情况下捕捉激波呢？
回复 支持 反对

使用道具 举报

 楼主| 发表于 2011-11-27 09:10:02 | 显示全部楼层

回复 9# zqb138 的帖子

同问
回复 支持 反对

使用道具 举报

lwd1981
发表于 2012-3-19 20:43:11 | 显示全部楼层
读读Frink的博士论文吧，你就明白了即使是基于线性重构二阶精度格式，有时并非一定要有限制器的，呵呵

[ 本帖最后由 lwd1981 于 2012-3-19 12:44 编辑 ]
回复 支持 反对

使用道具 举报

zqb138
发表于 2012-6-6 19:08:40 | 显示全部楼层
lwd1981,能否将你所说的文献贴一下，我没有找到啊！
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表