• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.69-76
• /
• 2007

High-frequency electromagnetic (EM) wave propagation associated with borehole ground-penetrating radar (GPR) is a complicated phenomenon. To improve the understanding of the governing physical processes, we employ a finite-difference time-domain solution of Maxwell's equations in cylindrical coordinates. This approach allows us to model the full EM wavefield associated with crosshole GPR surveys. Furthermore, the use of cylindrical coordinates is computationally efficient, correctly emulates the three-dimensional geometrical spreading characteristics of the wavefield, and is an effective way to discretise explicitly small-diameter boreholes. Numerical experiments show that the existence of a water-filled borehole can give rise to a strong waveguide effect which affects the transmitted waveform, and that excitation of this waveguide effect depends on the diameter of the borehole and the length of the antenna.

• Journal of the Korean GEO-environmental Society
• /
• v.9 no.6
• /
• pp.5-12
• /
• 2008

In this research, the plate load tests on sand which is reinforced by Tirecell mat were simulated by finite element method (FEM). Tirecell mat made by waste tires has the same function and similar shape to Geocell for soil reinforcement and it can also be used for civil engineering structure. The results were compared with those of field plate load tests for evaluation of suitability of modeling method. From the comparison of both results, it can be seen that the settlements by FEM were very similar to test results with small margin under the ultimate bearing capacity. For the ultimate bearing capacities of two results, difference was very small. After the confirmation of the modelling, reinforcing effects with variation of cover depth and number of reinforcement layers by Tirecell were analyzed additionally. Reinforcing effect decreases with increasing soil cover depth, and this is similar to previous test results by soil cover depth. As the number of reinforcing layers increased, reinforcing effect increased. However at more than 2 reinforcing layers, reinforcing effect was negligible. In conclusion, the modeling method in this research might be used for analysis of reinforced structures using Tirecell mat.

• Geomechanics and Engineering
• /
• v.17 no.1
• /
• pp.31-45
• /
• 2019

In this study the spatial variability of soils is substantiated physically and numerically by using random field theory. Heterogeneous samples are fabricated by combining nine homogeneous soil clusters that are assumed to be elements of an adopted random field. Homogeneous soils are prepared by mixing different percentages of kaolin and bentonite at water contents equivalent to their respective liquid limits. Comprehensive characteristic laboratory tests were carried out before embarking on direct shear experiments to deduce the basic correlations and properties of nine homogeneous soil clusters that serve to reconstitute the heterogeneous samples. The tests consist of Atterberg limits, and Oedometric and unconfined compression tests. The undrained shear strength of nine soil clusters were measured by the unconfined compression test data, and then correlations were made between the water content and the strength and stiffness of soil samples with different consistency limits. The direct shear strength of heterogeneous samples of different stochastic properties was then evaluated by physical and numerical modelling using FISH code programming in finite difference software of $FLAC^{3D}$. The results of the experimental and stochastic numerical analyses were then compared. The deviation of numerical simulations from direct shear load-displacement profiles taken from different sources were discussed, potential sources of error was introduced and elaborated. This study was primarily to explain the mathematical and physical procedures of sample preparation in stochastic soil mechanics. It can be extended to different problems and applications in geotechnical engineering discipline to take in to account the variability of strength and deformation parameters.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.1
• /
• pp.10-17
• /
• 1992

An elliptic finite difference numerical model is presented for the effective calculation of wave transformation on the water of slowly varying sea bed with currents. Elliptic governing equation has been used, which is derived by Joeng(1990) and Chae et al(1989). Numerical results are compared with analytic solution and laboratory data(Ito and Tanimoto, 1972), and show good agreement. This model is very efficient for modelling large area because of no numerical stable condition and high calculation speed. The capability and utility of the model are illustrated by application to current coexisting area.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2001

Wave-induced currents drive nearshore transport processes, and hence an accurate understanding of wave-current interaction is required for proper management of coastal zone. This paper presents details of an adaptive quadtree grid based numerical model of the coupled wave climate and depth-averaged current field. The model accounts for wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, mixing processes, bottom friction effects, and movement of land-water interface at the shoreline. The wave period- and depth-averaged governing equations arc discrctized explicitly by means of an Adarns¬Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Results from the numerical model are in reasonable agreement with the laboratory data of longshore current generated by oblique waves on a plane beach (Visser 1980, 1991).

• Computers and Concrete
• /
• v.15 no.1
• /
• pp.127-140
• /
• 2015

The process of chloride ingress in saturated concrete was presented by a previous study that used a mathematical model for the same as that concrete. This model is to be studied chloride ion diffusion which is considered as a chemical phenomenon and is to be represented the chloride diffusion process to be a nonlinear partial differential equation (PDE). In this paper, this nonlinear PDE is solved by the Kirchhoff transformation to render into a linear PDE. This linear PDE associated with initial and boundary conditions is also solved by the Laplace transformation to obtain an analytical solution. To verify the serviceability and reliability of this proposed method, the practical application should be supplied. The input parameters were cited from the previous study. The free chloride concentration profiles obtained by the analytical solution of mathematical model for saturated concretes after 24 and 120 hrs of exposure were compared with the previous study. The predicted results obtained from proposed method have a tendency with experimental results obtained by the previous study and trend toward numerical results approximated by finite difference technique.

• Journal of the Korean Solar Energy Society
• /
• v.21 no.1
• /
• pp.33-42
• /
• 2001

The paper discusses the modelling of the thermal storage roof with the air circulation system. In this system, the fully glazed absorber plate is put on the top of the conventional pitch roof made of massive concrete and acts as a solar air heater. Solar radiation collected into absorber is stored in the roof structure by radiation and convection so that it reduces the nighttime heating load through the roof. Another part of the energy is also transmitted to internal air drawn into the channel and is then introduced Into the room. To analysis the system, the energy balance equations are developed and are solved using a finite difference method. The calculation results show a good agreement with the measured ones obtained from our experiments. From the results, it is seen that the thermal storage roof with the air circulation system reduces significantly the conductive heat loss compared with that for the conventional roof and has the instantaneous solar collection efficiency of about 30%.

• Proceedings of the KIEE Conference
• /
• 1995.07a
• /
• pp.455-459
• /
• 1995

Catenary/pantograph system consists of overhead lines which have non-uniform elasticity and pantographs which move at high speed and give force to the lines, therefore happen to be failed in contacts between both from time to time. In this study, as the first step to develop a dynamic simulation program, the general theory is discussed for catenary/pantograph system and appropriate modelling. And comparison is conducted with the references after making a program which referred to the contact force equation algorithm. On this algorithm, the unknown contact force is computed by the equations which was induced as combining catenary and pantograph motion equations expressed in finite difference form. Another simulation program based on the assumed contact forces algorithm was developed. In this algorithm, numerical integraion of both the overhead line and pantograph equations, which without combining, are effected for two assumed values of contact force. The correct contact force is then obtained from these two sets of results by linear interpolation to satisfy the contact condition. Through the comparative review on the outputs from this program, it is verified that this algorithm is reliable.

• Geomechanics and Engineering
• /
• v.24 no.5
• /
• pp.443-456
• /
• 2021

Structural design of the vertical displacements and shear strains in the earth fill (EF) dams has great importance in the structural engineering problems. Moreover, far fault earthquakes have significant seismic effects on seismic damage performance of EF dams like the near fault earthquakes. For this reason, three dimensional (3D) earthquake damage performance of Oroville dam is assessed considering different far-fault ground motions in this study. Oroville Dam was built in United States of America-California and its height is 234.7 m (770 ft.). 3D model of Oroville dam is modelled using FLAC3D software based on finite difference approach. In order to represent interaction condition between discrete surfaces, special interface elements are used between dam body and foundation. Non-reflecting seismic boundary conditions (free field and quiet) are defined to the main surfaces of the dam for the nonlinear seismic analyses. 6 different far-fault ground motions are taken into account for the full reservoir condition of Oroville dam. According to nonlinear seismic analysis results, the effects of far-fault ground motions on the nonlinear seismic settlement and shear strain behaviour of Oroville EF dam are determined and evaluated in detail. It is clearly seen that far-fault earthquakes have very significant seismic effects on the settlement-shear strain behaviour of EF dams and these earthquakes create vital important seismic damages on the swelling behaviour of dam body surface. Moreover, it is proposed that far-fault ground motions should not be ignored while modelling EF dams.

• Journal of the Korean GEO-environmental Society
• /
• v.17 no.8
• /
• pp.29-37
• /
• 2016

In order to prevent leakage of contaminants in offshore landfill, vertical barrier should be installed. Vertical barrier should be installed at designed depth of seabed to prevent the horizontal transport of contaminant in the subsurface. In this study, the seepage and contaminant transport in the subsurface according to embedded depth of vertical barrier were analyzed by using 2-D finite element analysis program SEEP/W and 3-D finite difference analysis program Visual Modflow. Numerical modelling results show that seepage flux and contaminant transport in seabed was greatly reduced when vertical barrier was installed at certain depth of low permeable layer. Therefore, the determination of minimum embedded depth for preventing contaminant leakage is helpful to design the economical vertical barrier.