CFD中文网

LaTex用的人也挺多的

感觉自带的跟好用，键盘输入，我用mathtype编号

好吧，既然没有人回答，我就把我的笨办法贴出来
foamToVTK-->paraview-->打开算例文件和粒子Cloud文件-->File/Save Data-->得到每个时间步每个粒子信息.csv-->用python提取有用信息

如题，在一个lagrangian求解器中，设置一个面的粒子碰撞为escape

patchInteractionModel  standardWallInteraction; 
standardWallInteractionCoeffs
    {
        type            escape;
    }

那么怎么得到粒子escape时的速度？用ParaView能实现吗？或者用cloudFunctions?求知情人士告知，不甚感谢！

我觉得这本书不错Multiphase Flows with Droplets and Particles (Second Edition)-Clayton T. Crowe，网上有电子版，里面介绍离子的受力运动和粒子传热相变比较详细。

有没有用LES算一下，结果怎么样？

@kkguo http://120.25.245.130:4567/topic/450/有谁知道怎么在openfoam3-0-1上编译foamtotecplot360-求教/9

把我们办公室的另两位拉下水了，终于不是我一个苦逼的在学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

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

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

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

在一个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

@cfd-china 谢谢！我找到了，好像是主文件中的parcel.evolve()。

在lagrangian求解器中只能看到连续相的方程：

    tmp<fvVectorMatrix> UEqn
    (
         fvm::ddt(rho, U) + fvm::div(phi, U)
       + MRF.DDt(rho, U)
       + turbulence->divDevRhoReff(U)
      ==
         rho.dimensionedInternalField()*g
       + parcels.SU(U)
       + fvOptions(rho, U)
    );

请问离散相方程怎么找到？谢谢！

@OF九月九 好吧，有个sphereToCell

@wwzhao 你好，如果我想设置一个球形区域，应该在setFields中怎么设置？

@wwzhao 明白了，我思维定势了，谢谢！

在两相流中，如果设置成不同高度的粒子浓度不同，我在setFields是这样设置的

defaultFieldValues
(
    volScalarFieldValue alpha.air 1
    volScalarFieldValue alpha.particles 0
);

regions
(
    boxToCell
    {
        box (0 0 0) (35 35 1);
        fieldValues
        (
            volScalarFieldValue alpha.air 0.8
            volScalarFieldValue alpha.particles 0.2
        );
        box (0 0 1) (35 35 3);
        fieldValues
        (
            volScalarFieldValue alpha.air 0.9
            volScalarFieldValue alpha.particles 0.1
        );
    }
);

即设置了三层不同浓度的粒子，但是setFileds后

Reading setFieldsDict

Setting field default values
    Setting internal values of volScalarField alpha.air
    Setting internal values of volScalarField alpha.particles

Setting field region values
    Adding cells with center within boxes 1((0 0 1) (35 35 3))
    Setting internal values of volScalarField alpha.air
    Setting internal values of volScalarField alpha.particles

End

只默认设置了第二层的浓度

这是怎么回事？是不是这样设置不行？应该怎么办？求教？

@cfd-china :happy: THANK YOU!

如题所述，如果我运行

mpirun -np 4 ***Foam -parallel

在处理器设置为1或4的时候都没有出现错误。是否算例中processor的数目与虚拟机处理器的数量设置要一致？与主机的处理器数量又有什么关系？主机的处理器数量是否会限制虚拟机处理器数量的设置？求教！

@队长别开枪 果然如此，明白了，非常感谢！！！:cheeky: