@李东岳 感谢李老师！原来没有特意写出来的都是同一个:zoule:

在CFDEM中，颗粒间相互吸引的作用力是使用skjr模型来实现的吗，还是说通过力来实现。sjkr模型中内聚能密度又是什么呢？

@lwjetmann 多谢，我自己增加了一个solver，确实是可行的。

@李东岳

帮李老师补充一下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

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

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

自己写的求解器别人很难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能够直接计算的，还请各位大佬解惑