Numerical-simulation of low supersonic-flow around a sphere

The controlled equations defined in a physical plane are changed into those in a computational plane with coordinate transformations suitable for different Mach number M(infinity). The computational area is limited in the body surface and in the vicinities of detached shock wave and sonic line. Thus the area can be greatly cut down when the shock wave moves away from the body surface as M(infinity) --> 1. Highly accurate, total variation diminishing (TVD) finite-difference schemes are used to calculate the low supersonic flowfield around a sphere. The stand-off distance, location of sonic line, etc. are well comparable with experimental data. The long pending problem concerning a flow passing a sphere at 1.3 greater-than-or-equal-to M(infinity) > 1 has been settled, and some new results on M(infinity) = 1.05 have been presented.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

Domain decomposition hybrid method for numerical simulation of bluff body flows:theoretical model and application

The discrete vortex method is not capable of precisely predicting the bluff body flow separation and the fine structure of flow field in the vicinity of the body surface. In order to make a theoretical improvement over the method and to reduce the difficulty in finite-difference solution of N-S equations at high Reynolds number, in the present paper, we suggest a new numerical simulation model and a theoretical method for domain decomposition hybrid combination of finite-difference method and vortex method. Specifically, the full flow. field is decomposed into two domains. In the region of O(R) near the body surface (R is the characteristic dimension of body), we use the finite-difference method to solve the N-S equations and in the exterior domain, we take the Lagrange-Euler vortex method. The connection and coupling conditions for flow in the two domains are established. The specific numerical scheme of this theoretical model is given. As a preliminary application, some numerical simulations for flows at Re=100 and Re-1000 about a circular cylinder are made, and compared with the finite-difference solution of N-S equations for full flow field and experimental results, and the stability of the solution against the change of the interface between the two domains is examined. The results show that the method of the present paper has the advantage of finite-difference solution for N-S equations in precisely predicting the fine structure of flow field, as well as the advantage of vortex method in efficiently computing the global characteristics of the separated flow. It saves computer time and reduces the amount of computation, as compared with pure N-S equation solution. The present method can be used for numerical simulation of bluff body flow at high Reynolds number and would exhibit even greater merit in that case.

CFD结合降阶模型预测阵风响应

传统的阵风响应主要在频域内进行分析，气动载荷基于线性方法计算，不能考虑黏性和跨声速流动影响．飞机设计需考虑不同频率和不同形状阵风的响应，基于CFD的阵风响应预测由于计算工况太多，工作量巨大．本文发展了一种CFD结合非定常气动力ARMA（autoregressive—moving—averagemodel）降阶模型的阵风响应分析方法，CFD只要针对给定频率和形状的一种阵风响应进行计算，对获得的气动力时间历程运用线性最小二乘法参数辨识ARMA降阶模型的系数，则对任意频率和形状的阵风，代入降阶模型即可确定该阵风的响应，大大提高了计算效率．为效验发展的方法，先计算NACA0012翼型在低马赫数0．11的阵风响应，通过对比CFD、ARMA降阶模型及早期发展的不可压阵风响应预测方法的结果，验证了方法的有效性．再对比CFD、ARMA在跨声速马赫数0．8的阵风响应预测结果，证实所发展的方法对跨声速阵风响应预测亦是有效的．