In the numerical simulation of hypersonic chemical nonequilibrium flows, species density may become negative at the beginning of calculation, chemical stiffness is the main reason that affects the calculation’s stability and convergence. Based on the basic theory of source term’s linearazition, two linearazition methods were chosen. One is based on the physical rules of chemical reaction flows. The method’s advantage is that can improve the calculation’s stability and can avoid the species’ density to become negative during the calculation. The other method is the time preconditioner method, which has the function of diminishing the stiffness of the source term. The method was proved to be well approximates the partial differential equations. The calculation result shows that this method can effectively avoid negative density’s appearance and accelerate the calculation’s convergence. So the method effectively improves the calculation’s efficiency.
Thomas R A B, Rarll M M. A finite volume method for the calculation of compressible chemically reacting flows . AIAA-85-0331, New York:American Institute of Aeronautics and Astronautics, 1985
Kim S L, Jeung I S, Park Y H, et al. Approximate jacobian methods for efficient calculation of reactive flows . AIAA-2000-3384, Reston:American Institute of Aeronautics and Astronautics, 2000
Kim S S, Kim C, Rho O H. Multigrid algorithm for computing hypersonic chemically reacting flows[J].Journal of Spacecraft and Rockets.2001,38(6):865-874
���ø�,��ά�߳�������ճ������������ֵģ�� .����:������ѧ��ѧϵ,1991 Chen Rangfu. Numerical simulation of invisid steady hypersonic flow around three dimensional body .Beijing:School of Mathematical Sciences, Peking University, 1991(in Chinese)