dnsFoam初始化问题



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



  • @wzchen

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



  • :happy: 谢谢,在想问一下
    0_1461811236659_upload-a8c2479c-7dcc-49a1-bc7f-e32db3586c1d
    Ea k0的具体含义是什么,如何设置它们的值?有关于openfoam初始化方面的资料推荐一下吗:happy:



  • @wzchen

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

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



  • 好的,谢谢:happy:



  • @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;
    }
    


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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2