关于网格高度变化的动网格设置

xiaolin

@东岳
是的，李老师，在x和y方向沉淀都不一致，但是我对编程这块比较薄弱，所以好难啊
那如果是预定变化轨迹的网格高度的动网格变化是不是应该参照movingConeTopoFvMesh呢，还是有别的类型更优，才刚入手动网格这块，循序渐进吧

李东岳

那如果是预定变化轨迹的网格高度的动网格变化是不是应该参照movingConeTopoFvMesh呢

是的，这样简单的多

upc_ngh

@李东岳东岳老师，您好，目前我遇到了和上面那个同学类似问题@xiaolin ，具体问题如下：二维内，气固两相流动过程中，固态颗粒和壁面（半圆形壁面）碰撞，粘连形成沉积。目前的想法是当贴近壁面的第一层网格内，固相体积分数大于某值时，壁面的边界条件就向上移动一个单元，就是上面"alpha与point的关系“。目前发现相关动网格编程资料很少，盼望各位老师能给点建议或者参考资料学习一下，提前谢谢各位大神！
图片1.jpg

田畔的风

我写过类似的东西，网格可以根据当前的场值进行变动，用于研究风场作用下的积雪沉积/侵蚀。代码和算例可以参考[https://github.com/fightingxiaoxiao/driftScalarDyFoam]，如果有用的话，也可以引下我的论文[https://doi.org/10.3389/feart.2022.822140]。

大致的思路是这样的：

在dynamicMeshDict中定义速度相关的动网格：

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs (fvMotionSolvers);

motionSolver velocityLaplacian;   //displacementLaplacian;
diffusivity  quadratic inverseDistance 3(bottom.snow roof top);

定义pointMotionU边界条件为codedFixedValue

    bottom.snow
    {
        type            codedFixedValue;
        value          uniform (0 0 0);

        name         erosionDeposition;
    }

最后在codeDict嵌入代码，这里比较大的一个坑就是动网格使用节点值进行位移，但是实际计算拿到的都是面心值和体心值，需要做一个插值。另外，deltaH是我自己求解器里定义的场，需要替换成你自己的内容。

erosionDeposition
{
    code
    #{
        
        label patchIndex = patch().index();
        const volScalarField& deltaH = this->db().objectRegistry::lookupObject<volScalarField>("deltaH");

        primitivePatchInterpolation facePointInterp(deltaH.mesh().boundaryMesh()[patchIndex]);   //初始化插值类

        const scalarField& deltaHp = deltaH.boundaryField()[patchIndex];
        auto deltaHpp = facePointInterp.faceToPointInterpolate(deltaHp);                      //面心向节点插值
        Info << "Successfully interpolate "<< deltaHpp().size() << " points." << endl;
        Info << "Local points number = " << patch().localPoints().size() << endl;
        pointField pVec(deltaHpp().size());

        const scalar lowerLimit = 0.;

        forAll(pVec, i)
        {
            pVec[i][0] = 0;
            pVec[i][1] = 0;
            pVec[i][2] = deltaHpp()[i];

            if(patch().localPoints()[i][2] + pVec[i][2] * this->db().time().deltaTValue() < lowerLimit)
            {    
                pVec[i][2] = (lowerLimit - patch().localPoints()[i][2])/this->db().time().deltaTValue();
            }
        }

        (*this) == pVec;
    #};

    codeInclude
    #{
        #include "fvCFD.H"
        #include "primitivePatchInterpolation.H"
    #};

    codeOptions 
    #{
        -I$(LIB_SRC)/finiteVolume/lnInclude \
        -I$(LIB_SRC)/meshTools/lnInclude 
    #};

最终效果如下，一个方块周边的时变积雪

李东岳

@田畔的风大佬这个有点意思，感兴趣写个CFD界不

田畔的风

@李东岳可以的，李老师！不过做这部分内容的时候我对OF的认知还比较浅，很多细节可能还比较糙..

李东岳

@田畔的风大佬，参考这种风格，链接文本，如果你可以写一个比较容易让人理解的文章的话，可以发给我我在CFD界推一下

upc_ngh

@田畔的风非常感谢两位老师@李东岳，感觉久旱逢甘露，话不多说，我要开始汲取营养了，后面有相关问题和进展，再向老师们请教

疏影横斜水清浅

@田畔的风老师，你好，请问你在代码中求解壁面剪切力是采用的是volSymmTensorField Reff = turbulence->devReff(); ssp = (-Sfp/magSfp) & Reffp; scalarField UShear = sqrt(mag(ssp) / rhoAir);
这其中的Reff中有没有密度相？

田畔的风

@疏影横斜水清浅 OpenFOAM不可压求解器的密度是被约掉的，我这里的空气密度取的是1。当时的写法不太严谨，没检查量纲。

疏影横斜水清浅

@田畔的风好的，谢谢。我是看到你在后面又把空气密度除了一次，所以存在疑惑。如果你取空气密度为1，那么在数值上应该是没问题的。

CFD中文网

关于网格高度变化的动网格设置