NS 编程 单元的压力密度比为负值，求指点

ruby_nuaa

采用结构网格，有限体积法编程求解二维翼型NS方程，调程序的时候出现单元格的pressure/density即压力密度比为负值出错，出错地点为翼型的前缘，不知道错在何处，无从下手，求大神指点

ruby_nuaa

额  哪位大神帮帮忙啊？

sxhqjfl

有时候会因为各种原因导致密度为负，差分格式是重要原因之一（但不是唯一原因），所以有人专门做positivity-preserving scheme。如果密度为负的格点不多，最简单的方法就是将变成负值的密度取绝对值，看看结果怎么样再说。关于在计算中“保正”，是当前CFD的难点之一，我也在研究这方面的。