• International Journal of Ocean System Engineering
• /
• v.1 no.1
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

• Nuclear Engineering and Technology
• /
• v.51 no.5
• /
• pp.1261-1271
• /
• 2019

Droplet size and distribution are important parameters determining venturi scrubber performance. In this paper, we proposed physical models for a maximum stable droplet size prediction and upper limit log-normal (ULLN) distribution parameters. For the proposed maximum stable droplet size prediction model, a Eulerian-Lagrangian framework and a Reitz-Diwakar breakup model are solved simultaneously using CFD calculations to reflect the effect of multistage breakup and droplet acceleration. Then, two ULLN distribution parameters are suggested through best fitting the previously published experimental data. Results show that the proposed approach provides better predictions of maximum stable droplet diameter and Sauter mean diameter compared to existing simple empirical correlations including Boll, Nukiyama and Tanasawa. For more practical purpose, we developed the simple, one dimensional (1-D) calculation of Sauter mean diameter.

• Explosives and Blasting
• /
• v.37 no.3
• /
• pp.13-24
• /
• 2019

This paper was prepared to investigate the behavior of fragments in underwater torpedo explosion beneath a frigate or surface ship by using an explicit finite element analysis. In this study, a fluid-structure interaction (FSI) methodology, called the multi-material arbitrary Lagrangian-Eulerian (MM-ALE) approach in LS-DYNA, was employed to obtain the responses of the torpedo fragments and frigate hull to the explosion. The Euler models for the analysis were comprised of air, water, and explosive, while the Lagrange models consisted of the fragment and the hull. The focus of this modeling was to examine whether a worst-case fragment could penetrate the frigate hull located close (4.5 m) to the exploding torpedo. The simulation was performed in two separate steps. At first, with the assumption that the expanding skin of the torpedo had been torn apart by consuming 30% of the explosive energy, the initial velocity of the worst-case fragment was sought based on a well-known experimental result concerning the fragment velocity in underwater bomb explosion. Then, the terminal velocity of the worst-case fragment that is expected to occur before the fragment hit the frigate hull was sought in the second step. Under the given conditions, the possible initial velocities of the worst-case fragment were found to be very fast (400 and 1000 m/s). But, the velocity difference between the fragment and the hull was merely 4 m/s at the instant of collision. This result was likely to be due to both the tremendous drag force exerted by the water and the non-failure condition given to the frigate hull. Anyway, at least under the given conditions, it is thought that the worst-case fragment seldom penetrate the frigate hull because there is no significant velocity difference between them.

• Clean Technology
• /
• v.25 no.3
• /
• pp.206-222
• /
• 2019

Recently, various multiphase flows have been developed, and among them some models have been commercialized. However, most of them have been developed based on a pressure-based approach; therefore, various numerical difficulties were involved inherently. Accordingly, in order to overcome these numerical difficulties, a multiphase flow model, MultiPhaSe flow (MPS), following a fractional-flow based approach was developed. In this study, by combining a contaminant transport module describing an enhanced dissolution effect of a surfactant with MPS, a MultiPhaSe flow and TranSport (MPSTS) model was developed. The developed model was verified using the analytical solution of Clement. The MPSTS model can simulate the process of surfactant enhanced aquifer remediation including interphase mass transfer and contaminant transport in multiphase flow by using the coupled particle tracking method and Lagrangian-Eulerian method. In this study, a surfactant was used in a non aqueous phase liquid (NAPL) contaminated area, and the effect of hydro-geological heterogeneity in the layered media on remediation efficiency was studied using the developed model. According to the numerical simulation, when hydraulic conductivity in a lower layer is 10 times, 20 times, and 50 times larger than that in an upper layer, the concentration of dissolved diesel in the lower layer is much higher than that in the upper layer because the surfactant moves faster along the lower layer owing to preferential flow; thus, the surfactant enhances dissolution of residual non aqueous phase liquid in the lower layer.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.1
• /
• pp.10-15
• /
• 2012

Wave characteristics in the presence of energy damping are investigated using the linear shallow water equations. To get the phase and energy velocities, geometric optics approach is used and then these values are validated through numerical experiments. Energy damping affects wave height, phase and energy velocities which result in wave transformation. When the complex wavenumber is used by the Eulerian approach, it is found that the phase velocity decreases as the damping increases while the energy velocity increases showing higher values than the phase velocity. When the complex angular frequency is used by the Lagrangian approach, the energy-damping wave group is found to propagate in the energy velocity. The energy velocity is found to affect shoaling and refraction coefficient which is verified through numerical experiments for waves on a plane slope.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Journal of Korean Society for Atmospheric Environment
• /
• v.17 no.1
• /
• pp.1-17
• /
• 2001

The high-ozone episode in the Greater Seoul Area for the period of July 27 to August 1 1997 was modeled by the CIT(California Institute of Technology) three-dimensional photochemical model. Emission data were prepared by scaling the NIER(1994) data through and optimization method using VOC measurements in August 1997 and EKMA(Empirical Kinetic Modeling Approach). Two sets of meteorological data were prepared by the diagnostic routine. a part of the CIT model : one only utilized observations from the surface weather stations and the other also utilized observations from the automatic weather stations that were more densely distributed than those from the surface weather stations. The results showed that utilizing observations from the automatic weather stations could represent fine variations in the sind field such as those caused by topography. A better wind field gave better peak ozones and a more reasonable spatial distribution of ozone concentrations. Nevertheless, there were still many differences between predictions and observations particularly for primary pollutant such as NOx and CO. This was probably due to the inaccuracy of emission data that could not resolve both temporal and spatial variations.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.1023-1032
• /
• 1995

The multiple cooperating robot system plays an important role in the research of modern manufacturing system as the emphasis of production automation is more on the side of flexibility than before. While the kinematic and dynamic analysis of a single robot is performed as an open-loop chain, the dynamic formulation of robot in a multiple cooperating robot system differs from that of a single robot when the multiple cooperating robots form a closed kinematic chain holding an object simultaneously. The object may be any type from a rigid body to a multi-joint linkage. The mobility of the system depends on the kinematic configuration of the closed kinematic chain formed by robots and object, which also decides the number of independent input parameters. Since the mobility is not the same as the number of robot joints, proper constraint condition is sought. The constraints may be such that : the number of active robot joints is kept the same as mobility, all robot joints are active and have interrelations between each joint forces/torques, two robots have master-slave relation, or so on. The dynamic formulation of system is obtained. The formulation is based on recursive dual-number screw-calculus Newton-Eulerian approach which has been used for single robot analysis. This new scheme is recursive and compact symbolically and may facilitate the consideration of the object in real time.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.388-397
• /
• 2010

A two-dimensional time-domain, potential-theory-based fully nonlinear numerical wave tank (NWT) was developed by using boundary element method and the mixed Eulerian-Lagrangian (MEL) approach for free-surface node treatment. The NWT was applied to prediction of primary wave energy conversion efficiency of a bottom-mounted oscillating water column (OWC) wave power device. The nonlinear free-surface condition inside the chamber was specially devised to represent the pneumatic pressure due to airflow velocity and viscous energy loss at the chamber entrance due to wave column motion. The newly developed NWT technique was verified through comparison with given experimental results. The maximum energy extraction was estimated with various chamber-air duct volume ratios.

• 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.