求不稳定流动的稳态解

ice_flow

各位老师好，

文献中提到两种求解不稳定流动稳态的方法，一个是直接求解稳态方程；另一个是使用瞬态求解器求解，但是增大时间步长达到滤波的效果，也能达到求解稳态的目的。

我在计算不稳定流动的稳态解时尝试使用直接求解稳态的方案，可以较好的计算出对称结构的稳态解（例如simple求解Re=60的圆柱绕流）。但是使用瞬态求解时，pimple方法仍然受限于CFL数，很难增大时间步长。
我的网格最小尺寸大概在0.01，入口速度1，CFL限制下的deltaT大约在0.01

实际操作下确实最大时间步长不能超过0.008，过大就马上发散了。
但是pimple方法不是可以有很大的时间步长和库朗数吗，是不是我设置有问题，controlDict贴在下面。deltaT 0.05;是尝试达到的时间步长，但是是发散的。

application     pimpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         150;

deltaT          0.05;

writeControl    timeStep;

writeInterval   1000;

purgeWrite      50;

writeFormat     binary;

writePrecision  8;

writeCompression off;

timeFormat      general;

timePrecision   8;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           20;

还有pimple的设置

PIMPLE
{
    correctPhi          yes;
    nOuterCorrectors    50;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
	pRefCell 1001; 
    pRefValue 0; 
}
relaxationFactors
{
    fields
    {
    }
    equations
    {
        "U.*"           1;
    }
}

李东岳

我的网格最小尺寸大概在0.01，入口速度1，CFL限制下的deltaT大约在0.01

真肯定是不正常的。对于这种简单的算例，pimplefoam算到CFL在2-3都可以。

文献中提到两种求解不稳定流动稳态的方法，一个是直接求解稳态方程；另一个是使用瞬态求解器求解，但是增大时间步长达到滤波的效果，也能达到求解稳态的目的。

瞬态求解器为什么能获得稳态的结果？把sci发来我看一下

ice_flow

@李东岳
谢谢老师，那是我的计算设置有问题吗。
文章链接贴在下面
http://dx.doi.org/10.4236/ojfd.2015.52021

李东岳

他们这个文章需要重新写代码来植入。你不能直接用pimple来采用大CFL来进行计算。不过你的CFL=0.01肯定是有问题的。估计哪里没设置对。

ice_flow

@李东岳
基本理论是通过增加时间步长来尝试滤波，大致的时间步长对应的CFL在5左右，所以我想先用pimple试一下。
这里不是CFL=0.01，是实际应用时CFL只能在1以下，再提高CFL就会发散了。这个0.01是CFL=1对应的时间步长。另外想请问老师设置较大的CFL有什么要求吗，我是直接在controlDict里设置了maxCo。

李东岳

另外想请问老师设置较大的CFL有什么要求吗，我是直接在controlDict里设置了maxCo。

没有，自带的cavity，elbow算例之类的，你可以增加dt，co到5估计都没问题

CFD中文网

求不稳定流动的稳态解