计算出错，提示出现负温度，什么原因

shineway

>>> rho =   -901.23     at (     2,     18,      1)

>>> rho =  -0.16473E+18 at (     6,     22,      1)

>>> rho =  -0.22294E+17 at (     5,     23,      1)

>>> rho =  -0.17311E+19 at (     6,     23,      1)

>>> rho =  -0.23291E+19 at (     7,     23,      1)

>>> Additional negative density messages have been suppressed....


[WARNING]: (getpress) Negative or zero temperature encountered.
            Cell #        1027  has T =   -5146.13984678649     
            Reset to original value:    319.191444763736     

[WARNING]: (getpress) Negative or zero temperature encountered.
            Cell #        1255  has T =   -30.2951141472491     
            Reset to original value:    321.227044000696     

[WARNING]: (getpress) Negative or zero temperature encountered.
            Cell #        1262  has T =   -274.772726642636     
            Reset to original value:    321.255480582217     

[WARNING]: (getpress) Negative or zero temperature encountered.
            Cell #        1293  has T =   -40213.0815959950     
            Reset to original value:    320.923845300018     

[WARNING]: (getpress) Negative or zero temperature encountered.
            Cell #        1300  has T =   -381.028192249829     
            Reset to original value:    320.961144571770     

            (getpress) Additional negative-temperature messages suppressed...

    48     116   2.6676E+00   7.6823E-04   8.6874E+12   1.5051E-08   1.3555E-05   3.1621E+03
    49     116   2.5133E-01   3.5435E-05   1.2977E-04   2.6586E-06   2.0747E-05   3.1623E+03

>>> rho =  -0.81964E+31 at (     1,      1,      1)


[ERROR]: (upflim) Negative density at cell            1

  
--- Program Execution Terminated due to Error.
  

shineway

怎么没人回复，目前来看减小终止CFL number用点用，再算算看

nuaawubin

shineway 发表于 2014-11-19 16:18
怎么没人回复，目前来看减小终止CFL number用点用，再算算看

请稍安勿躁，大家都有自己的工作，不一定能按照你的想法及时回复你。
你给的信息只说明了负温产生，但是fastran产生负温基本上和设置不当有关，所以如果要解决问题，可能需要给出dtf文件（或者说明具体设置）

shineway

BC TYPE：inflow/outflow里面的static pressure是指的当地绝对压力吗？我给的101325pa。
三维机翼扰流，稳态SA模型，Ma=0.7; FDS first order;Final CFL=0.2;Time integration--implicit non-iterative。算了600步,残差还没下降一个数量级。应该修改什么加速收敛呐？

nuaawubin

shineway 发表于 2014-11-20 13:32
BC TYPE：inflow/outflow里面的static pressure是指的当地绝对压力吗？我给的101325pa。
三维机翼扰流，稳 ...

首先边界条件没有设置错误，但静压可能要根据实际情况换算出来，而不是直接给101325Pa。收敛性和收敛快慢可能与网格质量和cfl数有关，请关注一下这两个方面。

shineway

好，我改一下，谢谢指点

tigermeng11

请问楼主，你的负密度问题是如何解决的？我也出现的同样的情况。。。。然后我减小了初始CFL数和终止CFL数，倒是可以算了，
残差下降不明显，我觉得应该继续调整CFL数吧。请问楼主，还有做其他的修改吗？