An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptiv...An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptive time-step method is presented,which can improve the efficiency greatly.In the convection step,a third-order accurate weighted essentially non-oscillatory(WENO)method is adopted to reconstruct the solution in the unstructured grids.Numerical results show that our new method can capture the correct propagation speed of the detonation wave exactly even in coarse grids,while high order accuracy can be achieved in the smooth region.In addition,the proposed adaptive splitting method can reduce the computational cost greatly compared with the traditional splitting method.展开更多
Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is deve...Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.展开更多
Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the l...Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering (PLMS) technique. Two methods are developed to investigate the spatial distribu- tion and geometry of IWS: (1) the dominant mode extrac- tion (DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and (2) the phase tracking (PT) method, to identify the shape of IWS. The re- sults suggest that pressure perturbations account for the for- marion of IWS in the initial mixing region and the joint effect of dilatation and coherent vortices enhances IWS in the well- developed region. The large transverse (cross-flow) scale of the IWS and their relation to coherent vortices (CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the effective tools to study the geometrical features of wavy structures in compressible shear flows.展开更多
基金supported by the National Natural Science Foundation of China(Grants 51476152,11302213,and 11572336)
文摘An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptive time-step method is presented,which can improve the efficiency greatly.In the convection step,a third-order accurate weighted essentially non-oscillatory(WENO)method is adopted to reconstruct the solution in the unstructured grids.Numerical results show that our new method can capture the correct propagation speed of the detonation wave exactly even in coarse grids,while high order accuracy can be achieved in the smooth region.In addition,the proposed adaptive splitting method can reduce the computational cost greatly compared with the traditional splitting method.
基金supported by the National Natural Science Foundation of China (90716008,10572004 and 10921202)MOST 973 Project (2009CB724100) and CSSA
文摘Discontinuous Galerkin(DG) method is known to have several advantages for flow simulations,in particular,in fiexible accuracy management and adaptability to mesh refinement. In the present work,the DG method is developed for numerical simulations of both temporally and spatially developing mixing layers. For the temporally developing mixing layer,both the instantaneous fiow field and time evolution of momentum thickness agree very well with the previous results. Shocklets are observed at higher convective Mach numbers and the vortex paring manner is changed for high compressibility. For the spatially developing mixing layer,large-scale coherent structures and self-similar behavior for mean profiles are investigated. The instantaneous fiow field for a three-dimensional compressible mixing layer is also reported,which shows the development of largescale coherent structures in the streamwise direction. All numerical results suggest that the DG method is effective in performing accurate numerical simulations for compressible shear fiows.
基金supported by National Nature Science Foundation of China(90716008,10572004,and 11172006)by MOST 973 Project(2009CB724100)
文摘Semi-periodic structures namely inclined wavy structures (IWS) are experimentally observed in compressible mixing layers at two convective Mach numbers (Mc = 0.11 and 0.47). Flow structures are visualized by the laserinduced planar laser Mie scattering (PLMS) technique. Two methods are developed to investigate the spatial distribu- tion and geometry of IWS: (1) the dominant mode extrac- tion (DME) method, to extract the dominant modes of IWS from the streamwise gray-level fluctuation, and (2) the phase tracking (PT) method, to identify the shape of IWS. The re- sults suggest that pressure perturbations account for the for- marion of IWS in the initial mixing region and the joint effect of dilatation and coherent vortices enhances IWS in the well- developed region. The large transverse (cross-flow) scale of the IWS and their relation to coherent vortices (CV) indicate that the disturbance originated from CV in the mixing center propagates far into the free streams. The DME and the PT method are shown to be the effective tools to study the geometrical features of wavy structures in compressible shear flows.