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

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.391-401
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and the results are compared to the existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results show that as cross-sectional area decreases with the reduction in downstream pressure, flow rate increases and reaches the maximum when the speed index (mean velocity divided by wave speed) is near the unity at the point of minimum cross-section area, indicating the flow limitation or choking (flow speed equals wave speed) in one-dimensional studies. for further reductions in downstream pressure, flow rate decreases. The flow limitation or choking consist of the main reasons of waterfall effect which occurs in the airways, capillaries of lung, and other veins. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is 2, the area throat is located near the downstream end. As this ratio is increased to 3, the constriction moves to the upstream end of the tube.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.1
• /
• pp.11-19
• /
• 2015

We present the results of a multi-material numerical investigation of the propagation of a combustible gas mixture detonation in narrow metal tubes. We use an experimentally tuned one step Arrhenius chemical reaction and ideal gas equation of state (EOS) to describe stoichiometric $H_2-O_2$ and $C_2H_4-O_2$ detonations. The purely plastic deformations of copper and steel tubes are modeled using the Mie-Gruneisen EOS and Johnson-Cook strength model. To precisely track the interface motion between the detonating gas and the deforming wall, we use the hybrid particle level-sets within the ghost fluid framework. The calculated results are validated against the experimental data because the results explain the process of the generation and subsequent interaction of the expansion wave with the high-strain-rate deformation of the walls.

• Nuclear Engineering and Technology
• /
• v.56 no.5
• /
• pp.1712-1724
• /
• 2024

Spiral tubes are used in a wide range of applications and it is significant to understand the vibration introduced by two-phase flow in spiral tubes. In this paper, the numerical method is used to study the vibration induced by the gas-liquid two-phase flow in a spiral tube with different flow regimes. The pressure fluctuation characteristics at the pipe wall and the solid vibration response characteristics are obtained. The results show that the motion of small bubbles in bubbly flow leads to small pressure fluctuations with low-frequency broadband (0-50 Hz). The motion of the gas plug in the plug flow causes small amplitude periodic pressure fluctuation with a shortened low-frequency broadband (0-15 Hz) compared to the bubbly flow. The motion of the gas slug in the slug flow causes large periodic fluctuations in pressure with a significant dominant frequency (6-7 Hz). The wavy flow is very stable and has a distinct main frequency (1-2 Hz). The vibration regime in the bubbly flow and wave flow are close to the first-order mode, and the vertical vibrating component is dominant. The plug flow and slug flow excite higher-order vibration modes, and the lateral vibration component plays more important part in the vibration response.

• Journal of Aerospace System Engineering
• /
• v.12 no.2
• /
• pp.15-22
• /
• 2018

A compressible multiphase flow was investigated using a DIM consisting of seven equations, including the fifth-order MLP and a modified HLLC Riemann solver to achieve a precise interface structure of liquid and gas. The numerical methods were verified by comparing the flow structures of the high-pressure water and low-pressure air in the shock tube. A 2D air-helium shock-bubble interaction at the incident shock wave condition (Mach number 1.22) was numerically solved and verified using the experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.3
• /
• pp.315-324
• /
• 2012

This paper presents a numerical investigation of the deflagration-to-detonation transition (DDT) of flame acceleration by a shock wave filled with an ethylene/air mixture as the combustible gas, considering geometrical changes by using obstacles and bent tubes. The model used consists of the reactive compressible Navier-Stokes equations and the ghost fluid method (GFM) for complex boundary treatment. Simulations with a variety of bent tubes with obstacles show the generation of hot spots through flame and strong shock-wave interactions, and restrained or accelerated flame propagation due to geometrical effects. In addition, the simulation results show that the DDT occurs with a nearly constant chemical heat-release rate of 20 MJ/($g{\bullet}s$) in our numerical setup. Furthermore, the DDT triggering time can be delayed by the absence of unreacted material together with insufficient pressures and temperatures induced by different flame shapes, although hot spots are formed in the same positions.