• 한국가시화정보학회지
• /
• 제18권3호
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

• 천문학회지
• /
• 제40권4호
• /
• pp.161-164
• /
• 2007

We calculate the evolution of multiple supernova (SN) explosions inside a pre-exiting bubble blown up by winds from massive stars, using one-dimensional hydrodynamic simulations including radiative cooling and thermal conduction effects. First, the development of the wind bubble driven by collective winds from multiple stars during the main sequence is calculated. Then multiple SN explosion is loaded at the center of the bubble and the evolution of the SN remnant is followed for $10^6$ years. We find the size and mass of the SN-driven shell depend on the structure of the pre-existing wind bubble as well as the total SN explosion energy. Most of the explosion energy is lost via radiative cooling, while about 10% remains as kinetic energy and less than 10% as thermal energy of the expanding bubble shell. Thus the photoionization and heating by diffuse radiation emitted by the shock heated gas is the most dominant form of SN feedback into the surrounding interstellar medium.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2000년도 춘계학술대회논문집
• /
• pp.1185-1191
• /
• 2000

In general, the strength of hull structures can be estimated from stress evaluation considering static and hydro-dynamic load due to sea-wave. However, war ships such as submarine, have frequently experienced the underwater explosion and local structures of ship as well as hull girder can be damaged by the dynamic response excited from underwater non-contact explosion. When explosion happens at underwater, shock wave is radiated In early short time, then gas bubbles are generated, and expansion and contraction are repeated as they float to the surface. The shock wave causes the damage of equipment and its supporting structures, on the other hand, the hull girder strength can be lost by resonance between bubble pulsation and lowest ship natural vibration period. In this paper, the hydro-Impulse force due to bubble was calculated. Based on these results the hull girder strength of submarine was estimated from transient response analysis by using NASTRAN. Also, shock analysis for some equipment supporting structures was carried out by using DDAM. In order to evaluate the strength of these local structures due to shock wave.

• 한국전산유체공학회지
• /
• 제11권1호
• /
• pp.29-35
• /
• 2006

A code is developed to analyze a spherically symmetric underwater explosion. The arbitrary Lagrangian-Eulerian(ALE) Godunov scheme for two-phase flow is used to calculate numerical fluxes through moving control surfaces. For detonation gas of TNT and liquid water, the Jones-Wilkins-Lee(JWL) equation of states and the isentropic Tait relation are used respectively. It is suggested to use the Godunov variable to estimate the velocity of a material interface. The code is validated through comparisons with other results on the gas-water shock tube problem. It is shown that the code can handle generation of discontinuity and recovering of continuity in the normal velocity near the material interface during shock waves interact with the material interface. The developed code is applied to analyze a spherically symmetric underwater explosion. Repeated transmissions of shock waves are clearly captured. The calculated period and maximum radius of detonation gas bubble show good agreements with experimental and other numerical results.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.1201-1206
• /
• 2004

Sonoluminescence (SL) characteristics such as pulse shape, radiance and spectrum radiance from submicron bubbles were investigated. In this study, a set of analytical solutions of the Navier-Stokes equations for the gas inside bubble and equations obtained from mass, momentum and energy equations for the liquid layer adjacent the bubble wall were used to estimate the gas temperature and pressure at the collapse point, which are crucial parameters to determine the SL characteristics. Heat transfer inside the gas bubble as well as at the liquid boundary layer, which was not considered in the most of previous studies on the sonoluminescence was taken it into account in the calculation of the temperature distribution inside the bubble. It was found that bremsstrahlung is a very possible mechanism of the light emission from either micron or submicron bubbles. It was also found that the peak temperature exceeding $10^{6}$ K in the submicron bubble driven at 1 MHz and 4 atm may be due to the rapid change of the bubble wall acceleration near the collapse point rather than shock formation.

• 한국음향학회지
• /
• 제8권1호
• /
• pp.23-30
• /
• 1989

기포 혼합 유체 내에서의 압력파의 전파 현상을 수치해석으로 연구한다. 혼합 유체 영역을 지배하는 지배 방정식을 heuristic 한 방법으로 유도하고 기포 내부 영역에는 열전달 효과를 고려할 수 있도록 에너지 방정식을 도입한다. 기포 내부의 비등온 조건은 특히 기포가 고진폭을 가지고 진동할 때 매우 중요하다. 기포 역학 방정식으로서 Keller 방정식이 채택, 변형되어 기포 외부와 내부의 coupling을 맺어준다. 실제 문제로서 충격관내 충격파의 전파 현상을 수치해석 방법으로 해석한 결과가 Noordzij 및 van Wijngaarden 의 실험 결과와 거의 일치한다. 그러나 그들에 의해 설명된 충격파 구조의 변화 원인은 가스와 액체 간의 상대 운동인데 이는 본 모델에서 고려되지 않았기 때문에 가스와 액체 간의 열전달에 의해 충격파의 구조가 변화된다고 보는 것이 타당하다.

• 화약ㆍ발파
• /
• 제35권3호
• /
• pp.21-30
• /
• 2017

폭약이 수중에서 폭발하면 공기 중에서 폭발하였을 때와는 다르게 폭발 이후 발생하는 가스의 영향에 대한 고려가 필요하다. 수중 폭발 시에는 폭발압력의 전파속도가 공기 중에서 폭발했을 때에 비하여 빠르고, 발생하는 가스 또한 확산되어 에너지를 소실하기 전 물에 의하여 버블의 형태로 갇히게 된다. 이 때 버블은 팽창과 수축을 반복하며 충격파를 만들어낸다. 이러한 일련의 현상을 연구하기 위하여 내부를 관측할 수 있는 실린더형 철재 수조를 제작하고 폭발 실험을 수행하였다. 본 연구에서는 탄체가 없는 소량의 펜톨라이트를 수중에서 폭발시켰고, 이 때 발생하는 충격파를 계측하고 발생된 가스버블의 거동을 관측하여 그 결과를 관찰하였다.

• 대한조선학회논문집
• /
• 제57권5호
• /
• pp.271-277
• /
• 2020

Recoverability and vulnerability of navy ships under underwater explosion are critical verification factors in the acquisition phase of navy ships. This paper aims to establish numerical analysis techniques for the underwater explosion of navy ships. Doubly Asymptotic Approach (DAA) Equation of Motion (EOM) of primary shock wave and secondary bubble pulse proposed by Geers-Hunter was introduced. Assuming a non-compressive fluid, reference solution of the DAA EOM of Geers-Hunter using Runge-Kutta method was derived for the secondary bubble pulse phase with an assumed charge conditions. Convergence analyses to determine fluid element size were performed, suggesting that the minimum fluid element size for underwater explosion analysis was 0.1 m. The spherical and cylindrical fluid domains were found to be appropriate for the underwater explosion analyses from the fluid domain shape study. Because the element size of 0.1 m was too small to be applied to the actual navy ships, a very slender beam with the square solid section was selected for the study of fluid domain existence effect. The two underwater explosion models with/without fluid domain provided very similar results in terms of the displacement and stress processes.

• 대한조선학회논문집
• /
• 제42권6호
• /
• pp.631-636
• /
• 2005

One and three dimensional whipping response analyses of a naval surface combatant subjected to an underwater explosion bubble pulse were carried out to compare the efficiency and accuracy according to the modeling methods. In 1-D analysis, program UNDEXWHIP developed by KIMM was used, which is based on the thin-walled Timoshenko's beam theory and on the modal analysis method using wetted vibratory modes of the hull girder. In 3-D analysis, three finite element models were suggested using LS-DYNA/USA code, such as 3-D beam model considering geometric shape of wetted side shell, coarse and fine 3-D F.E. models. Through the comparison of results from the 1-D and 3-D analyses, it could be confirmed that 1-D analysis result is in good agreement with 3-D analysis ones, and that fine 3-D F.E. model, shock analysis one, is also used both in the shock response and whipping response analyses for the analyst effort and time savings.

• 한국수학교육학회지시리즈B:순수및응용수학
• /
• 제20권4호
• /
• pp.259-267
• /
• 2013

We present a high-order potential flow model for the motion of hydrodynamic unstable interfaces in cylindrical geometry. The asymptotic solutions of the bubbles in the gravity-induced instability and the shock-induced instability are obtained from the high-order model. We show that the model gives significant high-order corrections for the solution of the bubble.