• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.47-60
• /
• 2004

The permeability coefficients were calculated by the homogenization analysis method with sufficient consideration of fracture geometry dependent on aperture change. According to the results of aperture measurements using a confocal laser scanning microscope, apertures on each measuring point display different magnitudes, indicating that fracture walls can not be assumed as parallel feature. After construction of fracture model based on the aperture values measured on each pressure level, the homogenization analysis was conducted to compute permeability coefficients. The calculated permeability coefficients distribute in the ranges of $10^{-1}~10^{-3}cm/sec$. Most of the specimens show decreasing permeability coefficients with the increase of the applied pressure. However, the decreasing rates of permeability coefficients do not show a constant trend on each pressure level. This phenomenon is well matched to the observation results of Chae et al. (2003). It proves that aperture change strongly influences on permeability characteristics. Three sections of each specimen have all different values of permeability coefficient. It suggests that the variation of permeability coefficient depends sensitively on aperture magnitudes and characteristics of fracture geometry. It is very important to consider accurate fracture geometries for analysis of permeability characteristics in rock fractures bearing different aperture distribution. Therefore, it needs to consider sufficiently the fracture geometries for calculating the permeability coefficients of fractures.

• Geophysics and Geophysical Exploration
• /
• v.12 no.1
• /
• pp.99-104
• /
• 2009

We propose a waveform inversion method for SH-wave data obtained in a shallow seismic refraction survey, to determine a 2D inhomogeneous S-wave profile of shallow soils. In this method, a 2.5D equation is used to simulate SH-wave propagation in 2D media. The equation is solved with the staggered grid finite-difference approximation to the 4th-order in space and 2nd-order in time, to compute a synthetic wave. The misfit, defined using differences between calculated and observed waveforms, is minimised with a hybrid heuristic search method. We parameterise a 2D subsurface structural model with blocks with different depth boundaries, and S-wave velocities in each block. Numerical experiments were conducted using synthetic SH-wave data with white noise for a model having a blind layer and irregular interfaces. We could reconstruct a structure including a blind layer with reasonable computation time from surface seismic refraction data.

• Korean Journal of Materials Research
• /
• v.11 no.5
• /
• pp.372-377
• /
• 2001

A theory developed in Part I has been applied to calculate effective elastic and thermoelastic moduli of particle-strengthened, unidirectionally fiber-reinforced, and layered composites. For the unidirectional fiber composites the effect of fiber aspect ratio is taken into account. The analytical solutions obtained to the effective elastic moduli are compared with some of existing expressions and the following results are found. The effective bulk and shear moduli of the particle strengthened composites coincide with Korner's expressions, which correspond with the lower bounds of Hanshin and Shtrikman. The same expressions as the lower bounds of Hill and Hanshin are obtained for five independent moduli of the aligned continuous fiber composites, four of which coincide with Hanshin and Rosen's exact solutions for 'composite cylinder assemblage'.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.30 no.4A
• /
• pp.259-270
• /
• 2005

This paper proposes a mathematical model to analyze call-and packet-level performance of the TD-CDMA/TDD system which could serve a flexible radio resource management against multi-type heterogeneous and asymmetrical traffic conditions. On call-level analysis, the mathematical model based on queueing theory performs multi-dimensional operations using random vectors or matrices to consider multiple types of traffic and also deal with asymmetrical up- and down-direction transmissions separately. Employing the mathematical model, we obtain rail blocking probability for each type of traffic and also the optimum switching-point with the smallest call flocking probability. And on packet-level analysis, employing a non-prioritized queueing scheme between circuit and packet calls, we solve 2-dimensional random vector problem composed of the queue length for packets and the number of circuit calls being served. Finally, packet-level performance is analyzed in terms of the packet loss probability and the buffer size required under mixed-traffic conditions of multiple types of circuit and packet calls.

• Economic and Environmental Geology
• /
• v.53 no.6
• /
• pp.703-713
• /
• 2020

In many geological fields, including hydrogeology, it is of great importance to determine the heterogeneity of the subsurface media. This study briefly introduces the concept and theory of the method that can estimate the hydraulic properties of the media constituting the aquifer, which was recently introduced by Park (2020). After the introduction, the method was applied to the synthetic aquifer to demonstrate the practicality, from which various implications were drawn. The introduced technique uses a global optimization technique called the covariance matrix adaptation evolution strategy (CMA-ES). Conceptually, it is a methodology to characterize the aquifer heterogeneity by assimilating the groundwater level time-series data due to the imposed hydraulic stress. As a result of applying the developed technique to estimate the hydraulic conductivity of a hypothetical aquifer, it was confirmed that a total of 40000 unknown values were estimated in an affordable computational time. In addition, the results of the estimates showed a close numerical and structural similarity to the reference hydraulic conductivity field, confirming that the quality of the estimation by the proposed method is high. In this study, the developed method was applied to a limited case, but it is expected that it can be applied to a wider variety of cases through additional development of the method. The development technique has the potential to be applied not only to the field of hydrogeology, but also to various fields of geology and geophysics. Further development of the method is currently underway.

• Geotechnical Engineering
• /
• v.13 no.6
• /
• pp.147-164
• /
• 1997

This study has the assumption that scale effects in rock mass properties are atrributed to the discontinuous and inhomogeneous nature of rock masses. In order to escape the general equivalent material approach applied to the concept of representative volume element, this study presents the new method considering irregular i oink geometry and arbitrary numbers of i oink and arbitrary joint orientations. Based on the theoretical approach, this theory is applied to a real engineering project. Showing the property variations with size of rock mass element, various numerical experiments about scale effect are conducted. Particularly, to prove the adequacy of the verification process in scale effect with nomerical method, and to investigate the detailed source of scale effect, 4 models with increas ins number of joints are tested. On the basis of the experimental results, the test results of scale effects in 3-D rock mass are presented. From these experiments the effects of the mechanical properties of rock joints on the scale effects in rock mass strength and elastic constants are discussed. To verify the mechanism of scale effects in jointed rock mass, two models with different j oink geometries are studied.

• The Journal of Korean Society for Radiation Therapy
• /
• v.25 no.2
• /
• pp.137-143
• /
• 2013

Purpose: We are to find out the difference of calculated dose of treatment planning system (TPS) and measured dose in case of inhomogeneous organ structure. Materials and Methods: Inhomogeneous phantom is made with solid water phantom and cork plate. CT image of inhomogeneous phantom is acquired. Treatment plan is made with TPS (Pinnacle3 9.2. Royal Philips Electronics, Netherlands) and calculated dose of point of interest is acquired. Treatment plan was delivered in the inhomogeneous phantom by ARTISTE (Siemens AG, Germany) measured dose of each point of interest is obtained with Gafchromic EBT2 film (International Specialty Products, US) in the gap between solid water phantom or cork plate. To simulate lung cancer radiation treatment, artificial tumor target of paraffin is inserted in the cork volume of inhomogeneous phantom. Calculated dose and measured dose are acquired as above. Results: In case of inhomogeneous phantom experiment, dose difference of calculated dose and measured dose is about -8.5% at solid water phantom-cork gap and about -7% lower in measured dose at cork-solid water phantom gap. In case of inhomogeneous phantom inserted paraffin target experiment, dose difference is about 5% lower in measured dose at cork-paraffin gap. There is no significant difference at same material gap in both experiments. Conclusion: Radiation dose at the gap between two organs with different electron density is significantly lower than calculated dose with TPS. Therefore, we must be aware of dose calculation error in TPS and great care is suggested in case of radiation treatment planning on inhomogeneous organ structure.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Journal of Energy Engineering
• /
• v.22 no.2
• /
• pp.156-168
• /
• 2013

Drilling expedition #1 in 2007 and drilling expedition #2 in 2010 were performed for gas hydrate resources evaluation and optimal site selection of pilot test in Ulleung basin, East Sea, Korea. This study presents to build the 3D spatial distribution models using the estimated sedimentary facies, porosity, and gas hydrate saturation derived by well logs and core analysis data from UBGH1-4, UBGH1-9, UBGH1-10, UBGH1-14, UBGH2-2-1, UBGH2-2-2, UBGH2-6, UBGH2-9, UBGH2-10 and UBGH2-11. The objective of 3D spatial distribution modeling is to build a geological representation of the gas hydrate-bearing sediment that honors the heterogeneity in 3D grid scale. The facies modeling is populating sedimentary facies into a geological grid using sequential indicator simulation. The porosity and gas hydrate saturation modeling used sequential Gaussian simulation to populate properties stochastically into grid cells.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.40 no.3
• /
• pp.281-288
• /
• 2016

The ballistic performance data of composite materials is distributed due to material inhomogeneity. In this paper, the uncertainty in residual velocity is obtained experimentally, and a method of predicting it is established numerically for the high-speed impact of a bullet into laminated composites. First, the failure strain distribution was obtained by conducting a tensile test using 10 specimens. Next, a ballistic impact test was carried out for the impact of a fragment-simulating projectile (FSP) bullet with 4ply ([0/90]s) and 8ply ([0/90/0/90]s) glass fiber reinforced plastic (GFRP) plates. Eighteen shots were made at the same impact velocity and the residual velocities were obtained. Finally, simulations were conducted to predict the residual velocities by using the failure strain distributions that were obtained from the tensile test. For this simulation, two impact velocities were chosen at 411.7m/s (4ply) and 592.5m/s (8ply). The simulation results show that the predicted residual velocities are in close agreement with test results. Additionally, the modeling of a composite plate with layered solid elements requires less calculation time than modeling with solid elements.