@dzw05 @溯光 @李东岳 您好，我想请问一下这个interface里面的C表征的是流体对粒子一个什么样的力呀？我最近在验证波浪作用下一个粒子的输移运动，发现改变不同的C会使得粒子在水平和垂向位移上与stokes质量输移有较大的偏差

image.png

各位老师好，想请教一个问题，我现在在controlDict文件里写了一个function，里面需要读取一个字典中的值。字典中的值为一个标量，读取之后想改变这个值并重新写入字典中，这样下次读取之后就是新的值。
其中字典文件定义如下
1.png
function里面读取的代码如下
2.png
我尝试用write函数来改写，但是报错是这样
3.png
请问各位老师这个该怎么解决呢？谢谢！

@crispy 我去提bug后然后回复你

@Vortex 请问dynamicLagrangian湍流模型，flm和fmm的边界设置是什么样的

@Wayne 你好，请问一下这个问题怎么解决？

已解决 如果流场先动，颗粒后动，可观察颗粒运动时间节点t，将pre shift调整为-t。

@李东岳 谢谢师兄，因为模型有辐射模型，流固耦合，还有多物种传输，所以一直觉得是自己设置有问题，或者部分模型需要重新编译，忘了最基本的，刚才把所有收敛系数都调小了.可以跑了，应该就是前期不稳定。但是，也只是多跑了3步，我回头再看看

@Alvin 您好，您当时怎么解决这个问题的呢

openfoam7你就得按照主贴的方式试一下了

@李东岳 哦哦，懂了，那我都试试吧。谢谢李老师及时回复。:xiexie: 感恩的心:xinxin:

@李东岳 李老师，这是我自己编译的一个求解器，我在constant文件夹下将g从（0 0 -9.8）改为（0 0 0），后消失了，但是不符合正常要求呀

fvSolution里面的pFinal，relTol应该调到0吧。
矩阵最终求解的这些*Final一般都强制要求满足容差。
应该是这个原因导致你某一步没收敛，或者误差逐渐积累，导致发散。

@zhoushan 请问你解决了吗，具体在哪里切换驱动呢
18af7957-d664-4bae-9518-85ff95fa17ea-image.png

我以为是realizable本身的问题。看来需要详细debug一下OpenFOAM里面的这个湍流模型。

找到数值发散的原因了，是热流密度设置过大

谢谢啦，主要想找一个带相变的可压缩气液两相流的solver，我再找找试试~

好贴，谢谢李老师分享！！！

@小鲨鱼 在 OpenFOAM涡环模拟 中说：

数学上是有一个积分公式实现涡量场到速度场的转化的

写出来看看

如果可以的话，可以有偿指挥。

老哥，可以分享一下算例吗？我正在学习

ESI那面的OpenFOAM我没有仔细研究过，知识盲区了

这个就是发散了，目前看不出来什么原因

各位大神好，我在做一个关于可压缩流体的动网格算例时遇到了问题：我的模型是一个长方体，当我将x，y方向上的边界都设置为slip时，z方向的边界沿着z轴移动，从而实现压缩这个模型，这样是可以计算的。
但是当我想要实现沿x和z两个方向都压缩时则出现了问题，此时x和z的边界都沿着各自的轴移动，这样在x和z边界的接触的位置网格产生了畸变。想请教一下应该如何设置？

@李东岳 感谢东岳老师的回复。已解决。需要在system里有一个单独的forceCoeffs文件，然后在controlDict里的添加functions{#include "forceCoeffs"}，这样会在输出文件里输出每一步的Cd等相关信息。感谢老师！

@Dingcy 在 如何输出出口escape的颗粒数量随时间变化的曲线？ 中说：

如何输出出口escape的颗粒数量随时间变化的曲线？之前算单向流的时候我输出出口浓度随时间变化曲线是通过paraview那个plot selection over time, 但是现在这个颗粒好像不能那么整了....:136:
1c05a757-a68a-499f-bfc9-1582a3c49a26-image.png
前辈您好，请问这个操作怎么实现，得到下图:134:

我在无痛苦ns方程笔记里面增加了对这一行代码的解释：

fvScalarMatrix alpha1Eqn ( fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY??? - fvm::laplacian(alpha1alpha2f()*pPrimeByA_(), alpha1) ); alpha1Eqn.relax(); alpha1Eqn.solve();

尤其是其中的

fvm::ddt(alpha1) - fvc::ddt(alpha1) // WHY???

@Tens 恭喜恭喜！！！！隔壁师兄又要嫉妒死了！:mihu: :mihu: :mihu:

@李东岳 多谢李老师:146:

@李东岳 感谢东岳老师

这有一个已经做好的：

https://www.youtube.com/watch?v=fLddwYmp2Kg

A Closer Look At The Evaporation Process:
Illustration of evaporation process of stagnant water film of 4mm thickness at 25°C placed in an environment with dry air at 1 atm.
The interfacial mass transfer model is implemented in icoReactingMultiphaseInterFoam solver based on the work of Van P Carey [1] where evaporation and condensation fl
uxes are given based on the kinetic theory for a fl
at plate interface considering Hertz-Knudsen relation.

[1] Van P Carey. Liquid vapor phase change phenomena: an introduction to the thermophysics of vaporization and condensation processes in heat transfer equipment. CRC Press, 2018.

@pengdi LES这么添加是不对的。你那个公式是RANS的k，不是k_sgs。LES那面如果你用smagorinsky，不需要给k

文章出来的时候，烦请通知一下哈。

Distribute the map ... [cu01:53769] *** An error occurred in MPI_Bsend [cu01:53769] *** reported by process [140195840131073,17] [cu01:53769] *** on communicator MPI COMMUNICATOR 3 SPLIT FROM 0 [cu01:53769] *** MPI_ERR_BUFFER: invalid buffer pointer [cu01:53769] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort, [cu01:53769] *** and potentially your MPI job)

您好，请问您后来解决这个问题了吗？

@李东岳 谢谢李老师，我再回去看看是不是我公式的问题

我看了你的一些帖子，大部分都是基本的使用问题。这一类问题在学习用户手册后都可以被解决。 http://www.dyfluid.com/OpenFOAM用户指南-9.pdf

先处理出来涡量 ： http://dyfluid.com/code.html#id29 使用streamLine函数可以处理某个平面的速度，压力，涡量等值： http://dyfluid.com/code.html#id40

@hr Moodyn我没有试过，Moody是成功安装使用了的，后面阻力太大就放弃了。

@Stan同学 看你的图片，网格变形并不大，按理是不会出现负体积这个问题的，而且扭曲率也没有很离谱啊，为什么会错呢？openfoam其他的动网格方案都卡的比较死，我之前是放弃用openfoam实现较大网格变形的算例了。如果你真的要弄大变形，可能需要对Openfoam里的变形进行重新编程该写算法去做了。或者用FLuent的重构网格试试，也可以用Udf写新算法进行计算。

建立两个网格就行了。P一套网格。Q一套网格。

思路具体是需要对OpenFOAM中某一求解器中的头文件进行修改和变化，并修改其中其中一些项。是自己建立一个求解器然后将原来的自带求解器中代码粘贴到建立的求解器中修改，然后编译的过程吗？有什么注意事项吗？

@李东岳 谢谢东哥

@李东岳 李老师，在使用openFOAM-V2006版本的buoyantBoussinesqFoam求解器，出现如下报错：
--> FOAM FATAL ERROR:

request for basicThermo thermophysicalProperties from objectRegistry region0 failed available objects of type basicThermo are

0()

From const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::basicThermo] in file /home/lyc/OpenFOAM/OpenFOAM-v2006/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 463.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::exitOrAbort(int, bool) at ??:?
#2 Foam::basicThermo const& Foam::objectRegistry::lookupObjectFoam::basicThermo(Foam::word const&, bool) const at ??:?
#3 Foam::fv::fixedTemperatureConstraint::fixedTemperatureConstraint(Foam::word const&, Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#4 Foam::fv::option::adddictionaryConstructorToTableFoam::fv::fixedTemperatureConstraint::New(Foam::word const&, Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#5 Foam::fv::option::New(Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#6 Foam::fv::optionList::reset(Foam::dictionary const&) at ??:?
#7 Foam::fv::optionList::optionList(Foam::fvMesh const&, Foam::dictionary const&) at ??:?
#8 Foam::fv::options::options(Foam::fvMesh const&) at ??:?
#9 Foam::fv::options::New(Foam::fvMesh const&) at ??:?
#10 ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/buoyantBoussinesqPimpleFoam
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/buoyantBoussinesqPimpleFoam
请问是什么原因呢，我的constant 和0文件与tutorials一致，还是报错。

同问一下，我怎样能够提取某一截面中某一部分的速度，感谢

@李东岳 谢谢李老师

@heike256 hyper-v我都不知道是啥，大佬整的太高端了。我就知道虚拟机，wsl，真系统。那你试试真系统。

好了，没事了，我眼瞎

请问这个边界条件生成的流场脉动是线性变化的吗？

最后自己修改了mapFields把映射内部场的代码删掉了

@李东岳 在 OpenFOAM里面的dnsFoam为什么这么菜？ 中说：

@ZifeiYin

Hi,

请参考下述内容，把第一行第二句话自行删去：

鉴于CFD圈子的小众性，“CFD中国”执行最严厉的友善管理制度：所有用户请以专业的态度对待所有问题。不得在所有版块发表轻蔑性言论；
例如，下面的语句被认为是具有轻蔑性的：“楼主回去看看书吧”，“这种问题你也问？”，“不会用百度？”
“CFD中国”建议用户对每个问题严谨友善处理。例如上述内容可以这样表述：“建议楼主参考《数值传热学》第48页，有对xxx问题的详细表述”，“个人认为这个问题是一个非常普适性的问题，在大部分的教材中都有涉及，楼主可以参考xxx”，“这个问题和CFD的关系不大，建议楼主去百度。google去搜索相关内容”

哈哈哈哈，东岳老师干得漂亮

@李东岳 好的，明白！感谢老师！