CFD中文网

@李东岳 intensive variable是密集变量的意思么？图中那个截面的网格进行了网格加密，加密尺寸1mm。处理的变量是气相温度、组分质量分数这些。

@李东岳 谢谢李老师。请教下这样的表达式最后获得的应该是一个独立的数值，而不是流场内某一平面各个网格对应的数值?如果是这样的话，散点图就无法进行绘制啊，不知道理解对不对

@LiuWenchao 好的，非常感谢！

在使用PaSR模型计算燃烧流场时，根据一个截面上某一时刻的网格点出了瞬时混合分数、温度和CH2O的散点图。审稿人说网格分布不均匀（左图），这样的处理方式不合理，让用conditionally volume-averaged quantities方式进行后处理。请问这种方式如何后处理呢？我现在把数据点取为时均值，也就是平均的混合分数、温度和CH2O（右图），请问这种方法合理么？

@LiuWenchao 您好，请问可以分享下哪篇文章对甲烷层流火焰进行了验证，证明了δ/sf对当量比依赖性较低啊？谢谢！

@wangfei9088 对，是时间常数。您的意思是在PaSR模型中该常数与当量比无关？但是时间常数是由流动和化学反应特征时间定义的啊，火焰传播速度和厚度却又和当量比相关，感觉有点矛盾呢

PaSR模型中的模型常数（Cmix）取决于流动特征时间、化学反应特征时间，也就是和层流火焰速度和厚度相关。当流场内某个区域当量比变化非常剧烈的时候，比如从1.5降低至0.5，这时候层流火焰速度和厚度怎么确定的呢？

@Tens 前辈，请问您解决了么？我也在做输出每一个组分对应的扩散项，没弄出来

前期用sprayFoam计算了喷雾燃烧，审稿人提问计算过程中是否限制了parcel的最小数目和最小粒径，请问大家知道OpenFOAM中怎么实现的么？
当时并没有考虑到这方面的因素，只是根据文献（Srinivasan et al., Quality and reliability of LES, Springer, 2011）计算了下每个网格对应的parcel数目小于8，便开启流场计算了。目前的喷嘴设置模型如下图，破碎模型选用默认的ReitzDiwakar Model

@liuwenchao 好的，谢谢前辈！

@liuwenchao 请问前辈您有和输出焓值和反应速率相关的资源推荐么？我在文献上看到的也是您说的这个公式，一直没找到实现方法

@fkingdom 请问前辈找到方法了么？可以指导下么？我最近也在尝试后处理提取低温反应和高温反应的放热率，思路是每个反应的化学反应速率乘以其焓值，具体实现方法还在探索。

@星星星星晴 好的 我查下原始文献。谢谢前辈！

@星星星星晴 感谢前辈悉心指导！关于您说的问题，我今天查阅了一些公式，关于您提的四个问题，具体回复如下：

1. 根据图1中液滴阻力的公式，在100%和200%air中的阻力分别为3.7e-6N和2.3e-5N，也就是液滴受到的阻力变大了，是不是说明液滴更倾向于保持目前的球形状态呢？
2. 破碎模型用的是ReitzDiwakar model,具体公式如图中所示，计算得到100%和200%air中的rstable分别为18μm和1.2μm,也就是说后者的液滴粒径更小。这里您说的破碎模型如何处理质量是什么意思呢？是指公式中的密度么？
3. “Case2中的液滴数目（nParticles）”说错了，应该是流域内的parcels的数目，也就是图2中的8370.计算结果显示200%air case中的parcels数目比100%air case少了近一半。
4. Flowrateprofile constant已经改正。
   再次感谢前辈指点！
   
   

@星星星星晴 谢谢前辈！100%和200%代表计算域进口空气流量的相对值，比如100%是1kg/s的话，200%就是2kg/s。流场内空气速度增大了两倍左右，没有开启碰撞模型，sprayCloudProperties文件内容具体如下：

\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       dictionary;
    location    "constant";
    object      SprayCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;
    coupled         true;
    transient       yes;
    cellValueSourceCorrection on;
    maxCo           0.3;

    sourceTerms
    {
        schemes
        {
            rho             explicit 1;
            U               explicit 1;
            Yi              explicit 1;
            h               explicit 1;
            radiation       explicit 1;
        }
    }

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        thermo:mu       cell;
        T               cell;
        Cp              cell;
        kappa           cell;
        p               cell;
    }

    integrationSchemes
    {
        U               Euler;
        T               analytical;
    }
}


constantProperties
{
    T0              300;

    // place holders for rho0 and Cp0
    // - reset from liquid properties using T0
    rho0            744;
    Cp0             2090;
    Tvap            298;
    Tbp             489;

    constantVolume  false;
}


subModels
{
    particleForces
    {
        sphereDrag;
    }

    injectionModels
    {
        model1
        {
            type            coneNozzleInjection;
            SOI             0.029;
            massTotal       0.0044166;
            parcelBasisType mass;
            injectionMethod disc;
            flowType        constantVelocity;
            UMag            33.0;
            outerDiameter   0.0007;
            innerDiameter   0;
            duration        1.0;
            position        (0 0 -0.03584797);
            direction       (0 0 1);
            parcelsPerSecond 2000000;
            massFlowRate    0.0044166;
            flowRateProfile table
            (
                (0              0.0044166)
                (1              0.0044166)
                (2              0.0044166)
                (3              0.0044166)
                (4              0.0044166)
                (5              0.0044166)
                (6              0.0044166)
                (7              0.0044166)
                (8              0.0044166)
                (9              0.0044166)
                (10              0.0044166)
                
            );

            Cd              constant 0.9;

            thetaInner      constant 25;
            thetaOuter      constant 30;

            sizeDistribution
            {
                type        RosinRammler;

                RosinRammlerDistribution
                {
                    minValue        1e-05;
                    maxValue        6e-05;
                    d               4e-05;
                    n               3;
                }
            }
        }

        

    dispersionModel none;

    patchInteractionModel standardWallInteraction;

    heatTransferModel RanzMarshall;

    compositionModel singlePhaseMixture;

    phaseChangeModel liquidEvaporationBoil;

    surfaceFilmModel none;

    atomizationModel none;

    breakupModel    ReitzDiwakar; // ReitzKHRT;

    stochasticCollisionModel none;

    radiation       off;

    standardWallInteractionCoeffs
    {
        type            rebound;
    }

    RanzMarshallCoeffs
    {
        BirdCorrection  true;
    }

    singlePhaseMixtureCoeffs
    {
        phases
        (
            liquid
            {
                C12H26               1;
            }
        );
    }

    liquidEvaporationBoilCoeffs
    {
        enthalpyTransfer enthalpyDifference;

        activeLiquids    ( C12H26 );
    }

    ReitzDiwakarCoeffs
    {
        solveOscillationEq yes;
        Cbag            6;
        Cb              0.785;
        Cstrip          0.5;
        Cs              10;
    }

/*
    ReitzKHRTCoeffs
    {
        solveOscillationEq yes;
        B0              0.61;
        B1              40;
        Ctau            1;
        CRT             0.1;
        msLimit         0.2;
        WeberLimit      6;
    }
*/
    TABCoeffs
    {
        y0              0;
        yDot0           0;
        Cmu             10;
        Comega          8;
        WeCrit          12;
    }
}


cloudFunctions
{
    voidFraction
    {
        type   voidFraction;
    }
}



}

@星星星星晴 前辈，您好！我在使用sprayFoam喷射液滴的时候，在两个不同进气流量的case进行计算，case1为100%空气，case2为200%空气，所有的constant,system里的文件设置都一样，仅仅改变了进气流量。进行液滴后处理时发现，Case2的液滴数目（nParticles）比Case1少了近一半，请问您知道是什么原因么？液滴喷射模型设置具体如下：

model1
        {
            type            coneNozzleInjection;
            SOI             0.029;
            massTotal       0.0044166;
            parcelBasisType mass;
            injectionMethod disc;
            flowType        constantVelocity;
            UMag            33.0;
            outerDiameter   0.0007;
            innerDiameter   0;
            duration        1.0;
            position        (0 0 -0.03584797);
            direction       (0 0 1);
            parcelsPerSecond 2000000;
            massFlowRate    0.0044166;
            flowRateProfile table
            (
                (0              0.0044166)
                (1              0.0044166)
                (2              0.0044166)
                (3              0.0044166)
                (4              0.0044166)
                (5              0.0044166)
                (6              0.0044166)
                (7              0.0044166)
                (8              0.0044166)
                (9              0.0044166)
                (10              0.0044166)
                
            );

            Cd              constant 0.9;

            thetaInner      constant 25;
            thetaOuter      constant 30;

            sizeDistribution
            {
                type        RosinRammler;

                RosinRammlerDistribution
                {
                    minValue        1e-05;
                    maxValue        6e-05;
                    d               4e-05;
                    n               3;
                }
            }
        }

目前计算一个喷雾燃烧场时，想改变下燃油喷射位置，因此把几何模型的喷嘴结构(5mm圆孔)向上游移动了10mm, map前在system文件夹下创建了mapFieldsDict文件，具体信息如图1所示；然后在对应时刻文件夹下建立了和原始计算场相同时刻的文件（共125个）。在map时使用了1）mapFields /OF6/P25；2）mapFields /OF6/P25 -sourceTime 0.116555 -mapMethod mapNearest -targetRegion region0两种方式均出现phi未能成功map,同时PaSR_kappa和PaSR_Qdot，sprayCloud类文件除了sprayCloudTheta文件，其余均未能map过来，如图2所示。另外在现有map结果基础上试算了下，log和err文件如图3所示，在读取化学反应机理时浮点溢出了，尝试了foamChemistryReader和new format均未能解决。请问这俩问题大家有什么方法解决么？万分感谢！

@ir77 解决了，需要设置一个误差限，使if语句里二者的值小于误差限就默认二者相等了。

@星星星星晴 非常感谢！我找找看

@lts 您好！请问您是改成FoamChemsitryRead后成功了么？我也遇见您图里的错误，改成FoamChemsitryRead后依然报错，可否指教下？非常感谢！