气体的守恒变量平滑

对于气体，守恒变量ρ、ρV和ρεt基于在先前流体增量中计算的局部压力梯度进行平滑（在presgradient.f子程序中）。平滑在smoothshock.f例程中进行。平滑过程中的主要参数是shockcoef，其取值范围为0（无平滑）到2。该参数可以在输入文件的*STEP卡片上使用参数SHOCK SMOOTHING设置。如果求解发散（这通常表现为解中存在NaN（非数值），导致负时间步长或守恒变量到物理变量转换例程con2phys.f的发散），则计算以两倍于先前值的冲击平滑系数重新启动（如果先前的值为0，则为0.001）。一旦超过值2，程序将停止并显示错误消息。