• 한국추진공학회지
• /
• 제12권3호
• /
• pp.68-75
• /
• 2008

이전에 개발된 병렬 해석 코드를 이용하여 원형 관내에서의 삼차원 데토네이션 파 전파 특성을 관찰 하였다. 일련의 해석을 통하여 직경이 일정한 원형 관에서 반응 속도 상수 값 k에 의존하는 데토네이션 셀 생성 메커니즘을 알 수 있었다. 삼차원 유동에 대한 비정상 유동의 해석 결과는 two-, three-, four-cell 구조의 메커니즘을 잘 보여주었으며, 반응속도 상수 k가 증가함에 따라 생성되는 셀 수가 증가하고 셀의 폭과 길이는 감소하였다. 모든 다중-셀 모드에서 데토네이션 파의 구조와 그을음막 기록은 횡단파의 움직임에 의하여 만들어지며, 데토네이션 파면은 주기적으로 일정한 다각형 및 풍차 형상을 가지게 된다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2006년도 제32회 KOSCO SYMPOSIUM 논문집
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제24회 춘계학술대회논문집
• /
• pp.431-434
• /
• 2005

데토네이션 파의 셀 구조를 모사하기 위한 격자 해상도 요구 조건을 파악하기 위하여 여러 영역의 데토네이션 현상에 대하여 반응 속도 상수와 격자 해상도를 바꾸어가며 일련의 해석을 수행 하였다. 서로 다른 영역의 데토네이션 현상에 대하여 정상 상태 ZND 구조와 계산 결과를 비교하여 살펴봄으로써, 최소 격자 크기와 계산 영역의 크기에 대한 결론을 도출하였다.

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

A comprehensive numerical study was carried out to identify the on-set condition of the cell structures of oblique detonation waves (ODWs). Mach 7 incoming flow was considered with all other flow variables were fixed except the flow turning angles varying from 35 to 38. For a given flow conditions theoretical maximum turning angle is $38.2^{\circ}$ where the oblique detonation wave may be stabilized. The effects of grid resolution were tested using grids from $255{\times}100$ to $4,005{\times}1,600$. The numerical smoked foil records exhibits the detonation cell structures with dual triple points running opposite directions for the 36 to 38 turning angles. As the turning angle get closer to the maximum angle the cell structures gets finer and the oscillatory behavior of the primary triple point was observed. The thermal occlusion behind the oblique detonation wave was observed for the $38^{\circ}$ turning angle.

• Nuclear Engineering and Technology
• /
• 제49권7호
• /
• pp.1423-1430
• /
• 2017

Detonation cell width is an important parameter in hydrogen explosion assessments. The experimental data on gas detonation are statistically analyzed to establish a universal method to numerically predict detonation cell widths. It is commonly understood that detonation cell width, ${\lambda}$, is highly correlated with the characteristic reaction zone width, ${\delta}$. Classical parametric regression methods were widely applied in earlier research to build an explicit semiempirical correlation for the ratio of ${\lambda}/{\delta}$. The obtained correlations formulate the dependency of the ratio ${\lambda}/{\delta}$ on a dimensionless effective chemical activation energy and a dimensionless temperature of the gas mixture. In this paper, support vector regression (SVR), which is based on nonparametric machine learning, is applied to achieve functions with better fitness to experimental data and more accurate predictions. Furthermore, a third parameter, dimensionless pressure, is considered as an additional independent variable. It is found that three-parameter SVR can significantly improve the performance of the fitting function. Meanwhile, SVR also provides better adaptability and the model functions can be easily renewed when experimental database is updated or new regression parameters are considered.

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

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

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

In order to investigate the initiation and propagation processes of a spherical detonation wave induced by direct initiation, numerical simulations were carried out using two-dimensional compressible Euler equations with an axisymmetric assumption and a one-step reaction model based on Arrhenius kinetics with various levels of ignition energy. By varying the amount of ignition energy, three typical initiation behaviors, which were subcritical, supercritical and critical regimes, were observed. Then, the ignition energy of more than $137.5{\times}10^6$ in non-dimensional value was required for initiating a spherical detonation wave, and the minimum ignition energy(i.e., critical energy) was less than that of the one-dimensional simulation reported by a previous numerical work. When the ignition energy was less than the critical energy, the blast wave generated from an ignition source continued to attenuate due to the separation of the blast wave and a reaction front. Therefore, detonation was not initiated in the subcrtical regime. When the ignition energy was more than the minimum initiation energy, the blast wave developed into a multiheaded detonation wave propagating spherically at CJ velocity, and then a cellular pattern radiated regularly out from the ignition center in the supercritical regime. The influence on ignition energy was observed in the cell width near the ignition center, but the cell width on the fully developed detonation remained constant during the expanding of detonation wave due to the consecutive formation of new triple points, regardless of ignition energy. When the ignition energy was equal to the critical energy, the decoupling of the blast wave and a reaction front appeared, as occurred in the subcrtical regime. After that, the detonation bubble induced by the local explosion behind the blast wave expanded and developed into the multiheaded detonation wave in the critical regime. Although few triple points were observed in the vicinity of the ignition core, the regularly located cellular pattern was generated after the onset of the multiheaded detonation. Then, the average cell width on the fully developed detonation was almost to that in the supercritical regime. These numerical results qualitatively agreed with previous experimental works regarding the initiation and propagation processes.

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

This paper presents a numerical investigation on kerosene-air mixture detonation and behaviors of thermal elasto-plstic thin metal tube under detonation loading based on multi-material analysis. The detonation loading is modeled by the kerosene-air mixture detonation which is compared with CJ condition and experimental cell size. And the thermal softening effect on elasto-plstic model of metal tube is indicated by different dynamic response of detonation loaded tube in various temperature and tube thickness.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2005년도 제24회 춘계학술대회논문집
• /
• pp.343-346
• /
• 2005

전산 해석과 병렬처리를 이용하여 정사각 관 내부를 전파하는 데토네이션파의 삼차원 구조를 살펴보았다. 가변 비열비 공식과 간략화 된 일단계 Arrhenius 반응 모델과 연계된 비점성 유체 방정식을 MUSCL 기반 TVD 해법과 4단계 Runge-Kutta 적분 방법을 이용하여 해석하였다. 삼차원에서의 비정상 해석 결과로부터 그을음 막 기록(smoked-foil record)에서 같은 길이와 다른 폭을 가지는 수평 및 대각 방향 불안정에 의한 상세한 파면 구조를 파악할 수 있었다.