OpenFOAM

• U

  twoPhaseEulerFoam中的ITAE模型
  • upc_ngh  

  13
  13
  Posts
  319
  Views

  U

  晚上又测试了一下，在linux系统下是可以顺利的进行计算，谢谢东岳老师的帮助
• J

  Chang-Hsieh-Chen低雷诺数湍流模型边界条件设置
  • jinlinna  

  9
  9
  Posts
  185
  Views

  J

  @东岳 李老师，这样计算之后我发现计算结果不对，通过查找资料我发现这个边界条件可以用 来代替，所以我用codedfixedValue编写了这个边界条件，现在的问题是我好奇codedfixedValue里面的value值必须要设置成$internalField吗？谢谢！
• Y

  LES拉格朗日粒子速度时间平均
  • yfclark  

  4
  4
  Posts
  100
  Views

  Y

  速度平均还是不会做，不过发现湍流强度对计算结果影响很大，大家算两相流的注意一下
• 小

  水中比内能与温度的关系
  • 小考拉  

  3
  3
  Posts
  50
  Views

  

  比内能就是e啊，compressibleInterFoam顶层就有啊
• C

  反应速率，摩尔分数问题
  • CarmeloSun  

  1
  1
  Posts
  46
  Views

  C

  请问大家，oopenfoam里怎么输出多步反应的每一步的反应速率，还有怎么输出摩尔分数，谢谢啦
• 炎

  振荡来流的流场倒流问题
  • 炎冰灬灬  

  24
  24
  Posts
  227
  Views

  M

  @炎冰灬灬 检查网格看看呢，你也可以参考下这个帖子的设置 https://www.cfd-china.com/topic/664/cd计算不准的问题-les-re3900-三维圆柱绕流-pisofoam
• R

  Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam
  • random_ran  

  121
  121
  Posts
  78689
  Views

  炎

  @random_ran 你好。看了你的升阻力系数的图，我最近也出现了同样的问题。我算的是二维的刚性翼型的绕流问题，用的是SST k-w模型，pisoFoam求解器，算出来的升阻力系数的波动也是非常的大，不知道问题出在哪里，你能指点一二吗
• 

  p_rgh迭代次数过多
  • 同学博  

  6
  6
  Posts
  127
  Views

  

  @同学博 哈哈 好的
• 

  CFDEM结果用paraview后处理，流场与颗粒显示不同步
  • .J.  

  5
  5
  Posts
  1252
  Views

  J

  @linhan-ge 请问temporal shift scale 里pre shift， post shift，scale三个参数取值与算例的时间步和总时间的关系是？
• Z

  请教把OF数据导入MATLAB的方法
  • Zhixing18  

  3
  3
  Posts
  123
  Views

  Z

  @bestucan 先说声谢谢，我去仔细看看
• H

  openfoam中的omega初始化与fluent中的Omega初始值计算公式不同吗？
  • hery  

  3
  3
  Posts
  101
  Views

  H

  @东岳 用的komega湍流模型，有时在outlet处设置了intOut边界，value会设置成0，也尝试过zerogradient，没有value ，0文件夹的omega设置如下： internalField uniform 53.6; boundaryField { INLET { type turbulentMixingLengthFrequencyInlet; mixingLength 8.52e-4; value uniform 53.6; } OUTLET { type zeroGradient; } BOTWALL { type omegaWallFunction; value $internalField; } UPWALL { type omegaWallFunction; value $internalField; } HEATWALL { type omegaWallFunction; value $internalField; } }
• 

  关于OpenFOAM中K-Omega SST模型使用壁面函数时的Y+问题
  • tidedrinker  

  11
  11
  Posts
  2219
  Views

  H

  我后来的做法是将网格做到y+<1，然后不采用壁面函数，k方程壁面边界条件为kqRwallFuction, omega边界条件为omegaWallFuction, nut边界条件为nutkWallFunction 另外，我觉得将网格做到y+>30，然后指定壁面函数也是可以的 同样有疑问，指定了k w 和nut 的边界条件不就是采用了壁面函数吗？ 还是说壁面函数有什么另外的指定方法？ （还是我认知错误）
• 炎

  pisoFoam计算圆柱绕流出现拍现象问题
  • 炎冰灬灬  

  1
  1
  Posts
  70
  Views

  炎

  用pisoFoam计算了两百雷诺数的层流圆柱绕流问题。但是计算出来的升阻力系数出现了拍现象，有大神可以指点一下是什么问题吗。以下是升阻力系数图和算例设置
• F

  RANS动量方程中雷诺应力表达的问题
  • fangyuanaza  

  4
  4
  Posts
  101
  Views

  F

  @东岳 懂了~ 谢谢智慧的李老师~~
• S

  Tecplot无法显示内部流体
  • Stan同学  

  7
  7
  Posts
  158
  Views

  S

  @benqing 好的 谢谢您！ 我所用的求解器是olaDyMFlow。
• H

  有没有utility或者force，能使颗粒在RANS模拟中像LES那样发散开来
  • hotairballon1997  

  7
  7
  Posts
  173
  Views

  H

  抱歉没看到回复提醒，就是报错三楼贴出的问题，不过后来发现是我之前把DPMFoam里面的湍流模型调用修改了，导致出了问题。OpenFOAM里本来的DPMFoam是没有问题的~
• 炎

  pisoFoam和pimpleFoam计算出的升阻力系数波动过大问题
  • 炎冰灬灬  

  3
  3
  Posts
  89
  Views

  炎

  @东岳 刚性翼型用的pisoFoam，柔性翼型用的是fsiFoam。但是pisoFoam算出来的系数波动很大，而fsiFoam算出来系数波动就很小。不知道是为什么
• U

  topoSet设置半圆柱区域
  • upc_ngh  

  1
  1
  Posts
  79
  Views

  U

  各位前辈，如下计算区域模型，想要使用chtMultiRegionHeater求解壁面导热+内部空气水流动传热问题，怎么使用topoSet描述空气或水的半圆柱流道区域？对于空气和水的交界面应该给什么边界条件，还是compressible::turbulentTemperatureCoupledBaffleMixed这个边界条件吗?
• S

  关于cfdem-interFoam solver
  • shincoast  

  10
  10
  Posts
  4254
  Views

  D

  @zhuyuxuan 您好，请问这里“已有的例子只提供了耦合后的icoFoam和simpleFoam求解器”例子是在哪里可以找到？
• S

  关于sample fields
  • starDust  

  2
  2
  Posts
  95
  Views

  鲸

  现在解决了吗 我现在也是用toposet建立了cellzone 但是也是报错没有对应的cellzone
• K

  虚拟机openfoam安装出错，求助
  openfoam 安装 • • Kingkong  

  2
  2
  Posts
  73
  Views

  K

  已解决，虚拟机内存调到16g成功安装
• 鲸

  cellzone的问题
  • 鲸落  

  4
  4
  Posts
  127
  Views

  鲸

  @东岳 好的 谢谢老师 我现在知道利用topoSet来指定cellZone但是那个是静止的 不知道石英砂是动态的应该怎么指定
• 郑

  为什么二维方腔流用pisFoam算出的流线与下方没有涡
  • 郑澄远  

  6
  6
  Posts
  630
  Views

  

  @东岳 二维顶盖驱动方腔流无滑移时有涡，有滑移时不出现涡，这种情况可能是什么原因造成的？程序是fortran语言自编的
• S

  OpenFOAM-6的twoPhaseEulerFoam学习
  • Samuel-Tu  

  9
  9
  Posts
  348
  Views

  

  @Samuel-Tu kI项是雷诺应力项对照柯西应力项中正应力部分写出来的，而p正好是柯西应力中的正应力部分，所以两者放在一起
• 波

  请教：三维涡量识别方法
  • 波流力  

  6
  6
  Posts
  269
  Views

  好

  @波流力 可否推荐一下上交那篇植入其他涡量识别方法的文章。我目前也在搞相关的工作，想借鉴一下。谢谢！
• C

  OpenFOAM中随时间变化的属性设置
  • ChimneyFoam  

  8
  8
  Posts
  284
  Views

  王

  @Cr 十分感谢，问题已经解决了
• C

  源项类中的addSup函数是如何被调用的？
  • chaoscfd  

  4
  4
  Posts
  131
  Views

  C

  @Cr 原来是在这里啊，感谢感谢！
• 

  请教：第三代涡识别方法Liutex的代码在linux中运行的问题
  • yfleung  

  7
  7
  Posts
  217
  Views

  

  @yfleung 呃，那就看看$64^3$网格的时候用几个节点呀。一个节点基本就是一个主机了。大机器不是很重的任务分一个节点也差不多了。 看64分到几个节点，如果也是一个呢。 调试阶段还是时间步越小越好。确定正常工作了再调大
• H

  oscillatingFixedValue边界条件的问题
  • hellofoam  

  8
  8
  Posts
  2932
  Views

  L

  好的，谢谢您
• S

  olaDyMFlow动网格计算不收敛，求助
  • Stan同学  

  3
  3
  Posts
  94
  Views

  S

  @Cr 你好 我是把消波边界调整为2D消波边界了 然后可以计算了。 但是提取出来的坝前波面结果很奇怪，我的坝体振幅在三个方向均只有0.4m（初始水深设置的是50m），而水深却会突然变成0。提取出的波面时程曲线如下，暂时没有找到是什么原因导致的这种现象，还请各位老师指教。
• 搬

  OpenFOM里面的是RANS还是URANS
  • 搬运工不好当  

  7
  7
  Posts
  2192
  Views

  D

  @hotairballon1997 和steady Rans结果不同，因为每个时间步的结果不同。
• L

  模拟火箭发射问题
  • Lisiye  

  4
  4
  Posts
  182
  Views

  

  看来还是有只用动网格而不用自适应的。这种状况应该用在变形的铁板、高层建筑随风摆之类的在一定范围内变形的。这样网格即使没有自适应也没关系。一般用来做应力分析之类的。 查看了下网格发现网格变形为负了 你的这个初始网格就不好。在圆形边界上没有用O形网格。畸变本就很严重，然后一变形，还不是震荡变形，是一直一个方向走，网格又没有自适应，然后就负网格了。 我同时想在INLET（火箭喷管口）设置燃气射流环境，目前改了 0/p 0/U 0/T 下的边界条件，并在constant/thermophysicalProperties 设置了燃气参数，但不知道如何指定在INLET上不断喷射燃气。 最简单的把边界设定成固定压力梯度。要是还要计算燃烧过程就不懂了
• 小

  如何提取张量的某一个值
  • 小考拉  

  6
  6
  Posts
  129
  Views

  

  @小考拉 不忙不忙，就是不太懂 但是探索不懂的东西的方法都是相通的。你要解决的问题我没解决过，所以要解决花费的时间比较多。思路就是看代码。你也在看，但我们看的版本对不上，开源软件有个特点，太自由了，所以版本间的差距特别大。 在你的代码中找相关定义的方法： grep -R 'volTensorField' ./ 先用 cd 命令 切换到源代码所在的目录 ~/OpenFOAM/OpenFOAM-version/ 然后运行上面的那条命令。-R 是递归查找，包括子文件夹，'volTensorField' 是查找的内容，区别大小写， ./代表当前路径。 另外，在 ~/OpenFOAM/OpenFOAM-version/里查，查的太多。可以分别去这个文件夹下的src或者applications下面查。如果还是太多，就把输出重定向到一个文件再看。不过没颜色不太好看。 grep -R 'volTensorField' ./ > grep.log
• 小

  paraView-求展向平均问题
  • 小考拉  

  5
  5
  Posts
  175
  Views

  B

  @小考拉 我用的第二个 遍历网格单元的代码不知道你指的是什么意思 EXE = $(FOAM_USER_APPBIN)/spanwiseAvg这个就是编译一个程序，可以直接在终端调用，确实跟求解器的编译方法是一样的，其实求解器也是一个程序，类似于win里面的.exe的应用程序，只不过在这里需要通过终端调用
• 悬

  fluent总能比of得到更满意的结果？
  • 悬浮的猪  

  43
  43
  Posts
  3379
  Views

  J

  @东岳 谢谢李老师的回复！我确实修改了湍动能和湍流耗散率的初始值，得到了不同的结果，但我疑惑的是湍动能和湍流耗散率的初始值设置都是调试出来的吗？ 目前我找到的大家都说用湍动能经验公式估算这个值，但用这个公式计算出来的往往与文献中的结果不符。那如果是自己做研究，而不是验证文献，又该如何给定这个值呢？
• J

  传热撞击流模拟
  • jinlinna  

  2
  2
  Posts
  90
  Views

  J

  甚至我模拟的文献和李老师算例所用的文献都是同一位作者，且几何模型几乎一致，不知道是不是湍动能的初始值设置有问题
• K

  The influence of mesh on solution?
  • kimy  

  3
  3
  Posts
  114
  Views

  K

  Thanks. I created two meshes with the same body size but different type of meshes,namely structured and unstructured, the results are different. Note that both of two cases are grid independent. Do you have any idea about that ? since I saw a thread you opened in cfd online.
• B

  reactingFoam, 如何计算每个基元反应的反应速率
  • buxiangle2012  

  2
  2
  Posts
  1088
  Views

  C

  您好，请问您最后解决了嘛，怎么样输出每个基元反应的速率呀
• 悬

  调用新编译的湍流模型进行运算，只计算了一步，出现了这样错误提示，有老师碰到过吗？
  • 悬浮的猪  

  10
  10
  Posts
  291
  Views

  悬

  @King_RVM 没错的问题就在乱用了tmp指针，我已经解决了。关于tmp指针，我在网上搜索到过问题，祝好！
• 闻

  paraFOAM时的dlopen error
  • 闻久STU  

  7
  7
  Posts
  160
  Views

  F

  @东岳 是的~ 这样就不报错了 谢谢李老师~~
• V

  DPMFoam设置求解问题
  • veen  

  7
  7
  Posts
  195
  Views

  V

  @东岳 李老师，有没有可能在DPMFoam里边设置每步生成随机的点云啊？我看patch injection是根据面网格密度生成的点云，而且网格越密的地方粒子越多。目前我是用多个点云injection去替代的，我在想有没有方法可以实现，自定义点云的分布密度呢？
• 一

  关于调用 库函数的报错
  • 一颗鸭蛋  

  7
  7
  Posts
  190
  Views

  一

  我当时是由于湍流模型虽然编译通过但其实还是有问题的，具体就是H文件申明的虚函数divDevRhoReff，在C文件中没有具体的实现。你可以根据H文档最后几个虚函数（virtual函数）再好好对照一下看看C文档中是否都有对应的具体实现。其次我是2.3.1版本。
• L

  提取场内一个面的物理量 Function Obj 的方法
  • leilei  

  1
  1
  Posts
  63
  Views

  L

  请教大家： 提取场内一个面的物理量 Function Obj 的方法 谢谢 另外： 本人原帖 “湍流模型问题请教” 得各位相助， @闻久STU 君 "回访"，李博 @东岳 赐"西语句"，为不在原贴用回谢灌水，在此一并谢过
• 

  修改湍流模型，定义表达式的问题
  • King_RVM  

  17
  17
  Posts
  348
  Views

  

  @Cr 非常非常感谢你的分析和回复，我的编译通过了，编译通过之后运行算例有新的问题，原因是tmp这个内存管理小工具用得不好，我最新的通过编译的部分修改如下，供大家参考
• L

  湍流模型问题请教：
  • leilei  

  10
  10
  Posts
  329
  Views

  

  Just an English comment passing by
• 我

  CFL数如何计算？
  • 我是河滩  

  6
  6
  Posts
  385
  Views

  Y

  谢谢，学习到啦！
• 

  关于传热边界条件的问题
  • win  

  2
  2
  Posts
  132
  Views

  C

  这个求解器是不可压缩的求解器，他单独求解温度方程，没有用用到热物理模型库，所以就不支持这个边界条件，fluidThermo是热物理模型库的方法。 你可以去看看buoyantPimpleFoam或buoyantSimpleFoam下面的算例，这些算例的constant文件夹下都有thermophysicalProperties这个字典，这些求解器也都用到了热物理模型库求解焓方程，并且支持这个边界条件。
• 

  externalCoupledTemperature 怎么用在不可压缩solver中？
  • sibo  

  4
  4
  Posts
  73
  Views

  C

  或许可以将代码改为用不可压缩求解器求解的边界上的温度和导热率来获得，边界上的传热量。例如： if(db().foundObject<volScalarField>("T")) { const volScalarField& T = db().lookupObject<volScalarField>("T"); const fvPatchScalarField& Tp = T.boundaryField()[patchi]; qDot = kappa*Tp.snGrad(); } 其中kappa需要根据你的物性和边界上的温度场来计算获得
• 星

  jet 的边界条件如何选择？
  • 星星星星晴  

  13
  13
  Posts
  294
  Views

  星

  @hotairballon1997 算吧 不过不是map，我个人觉得不要轻易用map这个人东西。。因为网格密度不同，map出来的场可能是有问题的
• 

  记录一些自己用过的代码
  • 李东岳  

  6
  6
  Posts
  3380
  Views

  

  亲测实用