CFD中文网

采用interFoam修改的求解器计算算例，遇到网格太小带来的壁面处压力问题（图中为第一个输出时间步（约0.0001 s）的压力分布，出现了不正常最大和最小压力）。
算例的壁面（BOTTOM）为锯齿状壁面，图中所示部分为倾斜向上的一小部分，壁面附近的网格精度约为1微米，压力场初始为0，压力和速度壁面处边界条件如下：
0001.png
p_rgh

BOTTOM { type fixedFluxPressure; value uniform 0; }

U

BOTTOM { type fixedValue; value uniform (0 0 0); }

大家有没有遇到过这样的问题，有什么解决思路吗？

@wwzhao 在 液体与固体表面的接触角怎么定义和操作？ 中说：

接触角是指气-液界面的切线穿过液体与固-液交界线之间的夹角θ。θ小于90°代表固体表面亲水，大于90°代表固体表面疏水。

可以通过对wall类型的边界设置alpha边界条件来设置接触角。

参考constantAlphaContactAngle、dynamicAlphaContactAngle以及timeVaryingAlphaContactAngle。

想请教一下，如果像定义个类似于timeVaryingAlphaContactAngle的边界条件，只是把时间换成坐标x，y或者z，该如何创建新的边界条件呢？

因OpenFOAM现有格式不满足需求，需要实现新的插值格式，有没有熟悉的人说明一下思路，此外我不清楚MeshObjects文件夹下的 upwindCFCCellToFaceStencilObject和upwindCPCCellToFaceStencilObject的作用和区别是什么？谢谢

例如下面的两相界面处的jump condition，在表面张力系数σ不为常值函数时，沿界面方向的表面张力作用不能被忽略
[μ∇u · n] = ∇s σ
其中∇s =(I - n⊗n)∇，表示表面梯度算子
在interFoam中表面张力系数为常值函数，用CSF模型计算法向jump condition，想知道在多数CFD计算中，这种jump condition都会怎么处理?