【已附代码】求助一个激波捕捉问题，涉及物理平面->计算平面的边界转化。

chenye1717

参考李新亮老师的计算流体教程，我做了一些捕捉间断算法的算例，大多数集中在一维情况。最近想做一个涉及到坐标转换的二维绕圆激波捕捉问题（flow around the circular cylinder）。比较容易做的方法是用steger-warming算法，当我写好代码运行起来，得不到正确的数值解。思考了一些可能性的原因如下：
1）教程中的二维steger-warming算法的公式是否有错误？我对比了steger-warming原文文献的通量公式，没有差异。
2）我认为比较有可能出错的地方就是从物理平面向计算平面转换过程中的边界条件的转化，这个问题是比较基础性的，我把问题写成了一个PDF文件在此列出。希望各位朋友能够给我解惑！

谢谢。

[ 本帖最后由 chenye1717 于 2012-5-29 13:39 编辑 ]

Gerald135

楼主你好，我也非常想做个曲线坐标系下的激波问题，但是我坐标变换的代码弄的不是很清楚，好像坐标一变换，特征矩阵和通量分裂的项都要变，都要加上由于坐标变换引起的项，你能把你代码发上来大家看看吗？

通流

也许你应该看看边界条件是如何处理的。李老师的教程里不可能没有。

chenye1717

好的，我把代码传上来，大家一起看一看，我认为导致错误的问题在计算平面边界条件这块。

chenye1717

是的，我也认为问题出在边界转化上。我去李老师的教程中去找这块的知识补。谢谢斑竹的建议！

chenye1717

个人认为：坐标变换后，特征值（也就是波传递的特征速度）变了，所以特征矩阵和分裂公式都变。但是steger-warming通量公式是比较友好的，带入坐标变换系数后可以得到该方向上的通量，所以我认为主要的问题在 “边界条件转化” 上。

Gerald135

楼主，至于边界转换问题我也帮不了你，我还没学到那儿，你发的程序希望其他人能看到并帮到你，（PS：程序为什么不用FORTRAN写啊？）

chenye1717

好，多多交流。 我的Q：550567511

程序我一开始就用matlab写的，因为matlab的代码老是go to的 有点晕

chenye1717

好，多多交流。 我的Q：550567511

程序我一开始就用matlab写的，因为matlab的代码老是go to的 有点晕