@李东岳 感谢大佬

@李东岳 感谢大佬:146:

@wuyukai 您好，为什么我用这个命令，就提示报错呢

Build : 4.0-268bb07d15d8 Exec : surfaceConvert ibCylinder.stl ibCylinder.ftr Date : Jul 12 2021 Time : 19:24:51 Host : dyfluid PID : 4225 CtrlDict : "/home/dyfluid/foam/tutorials/immersedBoundary/movingCylinderInChannelIco/system/controlDict" Case : /home/dyfluid/foam/tutorials/immersedBoundary/movingCylinderInChannelIco nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 段错误 (核心已转储)

各位朋友大家好，有个问题想请教一下。

我在用OpenFOAM 6.0做圆管湍流传热的大涡模拟，入口速度设定为timeVaryingMappedFixedValue。速度数据是预先对一个短管进行进出口周期性计算，存储绝热条件下某一截面速度作为数据库，施加到加热区域的入口。为了缩短预计算时间，减小数据库容量，一个办法是将数据库循环读入到入口边界，也就是需要对timeVaryingMappedFixedValue这个边界条件进行一些改进。

偶然间发现CFD-online上有一个这样的程序，是针对OpenFOAM2.3写的。我把程序做微小改动后植入到OpenFOAM6，发现只能单核计算，并行计算出错。https://www.cfd-online.com/Forums/openfoam-programming-development/228505-modified-version-timevaryingmappedfixedvalue-bc-looping-through-data.html

于是我参照这个程序对OpenFOAM6原始的timeVaryingMappedFixedValue程序进行修改，命名为loopedTimeVaryingMappedFixedValue。程序编译能通多，但是测试时候发现无法计算，错误信息如下：
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam::loopedTimeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::updateCoeffs() at ??:?
#5 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
浮点数例外 (核心已转储)

我的OpenFOAM编程基础非常薄弱，自己折腾了两天，还是搞不定这个问题，请大家帮忙看看这个错误的原因。如果有编程基础好的朋友愿意花点时间把CFD-online帖子中针对OF2.3的程序改成OF6能用的，我也非常愿意付报酬，可以QQ359199350。

第一次发帖求助，打扰大家了，非常感谢！

@enthusiasm 在 如何设置初速度 中说：

本人目前在做一个物体斜射入水的问题，需要计算物体入水后的弹道和速度变化，使用的是重叠网格。
我的疑惑是如何定义入水时的初速度，众所周知，有些方法可以设置恒定的速度或角速度，可这与我的问题模型不符，因此诚恳请教各位如何设置初速度，这个初速度是不是应该在 6DOF 模型里设置呢

请问一下师兄怎么设置？我最近再算一个高速入水，之前用的商业软件。同样的困惑。如有解答十分感谢！

老哥你是哪个单位的？我也遇到了这个问题，一起交流交流？

@李东岳 在 interfoam控制方程的问题 中说：

@cfdngu http://dyfluid.com/interFoam.html 不是方程26么 这里面的推导跟nonNewtonianIcoFoam是一样的 那个更简单 你可以看看

好的，非常感谢东岳老师

@五好青年 你这个想法方法我个人想到了入口设置非均匀场解决，
如U文件入口场设置为非均匀，之后可以对每个网格节点单独指定速度，这个操作OpenFOAM一些自带算例里边就有。
个人意见，仅供参考，希望能对您有所帮助

@李东岳 谢谢老师，非常感谢:140:

@qingdong_wang Info << phi即可

@bestucan : :mihu: 任重道远啊我

@浪迹天大 谢谢带佬，我才知道原来github上也有机理:143:

@李东岳 谢谢东岳教授:140:

@shaneheee 这一块我涉及挺多的，希望可以帮到你

@qjh888 foam-extend中的滑移网格，不能并行计算吗？其他版本的滑移网格可以并行哎？你的为啥不能呢，能具体说下么

@winsway_zero 请问您的算例是否调好了呢？边界条件是否能分享一下呢？

@flyzx95 嗯 所以你直接忽略就行 另外OpenFOAM-v2106发布了

@李东岳谢谢东岳老师

谢谢@李东岳 和@zhouxu 两位回复，感觉这一项应该就是zhouxu所说的，对Cv进行相体积分数加权的调和平均了

@冠竹 你看下dnsFoam算例的tutorials，需要用boxTurb设置无三度初始场。如果你有更好的公式，用setFields也行