@李东岳

帮李老师补充一下OpenFOAM里面average的实现，对于vol场要先插值成surface场

template<class Type> tmp<GeometricField<Type, fvPatchField, volMesh>> average ( const GeometricField<Type, fvPatchField, volMesh>& vtf ) { return fvc::average(linearInterpolate(vtf)); }

所以|3|0|2|到average函数里面以后要先插值成

| 3 | 0 | 2 | 1.5 1

再基于

$v_{celli}=\frac{\sum (|\bfS_f| v_f)}{\sum |\bfS_f|}$
得到中间网格的值为1.25

@尚善若水 谢谢大佬，笔芯

@tidedrinker 非常感谢您分享您的经验！
我这个算例拿starccm跑过（二维网格，没有z轴厚度）结果并不理想，我也十分苦恼。

关于starccm网格导入openfoam，我个人的经验是先有starccm导出.ccm格式的三维网格，之后用ccm26ToFoam *.ccm操作就可以将三维网格导入openfoam了，值得注意的是，openfoam没有纯二维，都是通过设置empty边界条件实现二维的模拟，需要注意的是z方向网格数为1，这里会用到extrudeMesh的操作，具体内容见田东老师的视频
希望能够帮到你！

@李东岳 谢谢李老师，我按照您的方法已经解决了

自己写的求解器别人很难debug
有可能是边界条件没处理好

@李东岳 多谢老师解惑！:xinxin:

@guanguan0922 非常感谢！

各位大佬，我在看东子老师的势流算法文章，有些地方思考了很久还是不得解。问题①：出口边界条件为速度梯度=0，方程退化成方程（5），求解困难；假设出口足够远，又可以通过方程（5）来求解了，这俩表述看的我一脸懵62cfaf3b-86df-414e-974f-0351bf0a52b5-image.png 问题②：“方程（5）和方程（6）可以实现相同的效果”，是指啥效果呀？我手动算了一下(见图)，怎么保证两个c一样呢？∅_outlet的值如何确定呢？50cc50b2-4934-4b37-a443-2d507b865049-image.png 问题③：方程（10）我自己算了一遍，是个二次函数，和之前（5）、（6）都是一次函数相差太大了吧，不清楚（10）是如何能等于（6）8405b555-e186-4a2f-8ef1-9e5cc41fb6ec-image.png

@李东岳 好的老师，我这边自己再研究下，谢谢老师。

@李东岳 好的，谢谢老师。

@coolhhh 我去试下 谢谢老师指点

@尚善若水 大佬我试试，感谢感谢

@李东岳 谢谢李老师建议，我去试一下这种边界条件

@李东岳 谢谢老师

@haozhibin @youhaoyu @wangfei9088 感谢三位大佬

目前能看出来的，就是第二行dnew引入了volScalarField

增加了一种读取表格计算物性的方法，编译无报错但实际应用时显示下列报错信息eab5b1e7-fcb8-4b12-be0b-abecc65ae3fd-image.png

@AppleKiller 好的，谢谢

@tidedrinker 谢谢大佬非常感谢您

如题，后向台阶阻力表达式什么呢？我在网上查找到其阻力表达式分为两个部分一是摩擦阻力另一个是升压阻力，想减少边界层厚度，直接让阻力值最小是不是就能够达到目的？或者有其他更加直观的观测数据，最好是OpenFoam能够直接计算的，还请各位大佬解惑