• Journal of the Korean Geotechnical Society
• /
• v.24 no.5
• /
• pp.129-137
• /
• 2008

An active application of concrete track is being expected far the future constructions of Korean railroad. In comparison with the existing ballasted tract, a concrete track is very susceptible for the settlement, since its rehabilitation requires much time and cost. When a concrete track is constructed on fine-grained subgrade soil, excessive cumulative plastic settlements due to repetitive train road may occur. In this case, the settlement of the concrete track may be effectively reduced by installing a small number of small-diameter piles beneath the track. This paper presents the effect of pile installation on the reduction of cumulative plastic settlement of concrete track. A method combining experiential equation and numerical method is proposed. Using an existing experiential equation and the estimated earth pressure distribution, the cumulative plastic strain was calculated. From the results, it is verified that the effects of the pile installation is significant to effectively reduce the cumulative plastic settlement of concrete track. The reduction effects of the cumulative plastic settlement according to the pile number and pile arrangement are presented.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.22 no.1
• /
• pp.53-66
• /
• 1980

The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

• Nuclear Engineering and Technology
• /
• v.22 no.2
• /
• pp.116-127
• /
• 1990

The fission gas release model used In the SPEAR-BETA fuel performance code was modified by use of effective thermal conductivity for cracked fuel and by laking Into account axial fission-gas mixing between the fuel-clad gap and the plenum. With use of this modified model the fission gas release was analyzed under various power ramping conditions of P$_{max}$ and $\Delta$.fP. Effective fuel thermal conductivity that accounts for the effect of fuel tracking was used in calculation of the fuel temperature distribution and the Internal gas pressure under power ramping conditions. Mixing and dilution effects due to axial gas flow were also considered in computing the width and the thermal conductivity of the gap. The effect of axial gas flow w3s solved by the Crank-Nicholson method. The finite difference method was used to save running time in the calculation. The present modified fission-gas release model was validated by comparing its predicted results with experimental data from various lamping tests In the literature and calculated results with use of the models used In the SPEAR-BETA and FEMAXI-IV codes. Results obtained with use of the present modified model showed better agreement with experimental data reported in the literature than those results with use of the latter codes. The fuel centerline temperature calculated with introduction of effective thermal conductivity for centerline temperature calculated with Introduction of effective thermal conductivity for cracked fuel was 200 higher fission gas release predicted with use of the modified model was nearly 6% larger on the average than that calculated by use of the unmodified model used in the SPEAR-BETA code.e SPEAR-BETA code.e.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.3
• /
• pp.23-27
• /
• 2008

The new conjunctive surface-subsurface flow model at a large scale was developed by using a 1-D Diffusion Wave (DW) model for surface flow interacting with the 3-D Volume Averaged Soil-moisture Transport (VAST) model for subsurface flow for the comprehensive terrestrial water and energy predictions in Land Surface Models (LSMs). A selection of numerical implementation schemes is employed for each flow component. The 3-D VAST model is implemented using a time splitting scheme applying an explicit method for lateral flow after a fully implicit method for vertical flow. The 1-D DW model is then solved by MacCormack finite difference scheme. This new conjunctive flow model is substituted for the existing 1-D hydrologic scheme in Common Land Model (CLM), one of the state-of-the-art LSMs. The new conjunctive flow model coupled to CLM is tested for a study domain around the Ohio Valley. The simulation results show that the interaction between surface flow and subsurface flow associated with the flow routing scheme matches the runoff prediction with the observations more closely in the new coupled CLM simulations. This improved terrestrial hydrologic module will be coupled to the Climate extension of the next-generation Weather Research and Forecasting (CWRF) model for advanced regional, continental, and global hydroclimatological studies and the prevention of disasters caused by climate changes.

• Journal of the Korean Society for Railway
• /
• v.12 no.1
• /
• pp.31-38
• /
• 2009

The selection of the support system is an important design parameter in design and construction of the tunnel using the new Australian tunnel method. It is a common practice to select the support based on the rock mass grade, in which the rock mass is classified into five rock groups. The method is applicable if the characteristics of the rock mass are uniform in the direction of tunnel excavation. However, such case is seldom encountered in practice and not applicable when the properties vary along the longitudinal direction. This study performs comprehensive three dimensional finite difference analyses to investigate the ground deformation pattern for cases in which the rock mass properties change in the direction of the tunnel axis. The numerically calculated displacements at the tunnel crown show that the displacement is highly dependent on the stiffness contrast of the rock masses. The results strongly indicate the need to select the support type $0.5{\sim}1.0D$ before the rock mass boundary. The paper proposes a new guideline for selecting the support type based the results of the analyses.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.35 no.5
• /
• pp.291-298
• /
• 2015

An array type probe for monitoring metal tubes is proposed in this paper which utilizes peak value and peak time of a pulsed eddy current(PEC) signal. The probe consists of an array of encircling coils along a tube and the outside of coils is shielded by ferrite to prevent source magnetic fields from directly affecting sensor signals since it is the magnetic fields produced by eddy currents that reflect the condition of metal tubes. The positions of both exciter and sensor coils are consecutively moved automatically so that manual scanning is not necessary. At one position of send-receive coils, peak value and peak time are extracted from a sensor PEC signal and these data are accumulated for all positions to form an array type peak value signal and an array type peak time signal. Numerical simulation was performed using the backward difference method in time and the finite element method for spatial analysis. Simulation results showed that peak value increases and the peak appears earlier as the defect depth or length increases. The proposed array signals are shown to be excellent in reflecting the defect location as well as variations of defect depth and length within the array probe.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.1
• /
• pp.29-39
• /
• 2012

The maximum unit point resistance ($q_{max}$) of rock socketed drilled shafts subjected to axial loads was investigated by a numerical analysis. A 3D Finite Difference Method (FDM) analysis and a Distinct Element Method (DEM) analysis were performed with varying rock elastic modulus (E), discontinuity spacing ($S_j$), discontinuity dip angle ($i_j$), and pile diameter (D). Based on the results of obtained, it was found that the ultimate point resistance ($q_{max}$) increased as rock elastic modulus (E) and rock discontinuity spacing ($S_j$) increased. But, it was found that $q_{max}$ decreased as pile diameter (D) increased. As for the influence of the dip angle of rock discontinuity ($i_j$), it was shown that $q_{max}$ decreased up to 50% of maximum value within the range of $0^{\circ}$ < $i_j$ < $60^{\circ}$ due to the shear failure at rock discontinuities. Furthermore, it was found that if $20^{\circ}{\leq}i_j{\leq}40^{\circ}$, influence of $i_j$ should be taken into account because $q_{max}$ tended to approach a minimum value as $i_j$ approached a value near the friction angle of the discontinuity (${\phi}_j$).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.9 no.1
• /
• pp.9-23
• /
• 1997

A numerical study of a two-layer, stratified flow is investigated, using the implicit finite difference method in one dimension. The results of computational method have been tested and, in case of lock exchange flow, compared with the results of experimental data. The results of model experiments with various interfacial, bottom friction coefficients along with various time weighting factor of numerical scheme and dissipative interface are shown and discussed. Two-layer model experiment has been also carried out to investigate the generation and propagation characteristics of internal tidal wave over the steep bottom topography under stratified condition. The internal wave seems to well radiate through the downstream boundary under the experiments adopting radiation conditions both at two layers and only at upper layer, confirming the applicability of radiational boundary condition in stratified flows. It is also shown that the internal wave through the downstream boundary propagates more actively with increasing thickness of lower layer in the downstream. This implies that the potential tidal energy in the interface will depend upon the thickness of lower layer for the constant thickness of upper layer.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.2
• /
• pp.87-96
• /
• 1999

To analyze the dynamic behavior of structure, direct integration and mode superposition may be utilized in time domain analysis. As finite number of frequencies can give relatively exact solutions, mode superposition is preferable in analyzing structural behavior. In non-linear analysis, however, mode superposition is seldom used since time-varying element stiffness changes stiffness matrix, and the change of stiffness matrix leads to the change of essential constants - natural frequencies and mode shapes. In spite of these difficulties, there are some attempts to adopt mode superposition because of low cost compared to direct integration, but the result is not satisfactory. In this paper, a method using mode superposition in non-linear analysis is presented by separating local element stiffness from global stiffness matrix with the difference between linear and non-linear restoring forces to the external force vectors included. Moreover, the hysteresis model changing with the relative deformation in each floor makes it possible to analyze non-linear behavior of structure. The proposed algorithm is applied to shear beam model and the maximum displacement is compared with the result using direct integration method.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.139-147
• /
• 2006

We had conducted a numerical modeling to investigate seismic properties of gas hydrate with field parameters acquired over the East sea in 1998. We used a 2-D staggered grid finite difference method to generate synthetic elastic seismograms for multi-channel seismic survey, OBC (Ocean Bottom Cable) survey and VCS (Vertical Cable Seismic) survey. The results of this study showed that the method using staggered grid yielded stable results and could be used to seismic imaging. We could find out the high amplitude anomaly and the phase reversal phenomenon of reflection wave at interface between the gas hydrate layer and free gas layer such a BSR (Bottom Simulating Reflector) which is the evidence for existence of gas hydrate in seismic reflection data. And we computed the reflection coefficients at the incident angles corresponding to offset distance with the synthetic seismograms. The reflection coefficients acquired from the numerical modeling were nearly consistent with the reflection coefficient computed by Shuey's equation.