• Journal of the Korea Institute of Military Science and Technology
• /
• v.20 no.3
• /
• pp.384-392
• /
• 2017

A linear explosively formed penetrator(LEFP) is a modification of the explosively formed penetrator(EFP). An EFP is axisymmetric and has a dish-shaped liner while LEFP has a rectangular-shaped liner with curved cross section. Upon detonating LEFP forms laterally wide projectile like blade, leaving a long penetration hole on the target. On the other hand, a long-rod tungsten heavy alloy(WHA) penetrator is one of the major threats against most of the ground armored vehicles. In this paper, the feasibility of using an LEFP in protecting against a long-rod WHA penetrator by colliding LEFP into the threat was investigated through a set of numerical simulations. In this study, a scale-down WHA penetrator with length to diameter ratio(L/D) of 10.7 and 7.0 mm diameter was used to represent a long-rod WHA penetrator. LS-DYNA and Multi-Material ALE technique were employed for the simulation. For estimation of the protection effect by LEFP, residual penetration depths into RHA by the threat were compared according to various impact locations against the threat.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.5
• /
• pp.79-87
• /
• 2002

Acoustic responses to pressure oscillations in axisymmetric combustion chamber are numerically investigated to examine the qualitative trend of acoustic instability in liquid rocket engine. Chamber operating condition and excitation frequency of oscillating pressure are selected as exciting parameters of acoustic instability. Artificial perturbation is simulated by total-pressure oscillation with sine wave at chamber inlet. Many approximations and simplifications are introduced without losing the essence of acoustic pressure response. First, steady-state solution for each operating condition is obtained and next, transient analysis is conducted. Depending on operating condition and excitation frequency, the distinct response characteristics are brought. Weak-strength flames and high-frequency excitation tend to cause sensitive acoustic pressure response leading to unstable pressure field. These results are analyzed based on the correlation with acoustic pressure responses from the previous works adopting laminar flamelet model.

• International Journal of Highway Engineering
• /
• v.14 no.2
• /
• pp.11-17
• /
• 2012

This research describes how to predict a model of the tensile strain related to the fatigue cracking performance of several asphalt concrete structures through design of experiments(e.g., Response Surface Methodology) and harmony search(HS) algorithm. The axisymmetric analysis program of finite element method, which is the KICTPAVE, was used to determine the strain level at the interface layer between asphalt layer and lean concrete layer. Once the training database set of various strain levels was constructed under the several condition of layer stiffnesses and thicknesses in the asphalt concrete structures, the data set was trained through the HS algorithm in order to determine the regression coefficients defined based on a response surface methodology. Furthermore, the testing set, which was not used for the training procedure of HS algorithm, was also constructed in order to evaluate whether the regression coefficients of a prediction model can be appropriately applied for other cases in asphalt concrete structures.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.9
• /
• pp.615-622
• /
• 2017

A normal shock wave is initially formed in the shock tube that migrates towards the closed end of the tube, which, in turn, leads to the reflection of shock. Due to the interaction of the reflected shock with the boundary layer, bifurcation of the shock wave takes place. A shock train will be generated after the bifurcated shock wave approaches the contact surface. Until now, only a few studies have been conducted to investigate this shock train phenomenon inside the shock tube. For the present study, a CFD analysis has been performed on a two dimensional axisymmetric model of a shock tube using unsteady, compressible Navier-Stokes equations. In order to investigate the detailed characteristics of the shock train phenomenon, quantitative studies have been performed by varying shock tube length, diameter under fixed diaphragm, and pressure ratio inside a shock tube.

• Journal of the Korean Magnetics Society
• /
• v.6 no.2
• /
• pp.106-111
• /
• 1996

Electric machines such as motors which have rmving parts are designed for producing mechanical force or torque. The accurate calculations of electromagnetic force and torque are important in the design these machines. Electromagnetic force calculation method using the results of Finite Element Method(FEM) has been presented variously in 2-D problems. Typically the Maxwell's Stress Tensor method and the method of virtual work are used. The former calculates forces by integrating the surface force densities which can be expressed in terms of Maxwell Stress Tensor(MST), and the latter by differentiating the electromagnetic energy with respect to the virtual dis¬placement of rigid bodies of interest. In the problems including current source, magnetic vector potentials(MVP) have rmstly been used as unknown variables for field analysis by a numerical method; e. g. FEM. This paper, thus, introduces the two both methods using MVP in 3-D case. To verify the usefulness of presented methods, a solenoid model is chosen and analyzed by 3-D and axisymmetric FEM. It is found that the force calculation results are in good agreement for several mesh schemes.

• Geotechnical Engineering
• /
• v.13 no.4
• /
• pp.95-108
• /
• 1997

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation However, depending on the assumptions and analytical techniques, it could considerably vary even for a specific degree of consolidation. In this paper, a method is proposed to determine a consistent coefficient of consolidation over all ranges of degree of consolidation by applying the concept of the Optimum Design Technique. The initial excess pore pressure distribution is assumed to be obtainable by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of two dimensional linear-uncoupled axisymmetric consolidation analysis. The minimization of the differences between the measured and the predicted excess pore pressures was carried by BFGS unconstrained optimum design algorithm with one dimensional golden section search technique. By analyzing numerical and real field examples, it can be found that the adopted optimum technique gives a consistent and convergent results.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.1
• /
• pp.1-8
• /
• 1996

A potential-based panel method has been developed for numerical computation of wave resistance on a hybrid hydrofoil. Hybrid hydrofoil is composed of a main body, two struts and two hydrofoils. The main body, which is assumed to be an axisymmetric body for the present analysis, is normally used to support displacement of a body with its buoyancy. Normal dipoles and the sources are distributed on the body(main body, struts, hydrofoils) and the sources are distributed on the free surface. Linearized free surface and the radiation conditions are satisfied using the fourth order finite difference operator and the semi-linear pressure Kutta condition is used for the numerical computation of the hydrofoils. Poisson type free surface condition has been used for the numerical computation and hyperboloidal panel method has been used for better numerical accuracy. To verify this numeric method, model tests are performed in circulation water channel. From the comparison of experimental results with numeric ones, the present method can be used as a useful tool for the design of high speed vessels.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.133-138
• /
• 2002

A 3-D axisymmetric computer program is developed to predict the NO behavior in SNCR system for the stoker incinerator with the waste treatment capacity, 200ton/day. To this end a turbulent reacting flow field calculation is made using proper assumption and empiricism. The stoker bed is assumed to be a homogeneous waste-volatilized gaseous state. The initial composition or reactants are assumed based on the data of the ultimate analysis. Turbulent is resolved by k-e model and turbulent reaction is handled by eddy-breakup model harmonized with empirical chemistry data for gaseous combustion, NO and urea reaction. The liquid droplet is traced by Lagrangian method incorporated by aerodynamic drag, Coriolis and crntrifugal forces. Radiation is treated by sensible heat loss model. Calculation results are in good agreement with experimental data at the outlet of post combustion chamber in Daejon 4th industrial complex. The flue gas shows the temperature range of $900\sim1000^{\circ}C$, velocity of 5m/s and NO concentration of 140ppm at the exit while the measured temperature, flue gas velocity and NO concentration are $967^{\circ}C$, $3\sim4m/s$ and $100\sim200ppm$respectively. Using the developed computer program a parametric study has been made with the variation of heat content of waste, castable length and SNCR variables for the determination of proper injector location. In general, the calculated results are consistent and physically acceptable.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.729-732
• /
• 2010

The present study has been motivated by the development of a reliable numerical methodology for simulation of kerosene/LOx coaxial swirl injectors. To deal with thermodynamic non-ideality and anomalies of transport properties pronounced at supercritical pressures, a set of subroutine libraries has been constructed based on the cubic equations of state, and applied to an existing flamelet analysis code. For computational efficiency, two-dimensional axisymmetric RANS formulation with swirl was adopted and validated successfully against an isothermal coaxial swirling jet. For the actual problem with high pressure combustion, however, numerical results show that the RANS models yield excessive production of turbulence probably due to high density gradient magnitude in the vicinity of mixing layer of swirling film flow, and imply strongly further improvement of the turbulence models.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.2
• /
• pp.39-45
• /
• 2016

In a coaxial porous injector, a gas propellant is injected through the porous cylinder surface to the liquid jet which is encircled by a porous cylinder. In this study, to observe the differences in disintegration mechanisms between a shear coaxial injector and a coaxial porous injector, cold-flow tests and 2-D axisymmetric numerical analysis have been carried out. The shadowgraph images and Sauter mean diameters were compared in similar experimental conditions, and the effects of velocity distributions at the inner injector region on the disintegration of liquid jet were investigated through the numerical calculations. As a result, in high air mass flow rate condition, the disintegration performance of coaxial porous injector is better than shear coaxial injector, in spite of a lower velocity at the inner injector region.