• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.419-424
• /
• 2004

we performed the underwater explosion analysis for the liquefied oxygen tank - a kind of fuel tank of a mid-size submarine, and tried to verify the structural safety for this structure. First, we reviewed the theory and application of underwater explosion analysis using Structure-Fluid Interaction technique and its finite element modeling scheme. Next, we modeled the explosive and sea water as fluid elements, the LOX tank as structural elements and the interface between two regions as ALE scheme. The effect on shock pressure and impulse of fluid mesh size and shape are also investigated. As the analysis result, the shock pressure due explosion propagated into the water region and hit the structure region. The plastic deformation and the equivalent stress highly appeared at the web frame and the shock mount of LOX structure, but these values were acceptable for design criteria.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1022-1030
• /
• 1999

A numerical study was carried out to analyze the effect of flow distribution of stirred part and plug flow part on combustion efficiency at the coal gasification process in an entrained bed coal reactor. The model of computation was based on gas phase eulerian balance equations of mass and momentum. The solid phase was described by lagrangian equations of motion. The $k-{\varepsilon}$ model was used to calculate the turbulence flow and eddy dissipation model was used to describe the gas phase reaction rate. The radiation was solved using a Monte-Carlo method. One-step parallel two reaction model was employed for the devolatilization process of a high volatile bituminous Kideco coal. The computations agreed well with the experiments, but the flame front was closer to the burner than the measured one. The flow distribution of a stirred part and a plug flow part in a reactor was a function of the magnitude of recirculation zone resulted from the swirl. The combustion efficiency was enhanced with decreasing stirred part and the maximum value was found around S=1.2, having the minimum stirred part. The combustion efficiency resulted from not only the flow distribution but also the particle residence time through the hot reaction zone of the stirred part, in particular for the weak swirl without IRZ(internal recirculation zone) and the long lifted flame.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.6
• /
• pp.447-454
• /
• 2017

A three-dimensional two-phase large eddy simulation(LES) has been conducted to investigate the breakup and atomization of liquid jets such as a diesel jet in parallel flow and water jet in cross flow. Gas-liquid two-phase flow was solved by a combined model of Eulerian for gas flow and Lagrangian for a liquid jet. Two stochastic breakup models were implemented to simulate the liquid column and droplet breakup process. The penetration depth and SMD(Sauter Mean Diameter) were analyzed, which was comparable with the experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1335-1340
• /
• 2007

The media transport systems, such as printers, copy machines, facsimiles, ATMs, cameras, etc. have been widely used and being developed rapidly. In the development of those sheet-handling machineries, it is important to predict the static and dynamic behavior of the sheet with a high degree of reliability because the sheets are fed and stacked at such a high speed. Flexible media are very thin, light and flexible, so they behave in geometric nonlinearity with large displacement and large rotation but small strain. In the flexible media analysis, aerodynamic effect from the surrounding air must be included because any small force can make large deformation. In this paper, surrounding air was modeled by incompressible Navier-Stokes flow and an arbitrary Lagranigan-Eulerian(ALE) finite element method with automatic mesh-updating technique was formulated for large domain changes. In the numerical simulations, the results with consideration of the air fast decayed and converged into static results while the results without considering air oscillated continuously.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• Journal of the Korean Society of Combustion
• /
• v.4 no.1
• /
• pp.1-15
• /
• 1999

In order to analyze the sensitivity on the pulverized coal flames of the several variables, a numerical study was conducted at the gasification process. Eulerian approach is used for the gas phase, whereas lagrangian approach is used for the solid phase. Turbulence is modeled using the standard $k-{\varepsilon}$ model. The turbulent combustion incorporates eddy dissipation model. The radiation was solved using a Monte-Carlo method. One-step two-reaction model was employed for the devolatilization of Kideco coal. In pulverized flame of long liftoff height, the initial turbulent intensity seriously affects the position of flame front. The radiation heat transfer and wall heat loss ratio distort the temperature distributions along the reactor wall, but do not influence the reactor performance such as coal conversion, residence time and flame front position. The primary/secondary momentum ratio affects the position of flame front, but the coal burnout is only slightly influenced. The momentum ratio is a variable only associated with the flame stabilization such as flame front position. The addition of steam in the reactor has a detrimental effect on all the aspects, particularly reactor temperature and coal burnout.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.4
• /
• pp.36-41
• /
• 2011

In the present paper, a direct forcing/fictitious domain (DF/FD) level set method is proposed to simulate the FSI (fluid-solid interaction) in two-phase flow. The main idea is to combine the direct-forcing/fictitious domain (DF/FD) method with the level set method in the Cartesian coordinates. The DF/FD method is a non-Lagrange-multiplier version of a distributed Lagrange multiplier/fictitious domain (DLM/FD) method. This method does not sacrifice the accuracy and robustness by employing a discrete ${\delta}$ (Dirac delta) function to transfer quantities between the Eulerian nodes and Lagrangian points explicitly as the immersed boundary method. The advantages of this approach are the simple concept, easy implementation, and utilization of the original governing equation without modification. Simulations of various water-entry problems have been conducted to validate the capability and accuracy of the present method in solving the FSI in two-phase flow. Consequently, the present results are found to be in good agreement with those of previous studies.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.6
• /
• pp.12-16
• /
• 2009

In this study, the wave profiles and kinematics of highly nonlinear waves at various water depths were calculated using a 2D fully nonlinear Numerical Wave Tank (NWT). The NWT was developed based on the Boundary Element Method (BEM) with the potential theory and the mixed Eulerian-Lagrangian (MEL) time marching scheme by 4th-order Runge-Kutta time integration. The spatial variation of intermediate-depth waves along the direction of wave propagation was caused by the unintended generation of 2nd-order free waves, which were originally investigated both theoretically and experimentally by Goda (1998). These free waves were induced by the mismatch between the linear motion of wave maker and nonlinear displacement of water particles adjacent to the maker. When the 2nd-order wave maker motion was applied, the spatial modulation of the waves caused by the free waves was not observed. The respective magnitudes of the nonlinear wave components for various water depths were compared. It was found that the high-order wave components greatly increase as the water depth decreases. The wave kinematics at various locations were calculated and compared with the linear and the Stokes 2nd-order theories.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1368-1379
• /
• 2000

A numerical study on swirling pulverized coal combustion in an axisymmetric enclosure is carried out by applying the 2-phase weighted sum of gray gases model (WSGGM) approach with the discrete ordinate method (DOM) to model the radiative heat transfer equation. In the radiative transfer equation, the same polynomial equation and coefficients for weighting factors as those for gas are adopted for the coal/char particles as a function of partial pressure and particle temperature. The Eulerian balance equations for mass, momentum, energy, and species mass fractions are adopted with the standard and RNG k-${\varepsilon}$ turbulence model, whereas the Lagrangian approach is used for the particulate phase. The eddy-dissipation model is employed for the reaction rate for gaseous mixture, and the single-step and two-step first-order reaction model for the devolatilization process for coal. Special attention is given to establish the thermal boundary conditions on radiative transfer equation By comparing the numerical results with experimental ones, the radiation model used here is confirmed and found to provide an alternative for simulating the radiative transfer.