• Proceedings of the Korean Nuclear Society Conference
• /
• 2004.10a
• /
• pp.263-264
• /
• 2004

• Proceedings of the Korean Nuclear Society Conference
• /
• 2002.05a
• /
• pp.136-136
• /
• 2002

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.95-103
• /
• 2007

An analytical closed-form solution for wave propagation velocity and damping in saturated porous media is presented in this paper The fully coupled field model with compressible solid Brains and pore water were used to derive this solution. An engineering approach for the analysis of fully saturated porous media was adopted and closed-form solutions for one dimensional wave propagation in a homogeneous domain were derived. The solution is highly versatile in that it considers compression of the solid grains, compression of the pore water, deformation of the porous skeleton, and spatial damping and can be used to compute wavespeeds of first and second kind and damping coefficients in various geologic materials. This solution provides a means of analyzing the influence of material property variations on wavespeed and attenuation. In Part 2 of this work the theoretical solution is incorporated into the numerical code and the code is used in a parametric study on wave propagation velocity and damping.

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

• Journal of the Korean GEO-environmental Society
• /
• v.18 no.2
• /
• pp.53-58
• /
• 2017

In this paper, we propose a new method to debond between mixed finite elements for porous media in ABAQUS (2014). ABAQUS just provides debonding algorithm for the u-p model using cohesive elements in standard version. However, this approach has a drawback that it is hard to simulate complex debonding problems like element separation, rigid body motion, and contact between separated elements in standard version. ABAQUS-explicit can resolve these complex problems, but cohesive elements for the u-p model cannot be applied. We introduce a new algorithm for debonding for porous media instead of using cohesive elements. In this method, subroutines VUMAT to apply constitutive models and VDISP to separate elements in ABAQUS are used to simulate debonding problems. In addition, a simple 2-D example is demonstrated in the ABAQUS-explicit solver.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.5 no.4
• /
• pp.243-248
• /
• 1993

In this research, unsteady groundwater flow in unconfined and homogeneous three layer aquifers is studied theoretically and experimentally. Numerical solutions are obtained by Runge Kutta and Runge Kutta Gill method after transforming the governing nonlinear partial differential equations to nonlinear ordinary differential equations. Experimental apparatus includes a test section filled with fine, medium and coarse sands. Experimental results are compared with the numerical solutions and both experimental and numerical results correspond well with each other. This numerical approach may be also applied to the cases which have more aquifers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.9
• /
• pp.89-96
• /
• 2006

In this study, fluid flow and thermal characteristics in a catalyst bed for nitrous oxide catalytic decomposition which is introduced as a hybrid rocket ignition system for small satellites were theoretically considered. To analyze the thermal phenomena of the catalyst bed, a so-called porous medium approach has been opted for modeling the honeycomb geometry of the catalyst bed. Using a Brinkman-extended Darcy model for fluid flow and the one-equation model for heat transfer, the analytical solutions for both velocity and temperature distributions in the catalyst bed are obtained and compared with experimental data to validate the porous medium approach. Based on the analytical solutions, parameters of engineering importance are identified to be the porosity of the catalyst bed, effective volumetric ratio, the ratio of the radius of the catalyst bed to the radius of a pore, heat flux generated by a heater, and pumping power. Their effects on thermal phenomena of the catalyst bed are studied.

• Nuclear Engineering and Technology
• /
• v.45 no.2
• /
• pp.191-202
• /
• 2013

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy may be due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Nuclear Engineering and Technology
• /
• v.55 no.12
• /
• pp.4561-4569
• /
• 2023

A correct understanding of vibration-based degradation is crucial from the standpoint of maintenance for Steam Generators (SG) as crucial mechanical equipment in nuclear power plants. This study has established a novel approach to developing a model for investigating tube bundle degradation according to crack growth caused by two-phase Flow-Induced Vibration (FIV). An important step in the approach is to calculate the two-phase flow field parameters between the SG tube bundles in various zones using the porous media model to determine the velocity and vapor volume fraction. Afterward, to determine the vibration properties of the tube bundles, the Fluid-Solid Interaction (FSI) analysis is performed in eighteen thermal-hydraulic zones. Tube bundle degradation based on crack growth using the sixteen most probable initial cracks and within each SG thermal-hydraulic zone is performed to calculate useful lifetime. Large Eddy Simulation (LES) model, Paris law, and Wiener process model are considered to model the turbulent crossflow around the tube bundles, simulation of elliptical crack growth due to the vibration characteristics, and estimation of SG tube bundles degradation, respectively. The analysis shows that the tube deforms most noticeably in the zone with the highest velocity. As a result, cracks propagate more quickly in the tube with a higher height. In all simulations based on different initial crack sizes, it was observed that zone 16 experiences the greatest deformation and, subsequently, the fastest degradation, with a velocity and vapor volume fraction of 0.5 m/s and 0.4, respectively.

• Earthquakes and Structures
• /
• v.14 no.4
• /
• pp.313-322
• /
• 2018

The present study is dedicated to investigate the SH body-as well as Love-waves propagation effects in porous media with uncertain porosity and permeability. A unified formulation of the governing equations for one-dimensional (1-D) wave propagation in anisotropic porous layered media is presented deterministically. The uncertainties around the above two cited parameters are taken into account by random fields with the help of Monte Carlo Simulations (MCS). Random samples of the porosity and the permeability are generated according to the normal and lognormal distribution functions, respectively, with a mean value and a coefficient of variation for each one of the two parameters. After performing several thousands of samples, the mathematical expectation (mean) of the solution of the wave propagation equations in terms of amplification functions for SH waves and in terms of dispersion equation for Love-waves are obtained. The limits of the Love wave velocity in a porous soil layer overlaying a homogeneous half-space are obtained where it is found that random variations of porosity change the zeros of the wave equation. Also, the increase of uncertainties in the porosity (high coefficient of variation) decreases the mean amplification function amplitudes and shifts the fundamental frequencies. However, no effects are observed on both Love wave dispersion and amplification function for random variations of permeability. Lastly, the present approach is applied to a case study in the Adapazari town basin so that to estimate ground motion accelerations lacked in the fast-growing during the main shock of the damaging 1999 Kocaeli earthquake.