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

发布OpenCFD-EC2D; 更新 Ver 1.5.4 ; 三维软件使用协议

  [复制链接]
发表于 2011-4-9 14:27:45 | 显示全部楼层
谢谢,学习下!!
 楼主| 发表于 2011-4-11 20:01:46 | 显示全部楼层
2011-4-11 升级到了版本1.1, 支持多重网格(最多3重)。
程序、算例及更新版的理论手册可在本贴的主贴中下载。
下图为使用单重网格、二重网格及三重网格计算RAE2822翼型绕流问题的残差收敛情况。可以看出,多重网格对促进残差收敛很有帮助。

[ 本帖最后由 lixl-imech 于 2011-4-11 20:17 编辑 ]
export.jpg
 楼主| 发表于 2011-4-11 20:13:29 | 显示全部楼层
多块网格搜索算法的理论手册。
采用该算法,可以利用多块对接网格的坐标信息构建各块之间的连接关系。

多块网格连接信息搜索算法-ver1.4.pdf

304.2 KB, 下载次数: 1277

发表于 2011-4-26 23:11:38 | 显示全部楼层
李老师,请问您的软件有没有加入网格嵌套功能的想法?
 楼主| 发表于 2011-4-27 09:05:20 | 显示全部楼层
发表于 2011-4-27 09:16:07 | 显示全部楼层
李老师,这个程序会考虑增加湍流模型吗?
 楼主| 发表于 2011-4-27 15:03:52 | 显示全部楼层
原帖由 ch06 于 2011-4-27 09:16 发表
李老师,这个程序会考虑增加湍流模型吗?

以后会增加一个SST模型,计划支持BL, SA, SST三个模型。
 楼主| 发表于 2011-4-27 15:06:51 | 显示全部楼层

多块结构网格前处理程序: 利用几何信息,构建连接关系

发布我编写的一个小程序,可以利用多块对接网格的坐标信息构建各块之间的连接关系。
功能:根据多块对接结构网格的几何信息,通过搜索算法,构建网格连接信息。
要求:网格必须为多块对接网格(对嵌套网格无效);网格文件名为Mesh3d.dat,必须为PLOT3D格式。


使用方法:
使用时,将其与PLOT3D格式的网格文件Mesh3d.dat放在一个子目录里,编译并运行即可。
运行时,根据提示,输入如下信息
   1) 首先输入距离限度Dist_limit (阈值), 本程序使用距离搜索匹配连接点, 距离小于Dist_limit的两个点被认为是“连接点”。如果输入负值(如-1.),则软件自动确定该阈值 (选取全局最小网格间距的一半作为该阈值), 推荐输入-1., 让程序自动选取阈值。
2) 询问Mesh3d.dat是Formatted格式的文件还是Unformatted格式的文件。如是Unformatted格式的,则输入0;如果是Formatted格式的,则输入1.
3) 询问是否需要全局搜索匹配点,如未建立匹配信息,则需要搜索,需要输入 0。
输入这三个数后,程序开始运行,运行后产生BXCFD格式的连接信息文件 bcin。 bcin文件的格式与BXCFD的连接文件兼容,无任何连接点的子块,其连接号设置为-1 。
本程序只能识别内边界信息,无法区分其他边界(如固壁、远场、对称面等)。需要用户手工修改这些无连接号的子面,根据实际的边界条件,修改为-10 ,-20, -30等。

4) 1.4之后的版本附带“智能识别”功能,对于类似DLR-F4/F6之类的标模网格,可以“智能识别”固壁、远场及对称面。因为这类网格的计算区域通常为矩形,对称面为y=0的面,远场边界为该矩形的其他5个面,其余面为固壁。(仅)对于这种情况,可以使用“智能识别”功能。即,程序提示“Do you want to search the wall and Far field boundary (1 for yes, 0 for no ) ? ” 时,输入1即可。
其他情况下,“智能识别”出来的固壁、对称面及远场边界会出现问题。 请输入0。 然后手工编辑bcin文件,将标识为“-1”的边界,人工设定为固壁、远场或对称面。

[ 本帖最后由 lixl-imech 于 2011-4-29 11:23 编辑 ]

多块对接网格构建连接关系程序使用说明.doc

30 KB, 下载次数: 896

OpenCFD-EC理论手册-ver1.1.pdf

1.13 MB, 下载次数: 1022

get_bcin_v1.4.1.rar

8.59 KB, 下载次数: 1021

多块网格连接信息搜索算法-ver1.4.doc

1.52 MB, 下载次数: 1360

回复 支持 0 反对 1

使用道具 举报

 楼主| 发表于 2011-4-27 15:14:23 | 显示全部楼层
这个程序根据多块对接网格的几何信息,构建网格连接文件bcin. 该格式有些类似于CFL3D的inp文件,但语法有所改动。
    今后将对其进行修改,使其可生成兼容CFL3D的网格连接文件。
   
     算例可到我的网盘下载:
     http://cid-1cc0dcbff560c149.office.live.com/self.aspx/.Public/OpenCFD-EC/get%5E_bcin/getbcin.rar
    是从M6翼的网格文件Mesh3d.dat构建连接信息文件bcin的算例。需要说明的是,构建出来的网格连接文件无法区分固壁、对称面及远场,对于这些面的边界条件均设置为"-1", 需根据实际情况手工更改。

   用户可下载DLR-F4/F6之类的网格(只能支持多块对接网格),测试一下这个软件。

[ 本帖最后由 lixl-imech 于 2011-5-2 16:59 编辑 ]
发表于 2011-4-27 20:28:35 | 显示全部楼层
这个得仔细研读一遍
发表于 2011-4-28 01:18:43 | 显示全部楼层
李老师,会考虑加入预处理 以计算低马赫数问题吗
发表于 2011-5-3 10:34:48 | 显示全部楼层

Roe格式数值通量计算

李老师你好,我发现Roe格式的数值通量计算中似乎有个bug,在Harten型熵修正处,做if语句判断的时候,应该先给lamda(1~4)赋值为:
lamda(1)=abs(avu-ava)
lamda(2)=abs(avu)
lamda(3)=abs(avu)
lamda(4)=abs(avu+ava)
,否则,不管这四个特征值是多大,都将被进行熵修正,大大的加大的数值耗散。
另外,请问李老师,如果需要把多重网格改为W循环(1-2-3-2-3-2-1),我应该在哪里添加哪些语句?非常感谢。
发表于 2011-5-3 10:44:58 | 显示全部楼层
增加了多重网格后,计算速度提高了很多倍。李老师,是否考虑增加隐式推进格式LUSGS,LU-ADI等,进一步提高速度。
 楼主| 发表于 2011-5-3 17:25:32 | 显示全部楼层
原帖由 wingnuaa 于 2011-4-28 01:18 发表
李老师,会考虑加入预处理 以计算低马赫数问题吗

这个程序目前主要针对飞机的空气动力学计算。今后会增加预处理模块,以便快速地计算低速构型。
 楼主| 发表于 2011-5-4 09:58:52 | 显示全部楼层
原帖由 fanghongfox 于 2011-5-3 10:34 发表
李老师你好,我发现Roe格式的数值通量计算中似乎有个bug,在Harten型熵修正处,做if语句判断的时候,应该先给lamda(1~4)赋值为:
lamda(1)=abs(avu-ava)
lamda(2)=abs(avu)
lamda(3)=abs(avu)
lamda(4)=abs(avu+ ...

多谢你找出这个Bug.
如对多重网格的算法进行修改,仅需对sub_NS_multigrid.f90里面的子程序 NS_3stge_multigrid()进行修改即可。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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