CFD中文网

我是在多孔介质里面模拟气液两相流，为什么我的速度流线是这样的呢，蓝色气泡位置的流速很大
image.png

大家好，我目前在进行颗粒填充床干燥过程的CFD数值模拟，模拟对象为含湿颗粒（粒径约0.2mm~3mm）填充床热风干燥过程，涉及传热与传质的多物理场耦合问题。干燥床层可以视为多孔介质结构，现希望实现以下功能：

目标：
模拟颗粒填充床中水分的扩散/蒸发过程；

实现多孔介质区域内的非稳态水分传质（可包括对流+扩散）模型；

考虑温度场与水分场的耦合（温度影响扩散系数或蒸发速率）；

求解干燥过程中床层含湿量的空间分布与随时间的变化。

目前的思路/尝试：
将床层视为各向同性多孔介质，定义有效扩散系数、孔隙率、饱和度等参数；

使用UDF定义源项（如水分蒸发速率），或者通过Species Transport模块进行气-液相质量转移；

温度场通过Energy方程求解，热源项与水分蒸发挂钩；

遇到的问题：如何合理设置传质模型，特别是湿含量作为标量的定义方式不明确，Species模块中对水蒸气的处理较困难。

求助点：
是否可以通过UDF或者UDS自定义传质方程？有没有类似案例可以参考？

如果使用Species Transport模块，如何处理固相内的水分扩散与蒸发（不是自由流）？

有无建议的模型框架（例如某些文献中常用的多孔传质模型）可用于实现“含湿颗粒床的干燥”过程？

欢迎各位大佬分享经验或推荐文献、资料，非常感谢！

最近看《Holzmann T. Mathematics, numerics, derivations and OpenFOAM》，里边也提到大部分情况忽略k，对结果影响不大：
05a3730d-81ea-47f7-93d1-d3e7f8814605-image.png

@疏影横斜水清浅 你好，可以咨询一下这个的相关问题吗

另外，CFDEM如何实现壁面变形呢，像下面这样逐渐撑开流动通道
image.png

@Tens 请问下 如何在算例中调用 src\thermophysicalModels\thermophysicalProperties\liquidProperties的物性

http://dyfluid.com/datadrivenles.html 李老师方便提供完整编译的文件么？最近在编译数据驱动RANS模型，但是编译的时候，总是报错，libtorch与openfoam的typename冲突，不知道咋解决。

就是速度入口，柱子面积和入口面积的比×表观速度

您好，我最近也遇到了相同的问题，请问您是怎么解决的呢？我已经在snappyHexMeshDict中指定了面贴合参数

snapControls { nSmoothPatch 5; nSmoothInternal 5; tolerance 2; nSolveIter 300; nRelaxIter 5; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; }

@DY大世界 您好，我最近也遇到了相同的问题，请问您是怎么解决的呢？我已经在snappyHexMeshDict中指定了面贴合参数

snapControls { nSmoothPatch 5; nSmoothInternal 5; tolerance 2; nSolveIter 300; nRelaxIter 5; nFeatureSnapIter 10; implicitFeatureSnap false; explicitFeatureSnap true; multiRegionFeatureSnap false; }

各位老师好，有一个问题想请教一下，就是我在openfoam11中对fluid求解器进行了修改，重命名为efliud，并且编译成功，然后在controldict文件中输入：application foamrun; solvers efluid;计算的时候提示efluid solver为无效求解器，这是调用的问题么

各位老师号，有一个问题想请教一下，就是我在openfoam11中对fluid求解器进行了修改，重命名为efliud，并且编译成功，然后在controldict文件中输入：application foamrun; solvers efluid;计算的时候提示efluid solver为无效求解器，这是调用的问题么

这个算实际河道是不是很难用啊

@李东岳 在 driftFluxFoam适合用来做泥沙冲淤吗？ 中说：

新版OpenFOAM-10里面driftFluxFoam更新挺多的，加了几个新的粘度模型，你可以看下

这个模型是不是不能捕捉水面啊，那是不是意味着没办法算自由液面流呢

@wangfei9088 谢谢！只能慢慢试了

大家好。

我最近在学习cfdem，模拟了小球入水。流体网格尺寸0.003m，小球直径0.0254m，钢制小球。界面捕捉采用VOF方法，并采用浸没边界法IBM刻画小球形状。液体（水）表面张力0.07。

入水后小球位移模拟值与试验值对比良好。使用paraview-counter设置alpha.water=0.5提取液面，发现气液界面附近的空腔结构与试验照片对比较好，但是小球周围并没有出现空腔结构。

请问是算例设置相关问题还是paraview设置问题？

文献（https://doi.org/10.1016/j.oceaneng.2021.110198）
a72bdbf6-06ad-4f91-8cb0-05d6ada9f5cc-image.png

我的模拟结果小球周围没有出现空腔结构
2e351048-86c2-47d4-9f10-b76146ea1e8a-image.png

该书11.8节（中文版305页，英文版411页）中的例2，推导源项的表达式时，高阶格式减去一阶迎风格式，一阶格式的值取上上游的值φu，为什么这里不取φc呢？先谢谢了。
7a735b6e-01c9-4088-971a-359144da8189-image.png

@wangpc 有什么问题可以直接在这里问，集思广益嘛

@李东岳 driftFluxFoam可以捕捉自由液面吗

各位大佬好，我这边采用OpenFOAM5.x版本进行计算，在场2中存在一个计算变量A需要传递到场1的重叠区域中，目前两个计算域的网格如图片中所示存在重叠，且设置的边界保持一致。请问该如何保证数据传递的可靠性，我目前使用的指令如下：
mapFields $sourcedir$ -case $targetdir$ -noFunctionObjects -sourceTime '30'
指定map方法的三种（mapNearest, interpolate and cellPointInterpolate）我也试过，效果也并不是很好。各位大佬给点建议，有没有相类似的数据处理经验。感谢感谢！

14295afe-5391-425a-88dd-94e8f07b96ab-image.png