CFD中文网

@李东岳 文章还没见刊。

@Afr1yne 邮箱：tianshb230306@foxmail.com，你可以发个邮件，我邮件给你微信号

@zzzzz7777 大气边界层的雷诺应力应该怎么给，也有人说这个初始条件不重要，给的差不多后面也可以发展出正常的

@zzzzz7777 我给了一个相同的雷诺应力，不是廓线。

R      uniform (0.8 0 0 0.4 0 0.4);

@sxz0823 请问，

///////还有就是这一项代表了什么意思呢////////
surfaceScalarField phicForces
(
    fvc::flux(rAUccloudVolSUSu/rhoc) + rAUcf(g & mesh.Sf())
);

这个有答案了吗

@李东岳 我先将网格细化做一下测试，排除或者确认一下网格的原因。

@coolhhh 是的，yplus很大，三四千，这个是我在测试的，我是想做全尺寸的，所以想用的是IDDES，IDDES是否还要考虑yplus，我不是很清楚。测试算例，所以使用的网格很粗。

使用kOmegaSSTIDDES湍流模型计算大气边界层风场

simulationType      LES;

LES
{
    LESModel        kOmegaSSTIDDES;

    turbulence      on;

    printCoeffs     on;

    delta           IDDESDelta;

    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
    decayControl    true;//false;
    kInf            0;
    omegaInf        0;
    useSigma        false;
    kappa           0.41;
    CDESkom         0.82;
    CDESkeps        0.6;
    Ctrans          60;
    Cdt1            20;
    Cdt2            3;
    Cl              5;
    Ct              1.87;
    fe              true;
}

风速入口为对数廓线

            scalar kappa = 0.41;
            scalar Uref  = 2.0;
            scalar Zref  = 10.0;
            scalar z0    = 0.001;
            scalar zd    = 0.0;

            scalar Ustar = kappa*Uref/log((Zref+z0)/z0)；
            scalar z  = inletPatch.Cf()[i].z();
            Uf[i].x() = Ustar/kappa*log(((z-zd) + z0)/z0);
            Uf[i].y() = 0.0; 
            Uf[i].z() = 0.0;  

结果为

其中，9个位置分别为10m、200m、400m、500m、600m、700m、800m、900m、990m（计算域流向1000m，垂向100m，展向100m），网格划分(250 25 50) z是垂向，膨胀率 (1 1 5)。
请问，上部分的速度随高度基本不变，是什么原因？怎么才可以使得风速随高度呈对数形态。

请问在使用turbulentDFSEMInlet入口类型计算过程中，随着计算的进行，后半部分的风场慢慢的变得不是对数廓线了，请问是什么原因，有没有大佬知道。
我的初始边界信息：
计算域：xyz为1011（m）
U：

inlet
    {
        type            turbulentDFSEMInlet;
        delta           1;
        U
        {
            type        mappedFile;
            mapMethod   nearest;
        }
        R               uniform (0.8 0 0 0.4 0 0.4);
        L               uniform 0.004;
        d               1;
        nCellPerEddy    1;
        scale           1;
        value           $internalField;
    }
    outlet
    {
        type            advective;
        value           $internalField;
    }
    "(topWall|front|back)"
    {
        type            zeroGradient;
    }
    "(bottomWall|Model)"
    {
        type            noSlip;
    }

输入文件：
point：

(
( 0 0 0.01 )
( 0 0 0.05 )
( 0 0 0.1 )
( 0 0 0.15 )
( 0 0 0.2 )
( 0 0 0.25 )
( 0 0 0.3 )
( 0 0 0.35 )
( 0 0 0.4 )
( 0 0 0.45 )
( 0 0 0.5 )
( 0 0 0.55 )
( 0 0 0.6 )
( 0 0 0.65 )
( 0 0 0.7 )
( 0 0 0.75 )
( 0 0 0.8 )
( 0 0 0.85 )
( 0 0 0.9 )
( 0 0 0.95 )
( 0 0 0.99 )
)

U：

(
( 0 0 0 )
( 1.048097372 0 0 )
( 1.499488341 0 0 )
( 1.763535131 0 0 )
( 1.95087931 0 0 )
( 2.096194744 0 0 )
( 2.2149261 0 0 )
( 2.31531203 0 0 )
( 2.402270279 0 0 )
( 2.47897289 0 0 )
( 2.547585713 0 0 )
( 2.609653562 0 0 )
( 2.666317069 0 0 )
( 2.718442442 0 0 )
( 2.766702999 0 0 )
( 2.811632503 0 0 )
( 2.853661247 0 0 )
( 2.893141184 0 0 )
( 2.930363858 0 0 )
( 2.965573488 0 0 )
( 2.992431707 0 0 )
)

计算3600s后，平均风速廓线，U1U2U3U4U5U6U7U8U9的位置为0.1，2，4，5，6，7，8，9，9.9

你好，请问你编写的NSRFG，在使用时计算域是多大，网格尺度是多大，能否在整个计算域上保持低耗散甚至无耗散么

@脉动通风 你好，能详细说一下么，你是用的哪种颗粒注入方式，内部画网格是指生成面还是cellzone，你使用的方式是连续注入么

@Vortex 请问dynamicLagrangian湍流模型，flm和fmm的边界设置是什么样的

@SHUKK 是壁面函数的原因吗

@SHUKK 请问你这个问题解决了么

@李东岳 是的，我使用的OpenFOAM7，我安装个10试试。

@李东岳 我在controlDict中写入

cacheTemporaryObjects
(
    k
);

functions
{
    #includeFunc writeObjects("k")
}

结果出现错误，请问李老师，能知道是什么原因吗

@李东岳epsilon()函数要不要改

template<class BasicTurbulenceModel>
tmp<volScalarField> tshbSmagorinsky<BasicTurbulenceModel>::epsilon() const
{
   k_ = (this->k(fvc::grad(this->U_)));

   return volScalarField::New
   (
       IOobject::groupName("epsilon", this->alphaRhoPhi_.group()),
       this->Ce_*k_*sqrt(k_)/this->delta()
   );
}

@2019201300 请问，这个代码是能够将每个网格的体积存储在Vcell中吧？

volScalarField Vcell
(
    IOobject
    (
        "Vcell",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("", dimVolume, 0.0)
);

forAll(Vcell ,i)
{
    Vcell[i] = mesh.V()[i];
}

请问这个代码怎么修改，才可以使用网格的体积,

volScalarField muc
(
    IOobject
    (
        IOobject::groupName("mu", continuousPhaseName),
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    rhoc*laminarTransport.nu()
);

volScalarField Rep("Rep", dp/laminarTransport.nu()*mag(Up - U));
volScalarField CDRe("CDRe", 24/Rep+3.6*pow(Rep,-0.313));
scalar deltaTs = runTime.deltaTValue();

Up = Up + deltaTs*(3/4*muc/rhop/pow(dp,2.0)/mesh.V()*CDRe*(U - Up) + g);

其中的mesh.V()怎么修改才可以。

@ShaneHEEE 是的，和出口设置一样，就可以了