• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.34 no.1
• /
• pp.1-9
• /
• 2014

A two-dimensional numerical model based on the finite element method was built to simulate the wave propagation phenomena that occur during the ultrasonic time of flight diffraction (TOFD) process. First, longitudinal-wave TOFD was simulated, and the numerical results agreed well with the theoretical results. Shear-wave TOFD was also investigated because shear waves have higher intensity and resolution. The shear wave propagation was studied using three models with different boundary conditions, and the tip-diffracted shear-to-longitudinal wave was extracted from the A-scan signal difference between the cracked and non-cracked specimens. This signal showed very good agreement between the geometrical and numerical arrival times. The results of this study not only provide better understanding of the diffraction phenomena in TOFD, but also prove the potential of shear-wave TOFD for practical application.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2947-2950
• /
• 2007

A 3D human ventricular model is proposed to simulate an integrative analysis of heart physiology and blood hemodynamics. This consists of the models of electrophysiology of human cells, electric wave propagation of tissue, heart solid mechanics, and 3D blood hemodynamics. The 3D geometry of human heart is discretized to a finite element mesh for the simulation of electric wave propagation and mechanics of heart. In cellular level, excitations by action potential are simulated using the existing human model. Then the contraction mechanics of a whole cell is incorporated to the excitation model. The excitation propagation to ventricular cells are transiently computed in the 3D cardiac tissue using a mono-domain method of electric wave propagation in cardiac tissue. Blood hemodynamics in heart is also considered and incorporated with muscle contraction. We use a PISO type finite element method to simulate the blood hemodynmaics in the human ventricular model.

• Journal of Power System Engineering
• /
• v.14 no.5
• /
• pp.48-55
• /
• 2010

When a large ship is advancing in waves, ship undergoes the hydroelastic response, which has influences on structural stability and the fatigue destruction etc. of the ship. Therefore, to predict accurate hydroelastic response, it is necessary to analyze hydroelastic response including fluid-structure interaction. In this research, a ship is divided into many hull elements to calculate the fluid forces and wave exciting forces on each elements using three-dimensional source distribution method. The calculated fluid forces and wave exciting forces are assigned to nodes of hull elements. The neighbor nodes are connected with elastic beam elements. We analyzed hydroelastic responses, and those are formulated by using finite element method. Particularly, to estimate the influence of forward speed on the hydroelastic responses, we use two different methods : Full Hull Rotation Method(FHRM) and Sectional Hull Rotation Method(SHRM).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.174-186
• /
• 1998

A Galerkin's finite element model incorporating infinite elements for modeling of radiation condition at infinity has been developed, which is based on an extended mild-slope equation. To illustrate the validity and applicability of the present model, the example analyses were carried out for a resonance problem in the rectangular harbor of Ippen and Goda (1963) and for wave transformations over circular shoals of Sharp (1968) and Chandrasekera and Cheung (1997). Comparisons with the results obtained by hydraulic experiments and hybrid element method showed that the present model gives very good results in spite of the rapidly varying topography. Numerical experiments were also performed for wave transformations over a circular concave well which may be an alternative to conventional wave barriers.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.28 no.2
• /
• pp.199-204
• /
• 2008

The generation of axisymmetric Lamb waves and interaction with wall thinning (corrosion) defects in hollow cylinders are simulated using the finite element method. Guided wave interaction with defects in cylinders is challenged by the multi-mode dispersion and the mode conversion. In this paper, two longitudinal, axisymmetric modes are generated using the concept of a time-delay periodic ring arrays (TDPRA), which makes use of the constructive/destructive interference concept to achieve the unidirectional emission and reception of guided waves. The axisymmetric scattering by the wall thinning extending in full circumference of a cylinder is studied with a two-dimensional FE simulation. The effect of wall thinning depth, axial extension, and the edge shape on the reflections of guided waves is discussed.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.8
• /
• pp.957-964
• /
• 2011

This investigation presents a finite element method to obtain the transmission properties of bulk elastic waves in piezoelectric band gap structures(phonon crystals) for varying frequencies and modes. To this end, periodic boundary conditions are imposed on a three-dimensional model while both in-plane and out-of-plane modes are included. In particular, the mode decoupling characteristics between in-plane and out-of-plane modes are identified for each electric poling direction and the results are incorporated in the finite element modeling. Through numerical simulations, the proposed modeling method was found to be a useful, effective one for analyzing the wave characteristics of various types of piezoelectric phononic band gap structures.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.4
• /
• pp.574-581
• /
• 2001

Transient linear elastodynamic problems are numerically analyzed in a time-domain by the Finite Element Method, for which the variational formulation based upon the equations of motion in convolution integral is newly derived. This formulation is implicit and does not include the time derivative terms so that the computation procedure is simple and less assumptions are required comparing to the conventional time-domain dynamic numerical algorithms, being able to get the improved numerical accuracy and stability. That formulation is expanded using the semi-discrete approximation to obtain the finite element equations. In the temporal approximation, the time axis is divided equally and constant and linear time variations are assumed in those intervals. It is found that unconditionally stable numerical results are obtained in case of the constant time variation. Some numerical examples are given to show the versatility of the presented formulation.

• Transactions of the Korean Society of Pressure Vessels and Piping
• /
• v.14 no.1
• /
• pp.38-47
• /
• 2018

For safe operation of nuclear power plants, a loose-part monitoring system (LPMS) is used to detect and locate loose-parts within the reactor coolant system, and to estimate their mass and damage potential. There are several methods to estimate mass, such as the center frequency method based on the Hertz's impact theory, a frequency ratio method and so on, but it is known that these methods cannot provide accurate information on impact response for identifying the impact source. Thanks to increasing computing power, finite element analysis (FEA) method recently become an available option to calculate reliably impact response behavior. In this paper, a finite element analysis model to simulate the propagation behavior of the bending wave, generated by a metal ball impact, is validated by performing a series of impact tests and the corresponding finite element analyses for flat plate and shell structures. Also, a FEA-based metal sphere signal map is developed, and then blind tests are performed to verify the map. This study provides an accurate simulation method for predicting the metal impact behavior and for building a metal sphere signal map, which can be used to estimate the mass of loose-parts on site in nuclear power plants.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.137-139
• /
• 2007

To enhance capability on discerning local and regional seismic phases, such as, Pn, Pg, Sn, Rg, etc, within the crust, 2-D numerical forward modeling will be applied to the data obtained from local seismic stations by simulating almost all waves including not only body wave but also surface wave generated without having to explicitly include them under consideration of Q factor. In this study, after getting rid of instrumental response by deconvolution, pseudo-spectral method instead of relying on typical numerical methods, such as, FEM(Finite Element Method) and FDM(Finite Difference Method), will be implemented for 2-D numerical forward modeling by considering velocities of P-wave and S-wave, density, and Q factors. Ultimately, the Power of reaching the enhanced capability on discerning local and regional seismic phases will make it easier for us to identify the seismic source, whether it is originated from man-made explosion or pure earthquake.