lagrangian两相流中粒子质量负载的计算。

[OF九月九]

在一个lagrangian算例中，kinematicCloudProperties如下

solution
{
    active          true;
    coupled         true; //false;
    transient       yes;
    cellValueSourceCorrection on;

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        mu              cell;
    }

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
        schemes
        {
          U               semiImplicit 0.5;
        }
    }

}

constantProperties
{
    parcelTypeId    1;

    rhoMin          1e-15;   //parcel的最小密度
    minParticleMass 1e-15;   //parcel的最小质量

    rho0            1250;    //单个颗粒的密度
    youngsModulus   6e8;
    poissonsRatio   0.35;

    constantVolume  true;
}

subModels
{
    particleForces
    {
        sphereDrag;
        gravity;
    }
	
    injectionModels
    {
        model1
        {
            type            manualInjection;
            massTotal       0;
            parcelBasisType fixed;
            nParticle       1; //每个parcel的颗粒数，设成1就是DEM
            SOI             0;
            positionsFile   "kinematicCloudPositions";
            U0              ( 0 0 0 );
            sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   0.00125;
                }
            }
        }
    }

    dispersionModel none; //stochasticDispersionLES;

    patchInteractionModel standardWallInteraction;

    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel none;
	
    stochasticCollisionModel none;

    radiation off;

    standardWallInteractionCoeffs
    {
        //type            stick;
        type            rebound;
	     e               1;
        mu              0;
    }
}


cloudFunctions
{

}

kinematicCloudPositions如下：

50000
(
......
)

计算域为35m×35m×5m的长方体。Euler相密度为1.205kg/m3,则粒子的质量负载为

请问这样计算正确吗？谢谢？
@散漫守望2016

[chpjz0391]

你是想计算颗粒质量所占气体质量的百分比么？ 那你应该计算颗粒的体积而不是面积吧。 你这里计算的 π*(d/2)^2是面积不是体积。。看你的算例是不考虑颗粒之间的碰撞吧。。

[OF九月九]

首先，谢谢您的回复！的确，算例是不考虑颗粒之间的碰撞。按您说的修改成

这样算正确吗？计算结果质量分数太小了，总觉得和论文中的对不上。论文中是这样描述的：

难道是我的颗粒数50000设的太少了吗？而且如果我想要改变颗粒的质量分数，直接改变颗粒的密度就可以吗？谢谢！

[chpjz0391]

@OF九月九 你的颗粒数太少了。所以导致质量分数很小。你想一下35米乘以35米再乘5的空间。才5万个 1.25毫米的颗粒。 肯定不会大。 这样的话 体积分数更小了。在实际计算中可以忽略了体积分数了。并且这种情况下碰撞可以不用考虑。。 这篇文章可否共享下。

[OF九月九]

@chpjz0391 谢谢！应该怎么改进？增加kinematicCloudPositions吗？还是改变nParticle？
文章 http://xueshu.baidu.com/s?wd=02Sea+surface+drag+and+the+role+of+spray&rsv_bp=0&tn=SE_baiduxueshu_c1gjeupa&rsv_spt=3&ie=utf-8&f=8&rsv_sug2=1&sc_f_para=sc_tasktype%3D{firstSimpleSearch}&rsv_n=2

[chpjz0391]

@OF九月九 nParticel 是每个parcel里面包含的颗粒数。。 主要还是你的颗粒总数。关键看你要模拟什么了，这个数字要符合实际不是说瞎改的。。比如你实际系统中就5w个颗粒 。你怎么改体积分数和质量分数都很小。。