A third-order numerical scheme is presented to give approximate solutions to multi-dimensional hyperbolic conservation laws only using modified coefficients of an essentially non-oscillatory (MCENO) scheme without i...A third-order numerical scheme is presented to give approximate solutions to multi-dimensional hyperbolic conservation laws only using modified coefficients of an essentially non-oscillatory (MCENO) scheme without increasing the base points during construction of the scheme. The construction process shows that the modified coefficient approach preserves favourable properties inherent in the original essentially nonoscillatory (ENO) scheme for its essential non-oscillation, total variation bounded (TVB), etc. The new scheme improves accuracy by one order compared to the original one. The proposed MCENO scheme is applied to simulate two-dimensional Rayleigh-Taylor (RT) instability with densities 1:3 and 1:100, and solve the Lax shock-wave tube numerically. The ratio of CPU time used to implement MCENO, the .third-order ENO and fifth-order weighed ENO (WENO) schemes is 0.62:1:2.19. This indicates that MCENO improves accuracy in smooth regions and has higher accuracy and better efficiency compared to the original ENO scheme.展开更多
A dual-time method is introduced to calculate the unsteady flow in a certain vibrating flat cascade. An implicit lower-upper symmetric-gauss-seidel scheme(LU-SGS) is applied for time stepping in pseudo time domains,...A dual-time method is introduced to calculate the unsteady flow in a certain vibrating flat cascade. An implicit lower-upper symmetric-gauss-seidel scheme(LU-SGS) is applied for time stepping in pseudo time domains, and the convection items are discretized with the spatial three-order weighted non-oscillatory and non-free-parameter dissipation difference (WNND) scheme. The turbulence model adopts q-co low-Reynolds-number model. The frequency specmuns of lift coefficients and the unsteady pressure-difference coefficients at different spanwise heights as well as the entropy contours at blade tips on different vibrating instants, are obtained. By the analysis of frequency specmuns of lift coefficients at three spanwise heights, it is considered that there exist obvious non-linear perturbations in the flow induced by the vibrating, and the perturbation frequencies are higher than the basic frequency. The entropy contours at blade tips at different times display an intensively unsteady attribute of the flow under large amplitudes.展开更多
基金the National Natural Science Foundation of China(Nos.10676031,50675185)the Ph.D. Programs Foundation of Ministry of Education of China(No.20070530003)+1 种基金the Key Project of Chinese Ministry of Education(No.208093)the Scientific Research Fund of Hunan Provincial Education Department(No.07A068)
文摘A third-order numerical scheme is presented to give approximate solutions to multi-dimensional hyperbolic conservation laws only using modified coefficients of an essentially non-oscillatory (MCENO) scheme without increasing the base points during construction of the scheme. The construction process shows that the modified coefficient approach preserves favourable properties inherent in the original essentially nonoscillatory (ENO) scheme for its essential non-oscillation, total variation bounded (TVB), etc. The new scheme improves accuracy by one order compared to the original one. The proposed MCENO scheme is applied to simulate two-dimensional Rayleigh-Taylor (RT) instability with densities 1:3 and 1:100, and solve the Lax shock-wave tube numerically. The ratio of CPU time used to implement MCENO, the .third-order ENO and fifth-order weighed ENO (WENO) schemes is 0.62:1:2.19. This indicates that MCENO improves accuracy in smooth regions and has higher accuracy and better efficiency compared to the original ENO scheme.
基金This Project is supported by National Natural Science Foundation of China (No.50776056)National Hi-tech Research and Development Program of China (863 Program,No.2006AA05Z250).
文摘A dual-time method is introduced to calculate the unsteady flow in a certain vibrating flat cascade. An implicit lower-upper symmetric-gauss-seidel scheme(LU-SGS) is applied for time stepping in pseudo time domains, and the convection items are discretized with the spatial three-order weighted non-oscillatory and non-free-parameter dissipation difference (WNND) scheme. The turbulence model adopts q-co low-Reynolds-number model. The frequency specmuns of lift coefficients and the unsteady pressure-difference coefficients at different spanwise heights as well as the entropy contours at blade tips on different vibrating instants, are obtained. By the analysis of frequency specmuns of lift coefficients at three spanwise heights, it is considered that there exist obvious non-linear perturbations in the flow induced by the vibrating, and the perturbation frequencies are higher than the basic frequency. The entropy contours at blade tips at different times display an intensively unsteady attribute of the flow under large amplitudes.