CFD中文网

xzdbjx

谢谢您的分享

xzdbjx

@number44 我也遇到了相同的问题，试了您的方法，还是不行

xzdbjx

@陈冠江感谢您的回复，我修改下再看看结果。

xzdbjx

各位老师专家好，我目前在做ABL的LES湍流入口时也碰到一些问题想和大家交流学习。我之前一直用Fluent来生成LES的湍流入口，近半年开始对Openfoam感兴趣，于是开始尝试把Fluent的相关算法移植到Openfoam当中，我采用的算法是华南理工杨易老师课题组提出的NSRFG算法，参考文献：
[1] Yu Y, Yang Y, Xie Z. A new inflow turbulence generator for large eddy simulation evaluation of wind effects on a standard high-rise building[J]. Building and Environment, 2018, 138: 300-313.
先把Fluent和of的计算结果图贴两张，我主要对PSD和湍流强度与平均风剖面进行了验证。
计算域
计算域.png

Fluent流场图及验证图

P1 是入口处的测点，P2是模型处的测点
Openfoam流场图与验证图

这里貌似of的结果还行，但实际上我认为Fluent的结果更合理，因为数值耗散等因素的影响，湍动能一定会在计算过程中耗散，因此模型处的PSD一定会更加吻合目标值。所以我目前怀疑我的case设置有问题，目前我采用的LES是WALE模型，0文件夹下应该要设置k 和 nut两个文件，我的设置分别如下：

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0.0001;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    ground
    {
		type	        kLowReWallFunction;		//For high and low Re
        value           $internalField;
    }
    top
    {
        type            symmetry;
    }
    side1
    {
        type            symmetry;
    }
    side2
    {
        type            symmetry;
    }
    outlet
    {
        type            zeroGradient;

    }
}


/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 1e-8;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    ground
    {
		type            nutUSpaldingWallFunction;
        value           uniform 1e-8;
    }
    top
    {
        type            symmetry;
    }
    side1
    {
        type            symmetry;
    }
    side2
    {
        type            symmetry;
    }
    outlet
    {
        type            zeroGradient;

    }
}


// ************************************************************************* //

这里我不太确定两个文件中internalField参数应该如何取值？另外，这两个参数对结果影响是否很大？请各位老师专家解答。谢谢！（P.S. of计算速度是真的快）

xzdbjx

@xzdbjx

xzdbjx 发布的最新帖子