• Journal of the Korea Institute of Military Science and Technology
• /
• v.6 no.4
• /
• pp.89-99
• /
• 2003

Smoothed particle hydrodynamics, SPH, is a gridless Lagrangian technique which is a useful alternative numerical analysis method to simulate high velocity deformation problems as well as astrophysical and cosmological problems. The SPH method brings about some difficulties such as tensile Instability and stress oscillation. A new SPH method, so called normalized algorithm, was introduced to overcome these difficulties. In this paper we aimed to estimate this method and have developed an one-dimensional normalized SPH program. The high velocity impact model of an aluminum bar has been analysed by using the developed program and a commercial hydrocode, LS-DYNA. The obtained numerical results showed good agreement with the results of the same model in reference. The program also showed more stable results than those of LS-DYNA in stress oscillation. We hopefully expect that the developed one-dimensional normalized SPH program can be used to solve hydrodynamic problems especially for explosive detonation analysis.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.05a
• /
• pp.215-219
• /
• 2002

In this study, we show the numerical and the experimental results for fluid motions inside a rectangular container subjected to a background rotation added by a rotational oscillation. In the numerical computation, we used a parallel computer system of PC-cluster type. Attention is given to dependence of the flow patterns on the parameter change. It shows that the flow becomes in a periodic state at low Reynolds numbers and undergoes a transition to a chaotic motion at high Reynolds numbers. It also shows that the fluid motion tends to be depth-independent at ${\epsilon}$ up to 0.3 for Re lower than 5235.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.845-851
• /
• 2000

In this paper, we show the numerical and the experimental results of two-dimensional fluid motions inside a rectangular container subjected to a background rotation added by a rotational oscillation. In the PlY experiment we apply a new algorithm, new three step search(NTSS), to the velocity calculation. In the numerical computation, the linear Ekman-pumping model was used to take the bottom friction effect into account. It was found that it well produces the experimental results at low e number.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.4
• /
• pp.214-224
• /
• 2022

It is necessary to estimate manoeuvring characteristics of submerged bodies at the design stage in order to ensure the safe operations. In this study, added mass coefficients in the mathematical model of submerged bodies are estimated by captive model tests and numerical calculations. Two kinds of models, MARIN 'BB2'submarine model and AUV (Autonomous unmanned vehicle) model are utilized in the forced oscillation tests. Compared to BB2 submarine, AUV with cylindrical type hull form shows relatively small added masses in roll, pitch, and yaw directions. Next, numerical calculations based on potential theory are performed under the assumption that viscous effects on inertia forces are negligible. Added masses obtained by numerical calculations are in good agreements with forced oscillation test results. And if slow manoeuvres of submerged bodies are presumed, some of velocity coupled terms can be approximated by combinations of added mass coefficients.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1534-1539
• /
• 2004

The objective of this study is to simulate the etching characteristics with oscillation angle for the optimization of etching system. The etching characteristics were analyzed under different etching conditions. The spray characteristics were measured by Phase Doppler Anemometer (PDA). The correlation between the spray characteristics and the etching characteristics was investigated and used for fundamental data to simulate the etching characteristics with oscillation angle. The smaller coefficient of variation, the more uniform etching characteristic distribution became. It was found that numerical predictions of etching factor generally agreed well with the measured results with distance from nozzle tip. Oscillation leads to decrease of etching factor and increase of uniformity.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.282-286
• /
• 2004

Discharge current oscillation in a 20KHz range is a serious problem for Hall thruster performance, In our previous work, a possibility of controlling the oscillation amplitude by increasing the speed of neutral particles incoming to the ionization zone was predicted in our previous work. In this paper, the effects of diameter of anode orifice on the oscillation phenomena and the thrust performance were evaluated experimentally. The experimental results show that the measured amplitude of oscillation becomes smaller as the diameter of anode orifice. However, the larger orifice makes thrust performance lower. The results of numerical analysis of neutral particles show that these tendencies have much to do with neutral properties.

• 제어로봇시스템학회:학술대회논문집
• /
• 1996.10b
• /
• pp.612-615
• /
• 1996

In the continuous casting process, molten metal contacts the mold wall and the molten metal surface is subject to the mold oscillation. The mold oscillation results in the oscillation marks on the surface of solidified steel, which has undesirable effects on the quality of slabs. In order to reduce the oscillation marks by achieving soft contact of molten metal with the mold surface, alternating magnetic field is applied to the surface of molten metal. However, if the magnetic field strength becomes too strong, the melt flow induced by the magnetic field. causes the instability of the molten metal surface, which has also the bad influence on the slab quality. Therefore, it is very important to choose the optimal position of the inductor coil and the optimal level of electric power to minimize the surface defects. In the present work, as a first step toward the optimization problem of the process, numerical studies are performed to investigate the effects of coil position and the electric power level on the meniscus shape and the flow field. As numerical tools, the boundary integral equation method(BIEM) is used for the magnetic field analysis and the finite difference method (FDM) with orthogonal grid generation is used for the flow analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.115-119
• /
• 2007

In this paper, a flow visualization analysis has been carried out on an oscillating square section cylinder, numerically, using a commercially available code CFD-ACE. In this study, the square cylinder is forced to oscillate at different frequencies of excitation, viz., fe/fo=0.5, 1.0 and 2.0 (where, fe is the excitation frequency provided to the cylinder and fo is the natural vortex shedding frequency from the stationary cylinder at a particular Reynolds number (=5200). In all the cases, the peak-to-peak amplitude of oscillation is kept at 32% of the side dimension of the square cylinder. These studies are conducted to understand the influence of frequency of oscillation on the flow field features around the cylinder, particularly the mode of vortex shedding. Results indicate that, the flow field around a square cylinder is very much influenced by the excitation frequency, in particular the vortex shedding mode. It is also found that, the vortex street parameters are significantly influence by the oscillation frequency. Comparison with earlier reported experimental studies has also been attempted in this paper. In appears that, such a numerical exercise (as performed in this paper) is first of its kind. It is believed that, these studies would enable one to understand the mechanisms underlying the flow-induced vibrations of a square section cylinder.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.1
• /
• pp.116-124
• /
• 2002

The bubble model by Keller and Prosperetti is adapted to solve the nonlinear oscillation of a gas bubble. This formulation leads to accurate results since it introduces the energy equation instead of the polytropic assumption for the bubble interior. The numerical method used in this study is stable enough to handle large amplitude of bubble oscillation. The numerical results show some interesting nonlinear phenomena fur the bubble oscillator. The excitation changes the natural frequency of the bubble and makes some harmonic resonances at $f/f_0=1/2, 1/3$ and so on. The natural frequency of a bubble oscillator decreases compared with the linear case result, which means that the nonlinear bubble oscillation system is a "softening"system. In addition, the frequency response curve jumps up or down at a certain frequency. It is also found that there exist multi-valued regions in the frequency response curve depending on the initial conditions of bubble. The dependency of the bubble motion on the initial condition can generate extremely large pressure and temperature which might be the cause of the acoustic cavitation and the sonoluminescence.inescence.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.26-34
• /
• 2001

The flow past a circular cylinder forced to vibrate transversely is numerically simulated by solving the two-dimensional Navier-Stokes equations modified by the vibration velocity of a circular cylinder at a Reynolds number of 164. The higher-order finite difference scheme is employed for the spatial discretization along with the second order Adams-Bashforth and the first order backward-Euler time integration. The calculated cylinder vibration frequency is between 0.60 and 1.30 times of the natural vortex-shedding frequency. The calculated oscillation amplitude extends to 25% of the cylinder diameter and in the case of the lock-in region it is 60%. It is made clear that the cylinder oscillation has influence on the wake pattern, the time histories of the drag and lift forces, power spectral density and phase diagrams, etc. It is found that these results include both the periodic (lock-in) and the quasi-periodic (non-lock-in) state. The vortex shedding frequency equals the driving frequency in the lock-in region but is independent in the non-lock-in region. The mean drag and the maximum lift coefficient increase with the increase of the forcing amplitude in the lock-in state. The lock-in boundaries are also established from the present direct numerical simulation.