@赵一铭

感谢分享！最近安装低版本的OpenFOAM做原有的程序二次开发，用到了这个语句~~，

在此分享一个基于不同Linux系统的OpenFOAM不同版本安装总结，亲测OF211，OF231 再 Ubuntu 16.04 安装成功无Error！

https://openfoamwiki.net/index.php/Installation/Linux

再次感谢！

@东岳 谢谢老师，太久没看，给忘了

@Cp_Zhao 明白了，感谢指点！！

@gyzhangqm 是否有足够的数据支撑？有的话可以。

这是div(phi,T)与-div(phi,T)等的解释

微信图片_20200320103052.jpg

没有simpleR，http://dyfluid.com/simplefoam.html

可以直接在of中sample提取数据

@Samuel-Tu Lib那改到Foam_USER_Libbin

1.JPG
2.JPG

AdjointOptimisationFoam user Mannual 链接：
https://www.openfoam.com/documentation/files/adjointOptimisationFoamManual.pdf
建议使用之前看一下连续伴随方法中的相关公式推导、伴随变量的边界条件推导以及敏感度分析，目前仅有伴随S-A模型，适用于气动外形优化。目前OpenFOAM中实现连续伴随方法的主要难点在于动量伴随方程中ATC（adjoint transpose convection）项的处理，由于该项会导致矩阵刚度较大，如果处理不当会导致伴随方程求解发散。运用伴随方程可以方便的进行各种目标函数和约束的敏感度分析，包括表面敏度和体积敏度，分别对应形状优化和拓扑优化。欧洲相关学者已经研究得非常透彻了，未来应该会在openfoam新版本慢慢放出。已经商业化的是engys公司旗下的Helyx-adjoint产品。

@队长别开枪 队长您好，您可否帮我解答一下这个帖子里问的问题，就是为什么snappyHexMesh画网格始终无法移除固体，感谢队长！
http://www.cfd-china.com/topic/3408/snappyhexmesh为什么总是无法移除固体

就，就是普通的动网格，dynamicMesh

画网格过程中出现这个提示是什么意思
1.png

@东岳 好的，谢谢老师！

@东岳
我现在是对所有网格做遍历，遍历每个网格上的面，判断是内部的面还是边界的面，如果是边界的面的话，通量值就会和另一块processor上对应的网格上的值相关，目前就是不知道怎么去找这个网格。

@dzw05
不好意思, 刚仔细看了一遍看到了, 谢谢!

@东岳
好的我试一下，谢谢李老师

如下图，再同样的位置上，我的模拟的速度再0.01以后就远小于实验数据。
受力仅考虑drag和gravity。stochastic dispersion model。
但是也有实验室同事写的collision model（已经在他们的毕业论文中验证过了）

2-way coupling
flow和实验数据匹配。
injection 中的parcel 的size 和速度也和实验数据匹配。

问题就是在径向一段距离以后，parcel的动量损失很多。
不知道大佬们有没有什么相关经验？

4ebe1f80-272c-4107-ab8d-ae747f0e738f-image.png sudu 的

@东岳 谢谢东岳老师！

@羽之下 这貌似是一个类型转换函数，搜到了类似帖子https://bbs.csdn.net/topics/390537437

turbulenceProperties中object为RASProperties或turbulenceProperties有什么区别？

参考Hirsch 2007: 12.4.1 Basic Approach of Pressure Correction Methods

请问OpneFoam复杂模型的多孔介质流动是否有成功的案例或者相关资料供参考。
:xinlei: 心累，我啥也不懂，望各大神各显神通，感恩戴德。
邮箱：2685306635@qq.com。:zoule:

@东岳 谢谢东岳老师:ok3:

@东岳 谢谢东岳老师

大家好，请问各位老哥有使用过salome前处理的吗，最近打算研究一下，可以交流交流:chouchou:

@东岳 李老师，我试了一下，把出口和入口速度设置的速度一样，在开始的时候水池的水就不稳定，在入口处水面上升，在出口处水面下降。得好久才能接近一种平衡，而且导致液面晃荡。
1.png

2.png

哦哦，那应该就对网格没什么影响，可能是Pointwise这个软件里面才有的设置，谢谢李老师:chitang:

插值从表格里面做热物性吧？那估计是超临界啥的。算例本身算起来可能是会有点不好收敛。你层流验证的时候好收敛么，还是也很难收敛？如果层流很好收敛但湍流不收敛，那大概率是湍流的边界问题；如果层流本身也很难收敛，那就可以去找别的原因了。其实你可以试试，同样的设置下把热物性换成普通的（理想气体之类的）看看能不能算，这样也能排除一些问题。
单看log文件有两点我觉得值得考虑。1、你设置了pressureControl，限定了压力的最大最小值，你本身的工况是在这个压力范围之内么? 2、这个可能性比较低。你可以把SIMPLE的convergence criteria调小一些，看起来你的速度方程解了1步就收敛了，对线性求解器来说好像太容易了点。

零梯度边界，边界上的值一直等于相邻网格单元的值。

@刘雄国 请教一下, OpenFOAM 的结果怎么一步一步的打印出来那?

问题找到了，是newFunction1下的TrainingExp.H文件中TypeName(" ");的内容没对上。
表现为：可用function1类型中包含相似名称的类型，边界条件与该类型存在字符偏差。
相关链接：https://www.cfd-online.com/Forums/openfoam-programming-development/133466-unknown-patchfield-type.html

看一些边界条件，都会定义多个构造函数，每个构造函数的输入量不同，为什么需要多个构造函数呢？

@星星星星晴 在 不知道为什么会生成特别蛋疼的field 中说：

我的cell值都没啥事，但是point的值就很诡异了

那可能是边界条件的问题，

@东岳 今天刚好在看limitedLinear格式，后面的数字貌似指的就是这一串，max（min（2r/k，1），0）

给合作方写的文档，基本写完了，但需要跟合作方沟通各种事宜之后才能放出来，预计3月底。写成英文的，让老外也看看，然后去我的网站，看更多的中文内容，钓鱼:high:

目前看过的几个教程会拿最典型的案例进行操作，当我想要更详细地了解某个求解器（比如simpleFoam）时，虽然找到了tutorials下相应的案例，但每个案例都是特殊的，不只是单纯的blockMesh+simpleFoam。其他的文件如何使用是个问题，更谈不上选择和自己感兴趣的内容相近的案例进行改造。官网上似乎没有直接找到相关的案例说明，所以想请教一下看有没有具体的案例说明和操作。

@王慧博 看样子你是在计算超临界水方面的。请问这个问题您是如何解决的？我在计算圆管等热流加热算例，采用kwSST和kEpsilon模型，尝试了buoyantSimpleFoam和rhoSipleFoam都发散，下面是我用rhoSimpleFoam计算kEpsilon的错误信息，不太明白是什么导致的发散或错误。

#0 Foam::error::printStack(Foam::Ostream&)[40] #0 Foam::error::printStack(Foam::Ostream&)[39] #0 Foam::error::printStack(Foam::Ostream&)[38] #0 Foam: :error::printStack(Foam::Ostream&)[37] #0 Foam::error::printStack(Foam::Ostream&)[36] #0 Foam::error::printStack(Foam::Ostream&)[35] #0 Foam::error::printS tack(Foam::Ostream&)[34] #0 Foam::error::printStack(Foam::Ostream&)[33] #0 Foam::error::printStack(Foam::Ostream&)[32] #0 Foam::error::printStack(Foam::Ost ream&)[31] #0 Foam::error::printStack(Foam::Ostream&)[30] #0 [11] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[12] #0 Foam::error::printStack(Foam::Ostream&)[29] #0 [13] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[14] #0 Foam::error:: printStack(Foam::Ostream&)[28] #0 [15] #Foam::error::printStack(Foam::Ostream&)0 Foam::error::printStack(Foam::Ostream&)[16] #0 Foam::error::printStack(Foa m::Ostream&)[27] #0 Foam::error::printStack(Foam::Ostream&)[25] #0 [17] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[2 6] #0 Foam::error::printStack(Foam::Ostream&)[18] #0 Foam::error::printStack(Foam::Ostream&)[21] #0 Foam::error::printStack(Foam::Ostream&)[23] #0 Foam::e rror::printStack(Foam::Ostream&)[24] #0 Foam::error::printStack(Foam::Ostream&)[19] #0 Foam::error::printStack(Foam::Ostream&)[20] #0 Foam::error::printSta ck(Foam::Ostream&)[22] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [36] #1 Foam::sigFpe::sigHandler(int) at ??:? [41] #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libpthread.so.0" [29] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libpthread.so.0" [28] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? [37] #4 Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > co nst&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:? [36] #4 Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > co nst&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:? [35] #4 Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > Foam::operator*<Foam::volMesh>(Foam::tmp<Foam::DimensionedField<double, Foam::volMesh> > co nst&, Foam::DimensionedField<double, Foam::volMesh> const&) at ??:? [36] #5 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??: ? [37] #5 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??: ? 5 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct()Foam::RASMode ls::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct()[41] #5 Foam::RASModels::kEp silon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:? [30] #5 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??: ?

这里第一个错误应该是除数错误，下面的模型错误是什么导致，如果边界有问题的话也不会影响湍流模型？以及tmp（DimensionedField）是什么错误，也是边界错误吗？求助

矢量用黑体才对，我再查一下资料，这个公式确实是错的

我想到了，用fvc:;grad()&vector(0,1,0）求Y向的分量

:xinlei: :xinlei:

dear cfders。

遇到一个问题，我用piso瞬态算一个流场，当速度差不多稳定后，利用现在这个时刻的流场，稳态地求解一个浓度f场，这时又需要用simple稳态方法。simple完后，又要用piso重新计算流场，然后进行类似循环。
有没有可能把simple放入Piso中，先是算Piso，然后满足一定条件（用 if ?）后跳到simple中去。如下：

pisoLoop() { pisoFoam类似代码; if() { simpleFoam类似代码； } }

这样不知道fvSolution里面怎么设置。。

如果上述方法不行的话，我在想把稳态simple的浓度f场计算也改为piso，也就是说直接在Piso中加入浓度f场的微分方程就行了。这样是因为反正最后速度场是比较稳定的，那么受速度驱动的浓度f场最后也应该稳定（无内部源项，且inlet输入恒定），即使我用piso也没事。不知道这样想法可不可以。

请大家不吝赐教

只是一种压力降边界条件，see fanPressure

@东岳 谢谢

好像发现原因了，Newtonian没有从singlePhaseTransportMode继承，因此如果在singlePhaseTransportMode中定义虚函数是不行的

@Samuel-Tu 原来这样，简单的重构方法好像确实有一个代价就是网格要比较细，不然不准确。

加油加油

是的，都是绝对坐标系的

最近在家呆的闲，又没有实验室的资源，，想看看磁流体的流固耦合的代码和书最好中文的，，，求推荐啊～

滑移壁面没有剪切力