找回密码
 注册
楼主: lixl-imech

[新闻动态] 发布OpenCFD-SC 二维源代码 (更新:Mach 6 钝锥的完整算例及详细说明)

  [复制链接]
发表于 2011-12-4 17:11:17 | 显示全部楼层
下载了学习学习,谢谢李老师的无私贡献
发表于 2011-12-8 11:30:09 | 显示全部楼层

回复 24# merryfish 的帖子

面向对象的功能在fortran里纯粹是糟粕
发表于 2011-12-11 19:03:55 | 显示全部楼层
支持,下载
发表于 2012-2-22 14:41:36 | 显示全部楼层
感谢刘老师对国内CFD领域做出的贡献,开源需要很大的勇气!
发表于 2012-3-11 18:41:07 | 显示全部楼层
赞啊~~学习了
发表于 2012-3-22 21:03:45 | 显示全部楼层
李老师,请教柱坐标的奇异点怎么消除的?

我用有限体积写了个程序,压力在奇异点附近的控制体积是个三角形,解方程却发现不对.
发表于 2012-4-10 00:03:05 | 显示全部楼层
感谢李老师!
发表于 2012-6-6 15:42:51 | 显示全部楼层
国家应该大力支持李老师这种人才
 楼主| 发表于 2012-6-7 16:51:00 | 显示全部楼层
原帖由 hilllike 于 2012-3-22 21:03 发表
李老师,请教柱坐标的奇异点怎么消除的?

我用有限体积写了个程序,压力在奇异点附近的控制体积是个三角形,解方程却发现不对.


在奇性轴周围,利用物理量的旋转特性,给出奇点另一侧的值,从而在奇点周围也可以采用高精度差分法计算。
例如在柱坐标系下,r=0是奇点。 在r=0处不布置网格点(如rj=(j-1/2)*h  , j=1,2....N).  对应坐标线(seta=const), 利用柱坐标的性质,给出奇点另一侧的值,从而可以在j=1  点也可以采用高精度差分法计算。

当然,这种方法稳定性不是太好,奇性轴周围的流场容易变形。
最好的方法是在奇点采用有限体积法计算。

评分

1

查看全部评分

发表于 2012-6-19 21:40:33 | 显示全部楼层

回复 100# lixl-imech 的帖子

我这里再提供一种方案。对奇轴上的值可以取周围一圈临近点的值的平均值代替该奇轴点的值。这种方法的稳定性相对比较好一些。
发表于 2012-10-11 18:22:38 | 显示全部楼层

谢谢

中国CFD的发展需要开源和奉献,感谢李老师!
发表于 2012-10-14 19:20:59 | 显示全部楼层
不知道基于隐式时间离散方法的程序功能可否开源共享。
发表于 2012-11-27 16:31:13 | 显示全部楼层
感谢 无私学者
发表于 2013-1-17 22:46:36 | 显示全部楼层
发表于 2013-1-17 22:46:54 | 显示全部楼层
:loveliness: :loveliness:
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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