An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to captur...An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.展开更多
The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed o...The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed on the effect of chemical reactivity of mixture on the flame distortion and detonation initiation after the passage of the reflected shock wave. It is shown that the spatio-temporal characteristics of detonation initiation depend highly on the chemi- cal reactivity of the mixture. When the chemical reactivity enhances, the flame can be severely distorted to form a reactive shock bifurcation structure with detonations initiating at different three-dimensional spatial locations. Moreover, the detonation initiation would occur earlier in a mixture of more enhanced reactivity. The results reveal that the detona- tions arise from hot spots in the unburned region which are initiated by the shock-detonation-transition mechanism.展开更多
Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate ...Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.展开更多
Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS...Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.展开更多
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology (KFJJ09-13)
文摘An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.
基金supported by the National Natural Science Foundation of China (10972107)Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology (KFJJ12-4Y)Jiangsu Innovation Program for Graduate Education (CXLX11 0271)
文摘The interactions of a spherical flame with an incident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed on the effect of chemical reactivity of mixture on the flame distortion and detonation initiation after the passage of the reflected shock wave. It is shown that the spatio-temporal characteristics of detonation initiation depend highly on the chemi- cal reactivity of the mixture. When the chemical reactivity enhances, the flame can be severely distorted to form a reactive shock bifurcation structure with detonations initiating at different three-dimensional spatial locations. Moreover, the detonation initiation would occur earlier in a mixture of more enhanced reactivity. The results reveal that the detona- tions arise from hot spots in the unburned region which are initiated by the shock-detonation-transition mechanism.
基金supported by the National Natural Science Foundation of China (10872096)the Open Fund of State Key Laboratory of Explosion Science and Technology, Beijing University of Science and Technology (KFJJ09-13)
文摘Gaseous detonation propagating in a toroidal chamber was numerically studied for hydrogen/oxygen/nitrogen mixtures. The numerical method used is based on the three-dimensional Euler equations with detailed finiterate chemistry. The results show that the calculated streak picture is in qualitative agreement with the picture recorded by a high speed streak camera from published literature. The three-dimensional flow field induced by a continuously rotating detonation was visualized and distinctive features of the rotating detonations were clearly depicted. Owing to the unconfined character of detonation wavelet, a deficit of detonation parameters was observed. Due to the effects of wall geometries, the strength of the outside detonation front is stronger than that of the inside portion. The detonation thus propagates with a constant circular velocity. Numerical simulation also shows three-dimensional rotating detonation structures, which display specific feature of the detonation- shock combined wave. Discrete burning gas pockets are formed due to instability of the discontinuity. It is believed that the present study could give an insight into the interest- ing properties of the continuously rotating detonation, and is thus beneficial to the design of continuous detonation propulsion systems.
基金Project supported by the National Natural Science Foundation of China(Grant No.11672135)a Foundation for the Author of National Excellent Doctoral Dissertation of China(Grant No.201461)
文摘Based on the Fourier–Chebyshev spectral method, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with direct numerical simulation(DNS) methods for different Reynolds numbers. A formula is derived to express the relation between fluctuating velocities and the friction drag coefficient. With the application of electromagnetic force, the in-depth relations among the fluctuating velocities near the wall, Reynolds stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The fluctuating velocities along the streamwise and normal directions are suppressed significantly,while the fluctuating velocity along the spanwise direction is enhanced dramatically due to the spanwise electromagnetic force. However, the values of Reynolds stress depend on the fluctuating velocities along the streamwise and normal directions rather than that along the spanwise direction. Therefore, the significant effect of drag reduction is obtained. Moreover,the maximum drag reduction is weakened due to the decay of control effect for fluctuating velocities as the Reynolds number increases.