求助：【往密闭容器中充气】

td987

往密闭容器里充气，密闭容器的初始压强为一个大气压，本来想让它增加到1.2个大气压，但是结果才增加了1点多pa，我的初始条件是：density-based，unsteady，质量流量入口，ke湍流模型，idea-gas
容器的体积为0.1立方米，里面的初始压强为1个大气压，然后以流量：0.005立方米/秒往里充气4秒，本来容器内部要强变化应改为0.2个大气压左右，但是仿真完才变化增加1点多pa。
问题：是不是在设置质量流量入口时，Supersonic/Initial Gauge Pressure也要设置啊？

td987

知道的人请指导下啊~~~

通流

你应该提供更详细的内容。如果进口条件是什么？

xrs333

质量流量入口，为什么用体积流量单位？

td987

已知条件，已知一容器，如图 ，左端为进气口，进口直径为150mm，整个容器可容纳体积为0.9394 m3，入口质量流量为一周期T为30s的正弦函数变化的规律向容器中充气。质量流量用udf编写如下：
/**********************************************************************
   unsteady.c                                                         
   UDF for specifying a transient velocity profile boundary condition
***********************************************************************/

#include "udf.h"

DEFINE_PROFILE(unsteady_velocity, thread, position)
{
  face_t f;
  real t = CURRENT_TIME;

  begin_f_loop(f, thread)
    {   
      F_PROFILE(f, thread, position) =0.0241*sin(0.209*t);
    }
  end_f_loop(f, thread)
}
其中0.0241=ρAs，ρ为气体密度设为1.225kg/m3，A为正弦函数的幅值，经计算为1.11，s=πr2=0.0176，0.209=2π/T。
本想给容器充气15s，使容器里增入1/5的空气，但由于容器容量是固定的，所以会使压强增大为1.2个大气压，但是结果容器内部压强才增加1点多pa。
问题：1、我选的空气密度是idea-gas，在计算质量流量的时候我用的密度是1.225kg/m3这个有点问题吧，但我不知道怎么用。
2、容器内部刚开始时为一个大气压要初始化吗？要的话怎么初始化？我选的操作压强为默认的一个大气压，这个有问题吗？
3、质量流入口出的压强我本想设置为大气压，那Supersonic/Initial Gauge Pressure还要设置吗，默认为0不行吗？
4、在残差迭代计算过程中，我选择容器内部pressure/static pressure做监视，看到压强增加了上千pa，但是在读云图的时候看pressure/static pressure，压力增加只有1点多pa。
以上的问题可能我问的就有问题，请见谅~~~

[ 本帖最后由 td987 于 2011-10-8 15:29 编辑 ]

td987

我给整错了，不过我有一个问题，我选的空气的密度是idea-gas，那在计算质量流量时候的密度要去多少啊？

td987

图片如下

td987

5、和Courant Numbe 的取值有关，那这个只要怎么取呢？
6、和时间步长有关，我之间步长取的时0.01s，时间步数时1500次。是不是时间步长太小了？

xrs333

用了什么 进口边界条件 ?

td987

密度基，unsteady，质量流量进口，周期T为30s的正弦函数变化，函数为0.0241*sin(0.209*t)，Supersonic/Initial Gauge Pressure=0，水力直径=0.15

xrs333

这个正弦函数在一个周期内的平均值是零。

xrs333

mass flow specification method用的那个选项？

td987

我算的是半个周期。我现在改成压力基的情况好点了 不过还没算好

td987

我用的是Mass Flow Rate，多谢版主的关注~~谢谢

td987

压力基主要用于低速不可压的，密度基主要用于高速可压的，而我的这种情况入口是低速的但是到容器里又是可压的，像我的这种问题用哪一种方式好点啊？