dnsFoam初始化问题

[wzchen]

刚刚接触openfoam，问个简单的问题：对于tut里面的dnsFoam的例子，初始化时内部速度场的数据是怎么来呢

[CFD中文网]

@wzchen

你好，你可以把速度设置为uniform后，使用boxTurb程序来生成人工随机量。算例自带的例子只不过帮你完成了这个过程。:laughing:

[wzchen]

:happy: 谢谢，在想问一下

Ea k0的具体含义是什么，如何设置它们的值？有关于openfoam初始化方面的资料推荐一下吗:happy:

[李东岳]

@wzchen

我半个月前在写dnsFoam的资料，不过后来被别的任务打断了。dnsFoam模拟的是一盒子湍流，就是不带壁面的DNS模拟。因为没有壁面所以能量会被消耗掉导致无湍。dnsFoam在源项中添加了人工噪声，噪声通过快速傅里叶变换来生成力来保持湍流不被消耗掉。目前我只进行到这里。并没有对公式进行推导。最近事情非常非常的多。做一个类似的dnsFoam解析可能要2，3个月之后了。如果你感兴趣自己推导最好了！:cheeky:

另外，dsFoam的资料非常少，目前仅有一个查尔姆斯大学的tutorial，我建议你看你专业的DNS文献，里面提及了如何在源项添加力来保证有湍流。

[wzchen]

好的，谢谢:happy:

[TangShangyu]

@wzchen 请问您知道如何设置Ea和k0的值了吗？方便的话，请指教，谢谢。

[李东岳]

他们的计算在这个文件中Ek.H，目前如何指定的理论背景不详。

namespace Foam
{

inline tmp<scalarField> Ek
(
    const scalar Ea,
    const scalar k0,
    const scalarField& k
)
{
    tmp<scalarField> tEk = Ea*pow(k/k0, 4.0)*exp(-2.0*sqr(k/k0));

    /*
    scalarField& Ekf = tEk();

    label i;
    forAll(Ekf, i)
    {
        if (k[i] < 2 || k[i] > 10)
        {
            Ekf[i] = 0.0;
        }
    }
    */

    return tEk;
}

[Will_Wang]

@东岳 他这个是Von-Karman风速谱傅里叶变换的经验公式，李老师我想问一下这个模型算转捩的话扰动应该怎么给会好一些呢

[李东岳]

不太了解这个Von-Karman风速谱傅里叶变换公式。你哪有参考文献么？