各位大神，真心求教关于边界条件设置的问题。

minqi_lpps

本人最近在编写关于如下物理模型的CFD程序：控制方程为二维轴对称坐标系下的非定常欧拉方程组，初始条件为温度T=300K，压强P=5×10-4 Pa，两个方向的速度都为0。在程序编写过程中用有限差分法，时间上用双时间步的LU-SGS格式，空间上用Van-leer结合NND格式。边界上的情况如下图所示：

1

qint为流入气体的状态参数，流入气体在r<rlas的边界上进入计算域内。在程序中设置z=0处的边界条件是按如下方法：r<=rlas区域z=-1，-2虚网格上的值为流入气体的参数。r>rlas区域：z=-1，-2虚网格上的其他值与内侧的z=1，2上的值相等，只有z方向的速度(也就是图中的v) 与内侧z=1，2上的相反。然后照常使用数值格式计算和更新z=0上的数据，只不过每一步都做修正，使r>rlas区域z方向的速度为0。

在计算第一个时间步时就出现错误，具体错误为：最z=0的层上rlas+1这个点的内能密度[也就是压强/(比热容比-1)]为负值。请问各位这是否是因为我的边界条件给的不合适造成的？如果是应该怎么给？