• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.299-306
• /
• 2006

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.3
• /
• pp.474-481
• /
• 2002

Fouling deposits due to surface-coated particles have been calculated using CFD techniques. The sticking probabilities of the surface-coated particles are also evaluated on the basis of an energy balance. The sticking probabilities of the deposit surface are also included in the prediction of the deposition occurring through the multiple impaction of particles with the deposit surface. The sticking probability of an impacting particle is expressed in terms of such parameters as particle viscosity, surface tension, impact velocity, impact angle and the thickness of the sticky layer on a particle. Particulate behavior around a tube in cross flow was studied using the Lagrangian approach. Three important parameters i.e. impact velocity, impact angle, and particulate concentration, were used in the prediction of deposition rate. The computational predictions were found to be in good agreement with the experimental data.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.

• Structural Engineering and Mechanics
• /
• v.56 no.6
• /
• pp.917-938
• /
• 2015

The nonlinear numerical analysis of the impact response of reinforced concrete/mortar beam incorporated with the updated Lagrangian method, namely the Smoothed Particle Hydrodynamics (SPH) is carried out in this study. The analysis includes the simulation of the effects of high mass low velocity impact load falling on beam structures. Three material models to describe the localized failure of structural elements are: (1) linear pressure-sensitive yield criteria (Drucker-Prager type) in the pre-peak regime for the concrete/mortar meanwhile, the shear strain energy criterion (Von Mises) is applied for the steel reinforcement (2) nonlinear hardening law by means of modified linear Drucker-Prager envelope by employing the plane cap surface to simulate the irreversible plastic behavior of concrete/mortar (3) implementation of linear and nonlinear softening in tension and compression regions, respectively, to express the complex behavior of concrete material during short time loading condition. Validation upon existing experimental test results is conducted, from which the impact behavior of concrete beams are best described using the SPH model adopting an average velocity and erosion algorithm, where instability in terms of numerical fragmentation is reduced considerably.

• Proceedings of the Korea Institute of Fire Science and Engineering Conference
• /
• 2008.11a
• /
• pp.363-370
• /
• 2008

The present study has been performed to investigate the effect of statistical number of droplets on the simulation of the sprinkler spray using fire field model. In order to simulate the sprinkler spray characteristics, the present study uses NIST Fire Dynamics Simulator version 5.2. A group of Lagrangian particle with similar droplet characteristics, such as diameter, velocity, temperature and so on, is represented by parcel concept to decrease the total number of droplets tracked in the simulation. The present study introduces a new parameter to represent the ratio between real number of droplets and computational parcels. The dependency of the number of parcels on the fire suppression characteristics and spray patterns is quantitatively examined for different ratio between the real number of droplets and computational parcels.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2737-2742
• /
• 2007

Strong wind flow around a building complex was numerically studied by LES. The original motivation of this work stemmed from the efforts to develop a risk assessment technique for windstorm hazards. Lagrangian-averaged scale-invariant dynamic subgrid-scale model was used for turbulence modeling, and a log-law-based wall model was employed on all the solid surfaces including the ground and the surface of buildings to replace the no-slip condition. The shape of buildings was implemented on the Cartesian grid system by an immersed boundary method. Key flow quantities for the risk assessment such as mean and RMS values of pressure on the surface of the selected buildings are presented. In addition, characteristics of the velocity field at some selected locations vital to safety of human beings is also reported.

• Journal of the Korean Society of Industry Convergence
• /
• v.19 no.4
• /
• pp.206-215
• /
• 2016

Delta parallel robots are now widely used for high-speed applications. However, typical Delta robots, such as ABB Flexpicker suffer from rotating axis with passive prismatic joint subjected to critical speed and so requiring careful maintenance. In this paper, a novel 4-DOF high-speed parallel robot with four legs is presented, which consists of three legs with 90 degree arrangement for translational motions and one remaining leg with rack & pinion gears for rotational motion. The inverse kinematics, velocity, acceleration, statics, and inverse dynamics have been analyzed. From the workspace analysis and inverse dynamics simulation for 0.43 sec cycle time, the 4-axis parallel robot prototype with 12kg payload has been designed. In the future research, computed torque control methods will be developed for the prototype.

• International Journal of Ocean System Engineering
• /
• v.2 no.4
• /
• pp.200-208
• /
• 2012

Surge overflow may cause damage on earthen levees. Levee strengthened on the levee crest and landward-side slope can provide protection against the erosion damage induced by surge overflow. In this paper, surge overflow of a roller compacted concrete RCC strengthened levee was studied in a purely Lagrangian and meshless approach, the smoothed particle hydrodynamics (SPH) method. After verifying the developed model with analytical solution and comparing the results with full-scale experimental data, the roughness and erosion parameters were calibrated. The water thickness, flow velocity, and erosion depth at crest, landward-side slope and toe were calculated. The characteristics of flow hydraulics and erosion on the RCC strengthened levee are given. The results indicate that the RCC strengthened levee can resist erosion damage for a long period.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.