CFD中文网

zhe

        model1              
        {
            type            cellZoneInjection; //for injection at the patch
            SOI             14.15;   // Start time of injection
            parcelBasisType fixed;// parcel basis enumeration
            //patch       top;     
            cellZone       insandZone;
            //duration        100;           
            nParticle       1;   // number of particle per parcel
            numberDensity    5e7;//numberDensity multiplied by cellZoneVolume gives the number of particles to be injected
            U0              (0 0 0);
            sizeDistribution
            {
                type        normal;
                normalDistribution
                {
                    expectation 5e-4;
                    variance 1.25e-4;
                    minValue 2e-4;
                    maxValue 8e-4;
                }
            }
        }

你可以用这个，是可以的。

zhe

如题，在两小球碰撞过程中，绕接触点旋转，相对角速度是否可以直接用两小球自身角速度作差？

zhe

@wangtao121 之后发现，在openfoam中multiphase的solver中，有一个MPPICInterFoam的solver，把其中的MPPIC换成DPM就可以解决这个问题。实际上是在interfoam里加入了粒子效果。

zhe

@鲸落其实也不是只显示一层。我想你跟我当时一样，感觉摄入粒子很多，其实由于体积太小，所占的相分数其实不大。那么最后就会只有一层。如果有很多粒子的话，最后会形成很多层，因为只要满足你预设的最大相分数，它就不会再允许粒子进入，那么就会在其上方的网格中开始堆积。

zhe

如题，最近一直想不明白，在MPPIC中，怎么可以进行alpha和alphaPacked做差的。一个是fieldScalar，一个是scalar.不知道是不是哪里漏掉了，希望有人帮忙解惑，谢谢。

zhe

如题目所示，由于不想在高浓度的粒子群中存在碰撞，所以添加了一个体积分数theta的判定。当theta>0.5时，修正粒子速度为0。当theta<=0.5时，正常发生碰撞。下面是添加后的code，编译是通过的，但是出来的结果完全不对。theta的定义是第二个code。希望有经验的朋友帮忙解惑，谢谢拉。

    label nSubCycles = collision().nSubCycles();

    if (nSubCycles > 1)
    {
       Info<< "    " << nSubCycles << " move-collide subCycles" << endl;

       volScalarField::Internal& thetaIn = theta().ref();
       forAll(thetaIn, celli)
       {
           if( thetaIn[celli] > 0.5 )  //if volume fraction > 0.5,
           {
               nonMoveCollide(cloud, td, this->db().time().deltaTValue());
           }
           else
           {
               subCycleTime moveCollideSubCycle
               (
                   const_cast<Time&>(this->db().time()),
                   nSubCycles
               );

               while(!(++moveCollideSubCycle).end())
               {
                   moveCollide(cloud, td, this->db().time().deltaTValue());
               }

               moveCollideSubCycle.endSubCycle();
           }
       } 
    }    
    else
    {
         moveCollide(cloud, td, this->db().time().deltaTValue());
    }
}

这是theta的定义：

Foam::CollidingCloud<CloudType>::theta() const
{
    tmp<volScalarField> ttheta
    (
        new volScalarField
        (
            IOobject
            (
                this->name() + ":theta",
                this->db().time().timeName(),
                this->db(),
                IOobject::NO_READ,
                IOobject::NO_WRITE,
                false
            ),
            mesh_,
            dimensionedScalar(dimless, Zero),
            extrapolatedCalculatedFvPatchScalarField::typeName
        )
    );
    volScalarField& theta = ttheta.ref();
    forAllConstIter(typename CollidingCloud<CloudType>, *this, iter)
    {
        const parcelType& p = iter();
        const label celli = p.cell();

        theta[celli] += p.nParticle()*p.volume();
    }

    theta.primitiveFieldRef() /= mesh_.V();
    theta.correctBoundaryConditions();

    return ttheta;
}

zhe

在opemfoam中有一个内聚力，可以使粒子之间相互黏结。它是通过粒子之间的碰撞的能量转换实现的。其中通过一个叫cohesionEnergyDensity的系数来决定它的强度。这个密度越大，黏结越明显。希望我说的这个是你想问的。

zhe

麻烦帮忙改正一下。想要规定thetaBoundary为一个常数值，然后对theta（相分数）进行对比，总是实现不了。谢谢。

const Field<scalar>& theta();
const Field<scalar>& thetaBoundary();
thetaBoundary = {0.5};
if (theta < thetaBoundary)

zhe

@星星星星晴我是写在了cloud里面，因为在cloud里我会用到这个平均速度。但是感觉不对，一加进去这段，我的粒子就静止不动，连重力的影响都没有了。感觉好像出速度是0，然后就一直保持0了似的。

zhe

大家好，如题，想要得到网格内的平均速度。我编写的是类似这样的，但是得到的不是网格内的平均速度，麻烦有经验的朋友给更正一下，谢谢。

forAllIter(typename CloudType, this->owner(), iter)
{
    typename CloudType::parcelType& p = iter();
    const label cellI = p.cell();
    const fvMesh& mesh = this->owner().mesh();

    // average velocity of the lagrangian particle [m/s]
    volVectorField particleAverageVelocity
    (
        IOobject
        (
            "particleAverageVelocity",
            this->owner().db().time().timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedVector("zero", dimVelocity*dimMass/dimMass, vector::zero)
    );

    forAll(particleAverageVelocity,i)
    {
        particleAverageVelocity[i] +=p.nParticle()*p.mass()*p.U(); 
    }
    massTotal_ += p.nParticle()*p.mass();
    p.U() =(particleAverageVelocity[cellI] * (1.0/massTotal_));
}

zhe

@zhe

Latest posts made by zhe