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



  • 在一个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,则粒子的质量负载为
    0_1477575204116_QQ截图20161027213301.png
    请问这样计算正确吗?谢谢?
    @散漫守望2016



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



  • 首先,谢谢您的回复!的确,算例是不考虑颗粒之间的碰撞。按您说的修改成
    0_1477635757458_QQ截图20161028142219.png
    这样算正确吗?计算结果质量分数太小了,总觉得和论文中的对不上。论文中是这样描述的:0_1477636363869_QQ截图20161028143207.png

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



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





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


登录后回复
 

与 CFD 中国 的连接断开,我们正在尝试重连,请耐心等待