用间断有限元方法求解二维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. 现在还不能捕捉激波，因为没有加限制器。



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



哦 ^_^-Send a written discontinuous finite element program--- for solving two-dimensional Euler equations

The attachment is a C language program I wrote it myself. Discontinuous finite element method for solving two-dimensional Euler equations.

Example is the flow around a cylinder. Zip file which also contains the mesh generation process to write a MATLAB



Sequence, just do small changes can generate sparse different grid.



The program also some problems as follows:



Because there is no boundary, so the residuals of the calculation will not converge to a given decimal, such as



0.001 and so on. When the maximum residual is about 0.1, press the ESC key on the keyboard, the program calculates the end, will output



A file Result.dat, the use tecplot software to open this file, you can get some results, such as flow lines



And so on. Because the result does not converge, there will be two in the back of the cylinder vortices appear, if the calculation of longer, there will be



Th