采用interFoam修改的求解器计算算例，遇到网格太小带来的壁面处压力问题（图中为第一个输出时间步（约0.0001 s）的压力分布，出现了不正常最大和最小压力）。
算例的壁面（BOTTOM）为锯齿状壁面，图中所示部分为倾斜向上的一小部分，壁面附近的网格精度约为1微米，压力场初始为0，压力和速度壁面处边界条件如下：
p_rgh
BOTTOM
{
type fixedFluxPressure;
value uniform 0;
}
U
BOTTOM
{
type fixedValue;
value uniform (0 0 0);
}
大家有没有遇到过这样的问题，有什么解决思路吗？
@linhan-ge 测试过，另外可以参看一下Johan Roenby的说法："For you surface tension people, I believe you will also need an improved estimate of the curvature, right? This information is somehow contained in the difference in surface normal direction between adjacent surface cells. Hmm... still thinking about how to cook this up in a consistent and robust way for general meshes..." 来自https://www.cfd-online.com/Forums/openfoam-news-announcements-other/180439-isoadvector-release.html
@队长别开枪 好的，我研究一下，十分感谢！
哈哈，让我想起了 “Windows sucks！！！”
@linhan-ge 我顶多算刚入门，聊八卦有点多，没怎么聊正事
VOF-DEM不了解，不好说，有什么影响，难度如何，建议先试一试吧，试下就知道了。下次可以把参考论文post出来
@队长别开枪 非常感谢，很有帮助，抱歉回复得有点晚了
我起初参考论文[1]，通过PLIC重构，进而得到一个每个面上AOF(area of fraction)，来作为alphaEqn.H中的phiAlpha相可以显著提高界面精度。此外论文中实现了非结构化网格的PLIC，依据的是Maric T[2]的迭代方法，这部分可能和你的工作相关。
isoAdvector的实现确实很有参考意义，学习中，最近作者好像在测试AMR了，可能不久会发布。
之前Maric T说过等完成了会公开Code，见https://www.cfd-online.com/Forums/openfoam-programming-development/89713-plic-interfoam.html，后来似乎没有。猜测可能出书，办培训班的原因，看了一眼去年他的博士论文，里面涉及了PLIC和其他界面追踪方法，个人底子差，确实难懂...
关于我的课题，除了精确界面，还需要关注的是界面张力项，CSF模型有更大影响，所以在试着实现一个新的张力模型，在你给的论文中作者反复提到了“ without smearing, dispersing or wrinkling.”，大家都知道OpenFOAM原来的方法在smearing方面表现不好，其余wrinking方面是不是也表现不好，在2-D vortex测试中会出现不明锯齿状界面，那dispersing方面呢？
[1]Dianat M, Skarysz M, Garmory A. A Coupled Level Set and Volume of Fluid method for automotive exterior water management applications[J]. International Journal of Multiphase Flow, 2017, 91: 19-38.
[2]Maric T, Marschall H, Bothe D. voFoam-a geometrical volume of fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM[J]. arXiv preprint arXiv:1305.3417, 2013.
@linhan-ge 作者一直在这方面努力，好像最近有新的进展了，上周Yoube更新了新的AMR视频，之前说六月可能会发布 https://www.youtube.com/channel/UCt6Idpv4C8TTgz1iUX0prAA
@linhan-ge iso-advector很棒，不过考虑到自己添加自适应功能的难度，观望中...
@sjlouie91 linux家族庞大，HPC部署的应该是CentOS RH某版本吧？查一下具体版本
@队长别开枪 你好，因为研究内容界面张力主导，需要两相界面更精确一点，想要在OF中实现结构化网格的PLIC界面重构，找了一些资料，包括Fortran版的实现，几何计算实在看得不是很明白
请教一下有讲得好一点的资料推荐吗？计算几何的知识可以怎么补一下？感觉很多论文里网格层面操作很多啊
@浪迹天大 师兄可以加个微信吗？wallong95:laughing:
@浪迹天大 恩，师兄天大人？
@李东岳 李博士可以考虑写一个规范代码风格的脚本，像CSS那样的:cheeky:
@wwzhao 多谢指导:happy:
如题，请问大家有没有了解在多相流求解中，使用接触角的定义和操作方法？
因OpenFOAM现有格式不满足需求，需要实现新的插值格式，有没有熟悉的人说明一下思路，此外我不清楚MeshObjects文件夹下的 upwindCFCCellToFaceStencilObject和upwindCPCCellToFaceStencilObject的作用和区别是什么？谢谢
例如下面的两相界面处的jump condition，在表面张力系数σ不为常值函数时，沿界面方向的表面张力作用不能被忽略
[μ∇u · n] = ∇s σ
其中∇s =(I - n⊗n)∇，表示表面梯度算子
在interFoam中表面张力系数为常值函数，用CSF模型计算法向jump condition，想知道在多数CFD计算中，这种jump condition都会怎么处理?