CFD中文网

@上级 您好，这个可以在高斯点中开启高斯阵列，如下图

@上级 第一个问题的方法是
1、用foamToVTK命令把他转成VTK格式
2、用paraview打开这个VTK文件夹中的拉格朗日文件，并用高斯点表示
3、设置半径就可以了
第二个问题：确实就是用平均压力来表示该流化速度下的压降，但是这个流化速度下的压降随时间变化也不该这么剧烈，下面是一篇文献中的压降随时间变化的曲线

所以兄台你的压力波动变化大吗

用DPMFoam算二维喷动床流化现象，求解出颗粒相的运动好像一切正常，如下图

但是反应器的压降过于剧烈是什么情况？压力求解结果不对吗？

压力虽然波动剧烈，但是平均值都没问题，不知道是什么问题，换了几个入口速度都有类似的情况。按照文献压力波动不会这么大，±2000Pa以内比较正常。

各位老师，我之后在DPMFoam中验证这一槽道流型，不注射任何颗粒，发现怎么都验证不出来，依旧是入口在下，出口在上，讲道理速度剖面是个抛物线，出来却是个边界梯度很大的曲线，如下图，中心区域毫无梯度

用的是层流模型，想到秃头，边界条件设的以往一样，出口边界条件设的是zeroGradient，左右两个壁面是fixedValue，value设定为（0 0 0）
（其实我在coalChemistryFoam中不注射颗粒也跑了槽道，但是结果和DPMFoam一样，所以我怀疑是哪里设置出现了问题，导致跑不出抛物线来）
，收敛曲线如下

应该是收敛的把，虽然压力曲线在波动，但都是在收敛值以下

谢谢@东岳 老师回答，之后我解决了这个问题。

最近用icoFoam求解一个十分简单的二维槽道流动问题，入口速度设定为0.1m/s，不考虑重力，气流从下往上流。按照解析解来说，槽道出口截面速度分布应为抛物线形，抛物线最高点速度为0.15m/s。但是我用icoFoam的计算结果显示，出口截面速度在边界出现了匪夷所思的不连续现象。如图

但是速度分布在槽道内都是连续的抛物线，仅仅在出口界面的网格上不连续了，求问这是什么问题？有什么解释吗？我的出口边界条件设的是zeroGradient，左右两个壁面是fixedValue，value设定为（0 0 0）。

在OpenFOAM8.0和2.3中跑出来都是不连续的出口速度分布。

十分感谢！@东岳

奥第二个问题刚刚查代码搞清楚了，对于不可压缩流动的应力计算，密度都设为1了，所以速度方程代码的量纲没有问题。

谢谢@东岳 老师，但我依旧有两个后续问题想请教一下：
1、既然

- fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)

表示连续性方程和速度的乘积。
从物理意义上来说，连续性方程恒为0，那么将连续性方程放到速度方程里求解的意义是什么？
是为了更好的数值稳定性吗？还是说为了修正速度值？
2、我在湍流模型里确实找到了divDevTau的代码，如下

template<class BasicMomentumTransportModel>
Foam::tmp<Foam::fvVectorMatrix>
Foam::linearViscousStress<BasicMomentumTransportModel>::divDevTau
(
    const volScalarField& rho,
    volVectorField& U
) const
{
    return
    (
      - fvc::div((this->alpha_*rho*this->nuEff())*dev2(T(fvc::grad(U))))
      - fvm::laplacian(this->alpha_*rho*this->nuEff(), U)
    );
}

代码里面确实包含了密度rho和孔隙率alpha。由于我发现由于是不可压缩流体，速度方程代码两端都是除以密度rho表示的，然而应力求解出来却没有除以rho，这是不是与物理方程相违背

是不是速度方程应该改为

    fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
  - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
  + (1.0/rhoc)*continuousPhaseTurbulence->divDevTau(Uc)//这一项除以密度rho
 ==
    (1.0/rhoc)*cloudSU

比较正确？

最近看DPMFoam的求解器源代码，速度方程的代码一直没看懂，速度方程文件中的代码如下

fvVectorMatrix UcEqn
(
    fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)
  - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
  + continuousPhaseTurbulence->divDevTau(Uc)
 ==
    (1.0/rhoc)*cloudSU
);

UcEqn.relax();

fvOptions.constrain(UcEqn);

volScalarField rAUc(1.0/UcEqn.A());
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc));

surfaceScalarField phicForces
(
   fvc::flux(rAUc*cloudVolSUSu/rhoc) + rAUcf*(g & mesh.Sf())
);

if (pimple.momentumPredictor())
{
    solve
    (
        UcEqn
     ==
        fvc::reconstruct
        (
            phicForces/rAUcf - fvc::snGrad(p)*mesh.magSf()
        )
    );

    fvOptions.correct(Uc);
}

而实际的物理方程为

问题总结如下

1、代码

fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)

表示物理方程的左边两项好像比较清晰，但是代码

- fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)

表示的意义有点难以理解。其中

fvc::ddt(alphac) + fvc::div(alphaPhic)

貌似恰好是连续性方程？

2、应力项代码

+ continuousPhaseTurbulence->divDevTau(Uc)

为何与孔隙率alphac无关，好像和实际物理方程无法对应（左边第三项）？

求解！谢谢。

@马乔 非常感谢马乔老师的分享！我有一个问题想请教下老师，这里湍流模型需不需要用相分数alpha来修正？

老师您好，我想请教一个问题，我按照CFD Online 上关于 coalCollidingChemistryFoam的讨论，在OpenFOAM5.0 的coalChemistryFoam中添加碰撞模型，但是编译之后出现了下面这样的错误，是不是CFD Online上的方法不适用于OpenFOAM5.0？还是有其他需要修改的地方？

wmakeLnInclude: linking include files to ./lnInclude
could not open file surfaceWriter.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to No such file or directory 
could not open file makeParcelHeatTransferModels.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to Too many open files
could not open file makeReactingMultiphaseParcelCompositionModels.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to Too many open files
could not open file makeReactingParcelPhaseChangeModels.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to Too many open files
could not open file makeReactingMultiphaseParcelDevolatilisationModels.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to Too many open files
could not open file makeReactingMultiphaseParcelSurfaceReactionModels.H for source file coalParcel/makeCoalParcelSubmodelsAdd.C due to Too many open files
In file included from /public/software/apps/openfoam/OpenFOAM-5.0/src/lagrangian/intermediate/lnInclude/FacePostProcessing.H:182:0,
                 from /public/software/apps/openfoam/OpenFOAM-5.0/src/lagrangian/intermediate/lnInclude/makeParcelCloudFunctionObjects.H:31,
                 from coalParcel/makeCoalParcelSubmodelsAdd.C:27:
/public/software/apps/openfoam/OpenFOAM-5.0/src/lagrangian/intermediate/lnInclude/FacePostProcessing.C:29:27: fatal error: surfaceWriter.H: No such file or directory
 #include "surfaceWriter.H"
                           ^
compilation terminated.
make: *** [Make/linux64GccDPInt32Opt/coalParcel/makeCoalParcelSubmodelsAdd.o] Error 1