本文对传统的可压缩流动中的通量分裂方法进行了修正,使其可以用于扰动方程的计算。利用这一方法,对来流马赫数Ma为4.5、雷诺数Re为10^5的可压缩平板边界层中扰动的演化进行了数值模拟。对于给定的基本流,在计算域入口加入T-S波扰动,研究扰动的空间发展演化。对于小幅值的扰动,计算得到的扰动幅值演化和扰动法向分布与线性稳定性理论的结果符合的很好。对于有限幅值的扰动,结果表明,各次谐波无论幅值还是扰动分布都与完整的Navier-Stocks方程计算的结果符合的很好,但由扰动方程得到的平均流修正比N-S方程得到的大一些,这也许是由于在N-S方程计算时,基本流与零阶谐波一起计算而产生的误差所至。
Numerical method used in this paper is a amendment of the influx split method which is always used for compressible flow. It makes the method used in this paper fit for the calculation of perturbation equations. By this method, the evolution ofperturbance in compressible flat plate boundary layer with incoming flow of Ma 4.5 and Re 100000 is numerically simulated. For the given base flow, T-S wave at the entrance to research the evolution ofperturbance spatially is imposed. For the case of small amplitude perturbance, the result agrees well with the result gotten by LST. For the case of finite amplitude perturbance, the result is in consistant with the resuR gored by full N-S equation simulation in amplitude and disturbance distribution for all harmonic waves, but, for base flow amendment which is caused by perturbance, it is a little bit larger than the one gotted by full N-S equation simulation, the difference may due to the fact that we calculate base flow and zeroth harmonic together which would cause some error when do the full N-S equation simulation.