2010年中科院力学所《计算流体力学》课程（含课件和录像，共15讲）

lixl-imech

第10讲的PPT, 2010-5-24:15:00修改
第10讲PPT,2010-5-24:22:00修改；

[ 本帖最后由 lixl-imech 于 2010-5-24 13:54 编辑 ]

gutouqie

李老师，二维问题的steger-Warming 分裂，对于曲线坐标系下怎么处理？

lixl-imech

原帖由 gutouqie 于 2010-5-24 13:44 发表

                               
登录/注册后可看大图

李老师，二维问题的steger-Warming 分裂，对于曲线坐标系下怎么处理？

Steger-Warming分裂由直角坐标系推广到曲线坐标系中非常简单。 因为通常Steger-Warming 给出的都是针对组合通量：
                             F=alfa*f1+beta*f2
的分裂方法。 在直角坐标系中，令alfa=0, beta=1就得到了x方向的通量（分裂），令alfa=1,beta=0即y方向的通量（分裂）。
    在曲线坐标系中，通量为：
                             E=kx/J * f1 + ky/J*f2; F=ix/J*f1+iy/J*f2
                              （J 为坐标变换的Jocabian行列式，kx,ky,ix,iy为四个坐标变换系数）
令alfa=kx/J, beta=ky/J  (或alfa=ix/J, beta=iy/J) ，带入Steger-Warming分裂的表达式，就得到了曲线坐标系下的分离形式。

[ 本帖最后由 lixl-imech 于 2010-5-24 14:01 编辑 ]

守望者811

虽然我看不懂，但是我还是学

gutouqie

原帖由 lixl-imech 于 2010-5-24 22:00 发表

                               
登录/注册后可看大图

Steger-Warming分裂由直角坐标系推广到曲线坐标系中非常简单。 因为通常Steger-Warming 给出的都是针对组合通量：
                             F=alfa*f1+beta*f2
的分裂方法。 在直角坐标系中，令alfa=0, be ...

那么， 就是E（T)=kx/J * f1 (T)+ ky/J*f2 (T); F (T)=ix/J*f1 (T)+iy/J*f2 (T)吗？T为特征值

lixl-imech

原帖由 gutouqie 于 2010-5-25 00:29 发表

                               
登录/注册后可看大图

那么， 就是E（T)=kx/J * f1 (T)+ ky/J*f2 (T); F (T)=ix/J*f1 (T)+iy/J*f2 (T)吗？T为特征值

是的。 具体可参考： 《计算空气动力学》(傅德薰 主编） 第30-31页，第148-150页及160-162页。

lixl-imech

第10 讲录像已传至网盘
http://cid-1cc0dcbff560c149.skydrive.live.com/browse.aspx/.Public

同时上传到了力学所的FTP，位置为：   /Incoming/document/李新亮-CFD-讲课录像
所内用户可以从该FTP下载，速度会快些。

lixl-imech

程序Sod-HLL.f90 : 使用HLL及HLLC通量计算Sod问题
程序刚编制完，用时一个小时左右，比较仓促，估计里面有Bug, 因为中心区压力算得不好。 程序在附件中，欢迎大家下载参考。

明天有时间我再修改一下，推出新版的程序。
调试完成，发现最后绘图时，把压力算错了。 新的程序（Ver 1.1）中已改正过来。 目前的结果是正确的。

本版程序，欢迎下载使用及传播。 但请别忘记标注 （在论文或报告中引用作者文献或致谢均可）。

[ 本帖最后由 lixl-imech 于 2010-5-31 03:27 编辑 ]

sariel

感谢老师！！！呵呵！！！！

lixl-imech

测试算例1： 一维Sod激波管问题， t=0.14时刻的数值解。
测试算例2： Shu-Osher 激波-密度扰动波干扰问题，t=1.8时刻的数值解。

方法1： 差分法，Steger-Warming FVS,  5阶WENO差分格式；
方法2： 有限体积法， 5阶WENO重构， 使用HLLC型近似Riemann解构造通量；
方法3： 有限体积法， 5阶WENO重构， 使用HLL型近似Riemann解构造通量。
结果见下图。
分析：
       HLL与HLLC 解效果相差不大， HLLC比HLL在接触间断处分辨率更好些。
       针对Sod问题，HLL/HLLC型FDS方法效果优于Steger-Warming型FVS方法。
       对于Shu-Osher问题， HLL/HLLC及Steger-Warming FVS方法效果接近。
结论：
     与Steger-Warming FVS方法相比，HLL/HLLC更适合分辨强间断，其中HLLC更好些。
       在光滑区，分辨率主要受重构方式（即差分格式）决定，三种通量形式带来的差别比较小。

[ 本帖最后由 lixl-imech 于 2010-5-26 08:58 编辑 ]

lixl-imech

第11讲 （代数方程组的求解， 网格生成 ) PPT
2010-5-31 13:40 修改

[ 本帖最后由 lixl-imech 于 2010-5-31 05:30 编辑 ]

lixl-imech

使用守恒变量与原始变量重构

lixl-imech

第11讲的录像已传至网盘 http://cid-1cc0dcbff560c149.skydrive.live.com/browse.aspx/.Public

zinsser

个人比较喜欢原始变量，守恒变量重构时很难保证正压力和正密度。另外也可以利用特征变量重构，在做DG等高阶方法时一般需要采用特征变量才能得到好的高阶效果。

通流

守恒变量重构时很难保证正压力和正密度
------------------------------------------------------

有这事？