• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2012년도 제45회 KOSCO SYMPOSIUM 초록집
• /
• pp.359-362
• /
• 2012

The results are shown highly unstable mode's detonation dynamics by compared with weakly unstable mode. And we investigate the difference and similar features of 2D and 3D results. By using PSD via FFT, the effects of pre-exponential factor difference and of unstable mode were investigated in this study. The result of PSD is shown pulsed features in weakly unstable mode, but noselike in highly unstable mode. By compared between Sheliren image and overlaid slice image, the irregular feature of detonation waves structure was discussed in highly unstable mode.

• 한국추진공학회지
• /
• 제15권4호
• /
• pp.57-64
• /
• 2011

펄스 데토네이션 엔진에서와 같이 탄화수소를 연료로 하는 데토네이션 파는 강한 불안정성을 가지며 난류 연소 효과를 고려한 연구를 수행하여야 함이 제시된 바 있다. 본 연구에서는 강한 불안정성을 가지는 데토네이션 파의 구조를 이해하기 위하여 비점성 해석, 점성 해석, 난류 모델 및 간단한 난류 연소 모델을 고려한 수치 해석 연구를 수행하였다. 모델링 수준에 따른 연구를 통하여 점성 및 난류는 저주파 특성에는 거의 영향이 없으나, 고주파 특성을 강화하는 경향이 있는 것으로 보인다. 한편, 데토네이션 연구를 위한 난류-연소 상호 작용 모델에는 활성화 에너지의 영향이 고려되어야 하는 것으로 여겨진다.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2007년도 제29회 추계학술대회논문집
• /
• pp.367-370
• /
• 2007

Present study examines detonation wave propagation characteristics in annular channel. A normalized value of channel width to the annular radius was considered as a geometric parameter. A parametric study was carried out for a various regimes of detonation waves from weakly unstable to highly unstable detonation waves. Numerical approaches that used in the previous study of numerical requirements of the simulation of detonation wave propagations in 2D and 3D channel were used also for the present study with OpenMP parallization for multi-core SMP machines. The major effect of the curved geometry on the detonation wave propagation seems to be a flow compression effect, regardless of the detonation regimes. The flow compression behind the detonation wave by the curved geometry of the circular channel pushes the detonation wave front and results in the overdriven detonation waves with increased detonation speed beyond the Chapmann-Jouguet speed. This effect gets stronger as the normalized radius smaller, as expected. The effect seems to be negligible beyond the normalized radius of 10.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2006년도 제26회 춘계학술대회논문집
• /
• pp.314-318
• /
• 2006

연소파를 이용한 로켓 점화장치 등에의 응용으로 관내에서 전파되는 데토네이션에 대한 관심이 증대되고 있으나 그 특성에 관하여 알려진 바는 많지 않다. 본 실험적 연구에서는 프로판-산소 혼합기로 채워진 직경 5mm, 길이 30m의 투명한 관을 이용하여 데토네이션 파의 진행을 고속촬영을 통하여 가시화하였다. 평균 속도를 측정한 결과, 압력이 감소함에 따라 정상적인 Chapman-Jouguet(CJ) 모드에서 약 $0.5V_{CJ}$에 이르는 저속 모드로 평균 속도가 감소하는 천이 영역이 존재한다. 관 길이 전체에 걸쳐 데토네이션 전파의 동적 특성을 관찰한 결과, 천이 영역에서는 데토네이션이 매우 불안정하여 주기적이거나 간헐적인 속도의 출렁임이 발생한다.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2014년도 제49회 KOSCO SYMPOSIUM 초록집
• /
• pp.169-174
• /
• 2014

Detonation propagation studies is recently getting more attention in these days for its feasibility in aerospace application. Another motivation for this study is the safety concern in industries, since the detonation can cause failure to the mechanical components particularly when the flame accelerates within a pipe or tubes. In this study we numerically simulated a Moderately unstable detonation case with various grid systems and fluid dynamic length scales and have compared in the contents. Moderately Unstable detonation case was selected for this study and detailed Hydrogen-Air Reaction Mechanisms proposed by Jachimowski was used in this study with N2 as inert species.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 춘계 학술대회논문집
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• 한국추진공학회지
• /
• 제10권2호
• /
• pp.1-14
• /
• 2006

본 논문은 약한 불안정 데토네이션 영역부터 강한 불안정 데토네이션 영역까지 여러 영역에 걸친 데토네이션 파 셀 구조 모사에 대한 수치적 문제점들을 살펴보았다. 비열 비 값이 변하는 점성 유체 역학 방정식 및 1단계 Arrhenius 반응 모델 해석을 위하여 MUSCL-type TVD 기법을 이용한 공간 차분과 4차 정확도의 Runge-Kutta 시간 적분을 이용하였다. 일련의 수치해석 연구는 여러 반응 상수 및 격자 해상도에 따른 데토네이션 셀 구조를 해석하기 위하여 요구되는 계산 조건을 구하기 위하여 다양한 데토네이션 현상 영역에서 수행되었다. 다른 영역의 데토네이션 현상에서 셀 구조를 포착하기 위한 계산 영역의 크기와 최소 격자 해상도를 찾아내기 위하여 정상 1차원 ZND 해석 결과와 전산 해석 결과를 비교 검토하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.341-345
• /
• 2006

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.364-370
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable and unstable pitch modes for the lower and higher activation energies, respectively. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of two modes. The maximum pressure history in the stable pitch remained nearly constant, and the single Mach leg existing on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the unstable pitch due to the generation and decay of complex Mach interaction on the shock front shape. The high frequency oscillation was self-induced because the intensity of the transverse wave was changed during propagation in one cycle. The high frequency behavior was not always the same for each cycle, and therefore the low frequency oscillation was also induced in the pressure history.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.