好的，我认真看看，非常感谢您的回答，谢谢

谢谢东岳老师。
我的研究目的是看使用毛巾往复擦桌子过程中，由于毛巾的往复运动，引起的周围气流变化，毛巾速度是0.2m/s，如下图，先向左前进1s，再向右运动1s，回来。d49e6d5e-44a9-45c2-add6-7afe01a3bee3-image.png
这种匀速的往复匀速，速度边界条件应该如何给定？

谢谢您。

@zhe 直接设置水体区域，没有另外处理

dyfluid@dyfluid:~/twoWayHyQMOMFoam/tutorials/PTC_2D/testCase2$ oneWayHyQMOMFoam oneWayHyQMOMFoam: symbol lookup error: oneWayHyQMOMFoam: undefined symbol: _ZN4Foam6HyQMOMC1ERKNS_7PtrListINS1_INS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEEEE

原因：当时调整了src下的的构造函数，紧跟着solver那面的构造函数也加了一个参数进去。然后重新编译了src，忘记了编译solver，出错。
方案：重新编译solver后解决。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"

原因：在class里面增加了一个私有成员，编译库之后，求解器可以对类进行初始化，随即报错。报错只有3行，没有任何其他信息
方案：编译class之后，尤其这种增加私有成员的时候，记得重新编译solver一下就可以解决

@东岳 ok:huahua:

嗯嗯，好的，谢谢李老师指点:chitang:

边界条件、网格、格式

不可压缩LES模型算法非常成熟了，一般都是设置不正确带来的发散

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn ( fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1) - fvm::Sp(contErr1, he1) + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1) - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + p*fvc::ddt(alpha1) : -alpha1*dpdt ) - fvm::laplacian ( fvc::interpolate(alpha1) *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())), he1 ) ); E1Eqn.relax(); E1Eqn -= ( Kh*(thermo2.T() - thermo1.T()) + Kh*he1/Cpv1 - fvm::Sp(Kh/Cpv1, he1) + alpha1*rho1*(U1&g) + fvOptions(alpha1, rho1, he1) );

a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？

请问：固体颗粒的turbulent viscosity 怎么编程让其等于流体的？

请教一下大家，我现在取一小的truncated domain作为LES计算域，初始RANS结果是在20倍的弦长上计算的，截断区域是法向距离为0.5倍弦长，尾迹区域为1倍的弦长，如图所示。
truncated domain.jpg
入口设置为RANS结果插值速度和压力，通常OF入口压力都设为ZG，我不清楚这种插值速度设置是否可行，大家如果有知道的是否可以提供一些建议参考？
出口速度设置为advective，请问此时压力该如何设置？我设置为ZG, 计算结果直接从出口开始发散并向内部区域传递。

我在使用solidDisplacementFoam求解器进行热应力计算时，需要对位移向量场重新设置边界条件：其中的x,y方向为zeroGradient边界，z方向为fixedValue边界。请问OpenFOAM里有类似的边界条件或者知道该怎么设置的吗？？

李老师，我之前在of6和of4.1中的 controlDict里面设置为 边计算边输出空穴体积变化，现在2.3中我目前没有找到 边计算边监控体积变化的方法。。。

老师，您可否指点一个实现上述过程的思路，我去试试:xinxin:

貌似需要两边都懂才能讲出区别... 就像你明白你那面，我明白openfoam这面，但好像联系不起来，等待大佬分享..

建立了一个长宽高（x-y-z）为3015的立方体，目的是建立水槽。
1.jpg
使用setFieldDict字典设置水域范围，目标水域高度3米

FoamFile { version 2.0; format ascii; class dictionary; location "system"; object setFieldsDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // defaultFieldValues ( volScalarFieldValue alpha.water 0 ); regions ( boxToCell { box (-1.0 2.0 -1.0) (30.0 -2.0 3.0); fieldValues ( volScalarFieldValue alpha.water 1 ); } );

只计算到第5个时间步电脑就不运行了，有的时候可以一直运行，但是结果是全流域全部都是空气没有set的water field，而且我的建的流域模型X方向上的长度变成了15m，被砍掉了一半。
21.jpg
请老师们指教！感谢！

先对路过的大佬表示下感谢！
我看到不同算例中使用约束有constraints和restraints两种不同关键字，它们有区别吗？其中fixedPoint和fixedLine能理解，但fixedAxis代表什么呢？另外，我模拟的是浮体在水面上自由漂浮，不受约束，设置时直接将这个关键字删除了，但是运行过程发现浮体不动了，这是怎么回事呢？

5b9ac92e-ef00-485c-92bc-983ef6807c88-image.png
我想在求解浓度标量传输方程的时候将浓度的边界条件设置成图片的形式，其中速度是根据我0文件下的速度进行提取的，这个边界条件openfoam里可以直接设定吗，还是要我自己改变边界条件？

用tecplot得到圆管中切片x=15D的速度流场
small.JPG
希望能够得到如下图的速度曲线
案例.png
目前导出所有数据，摸索用matlab取一些半径上的平均速度得到的折线太丑了:xinlei: ；(https://www.cnblogs.com/code-saturne/p/10491817.html)的方法没怎么搞懂。请问用tecplot或paraview该怎样进行处理，恳求各位老师指点一下:xiexie:

Hello foamers:
如题
比如我想画水体的流场矢量图
那么应该怎样只把液相（体积分数>0.5）的速度场数据提取出来呢？
谢谢！

@范准
兄弟
这个问题解决了吗:chouchou:

@东岳 谢谢东岳老师，就是下面A beta Ta三个参数求化学反应速率常数，化学反应速率中反应物质浓度次幂直接在方程里面表现出来了哇，没有是默认为0次吗？