• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.3
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• Journal of Radiation Protection and Research
• /
• v.28 no.2
• /
• pp.97-108
• /
• 2003

A three-dimensional wind field model based on the variational technique has been developed for estimating the overall wind patterns over a complex terrain. The three-dimensional elliptic partial differential equations on Cartesian and terrain-following coordinates have been established to obtain the Lagrangian multiplier and the adjusted wind velocity. The simulations were performed to evaluate the variations of the velocity vectors on the hemisphere, half-cylinder, and saddle type obstacles. Also, the wind field model in the terrain-following coordinate has been applied for evaluating the characteristics of wind patterns according to the variations of Gauss precision moduli on the hemispheric topography. The results showed that the horizontal and vertical wind components were strongly governed by the selection of the values of Gauss precision moduli.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.959-967
• /
• 1988

In this paper, 2-DOF vehicle suspension system with a semiactive on-off damper was studied for improving the ride comfort. It is known that a nonlinear hydraulic damper, which generates force proportional to the square of the relative velocity, can describe the actual fluid resisting type damper more properly than the traditional viscous damping model. On the other hand, hydraulic damper adoption in analysis makes the system nonlinear and causes difficulties to get the system response. In this work, time domain direct integration method was used to calculate system displacement and acceleration. first of all, the response of the suspension system experiencing a given road profile was optimized by Lagrangian multiplier method within the range of given damping coefficients. The appropriate on-loaf damping values were obtained by averaging the already calculated optimum damping coefficients from Lagrangian techniques. The criterion to control the on-off mechanism was determined by examining the suspension efficiency. It was found that the best out of practically applicable criteria is following the sign (positive and negative) of the multiplication of relative displacement and velocity. Judging from the theoretical calculations, it was proved that the semiactive on-off damper can increase suspension efficiency as much as 8-11% in object function.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• v.1
• /
• pp.181-186
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. The present model of the fluid motion is based on the Navier-Stokes equations incorporating velocity-pressure formulation because of need to model the nonlinear wave interaction with spatially non-uniform current field. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. As an intermediate stage of development, solution procedure and characteristics aspects of the present modeling and numerical method are addressed in detail, and preliminary numerical results prove its accuracy and convergence.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• Ocean Systems Engineering
• /
• v.7 no.3
• /
• pp.299-318
• /
• 2017

The paper aims at illustrating several key issues and ongoing efforts for development of a reliable fully-Lagrangian particle-based solver for simulation of hydroelastic slamming. Fluid model is founded on the solution of Navier-Stokes along with continuity equations via an enhanced version of a projection-based particle method, namely, Moving Particle Semi-implicit (MPS) method. The fluid model is carefully coupled with a structure model on the basis of conservation of linear and angular momenta for an elastic solid. The developed coupled FSI (Fluid-Structure Interaction) solver is applied to simulations of high velocity impact of an elastic aluminum wedge and hydroelastic slammings of marine panels. Validations are made both qualitatively and quantitatively in terms of reproduced pressure as well as structure deformation. Several remaining challenges as well as important key issues are highlighted. At last, a recently developed multi-scale MPS method is incorporated in the developed FSI solver towards enhancement of its adaptivity.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.6
• /
• pp.38-46
• /
• 2011

In order to simulate the ignition-gas injection in the interior ballistics, a two-dimensional analytic code for two-phase flows has been developed. The Eulerian-Lagrangian approach and the low-dissipation simple high-resolution upwind scheme(LSHUS) have been adopted in the numerical code for the propellant combustion of the gun propelling charges. The ghost-cell extrapolation method has been used for the moving boundary in the chamber with the projectile movement. The calculation results of the developed code have been compared and verified through those of the dimensionless IBHVG2 code and the previous one-dimensional code. In comparison with the two-phase flows according to the drag models, the numerical analysis of the muzzle velocity has been affected by the drag model.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.2
• /
• pp.1-9
• /
• 2010

LES(Large eddy simulation) of breakup and droplet atomization of a liquid jet into cross turbulent flow was performed. Two phase flow of gas and liquid phases were modeled by the mixed numerical scheme of both Eulerian and Lagrangian methods for gas and liquid droplet respectively. The breakup process of a liquid column and droplets was observed by implementing the blob-KH wave breakup model. The penetration depth into cross flow was comparable with experimental data for several variants of the liquid-gas momentum flux ratio by varying liquid injection velocity. SMD(Sauter Mean Diameter) distribution downstream of jet was analyzed.