• Coupled systems mechanics
• /
• v.2 no.1
• /
• pp.85-110
• /
• 2013

Generally, mass concrete structural behavior is governed by the strain components. However, relevant guidelines in dam engineering evaluate the structural behavior of concrete dams using stress-based criteria. In the present study, strain-based criteria are proposed for the first time in a professional manner and their applicability in seismic failure evaluation of an arch dam are investigated. Numerical model of the dam is provided using NSAD-DRI finite element code and the foundation is modeled to be massed using infinite elements at its far-end boundaries. The coupled dam-reservoir-foundation system is solved in Lagrangian-Eulerian domain using Newmark-${\beta}$ time integration method. Seismic performance of the dam is investigated using parameters such as the demand-capacity ratio, the cumulative inelastic duration and the extension of the overstressed/overstrained areas. Real crack profile of the dam based on the damage mechanics approach is compared with those obtained from stress-based and strain-based approaches. It is found that using stress-based criteria leads to conservative results for arch action while seismic safety evaluation using the proposed strain-based criteria leads to conservative cantilever action.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.1-9
• /
• 2008

The free-surface motions interacting with structures are investigated numerically using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka et al. (1996) for solving incompressible flow. In the method, Lagrangian moving particles are used instead of Eulerian approach using grid system. Therefore the terms of time derivatives in Navier-Stokes equation can be directly calculated without any numerical diffusion or instabilities due to the fully Lagrangian treatment of fluid particles and topological failure never occur. The MPS method is applied to the numerical study on the fluid impact loads for wet-drop tests in a LNG tank, and the results are compared with experimental ones.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.5
• /
• pp.424-432
• /
• 2002

Numerical analysis has been conducted to characterize deposition rates of aerosol particles onto a heated, rotating disk with electrostatic effect under the laminar flow field. The particle transport mechanisms considered were convection, Brownian diffusion, gravitational settling, thermophoresis and electrophoresis. The aerosol particles were assumed to have a Boltzmann charge distribution. The electric potential distribution needed to calculate local electric fields around the disk was calculated from the Laplace equation. The Coulomb, the image, the dielectrophoretic and the dipole-dipole forces acting on a charged particle near the conducting rotating disk were included in the analysis. The averaged particle deposition vetocities and their radial distributions on the upper surface of the disk were calculated from the particle concentration equation in a Eulerian frame of reference, along with a rotation speed of 0∼1,000rpm, a temperature difference of 0∼5K and a charged disk voltage of 0∼1000V.Finally, an approximate deposition velocity model for the rotating disk was suggested. The present numerical results showed relatively good agreement with the results of the present approximate model and the available experimental data.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.26 no.1
• /
• pp.42-47
• /
• 2014

Adequate relative humidity is needed in the common living space, for human health, and industrial space, for maintenance and efficiency. Evaporative humidifier systems generally have slow response tendencies of air loads, but relatively low initial investment and maintenance costs. Humidification efficiency in the evaporative humidifier is dependent mostly on the shapes and arrays of humidifying modules. So, we tried to apply the computational fluid dynamics to the evaporative humidifier systems, and studied the humidification efficiency of evaporative humidifiers, through the outlet air conditions.

• Journal of The Korean Astronomical Society
• /
• v.40 no.4
• /
• pp.165-169
• /
• 2007

The numerical simulations associated with the interaction of High Velocity Clouds (HVC) with the Magnetized Galactic Interstellar Medium (ISM) are a powerful tool to describe the evolution of the interaction of these objects in our Galaxy. In this work we present a new project referred to as Theoretical Virtual i Observatories. It is oriented toward to perform numerical simulations in real time through a Web page. This is a powerful astrophysical computational tool that consists of an intuitive graphical user interface (GUI) and a database produced by numerical calculations. In this Website the user can make use of the existing numerical simulations from the database or run a new simulation introducing initial conditions such as temperatures, densities, velocities, and magnetic field intensities for both the ISM and HVC. The prototype is programmed using Linux, Apache, MySQL, and PHP (LAMP), based on the open source philosophy. All simulations were performed with the MHD code ZEUS-3D, which solves the ideal MHD equations by finite differences on a fixed Eulerian mesh. Finally, we present typical results that can be obtained with this tool.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.177-183
• /
• 2018

This study simulated ice load and the motion response of a moored semi-submersible rig in pack-ice conditions using a finite element method. Ice flows of random size and shape were modeled, and interactions for ice-sea, ice-structure, ice-ice were simulated using a simplified method. Parameters for the simplified method such as drag force coefficient and the pressure-penetration relation were obtained based on the result of detailed analysis using the coupled Eulerian-Lagrangian method. The mooring lines were modeled by spring elements based on their stiffness. As a result of the simulation over 1,400 seconds, the force and motion response of the rig were obtained and validated using discrete elements and compared with the results found by the Krylov State Research Centre.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.4
• /
• pp.245-252
• /
• 2020

Wave run-up is an important phenomenon that should be considered in ocean structure design. In this study, the wave run-up of a surface-piercing circular cylinder was calculated in the time domain using the three-dimensional linear and fully nonlinear numerical wave tank (NWT) techniques. The NWT was based on the boundary element method and the mixed Eulerian and Lagrangian method. Stokes second-order waves were applied to evaluate the effect of the nonlinear waves on wave run-up, and an artificial damping zone was adopted to reduce the amount of reflected and re-reflected waves from the sidewall of the NWT. Parametric studies were conducted to determine the effect of wavelength, wave steepness, and the draft of the cylinder on the wave run-up of the cylinder. The maximum wave run-up value occurred at 0°, which was in front of the cylinder, and the minimum value occurred near the circumferential angle of 135°. As the diffraction parameter increased, the wave run-up increased up to 1.7 times the wave height. Furthermore, the wave run-up was 4% higher than the linear wave when the wave steepness was 1/35. In particular, the crest height of the wave run-up increased by 8%.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.21-38
• /
• 2006

Programmed grammars, one of the most important and well investigated classes of grammars with context-free rules and a mechanism controlling the application of the rules, can be described by graphs. We investigate whether or not the restriction to special classes of graphs restricts the generative power of programmed grammars with erasing rules and without appearance checking, too. We obtain that Eulerian, Hamiltonian, planar and bipartite graphs and regular graphs of degree at least three are pr-universal in that sense that any language which can be generated by programmed grammars (with erasing rules and without appearance checking) can be obtained by programmed grammars where the underlying graph belongs to the given special class of graphs, whereas complete graphs, regular graphs of degree 2 and backbone graphs lead to proper subfamilies of the family of programmed languages.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

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