fsiFoam自带算例的flowModel都是consistentIcoFlow，是计算层流的。修改flowModel为pisoFlow后可以计算湍流，请问有没有大神知道改为pisoFlow后，fvSchemes和fvSolution文件中，是否要做一些相应的修改呢。目前我使用pisoFlow后的fvSchemes和fvSolution文件如下，：

@UPC_王骅钟 链接：https://pan.baidu.com/s/1me5Nz_C36_6HxjWbC7Prcg 提取码：8o40

请问你解决了吗，我最近也遇到了同样的问题。我不是刚开始出现的，我是计算到三分之一的时候出现的这个问题

我是初学openfoam的研究生，目前研究的方向是 管内湍流拟序结构的研究。一开始一直用的是ansys.现在在运用周期性边界条件时，出现错误。能帮我看看我是不是constant设置错误。非常感谢！这次网格是使用icem非结构网格然后转换，管内无扰流器（光管）

我猜一下。 如果我写伪2D的求解器，我就把密度设成体积相关。这样就保证质量守恒了。这样体积怎么变，单位网格的质量不变。质量不变，力不变，单位时间内冲量、动量效果就一样。 比如我认为读入的密度是当厚度为1的密度。具体计算中设定厚度是多少，我再换算成特定厚度时的密度。 而且这种变化只影响与坐标相关的，比如扩散，但是这种变化只在z轴上，z轴上的变化又是不必关心的。 兴许扣扣那个求解器的密度当2D的时候怎么定义的能找出答案

interFoam是直接模拟，不能模化升力。 双流体模型是模型化后的模型，可以模化升力。 因此 升力的方程表达式中出现气相与液相速度，这两者速度与interFoam控制方程中的速度有什么关联呢？ 二者毫无关系

顶一下

代码导进去一都是错误怎么回事大家有遇到这个问题吗

非常感谢大佬们的对话，受益匪浅啊

@luofq-sysu 目前版本也没法用，应该也得自己编译这个网格转换命令。

请问楼主这个问题解决了吗?

用fsiFoam计算湍流，将流体求解器由constentIcoFlow改为pisoFlow，再添加湍流相关文件后开始计算。计算时间设置为6秒，时间步为0.00005，刚开始计算时很容易收敛，当计算到2秒的时候，每一个时间步都要迭代到最大迭代步数才会进行下一个时间步的计算，有什么方法可以改进吗？以下是计算过程迭代截图、流体求解器设置、耦合求解器设置以及固体属性截图。

用了两年OpenFOAM了，感觉至今感觉还没太看懂icoFoam中的东西，慌得一批

最近在做一个高分子溶胀-溶解的项目，里面涉及到固体的溶胀、固体的溶解这两个过程，万能的CFD中文网，有没有做这方面研究的CFDer？ 在知网上搜了一下模拟糖块在水中溶解的文章，发现竟然没有，是这个太简单了没人整么？

@bestucan 好的，谢谢

看了看我的程序 ，就是你说的样子，气相是wall，不过给了颗粒一个Y向速度，还有Y的进口体积分数，不知道算不算是

@yxg2018 谢谢，我将网格按照单元类型分离并将边界层网格也设定为刚体运动区域，计算依旧会报错，提示Divergence detected in AMG solver: x-momentum，我将迭代次数增加并改小了动量的松弛因子，计算还是发散。

你这是在耦合计算中没有找到transportProperties这个字典文件。可以仔细看一下，耦合计算和单场计算时，流体域是不是用的相同的求解器。

大家好，遇到一个动网格入口速度边界的问题，希望大家帮帮忙。 简述： 在一个做六自由度运动的网格的一个边界上定义随时间及高度变化的速度入口。 平台：OpenFOAM 2.4.0, interDyMFoam 物理背景： 做俯仰运动的明渠中的流动，在左边界施加随深度变化的平行于底床的入流速度，考虑到明渠本身也在做俯仰运动，所以入口边界的绝对速度还应该加上俯仰运动所产生的速度。所以就是一个随时间及高度变化的速度。 描述：网格的运动通过solidBodyMotionFvMesh 实现，然后通过6DOF的utility来指定运动方式。对于时变或随位置变化的入口边界条件，在论坛里面找到了一些帖子，然后用的codedFixedValue来实现这个，但是在实现的时候遇到错误，自己编程不是很熟悉，想请大家帮忙看看代码的编写上的问题。同时对于在动网格中的入口速度边界的指定，不知道我目前采用的这种方法是否合理？是否还有其他更好的方法？谢谢大家。 inlet { type codedFixedValue; value $internalField; redirectType abcde; code #{ const fvPatch& boundaryPatch = this->patch(); vectorField v= boundaryPatch.Cf(); scalar ymin=min(this->patch().Cf().component(1)); Info<<"ymin="<<ymin<<nl; scalar thetaamp=0.056119; scalar omegaa=1.4564; scalar pitchphase=2.0125828; scalar tt = this->db().time().value(); scalar udeck=-thetaamp*omegaa*50*cos(omegaa*tt+pitchphase); scalar uslope=0.09769*sin(1.488*tt-2.748)-0.00532; scalar ucut=0.2045*sin(1.692*tt-3.203)-0.0601; scalar theta=thetaamp*sin(omegaa*tt+pitchphase); forAll(boundaryPatch,faceI) { scalar y = this->patch().Cf().component(1); scalar ux=-udeck*sin(theta)+(uslope*(y-ymin)+ucut)*cos(theta); scalar uy=udeck*cos(theta)+(uslope*(y-ymin)+ucut)*sin(theta); v[faceI] = vector(ux,uy,0); operator==(v); } #} } 错误信息如下： // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solid-body motion function tabulated6DoFMotion Applying solid body motion to entire mesh PIMPLE: Operating solver in PISO mode Reading field p_rgh Reading field U Using dynamicCode for patch inlet on field U at line 26 in "/home/min/OpenFOAM/min-2.4.0/run/pitchDamBreak/pitchCases/pitch14-8update3-mapped/0/U.boundaryField.inlet" Creating new library in "dynamicCode/pitchdam/platforms/linuxGccDPOpt/lib/libpitchdam_89672c10695eb4de1d809896f609f8e1131d9634.so" Invoking "wmake -s libso /home/min/OpenFOAM/min-2.4.0/run/pitchDamBreak/pitchCases/pitch14-8update3-mapped/dynamicCode/pitchdam" wmakeLnInclude: linking include files to ./lnInclude Making dependency list for source file fixedValueFvPatchFieldTemplate.C /home/min/OpenFOAM/min-2.4.0/run/pitchDamBreak/pitchCases/pitch14-8update3-mapped/0/U.boundaryField.inlet: In member function ‘virtual void Foam::pitchdamFixedValueFvPatchVectorField::updateCoeffs()’: /home/min/OpenFOAM/min-2.4.0/run/pitchDamBreak/pitchCases/pitch14-8update3-mapped/0/U.boundaryField.inlet:53:43: error: cannot convert ‘Foam::tmp<Foam::Field<double> >’ to ‘Foam::scalar {aka double}’ in initialization make: *** [Make/linuxGccDPOpt/fixedValueFvPatchFieldTemplate.o] Error 1 --> FOAM FATAL IO ERROR: Failed wmake "dynamicCode/pitchdam/platforms/linuxGccDPOpt/lib/libpitchdam_89672c10695eb4de1d809896f609f8e1131d9634.so" file: /home/min/OpenFOAM/min-2.4.0/run/pitchDamBreak/pitchCases/pitch14-8update3-mapped/0/U.boundaryField.inlet from line 26 to line 68. From function codedBase::createLibrary(..) in file db/dynamicLibrary/codedBase/codedBase.C at line 213. FOAM exiting 谢谢了。

测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试 测试

这几天发出去几十个邀请码，大部分都是通过1.1、1.2条的方式获取的，好多老师的研究都很高端，学习！

@东岳 感谢李老师！我去找一下相关内容！

@winsway_zero 多谢！

或许CATIA的多截面曲面可以实现这个的建模。

你可以试试pointwise，也是一款优秀的CFD网格划分软件，结构化网格需要自己从点线面体一层层的实现，质量与个人水平有关，祝好运

我好像发现什么了，这个结构是一样的，不过tutorial case里obj生成polyMesh中，有分组。。。 ( aerofoil { type wall; inGroups List<word> 1(wall); nFaces 120; startFace 31760; } inlet { type patch; inGroups List<word> 1(freestream); nFaces 200; startFace 31880; } outlet { type patch; inGroups List<word> 1(freestream); nFaces 160; startFace 32080; } back { type empty; inGroups List<word> 1(empty); nFaces 16000; startFace 32240; } front { type empty; inGroups List<word> 1(empty); nFaces 16000; startFace 48240; } ) 所以，/0下文件头上freestream不是类型名，是边界（组）名。。。 感觉这个地方有误导性，导致我以为文件结构不同。。。 建议改个组名！

@星星星星晴 前辈您好，我再问一下，您的cloudFunction是怎么写的么？最好是截个图，谢谢

@浪迹天大 您好老师，有没有关于LeMOS的参数设置的计算方法，或者文献，感谢老师

@BlookCFD 非常感谢

各位大佬，有可以分享一下有关模拟平板作正弦柔性波动的udf吗？需要用到define_grid_motion 这个宏，初学udf，可参考的案例实在太少了。在此表示不胜感激，小弟邮箱为：815000514@qq.com 再次表示感谢！

@程迪 大佬可以详细的解释一下cyclic的计算过程中吗？比如这个CL和CU到底对应哪一项，initMatrixInterfaces怎么计算了Cl*x0.

@BlookCFD 谢谢，这个问题已解决，这个主要涉及到interFoam系列算法子循环求解相连续性方程并修正质量流率，可以将surfaceScalarField rhoPhi在循环外声明创建，但是循环内部每次循环初始值要赋0.

@东岳 好的，谢谢

自己回复一下，这里能选的都是标量，而速度本身是矢量，所以没有，想显示出速度值，就是通过计算器，自己计算出某个方向的速度（X，Y，Z），之后才可选~

奇怪，我在tutorials里面allrun了一个例子，paraview的时候就打开不了，会有error，说是polymesh/faces里面符号错误。。 按理说，tutorials里面的应该没问题啊。。我找了很久的原因，发现可能是二进制的问题，所以我把controlDict里面的输出从Binary改成了ASCII，结果allrun之后就可以paraview打开看见了。。 好奇怪，这是什么原因

@Samuel-Tu 具体忘了，网上搜下吧，可以设置每个软件调用独显还是集显

@xpqiu 感谢

测试编辑

我把代码移植到OF-7上了，并附上了我的算例以供参考： https://github.com/ZhangYanTJU/boundaryConditions

@bestucan 感谢回复！已经Allclean了 ，另外，不更改重新跑了一遍又没问题了，就是对出现的问题感觉比较迷。。。

速度的物质倒数两端乘以密度，再化简，就是时间项和对流项加在一起。 密度就是单位体积的质量，速度的物质导数就是加速度。 所以时间项和对流项加在一起就是单位体积内的$F=ma$，即$F=\rho a$。这就是惯性力吧？ 力分为体积力和表面力，本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了（我猜的），力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力，右边是流体微元真实受到的力，两者相等，就推出来流体微元怎么变化的了。 物理意义有助于了解方程怎么来的，但是后来方程形式因为代换化简早没有了当初的模样，只能从大概上理解。像广义的源项，这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程，运输的物理量分别是密度、速度、熵还是啥。源项，就像一个无源封闭体内进来多少就出去多少，但是有了源，就兴许进来的少出去的多。就是对运输产生了扰动。 要是搁到动量方程上，比如加个源项，磁流体在磁场中受电磁力，或者离心机里受额外的等效重力。 另外粘性造成的应力是和速度有关系的吧？叫本构关系，用来区别牛顿和非牛顿流体的。系数小不打紧，只要湍流涡够多，速度够快。耗散能量的能力就很强。

@cccrrryyy 非常感谢

坚果云竟然挂了？？？？

@xiaolin 我是因爲安裝了不同版本的OpenFOAM而導致的問題.比如說OpenFOAM6和OpenFOAM18.06,這種不同編號的版本是不同組織創建的.所以裏面有一些東西不兼容.你需要下載相應的版本.

@bestucan 谢谢你的回复，有所帮助的，不是猜的，您谦虚了啊。不过，我这个方法是比较特殊的MPPIC，跟普通的粒子追踪有些差别。我又翻了两天，还是发现，那个文献里提到的particle distribution function是不被OpenFOAM进行求解的。只是一个描述方程。欢迎您多多留言啊。很有帮助！

@马乔 您好，我参照你的方法将transportProperties代码写进去，但是很奇怪的是，报错显示，request for dictionary transportPropertie ,而可获得的字典里是transportProperties,两者总是存在单复数的差别（就像你贴的代码里一样）。本来想着直接用令锑度边界条件绕开这个错误，但是最近发现还是得用alphat壁面函数。真是很奇怪呀。

请问有人知道重叠网格里嵌套网格和背景网格二者之间的尺寸关系需要满足什么条件吗？另外，如果运动物体有薄壳结构（厚度0.02m），那在划分背景网格时尺寸是不是一定要小于0.02m？恳请回复！

很有可能，cyclicAMI有可能把粒子算丢了

测试出来了QAQ： 𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑠𝑜𝑢𝑟𝑐𝑒=(𝐵𝑖𝑛 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛 𝑐ℎ𝑎𝑛𝑔𝑒×𝐶𝑢𝑟𝑟𝑒𝑛𝑡 𝑝ℎ𝑎𝑠𝑒 𝑣𝑜𝑙𝑢𝑚𝑒 𝑓𝑟𝑎𝑐𝑡𝑖𝑜𝑛)/(𝐶𝑒𝑙𝑙 𝑣𝑜𝑙𝑢𝑚𝑒 ×𝑇𝑖𝑚𝑒 𝑠𝑡𝑒𝑝 𝑠𝑖𝑧𝑒) Eg: Bin 1 bin fraction decrease 2% in 0.1s, cell volume is 0.001 m3, inclusion phase volume fraction is 1% 2 is the number that we give fluent, 0.02 is actual bin fraction change.