@东岳 好的谢谢，我先去学习下。

有没有使用过 CFDEM®project 的大佬，我最近在调研一次雾化算法方面的文章，想改进一下OpenFoam的拉格朗日库，偶然看到 CFDEM®project，不知道能否解决这个问题，关于这部分的库，是开放的还是需要另外购买

我个人推荐一般需要我认可。但我学习OpenFOAM程序都是自学，只看过《OpenFOAM用户指南》，《OpenFOAM编程指南》，并且我断断续续看了一年《C++ primer plus》，除了这三个因此没有文章推荐给你因为我也没看过不好评价

学习CFD这面倒是有几本书推荐给你，比如《数值传热学》，《Computational Methods for Fluids Dynamics》，应该会有人推荐别的书，比如OpenFOAM Matlab那个、不过我个人觉得没上面这俩本好

另外可以看看一些OpenFOAM的博士论文，比如Jasak的博士论文是个很好的资料

暂时只想到这么多

@星星星星晴 感谢

@wwzhao

谢谢老师！问题已经解决。果然是 输出文件的路径不对！应该是 放到$FOAM_USER_APPBIN下面:chitang:

@东岳 谢谢，速度方向没有问题，网格就是用blockMesh画的结构化网格，刚开始算，网格没有画太多，比较稀疏

@散漫守望2016 多谢您的指导，用自己的foam算过软件自带的例子，结果和标准的例子结果比较吻合，应该不是程序的问题。根据您的分析，我觉得很可能还是我模型结构的问题。十分感谢。

显示的柯朗数也很小，10-4级别，应该有空余增大时间步长，但是增大一丁点就会发散

我是从事水利行业

需要动网格实现的功能是，网格本身位置不动，当水深小到某个临界值时，这个网格由湿的转变成干的，也就是变成固体边界，水深超过某个临界值时，干网格又恢复成湿网格，这个不知道怎么实现，有没有相关算例可供参考

不需要

我也没看懂，楼主写得太复杂了，别人没功夫仔细看你的公式。我建议你自己用matlab算好每个mesh上的数值，然后写成csv格式用swak4Foam导入。

运动物体计算，浸入边界。

@dzw05 好的谢谢了！！！

谢谢李老师~

编译C文件的时候，我总是使用注释来找出错误的地方，然后保存，再在终端wmake一下。但感觉这样的编程效率太低，OF里能否断点编译？能的话需要使用到哪些软件工具，大家有什么推荐吗？

我用OpenFOAM计算外流问题，当来流速度设置为0时可以正常计算，但当来流速度大于零时最多计算一步流道就堵塞了，我怀疑是边界条件有问题。我的边界条件设置如下，请教大家这样搭配可以吗？如果不行应该是怎样呢？
U

FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform (0 0 50); } outlet { type zeroGradient; } symm { type slip; } propellerW-f { type movingWallVelocity; value uniform (0 0 0); } propellerW-s { type movingWallVelocity; value uniform (0 0 0); } propellerW-h { type fixedValue; value uniform (0 0 0); } propellerW-t { type fixedValue; value uniform (0 0 0); } }

p

dimensions [1 -1 -2 0 0 0 0]; internalField uniform 101325; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type zeroGradient; } outlet { type fixedValue; value uniform 101325; } symm { type slip; } wall { type zeroGradient; } }

T

dimensions [0 0 0 1 0 0 0]; internalField uniform 293; boundaryField { inlet { type fixedValue; value uniform 293; } outlet { type inletOutlet; inletValue uniform 293; value uniform 293; } symm { type slip; } "inter.*" { type fixedValue; value $internalField; } wall { type zeroGradient; } #includeEtc "caseDicts/setConstraintTypes" }

nuTilda

dimensions [0 2 -1 0 0 0 0]; internalField uniform 1e-5; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform 1e-5; } outlet { type inletOutlet; inletValue uniform 1e-5; value uniform 1e-5; } symm { type slip; } wall { type fixedValue; value uniform 0; } }

nut

dimensions [0 2 -1 0 0 0 0]; internalField uniform 8.58e-6; boundaryField { //- Set patchGroups for constraint patches #includeEtc "caseDicts/setConstraintTypes" inlet { type fixedValue; value uniform 8.58e-6; } outlet { type inletOutlet; inletValue uniform 8.58e-06; value uniform 8.58e-06; } symm { type slip; } wall { type nutUSpaldingWallFunction; value uniform 0; } }

alphat

dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 0; } outlet { type calculated; value uniform 0; } symm { type slip; } "inter.*" { type fixedValue; value $internalField; } wall { type compressible::alphatWallFunction; Prt 0.85; value uniform 0; } #includeEtc "caseDicts/setConstraintTypes" }

其中inlet，outlet，分别是进口和出口，几何边界条件为patch；symm为四周的自由边界，为柱状，几何边界条件为patch，propeller是计算的螺旋桨，几何边界条件为Wall；计算所用湍流模型为SA模型，所用求解器为rhoSimpleFoam,设置了两块MRF区域，MRF设置如下：

FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object MRFProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // MRF1 { cellZone sector_part1; active yes; Fixed patches (propellerW-f); nonRotatingPatches (); origin (0 0 0); axis (0 0 1); omega 112.5737; // rad/s } MRF2 { cellZone sector_part2; active yes; Fixed patches (propellerW-s); nonRotatingPatches (); origin (0 0 0); axis (0 0 1); omega -112.5737; // rad/s }

@Calf-Z-DNS delta Y+的计算公式你已经写出来了，和delta X+是一样的。平板边界层的估算算法：https://www.cfd-online.com/Wiki/Y_plus_wall_distance_estimation

算例为土柱冻结算例，起始土柱内部温度为6.7℃，土柱顶部施以温度为-6℃的循环液流(the circulating fluid)，使土柱从上到下逐步冻结，下图为示意图，
捕获.JPG
温度边界条件：
变热通量上边界条件/a variable heat flux upper boundary condition：
qh=hc(TCoolant-TTop), hc=28 W m-2 K-1
其中：TCoolant是循环液流的温度值(-6 ℃)，而TTop是土柱顶部温度，我理解的是每一次时间步长计算后都会更新，所以顶部通量是一个动态更新的过程。
1、之前有看过关于fixed heat flux boundary的帖子，采用fixedGradient边界条件代替，但没找到关于 variably heat flux boundary condition
2、 采用swak4Foam？
希望能给我一些建议，共同探讨下。非常感谢！

大家好，我最近想做一个这样一个多孔介质：在多孔介质的不同部位，接触角也是不同的。即通过指定两点坐标，来设置两点所包含几何体的接触角，类似于setFields中设置水相分布的功能。这能否通过在openfoam中编译程序来实现？谢谢大家了。

@东岳 我在做单液滴破碎时没有发散，但是用同样的物性和离散格式等做液滴碰撞case就会出现发散