@李东岳

@李东岳 谢谢李老师！0.5暂时还没有找到足够有说服力的依据，只是看到有些文献里直接用的μ=0.5pfsin(phi)；
但我觉得wachem2000里面公式不太对，因为他引用的JJ的文章是这样定义的：

请教各位前辈一个问题：

phi是一个surfaceScalarField，和volScalarField不矛盾吗？

此外，现在需要对速度散度判断，以确定n的取值：

n是一个volScalarField或者scalar也可以，要计算：

写循环的话，应该用div（U），还是div（phi）？（ps和pc都是volScalarField）

@Yongbo 谢谢

@李东岳 李老师，请问发散后继续计算又不发散了，比如如下报错：

或者报错与solver有关（gauss-seidel），这证明植入的模型还是不容易收敛是吧？

@veen 还在整理研究结果 ，但是也有了一点新认识：
1）沉降效果与曳力模型和颗粒粒径有关；
2）砂堤回缩与边界条件和固相应力模型（猜测，用JJ比JJShaffer要好，且正在植入其他模型）有关；
3）达不到alphaMax我猜想可能也和固相应力模型有关。
争取早日把成果研究整理分享出来

请教各位前辈老师：


JJ模型中的0.5是哪来的呢？

请教各位前辈：

针对成员函数：

定义：

返回值：

报错：

但直接返回：

计算没问题，不知问题出在哪里

在dargmodel.C这个衍生类中没有mesh这个成员对象，不知道该如何引入

如果改变写法：

Foam::dragModels::wallEffect::wallEffect
(
    const dictionary& dict,
    const phasePair& pair,
    const bool registerObject
   
)
:
    dragModel(dict, pair, registerObject),
    residualRe_("residualRe", dimless, dict),
    y_(nullptr)
       
{}

Foam::dragModels::wallEffect::wallEffect
(
    const dictionary& dict,
    const phasePair& pair,
    const bool registerObject,
    const fvMesh& mesh
   
)
:
    dragModel(dict, pair, registerObject),
    residualRe_("residualRe", dimless, dict),
    y_(&wallDist::New(mesh).y())
       
{}

计算的时候会报错：

@李东岳 也烦请东岳老师帮忙看看这个问题如何实现~ 三月再您的课上还问过您这个问题，但是下来自己搞没搞通