CFD中文网

@liuzhuang
非常感谢您的回答。
您的意思也就是说TRrf取300就可以了是吗？我自己的case是一个带传热的喷射流，这样的话入口和出口的p_rgh该怎么取呢？
另外想问您一个问题，system文件夹下的setFields有什么作用，该如何填写？

最近老师让我用buoyantBoussinesqSimpleFoam模拟一个带传热的case，我参考hotRoom这个算例进行设置，碰到几个问题不太懂，恳求各位老师指点。

1，

constant文件夹下的transProperties里面的Reference temperature是什么意思，Pr和Prt该怎么取？

2，

0文件夹下的p_rgh是什么意思，各个位置的值该怎么取？

3，

0文件夹下的alphat是什么意思，各个位置的值该怎么取？

4,

为什么这个算例里面0文件夹下p的壁面不取zeroGradient？

各位大家好，最近我老师让我用C++实现如上图所示的空腔内的二维稳态层流的流场。我自己弄了几天，但是始终没办法得到正确的结果。各位有相关的参考资料吗？

@hy1112006 非常感谢您的回答，这篇文章解决了我的疑惑。请问您对turbulentInlet这个入口边界条件是如何理解的呢？我在网上找到的都是一些相当模糊的解释，例如“通过添加场的随机分量来表现入口边界的波动”，但是具体如何获取fluctuationScale却都没解释清楚。

各位老师，大家好。我在使用LES进行模拟时遇到了两个困惑，请各位老师指点。

1.LES的边界条件设置时，U里面的fluctuationScale的具体意义,该怎么设置？

2.fieldAverage这个函数是在多长时间上取平均值？

@东岳 老师，那为什么pitzDaily使用二维的LES和RANS的结果吻合的很好呢？

@东岳 不，这个case是二维的。

@东岳 非常感谢东岳老师的帮助。另外我还有一个问题，我计算时取不同的时间步长会得到相差非常大的结果，有的时间步长还会报错。目前我的尝试的结果是，0.1s的步长比0.01s和0.001s的步长计算结果准确，但是和实验结果相差还是挺大的，没有k-epsilon吻合的好。老师，您能帮我分析一下吗？

各位老师大家好，

我之前使用k-epsilon模型计算过上图所示的空腔内流场，结果和实验数据吻合的很好。现在想使用LES模型来计算上图所示的空腔内的流场。并且将结果和实验数据、k-epsilon模型计算结果进行对比。已知的条件是入口处速度场为（0.455 0 0）m/s，假定湍流强度为4%
我模仿tutorials里面pitzDaily算例进行计算，选择了dynamicKEqn模型，求解器使用pisoFoam。现在碰到了以下几个疑问。

1、LES模拟中如何确定0文件夹下k的初值，是否像k-epsilon一样存在计算k初值的公式？
2、为什么pitzDaily用LES模型计算到0.1s的结果取时间平均得到的流场和使用RANS模型计算到297s得到的结果是一致的？
3、为什么我所计算的空腔使用LES模型计算到0.1s的结果取时间平均得到的流场并未充分发展？
4、k-epsilon和LES结果进行对比时，应该取LES什么时候的结果进行对比？

恳求各位老师赐教，不胜感激！

@东岳 老师，是用filedAverage这个函数来获取LES流场的平均p、U、k吗？但是用LES计算时要求时间步很小，计算到RANS稳态结果的时间计算量非常大

@东岳 东岳老师，“去LES/DNS流场的平均结构，可以和稳态的对比”这句话是什么意思呢，能解释详细点吗?

@东岳 老师，LES可以得到pitzDaily稳定下来的结果吗？

@东岳 老师，tutorials里面simpleFoam里的pitzDaily是用RAS（k-epsilon）算的，我需要用LES模型计算pitzDaily，得到稳定的结果，然后和RAS的结果进行对比。

各位老师大家好，
我想分别使用RAS和LES模型计算pitzDaily这个算例，并进行结果对比。求解器用的是simpleFoam。tutorials/incompressible/pisoFoam/les/pitzDaily里面pitzDaily用LES模型计算时用的是pisoFoam，我对fvsolution和controlDict文件进行了一点修改，想用simpleFoam计算得到稳定之后的结果。但是运行simpleFoam后计算一直无法无法稳定下来（无论计算时间设置多长）。以下是我修改的内容:

fvsolution

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-06;
        relTol          0.1;
        smoother        GaussSeidel;
    }

    pFinal
    {
        $p;
        smoother        DICGaussSeidel;
        tolerance       1e-06;
        relTol          0;
    }

    "(U|k|B|nuTilda|s)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-05;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-2;
        "(pFinal|U|B|nuTilda|s)" 1e-3;
    }
}

relaxationFactors
{
    equations
    {
        U               0.9; // 0.9 is more stable but 0.95 more convergent
        ".*"            0.9; // 0.9 is more stable but 0.95 more convergent
    }
}

controlDict

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         2000;

deltaT          1;

writeControl    timeStep;

writeInterval   100;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

谢谢各位指教

在这个算例里面，我可以找到更密和更疏的网格使它不发生这样的错误。为什么特定的网格会出现这样的结果呢

请教各位老师，使用simpleFoam时出现了这个错误是什么意思？

@东岳 老师，是这个volScalarField吗？该怎么声明呢？

我在使用simpleFoam运行这个模型时出现了这样的错误，这是为什么？

--> FOAM FATAL ERROR: 

    request for volScalarField::Internal LamBremhorstKE:G from objectRegistry region0 failed
    available objects of type volScalarField::Internal are

8
(
nut
yWall
k
nu
p
(sqr((1-exp((-0.0165*((sqrt(k)*yWall)|nu)))))*(1+(20.5|((sqr(k)|(nu*epsilon))+1e-15))))
(sqr(k)|(nu*epsilon))
epsilon
)


    From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::DimensionedField<double, Foam::volMesh>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::DimensionedField<double, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::DimensionedField<double, Foam::volMesh> >(Foam::word const&) const at ??:?
#3  Foam::epsilonLowReWallFunctionFvPatchScalarField::calculate(Foam::turbulenceModel const&, Foam::List<double> const&, Foam::fvPatch const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#4  Foam::epsilonWallFunctionFvPatchScalarField::calculateTurbulenceFields(Foam::turbulenceModel const&, Foam::Field<double>&, Foam::Field<double>&) at ??:?
#5  Foam::epsilonWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#6  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() at ??:?
#7  Foam::incompressible::RASModels::LamBremhorstKE::correctNut(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8  Foam::incompressible::RASModels::LamBremhorstKE::correctNut() at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?

@dzw05 就是说不是所有的文件都有用，只需要填写自己使用模型的变量的边界条件就可以是吗（在LB模型中只填写nut，k，epsilon，p，u）