找回密码
 注册
流体中文网»论坛 › 教学科研 › CFD基础理论 › 发布OpenCFD-EC2D; 更新 Ver 1.5.4 ； 三维软件使用协议
1 ...234567891011... 24下一页
返回列表 发新帖
楼主: lixl-imech

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

  [复制链接]
 楼主| 发表于 2011-5-4 09:59:27 | 显示全部楼层
原帖由 fanghongfox 于 2011-5-3 10:44 发表

                               
登录/注册后可看大图

增加了多重网格后，计算速度提高了很多倍。李老师，是否考虑增加隐式推进格式LUSGS，LU－ADI等，进一步提高速度。

今后的版本应当会增加隐格式功能。
回复 支持 反对

使用道具 举报

ch06
发表于 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 | 显示全部楼层
谢谢李老师了。。。我们很期待你的成果。。。加油。。。
回复 支持 反对

使用道具 举报

nwpuzww
发表于 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型多重网格（只显示最细网格上的残差）。
回复 支持 反对

使用道具 举报

liurui
发表于 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），还请李老师点拨下！
回复 支持 反对

使用道具 举报

下一页 »
1 ...234567891011... 24下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

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