• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.5
• /
• pp.1158-1168
• /
• 1994

For BEM analyses of the impact problems of dissimilar materials, the connected multi-region method using perfect bonded conditions on the interface boundaries was added to two-dimensional Laplace transformed-domain BEM program for a single region analysis. It was confirmed that the BEM results of impact problems of a single-region and multi-regions for a homogeneous isotropic material are agreed well. The two-dimensional Laplace transformed-domain BEM program combined with connected multi-region method was applied to analyse several impact problems of dissimilar materials. Also the feasibility of BEM impact analyses was investigated for dissimilar materials by the analysis of the BEM results for impact problems of dissimilar materials in terms of physical aspects. As for an application, the two-dimensional Laplace transformed BEM concerning impact problems of cracks at the interface of dissimilar materials and the determinating process of the dynamic stress intensity factors by extrapolation method are presented in this paper.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.1
• /
• pp.23-32
• /
• 2011

Numerical analysis of hydrodynamic performance of a movable submerged breakwater was carried out as an eco-friendly marine structure for coastal and harbor protection. Using boundary elements method with two-dimensional frequency-domain reflection and transmission coefficients and wave forces acting on the submerged flat plate were calculated with various submerged depths and respective motion allowable modes. The movable breakwater was found to be more efficient in wave-blocking than the fixed structure. Variation of reflection coefficients was significantly influenced by vertical motion of the body.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.37-49
• /
• 2000

This paper is for the first essential study on the development of unified finite element formulations for solving problems related to the dynamics/thermoelastics behavior of solids. In the first part of formulations, the finite element method is based on the introduction of a new quantity defined as heat displacement, which allows the heat conduction equations to be written in a form equivalent to the equation of motion, and the equations of coupled thermoelasticity to be written in a unified form. The equations obtained are used to express a variational formulation which, together with the concept of generalized coordinates, yields a set of differential equations with the time as an independent variable. Using the Laplace transform, the resulting finite element equations are described in the transform domain. In the second, the Laplace transform is applied to both the equation of heat conduction derived in the first part and the equations of motions and their corresponding boundary conditions, which is referred to the transformed equation. Selections of interpolation functions dependent on only the space variable and an application of the weighted residual method to the coupled equation result in the necessary finite element matrices in the transformed domain. Finally, to prove the validity of two approaches, a comparison with one finite element equation and the other is made term by term.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.375-387
• /
• 2013

Soils consist of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled with water and air. Hence, soil behaviour must be analyzed by incorporating the effects of the transient flow of the pore-fluid through the voids, and therefore requires a two-phase continuum formulation for saturated porous media. The present paper presents briefly the Biot's basic theory of dynamics of saturated porous media with u-P formulation to determine the responses of pore fluid and soil skeleton during cyclic loading. Kelvin elements are attached to transmitting boundary. The Pastor-Zienkiewicz-Chan model has been used to describe the inelastic behavior of soils under isotropic cyclic loadings. Newmark-Beta method is employed to discretize the time domain. The response of fluid-saturated porous media which are subjected to time dependent loads has been simulated numerically to predict the liquefaction potential of a semi-infinite saturated sandy layer using finite-infinite elements. A settlement of 17.1 cm is observed at top surface. It is also noticed that liquefaction occurs at shallow depth. The mathematical advantage of the coupled finite element analysis is that the excess pore pressure and displacement can be evaluated simultaneously without using any empirical relationship.

• Proceedings of the KIEE Conference
• /
• 1988.07a
• /
• pp.370-374
• /
• 1988

This paper is concerned with the investigation of the impurity redistribution process in a two step diffusion. In integrated circuit technology, two step boron diffusion involving a deposition step followed by a drive-in step in commonly encounted. The drive-in process is usually performed in oxidizing atmosphere resulting in redistribution of impurity (boron) within the semiconductor. This paper proposes a new numerical analysis method; Bounary Element Method to determine impurity profile at the arbitrary point in domain by its coordinate and boundary value.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.3
• /
• pp.562-577
• /
• 2014

A reliable steady/transient hydro-elastic analysis is developed for flexible (composite) marine propeller blade design which deforms according to its environmental load (ship speed, revolution speed, wake distribution, etc.) Hydro-elastic analysis based on CFD and FEM has been widely used in the engineering field because of its accurate results however it takes large computation time to apply early propeller design stage. Therefore the analysis based on a boundary element method-Finite Element Method (BEM-FEM) Fluid-Structure Interaction (FSI) is introduced for computational efficiency and accuracy. The steady FSI analysis, and its application to reverse engineering, is designed for use regarding optimum geometry and ply stack design. A time domain two-way coupled transient FSI analysis is developed by considering the hydrodynamic damping ffects of added mass due to fluid around the propeller blade. The analysis makes possible to evaluate blade strength and also enable to do risk assessment by estimating the change in performance and the deformation depending on blade position in the ship's wake. To validate this hydro-elastic analysis methodology, published model test results of P5479 and P5475 are applied to verify the steady and the transient FSI analysis, respectively. As the results, the proposed steady and unsteady analysis methodology gives sufficient accuracy to apply flexible marine propeller design.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.513-528
• /
• 2013

A floating Oscillating Water Column (OWC) wave energy converter, a Backward Bent Duct Buoy (BBDB), was simulated using a state-of-the-art, two-dimensional, fully-nonlinear Numerical Wave Tank (NWT) technique. The hydrodynamic performance of the floating OWC device was evaluated in the time domain. The acceleration potential method, with a full-updated kernel matrix calculation associated with a mode decomposition scheme, was implemented to obtain accurate estimates of the hydrodynamic force and displacement of a freely floating BBDB. The developed NWT was based on the potential theory and the boundary element method with constant panels on the boundaries. The mixed Eulerian-Lagrangian (MEL) approach was employed to capture the nonlinear free surfaces inside the chamber that interacted with a pneumatic pressure, induced by the time-varying airflow velocity at the air duct. A special viscous damping was applied to the chamber free surface to represent the viscous energy loss due to the BBDB's shape and motions. The viscous damping coefficient was properly selected using a comparison of the experimental data. The calculated surface elevation, inside and outside the chamber, with a tuned viscous damping correlated reasonably well with the experimental data for various incident wave conditions. The conservation of the total wave energy in the computational domain was confirmed over the entire range of wave frequencies.

• Smart Structures and Systems
• /
• v.19 no.1
• /
• pp.39-46
• /
• 2017

Finite element model updating is very effective procedure to determine the uncertainty parameters in structural model and minimize the differences between experimentally and numerically identified dynamic characteristics. This procedure can be practiced with manual and automatic model updating procedures. The manual model updating involves manual changes of geometry and analyses parameters by trial and error, guided by engineering judgement. Besides, the automated updating is performed by constructing a series of loops based on optimization procedures. This paper addresses the ambient vibration based finite element model updating of long span reinforced concrete highway bridges using manual model updating procedure. Birecik Highway Bridge located on the $81^{st}km$ of Şanliurfa-Gaziantep state highway over Firat River in Turkey is selected as a case study. The structural carrier system of the bridge consists of two main parts: Arch and Beam Compartments. In this part of the paper, the arch compartment is investigated. Three dimensional finite element model of the arch compartment of the bridge is constructed using SAP2000 software to determine the dynamic characteristics, numerically. Operational Modal Analysis method is used to extract dynamic characteristics using Enhanced Frequency Domain Decomposition method. Numerically and experimentally identified dynamic characteristics are compared with each other and finite element model of the arch compartment of the bridge is updated manually by changing some uncertain parameters such as section properties, damages, boundary conditions and material properties to reduce the difference between the results. It is demonstrated that the ambient vibration measurements are enough to identify the most significant modes of long span highway bridges. Maximum differences between the natural frequencies are reduced averagely from %49.1 to %0.6 by model updating. Also, a good harmony is found between mode shapes after finite element model updating.

• Ocean Systems Engineering
• /
• v.2 no.2
• /
• pp.147-158
• /
• 2012

A time-domain simulation of a land-based Oscillating Water Column (OWC) with various irregular waves as a form of PM spectrum is performed by using a two-dimensional fully nonlinear numerical wave tank (NWT) based on the potential theory, mixed Eulerian-Lagrangian (MEL) approach, and boundary element method. The nonlinear free-surface condition inside the OWC chamber was specially devised to describe both the pneumatic effect of the time-varying pressure and the viscous energy loss due to water column motions. The quadratic models for pneumatic pressure and viscous loss are applied to the air and free surface inside the chamber, and their numerical results are compared with those with equivalent linear ones. Various wave spectra are applied to the OWC system to predict the efficiency of wave-energy take-off for various wave conditions. The cases of regular and irregular waves are also compared.