CFD中文网

@逆变张量 是的速度入口不合适，给的总压

@ibelief 解决了呢，没问题。

自带速度型入口算例中，经过编译通过了，在速度入口边界条件中选择，计算发现仅在第一层网格内执行了，第二层网格已经是来流速度了，单位、网格均没有问题。
udf：

DEFINE_PROFILE(x_velocity,t,i)
 {
    real y, del, h, x[ND_ND], ufree;  /* variable declarations */
    face_t f;

    h = YMAX - YMIN;
    del = DELOVRH*h;
    ufree = UMEAN*(B+1.);

    begin_f_loop(f,t)
      {
         F_CENTROID(x,f,t);
         y = x[1];
         if (y <= del)
           F_PROFILE(f,t,i) = ufree*pow(y/del,B);
         else
           F_PROFILE(f,t,i) = ufree*pow((h-y)/del,B);
      }
    end_f_loop(f,t)

@zxyf1 flux 在哪里改呢请教

@static还没呢，已经对这个求解器失去信心

什么求解器用的，sonicfoam可以用wavetransmissive，一般的直接给速度入口压力固定值出口。

@d-benjamin 好的。

如题，请问是具体在哪里设置，怎么使用？

附上温度限制：

temperatureLimit
{
    type limitTemperature;
    active true;

    limitTemperatureCoeffs
    {
        selectionMode all;
        min <Tmin>;
        max <Tmax>;
    }
}

div改成upwind有所改善，但不能根本上解决问题。
如果温度使用uniformFixedValue将压力从小到大随时间增大，为何实际算的时候也不起作用？
初始压力极小如绝对压力为1，即使温度也很低，则能正常计算，压力和温度同事增大或减小才可以计算，这是什么原因？

问题描述：
以openfoam自带算例斜劈（或其他算例均可）为例，修改如下：
M=6或10等，热物理参数使用pefectgas, 入口温度为60K，静压1270pa，模拟试验状态，因为总温有限，用等熵公式算出静温一般较低；边界条件不变
问题是，不管Dt和CFL数怎么调，很容易在计算若干步出现初始T0为负值的情况；
个人理解：M较大时，流场参数的梯度较大，过低的初始温度和静压加剧了计算不稳定，若改成自然来流，则能正常计算，通过fvOptions对温度进行限制根本不起作用，请大家指导一二