CFD中文网

@capillaryFix 好的，太感谢了，希望有机会跟您进行广泛交流

@capillaryFix 您好，我刚刚认真看了您的两篇文章，打算参考您的光滑方式结合半拉格朗日方法试试看？日本小组的文章中结果确实很好，我个人能力有限，当时也没有复现出来，包括解析方法的表面张力的引入。后面我使用别的方法考虑接触角的准确性，暂时忽略了表面张力对颗粒的作用，也就暂时停止了这种外推的计算。

@capillaryFix 太好了，感谢感谢，我会好好看看。我目前一直使用的消除寄生流的方法是基于sharpen surface model来的，这个效果也是很好。

@EricLiu 不是大佬。这个讨论顶部文章附带一些代码，或者直接github搜索clsvof也能有很多收获。

@Sk 我的问题是与phiIB计算时的局部参考点有关。由于默认设置为0号网格，对每个核来说，都会在角点位置出现一个异常的结果。我采用的方式是对phiIB这个量的计算避免使用非迪利克雷边界，进而避免引入参考点场值。

@郑学习 感谢老师！我有相关进展一定与您分享讨论

@chapaofan 没啥没啥，下次如果自己解决了也可以自己贴出来奥

hard_particles yes是单独的命令，不是跟fix混在一起的。在使用设置物性前使用。

@五好青年 老师，请问您现在找到相关资料了吗。

@队长别开枪 好的，感谢老师！

@chapaofan 哦哦，突然发现之前那个回复不是很全。因为想到大家的输出精度不一样。reconstructParMesh报错一般确实是mergeTol不匹配，不过他会给一个提示的说明这个选项问题，于是我之前的回复中忽略了这个。

@chapaofan 好像是因为开了自适应呢

您好，我想请教一下，您在原本的矩形网格上是否有尝试过削弱这种铺展的行为呢？

@队长别开枪 老师，我想问一个问题。如果我手动移动某个内部点的坐标，这样会导致该点参与的面上的四个点（比如原本是笛卡尔网格）不共面。这种情况下，OpenFoam能够正常计算和处理吗？一般来说是如何处理呢？

推荐一本CFD-DEM基础书籍，感觉其中的unresolved CFD-DEM的理论公式说明比较清晰。
Norouzi 等 - 2016 - Coupled CFD‐DEM Modeling: Formulation, Implementation and Application to Multiphase Flows

有时候重新计算案例，上一次的东西没有删除干净，直接计算后重组似乎也会有这个问题。

@李东岳 李老师，这个是浸入边界法中的锐利边界方法。每次计算的时候是将固相和液相通过切割完全分开以提高计算的准确性。那个球实际上就是固体，球体内部是有网格的。这里只是把固体内部网格抠掉显示了。

@队长别开枪 感谢老师，我再找找资料看看。说不定边界单元变形也能一定程度解决我的问题。

想请教一下大家，可能与IBM更相关。在OpenFoam中是否能够实现手动的网格切割，或者说有什么实现思路与学习方法？主要是想将立方体网格切割，使固体边缘更加光滑。如下图所示。
我没有进行过相关的网格切割测试，很担心在拓扑性那一块就过不去了。个人想法的话是直接阅读OpenFoam中的dynamicRefineFvMesh 或 snappyHexMesh？希望各位老师提一些建议。或者说是否有一些基础的代码能够学习。

谢谢大家

参考文献：
Schneiders et al., 2016, An efficient conservative cut-cell method for rigid bodies interacting with viscous compressible flows.
Schneiders et al., 2015, Sharp resolution of complex moving geometries using a multi-cut-cell viscous flow solver

@李东岳
@coolhhh
感谢两位老师，我提出这个问题的主要原因是在现有的CFD-DEM流固耦合计算中，一种基于FD/IBM的细网格方法是需要对速度场U进行改变，并针对连续性方程进行修正的。而这个连续性修正的过程就需要再次显式构建连续性条件：div(U)=0。当然这是理论上的条件。
$\nabla \cdot (\nabla \phi) = \nabla \cdot U$
而下面代码中的连续性修正因子的phiIB理论上只在颗粒存在的条件下才有意义，同样才有值。但是如果div(U)在颗粒不存在得情况就不为0了，那么这样的修正就是错误的。虽然当前这个方面的主要算法就是下面展示的这样。

目前看来在PISO外部，U和phi是分离的，但是CFD-DEM的修正是针对U的，这样操作应该是有问题的，误差大小有待考察。而我原本的目标是希望能够在PISO外，得到满足连续性的U或者phi，不过目前看来是很难实现的。

上述过程具体代码如下：

void Foam::cfdemCloud::calcVelocityCorrection
(
    volScalarField& p,
    volVectorField& U,
    volScalarField& phiIB,
    volScalarField& voidfraction//颗粒体积分数（颗粒内部为1
)
{

    
    void Foam::cfdemCloudIB::calcVelocityCorrection
    
    (
    
        volScalarField& p,
        volVectorField& U,
        volScalarField& phiIB,
        volScalarField& voidfraction
    
    )
    {
        setParticleVelocity(U);//改变颗粒所在区域流体的速度
        // 修改速度会使其不满足连续性，因此使用一个phiIB进行修正
    
        fvScalarMatrix phiIBEqn 
        (  
            fvm::laplacian(phiIB) == fvc::div(U) + fvc::ddt(voidfraction)
        );
    
        if(phiIB.needReference()) 
    
        {
    
             phiIBEqn.setReference(pRefCell_, pRefValue_);
    
        }
    
        phiIBEqn.solve();

        U=U-fvc::grad(phiIB);
    
        U.correctBoundaryConditions();

        // correct the pressure as well
    
        p=p+phiIB/U.mesh().time().deltaT(); 
    
        p.correctBoundaryConditions();  
    }