CFD中文网

@coolhhh 感谢

挖老坟 我今天下午看fluent 2023R1 theory guide 的相间传质 14.7.7节 你要是多相多组分 不能使用VOF 应该使用mixture 或者双欧拉模型 才能进行多相多组分之间的传质计算 我使用VOF计算了好久好久 看fluent 手册 才知道这

@J 9年过去了...

@ibelief 资料能发我一下吗，谢谢，邮箱2682577721@qq.com

参考：https://zhuanlan.zhihu.com/p/427686653
另外注意网格的尺寸，of里面是米的单位，一般需要使用transformPoints从mm转到m

@李东岳 李老师好，这里我请教一个问题->链接

为了简化算例，我把网格三维转化为二维，简化了结构形式，以便排错。

【算例云图】
6101495a-adbe-455d-ab41-45e2b1b7b415-image.png
bb59f1c1-9dde-4169-83d9-cf00c0daf433-image.png

【简化版算例下载】
functionalTest_NoPZRefine2_2D.zip
生成网格createMesh，
跑算例runCaseOnly，
全流程runCase

感谢~

我已经忘记我2020年写的这个'OpenFOAM拉格朗日粒子包内嵌的算法实现（每个网格一个粒子）'是咋实现的了。不过我现在记得有一个潜在的问题，就是这种方法是一种反向找粒子的方法，我记得我当时边界条件好像没太想明白。只能想起这么多了。

想请问怎么统计一段时间内欧拉多相多组分出口处某一物种的质量呢，使用intergral对出口处积分，specise的变量只有质量浓度，找不到质量流率，所以是不是需要自己使用表达式
啊？f094ca30-ada6-40d5-af75-b330c0b43f2c-3a2845d63bf7d22609216510d354461.png

@hy1112006 您好，如果方便的话可以有您的邮箱或者联系方式交流一下使用pisoFOAMopenfast吗？

想尝试与一篇文献中实验结果作对比，但根据文中所提供的条件算出的入口速度很小（0.18m/s），以此作为入口的边界条件；同时入口边界的alpha.water也固定为1，计算一段时间后如图所示，发现入口处几乎没有水相的存在了。这是为什么？

ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

入口处的边界条件如下：

U { type fixedValue; value uniform (0.18 0 0); } p_rgh { type zeroGradient; } alpha.water { type fixedValue; value uniform 1; } alpha.air { type zeroGradient; }

多谢老师

我觉得这个是一个自圆其说的东西。如果看sci，肯定有人对比LES和DNS的结果，证明他们的LES或者DNS的结果更好。

但我认为，本身欧拉拉格朗日里面存在了很多模化，在粒子受力的角度，存在更多的模化，尤其是湍流分散力那面，基本都是参数。本身DPMFoam和DNS就是冲突的东西。所以如果做这方面研究，自圆其说是最重要的。

AI-CFD有时候可以实现CFD完全没法实现的问题。这样就比较好。但是有很多AI-CFD跟经典CFD是完全重复的，在这种情况下，目前来看还有很大提升。比如PINN的训练效率问题。但如果用PINN做隐藏流体力学，这就好多了。改天给你发个公众号问问。哈哈。看看评论区怎么说。

https://cfd-china.com/topic/8018

能做的话可以上传到这里我给你看看

@余正东 这个我关注了，未来可期，目前只是初步移植了一些

@先先贝 我的方法都在之前写的帖子里面分享出来了。算几个时间步很快的。而且串行计算三个时间步就差不多了，你设置的自适应加密等级也基本到位了。这个时候再运行一次setFields，给加密后的较密的网格赋值，再去分块。
分块的话也得想想加密的部分会不会在你流场里面移动，你的不同核之间的负载要基本平衡，可能还需要你手动去划分区域。

OpenFOAM中使用multiDimAMR编译和运行正常，使用reconstructParMesh也正常，但是之后再使用reconstructParMesh会报错。
第一个时刻还是正常的，第二个时刻发生报错，还有一个算例是第5个时刻开始报错，报错内容一样。使用paraview中的decomposePar也能查看结果，为什么会出现这个问题呢~
Reconstructing fields for mesh region0
Time = 0.01
Reconstructing FV fields
Reconstructing volScalarFields
alpha.water
alpha.water_0
p
p_rgh
rAU
Reconstructing volVectorFields
U
Reconstructing surfaceScalarFields
alphaPhi0.water
phi
Reconstructing surfaceVectorFields
Uf
Reconstructing point fields
No point fields
No lagrangian fields
No FA fields

Time = 0.02

--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.
--> FOAM Warning :
From virtual Foam::polyMesh::readUpdateState Foam::polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 210
Number of patches has changed. This may have unexpected consequences. Proceed with care.

**--> FOAM FATAL ERROR: (openfoam-2012)
Processor 2 has a different polyMesh at time 0.02 compared to any previous processors.
Please check time 0.02 directories on all processors for consistent mesh files.

From Foam::polyMesh::readUpdateState Foam::processorMeshes::readUpdate() in file processorMeshes.C at line 192.**

FOAM exiting

刚好在用kinematicParcel ，顶大佬一波！！

@学流体的小明 是的 想通过设置入口湍流速度边界来发展出槽道内湍流。刚接触这一块，觉得是不是可以 用贴子这个初始扰动的方法，设置成入口边界，结果试了一下单纯加上去不太行。就想问问有没有类似的设置入口边界的方法。

大佬们，这个算例用rhoCentralFoam能跑通。因为后面想加入颗粒计算，现在在同样的边界条件下我用rhoPimpleFoam计算时，却一直卡在这里，入口的压力也不对，总压给定的是3e6。这是什么问题呢？
image.png

image.png

image.png