A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avo...A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.展开更多
With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compr...With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compressible viscous flow around rigid airfoils in arbitrary unsteady motion. The selection of physical time step is not restricted by stability condition any more, and most of the successful acceleration techniques used in steady calculations can be implemented to increase the computation efficiency.展开更多
基金Supported by the National Basic Research Program of China(2012CB825601)the National Natural Science Foundationof China(41031066,41231068,41274192,41074121,41204127)+1 种基金the Knowledge Innovation Program of the ChineseAcademy of Sciences(KZZD-EW-01-4)the Specialized Research Fund for State Key Laboratories
文摘A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.
文摘With the cell vertex finite volume discretization in space and second order backward implicit discretization in time, 2D unsteady Navier Stokes equations are solved by a dual time stepping method to simulate compressible viscous flow around rigid airfoils in arbitrary unsteady motion. The selection of physical time step is not restricted by stability condition any more, and most of the successful acceleration techniques used in steady calculations can be implemented to increase the computation efficiency.
