发一个自己写的间断有限元程序---求解二维Euler方程

glandetian

附件是我自己写的一个c语言程序。用间断有限元方法求解二维Euler方程。
算例是一个圆柱绕流问题。压缩文件里面还包含了一个用MATLAB写的网格生成程序，只需要做很小的改动就可以生成稀疏程度不同的网格。

程序还有些问题如下：

1.  因为没有对边界处作处理，所以计算的残差不会收敛到一个给定的小数，如0.001等。当最大残差为0.1左右时，按下键盘上的ESC键，程序计算结束，会输出一个文件Result.dat，用tecplot软件打开这个文件，可以得到一些结果，如流线等。因为结果不收敛，在圆柱的后面会有两个涡出现，如果计算时间更长，会有涡脱落的现象出现。如果想计算收敛，则需要对边界作特殊处理。处理方法见2006年的文献High-order accurate implementation of solid wall boundary conditions in curved geometries 。或者使用高阶的单元，如6节点三角形单元。

2.  现在还不能捕捉激波，因为没有加限制器。

如果您在这个程序的基础上解决了上面提到的问题，请联系我，也共享一下程序哦 ^_^


有同学说下载不了，我也不知道是怎么回事，阅读权限设为0 表示不限制 已经设为0了。

[ 本帖最后由 glandetian 于 2012-6-7 09:24 编辑 ]

zqb138

怎么下载不了？阅读权限255啊

glandetian

重新发一次，权限设为100

glandetian

重新发一次

030301029

楼主强大

onesupeng

也就是说程序不对了啊

glandetian

程序的主体框架是对的，只需要作2处改动就行了，原文已经说明白了。

ch06

阅读权限设置为100是只给版主看的吗？

glandetian

上传的时候，不知道权限是干什么用的，随便设了一个，我到现在也没有搞清楚有什么用。

eryun05

什么意思，下载不了的

eryun05

权限不够呀

eryun05

楼主能否把权限设置的低一点，好让更多人学习呀

铁洲宇

求助，我也用的间断有限元，发现有个问题，质量矩阵是一定的，那如果，我改变局部索引对应全局索引的顺序，总的质量矩阵会变。。比如局部索引1,2,3，对应全局索引是17,18,19,此时a(1,2)=a(17,18)，换成18,19,17就变成a(1,2)=a(18,19)。。。还是说这样都没关系？？？

铁洲宇

求助，我也用的间断有限元，发现有个问题，质量矩阵是一定的，那如果，我改变局部索引对应全局索引的顺序，总的质量矩阵会变。。比如局部索引1,2,3，对应全局索引是17,18,19,此时a(1,2)=a(17,18)，换成18,19,17就变成a(1,2)=a(18,19)。。。还是说这样都没关系？？？

glandetian

每个人编程都不一样。
我用的c++ 用类封装单元 每个单元对应一个类的对象，并且唯一对应一个单元的编号。然后定义该类的数组，计算时，依次计算各个单元就可以了。没有你所说的全局索引什么的。