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

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

  [复制链接]
 楼主| 发表于 2011-5-4 09:59:27 | 显示全部楼层

今后的版本应当会增加隐格式功能。
发表于 2011-5-11 19:35:48 | 显示全部楼层
李老师,这个程序的前处理会改进成与商业软件结合吗? 比如读入ICEM或gambit的多块网格
 楼主| 发表于 2011-5-14 20:30:52 | 显示全部楼层
原帖由 ch06 于 2011-5-11 19:35 发表
李老师,这个程序的前处理会改进成与商业软件结合吗? 比如读入ICEM或gambit的多块网格

目前还不能从ICEM或gambit直接读入网格。今后会向这个方向发展。
发表于 2011-5-18 09:20:27 | 显示全部楼层

关于OpenCFD-EC2D-1.1程序的问题

C:\Documents and Settings\Administrator\桌面\未命名.PNG你好,李老师!你的OpenCFD-EC2D-1.1程序在我的机器计算RAE2822流场时出现了如下图的问题,不知道是什么原因?期待你的解答,谢谢!

[ 本帖最后由 chenhuanlong 于 2011-5-18 01:23 编辑 ]
未命名.PNG
 楼主| 发表于 2011-5-22 19:54:30 | 显示全部楼层
原帖由 chenhuanlong 于 2011-5-18 09:20 发表
C:\Documents and Settings\Administrator\桌面\未命名.PNG你好,李老师!你的OpenCFD-EC2D-1.1程序在我的机器计算RAE2822流场时出现了如下图的问题,不知道是什么原因?期待你的解答,谢谢!


哦,也许是程序中有Bug.
附件中的文档是OpenCFD-EC2D ver 1.12, 修正了ver 1.1中的几个Bug. 我刚调试过,应当可以运行。

运行方法:把附件中的文档解压开,仅需编译运行opencfd_ec2d_v1.12.f90即可。

OpenCFD-EC2D-1.12.rar

263.54 KB, 下载次数: 670

发表于 2011-5-23 15:43:52 | 显示全部楼层
有木有三维的圆钝体的高超飞行模拟啊。。。。球李老师给我发一份吧  我是个新手菜鸟  想看看程序。。。在论坛转了好久。。。我都不知道从何入手了。。。前几天老关注FASTRAN  软件。。发现在7组元和11组元的问题上不是很理想。。二维的模拟。。。公式中e 的符号没有办法敲入。。好像他的软件只针对5组元。。。2组元分析。。求李老师帮助啊   。。。我的邮箱123191492@qq.com。。期待ing。。
 楼主| 发表于 2011-5-23 16:51:42 | 显示全部楼层
原帖由 hua123191492 于 2011-5-23 15:43 发表
有木有三维的圆钝体的高超飞行模拟啊。。。。球李老师给我发一份吧  我是个新手菜鸟  想看看程序。。。在论坛转了好久。。。我都不知道从何入手了。。。前几天老关注FASTRAN  软件。。发现在7组元和11组元的问题上不 ...


  带真实气体效应及化学反应的高速流动程序我们也刚刚开发。目前由我的学生在编制程序,估计得一段时间后才能调试通过。
发表于 2011-5-24 00:00:59 | 显示全部楼层
谢谢李老师了。。。我们很期待你的成果。。。加油。。。
发表于 2011-5-25 17:22:27 | 显示全部楼层
李老师您好!我想请教一下您EC版程序里面插值重构的WENO格式和SC版里的WENO格式区别在哪里呢?刚入门CFD,还不太明白,谢谢您!
发表于 2011-5-25 19:18:43 | 显示全部楼层
我用1。1倒是能算下去。
从mesh3,算到2,再1,再1。最后只是在mesh 1上算了。我还以为多重网格要3-2-1-2-3-2-1。。。的轮着算呢。

1和1之间为什么有个  Initialize?

变mesh的时候,t重新开始了。mesh 1到1间的t不知道怎么变的。例子是稳态的还是瞬态的呀?

  1. ----------------- OpenCFD-EC2D ver 1.1---------------------------------
  2.          Copyright by Li Xinliang, lixl@imech.ac.cn                     
  3.          Programming by Li Xinliang & Leng Yan, 2011-04-11              
  4. -----------------------------------------------------------------------
  5. Read Mesh2d.dat, Check if Multi-Grid can be used ...
  6.   Finished check Mesh2d.dat,  Most stage is            5
  7. Check bc2d.in ...
  8. Check multigrid OK
  9. Total Block number is            3 Total Cell number is        23552
  10. Most stage number of multi-grid is            5
  11. --------------------------------------------------------------
  12.            3   stage grids is used !
  13. -------------------------------------
  14. read Mesh2d.dat
  15. read Mesh2d.dat OK
  16. read bc2d.in ......
  17. read bc2d.in OK
  18. Creat Mesh ......           2
  19. Creat Mesh            2  OK Total Cell number is        5888
  20. Creat Mesh ......           3
  21. Creat Mesh            3  OK Total Cell number is        1472
  22. Kstep, t=           1  4.99999999999999977E-007
  23. ----------The Max Residuals are--------  ---Mesh---           3
  24.     0.7753608972E-02    0.7348449463E-02    0.5135622696E-03    0.4041801378E-01
  25.    The R.M.S Residuals are
  26.     0.8649557237E-03    0.8714159899E-03    0.6533435988E-04    0.4493792359E-02
  27. Kstep, t=           2  9.99999999999999955E-007
  28. ----------The Max Residuals are--------  ---Mesh---           3
  29.     0.7582504185E-02    0.6284729332E-02    0.1530098442E-02    0.3973969432E-01
  30.    The R.M.S Residuals are
  31.     0.8041464468E-03    0.8072630355E-03    0.2231385743E-03    0.4162616881E-02
  32. write data file ...flow2d-3.dat                                      
  33.   Prolong  to mesh            2    OK
  34. Kstep, t=           1  4.99999999999999977E-007
  35. ----------The Max Residuals are--------  ---Mesh---           2
  36.     0.4176128558E-02    0.3605486775E-02    0.2581921001E-02    0.2356941347E-01
  37.    The R.M.S Residuals are
  38.     0.2694712456E-03    0.3139314298E-03    0.2993138971E-03    0.1434276842E-02
  39. Kstep, t=           2  9.99999999999999955E-007
  40. ----------The Max Residuals are--------  ---Mesh---           2
  41.     0.3758138371E-02    0.2319293129E-02    0.2184189226E-02    0.2114269795E-01
  42.    The R.M.S Residuals are
  43.     0.2504025714E-03    0.2578244080E-03    0.2659229688E-03    0.1320697079E-02
  44. write data file ...flow2d-2.dat                                      
  45.   Prolong  to mesh            1    OK
  46. Kstep, t=           1  4.99999999999999977E-007
  47. ----------The Max Residuals are--------  ---Mesh---           1
  48.     0.1370238002E-02    0.6537118620E-03    0.8272982190E-03    0.8397708052E-02
  49.    The R.M.S Residuals are
  50.     0.6498282974E-04    0.4474184357E-04    0.8008293440E-04    0.3495755542E-03
  51. Kstep, t=           2  9.99999999999999955E-007
  52. ----------The Max Residuals are--------  ---Mesh---           1
  53.     0.1227673672E-02    0.8052771769E-03    0.8388357098E-03    0.7570845066E-02
  54.    The R.M.S Residuals are
  55.     0.6371551553E-04    0.4327602958E-04    0.7724725808E-04    0.3394780854E-03
  56. write data file ...flow2d.dat                                       
  57.   Initialize OK ......
  58.   Start ......
  59. Kstep, t=           3  1.50000000000000004E-006
  60. ----------The Max Residuals are--------  ---Mesh---           1
  61.     0.1009210805E-02    0.1033928020E-02    0.9534873513E-03    0.6233154383E-02
  62.    The R.M.S Residuals are
  63.     0.6171766215E-04    0.4463595643E-04    0.7576167456E-04    0.3249998190E-03
  64. Kstep, t=           4  1.99999999999999991E-006
  65. ----------The Max Residuals are--------  ---Mesh---           1
  66.     0.5857027438E-03    0.1411333515E-02    0.8428730499E-03    0.3576803841E-02
  67.    The R.M.S Residuals are
  68.     0.6137168812E-04    0.5281255867E-04    0.7364619711E-04    0.3235287821E-03
复制代码
 楼主| 发表于 2011-5-27 09:02:46 | 显示全部楼层
原帖由 shirazbj 于 2011-5-25 19:18 发表
我用1。1倒是能算下去。
从mesh3,算到2,再1,再1。最后只是在mesh 1上算了。我还以为多重网格要3-2-1-2-3-2-1。。。的轮着算呢。

1和1之间为什么有个  Initialize?

变mesh的时候,t重新开始了。mesh 1到1间 ...


这个算例是稳态计算。
最开始的初始化过程是:最粗的网格(网格3) 先算一定步数(如100步),然后插值到网格2, 算一定步数(如50步);然后再插值到网格1, 然后网格1算一定步数(如10步)。

初始化结束后,就进行多重网格计算,为V型多重网格(只显示最细网格上的残差)。
发表于 2011-5-27 09:43:24 | 显示全部楼层
李老师。J.Blazek书上边界条件的关于subsonic inlet
知道总温总压推出入口流动参数,可不可以给解释下
发表于 2011-5-27 16:50:43 | 显示全部楼层
原帖由 lixl-imech 于 2011-5-22 11:54 发表


哦,也许是程序中有Bug.
附件中的文档是OpenCFD-EC2D ver 1.12, 修正了ver 1.1中的几个Bug. 我刚调试过,应当可以运行。

运行方法:把附件中的文档解压开,仅需编译运行opencfd_ec2d_v1.12.f90即可。

李老师,谢谢你的解答!你的程序OpenCFD-EC2D ver 1.12我在实验室的几台机器上运行依然出现94#同样的问题(我用的是Compaq Visual Fortran Professional Edition 6.6.0,编译没有问题,计算一段时间就出问题),麻烦你再帮我找找原因,谢谢!

[ 本帖最后由 chenhuanlong 于 2011-5-27 08:52 编辑 ]
发表于 2011-6-3 19:41:44 | 显示全部楼层

我对OpenCFD-EC的几点认识

1,在modual Global_Var中声明全局变量时,应该要显式地加上common或者是save这样的修饰词,这样才能保证各个子函数都能准确地可见,而且不会随不同的编译器出现不同的问题。(参考Fortran95教材)
2,是个bug,在output_Res子函数中,有句write(99,"(I8,9E20.10)"),应为write(99,"(I8,8E20.10)")
3,对于程序中处理角点处的值,使用U_average_conner的原理是什么,请李老师给予说明,同时,再问下,若在三维下,U(0,0,0)这样的极端角点,又应该如何处理?
4,为什么不把几何信息在迭代计算之前,保存起来,这样可以更快地计算,虽然牺牲了些内存,但是现在的电脑的内存一般都能满足中等量得计算。这是我一般的做法,请问李老师之前是从哪个角度考虑的,处理内存与计算速度这对矛盾体的?想知道同样做CFD的其他人的思考方式。
发表于 2011-6-4 10:32:38 | 显示全部楼层

OpenCFD-EC 与 LES

李老师,您好,目前,我已经把您的OpenCFD-EC的程序看了一遍,想在OpenCFD-EC中加入LES子程序,这样可以用LES计算流体流动。目前想先试算下channel flow,之前看过您用DNS 算过channel flow,所以想就channel flow的具体计算的实施上问几个问题。诸如:初始条件如何给定和保证流动处于充分发展段,源项的问题?(是可压缩的channel flow),还请李老师点拨下!
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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