In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the...In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.展开更多
A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated b...A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated by utilizing an eigenvalue analysis. Two test problems of wave propagation with initial disturbance consisting of a Gaussian profile or rectangular pulse are performed. And the performance of the schemes in short,intermediate,and long waves is evaluated. Moreover,numerical results between the nodal discontinuous Galerkin method and finite difference type schemes are compared,which indicate that the numerical solution obtained using nodal discontinuous Galerkin method with a pure central flux has obviously high frequency oscillations for initial disturbance consisting of a rectangular pulse,which is the same as those obtained using finite difference type schemes without artificial selective damping. When an upwind flux is adopted,spurious waves are eliminated effectively except for the location of discontinuities. When a limiter is used,the spurious short waves are almost completely removed. Therefore,the quality of the computed solution has improved.展开更多
A second-order optimized monotonicity-preserving MUSCL scheme(OMUSCL2) is developed based on the dispersion and dissipation optimization and monotonicity-preserving technique.The new scheme(OMUSCL2) is simple in expre...A second-order optimized monotonicity-preserving MUSCL scheme(OMUSCL2) is developed based on the dispersion and dissipation optimization and monotonicity-preserving technique.The new scheme(OMUSCL2) is simple in expression and is easy for use in CFD codes.Compared with the original second-order or third-order MUSCL scheme,the new scheme shows nearly the same CPU cost and higher resolution to shockwaves and small-scale waves.This new scheme has been tested through a set of one-dimensional and two-dimensional tests,including the Shu-Osher problem,the Sod problem,the Lax problem,the two-dimensional double Mach reflection and the RAE2822 transonic airfoil test.All numerical tests show that,compared with the original MUSCL schemes,the new scheme causes fewer dispersion and dissipation errors and produces higher resolution.展开更多
基金NKBRSF CG 1990 3 2 80 5 National Natural Science F oundation of China !( No.5 98760 0 2 )
文摘In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.
基金Supported by the National Natural Science Foundation of China(51106099,50976072)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)
文摘A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated by utilizing an eigenvalue analysis. Two test problems of wave propagation with initial disturbance consisting of a Gaussian profile or rectangular pulse are performed. And the performance of the schemes in short,intermediate,and long waves is evaluated. Moreover,numerical results between the nodal discontinuous Galerkin method and finite difference type schemes are compared,which indicate that the numerical solution obtained using nodal discontinuous Galerkin method with a pure central flux has obviously high frequency oscillations for initial disturbance consisting of a rectangular pulse,which is the same as those obtained using finite difference type schemes without artificial selective damping. When an upwind flux is adopted,spurious waves are eliminated effectively except for the location of discontinuities. When a limiter is used,the spurious short waves are almost completely removed. Therefore,the quality of the computed solution has improved.
基金supported by the National Natural Science Foundation of China (Grant Nos.10632050,10872205,11072248)the National Basic Research Program of China (Grant No.2009CB724100)+1 种基金the National High Technology Research and Development Program of China (Grant No.2009AA010A139)the Chinese Academy Sciences Program (Grant No.KJCX 2-EW-J01)
文摘A second-order optimized monotonicity-preserving MUSCL scheme(OMUSCL2) is developed based on the dispersion and dissipation optimization and monotonicity-preserving technique.The new scheme(OMUSCL2) is simple in expression and is easy for use in CFD codes.Compared with the original second-order or third-order MUSCL scheme,the new scheme shows nearly the same CPU cost and higher resolution to shockwaves and small-scale waves.This new scheme has been tested through a set of one-dimensional and two-dimensional tests,including the Shu-Osher problem,the Sod problem,the Lax problem,the two-dimensional double Mach reflection and the RAE2822 transonic airfoil test.All numerical tests show that,compared with the original MUSCL schemes,the new scheme causes fewer dispersion and dissipation errors and produces higher resolution.
