• Proceedings of the IEEK Conference
• /
• 2002.07a
• /
• pp.650-653
• /
• 2002

The increaseing speeds are accompanied by decreases in pulse rise and fall time in VLSI circuits. These accenturate the high frequency spectral contents of the signals and cause the frequency dependent loss of the conductors which interconnect the various sub-circuits composing of VLSI circuit. The lossy effect is approximated by the square root of frequency dependence of the per unit length resistance. In the practical applications, several problems may arise along with this approximation, so we extend our investigation of the lossy effect by numerical Laplace inversion method.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.25 no.3
• /
• pp.252-261
• /
• 1988

Employed the new invese scattering scheme based on the moment mehtod, which was presented in the Part I of these companion papers, numerical simulations are performed to investigate the effect of measurement errors and noise contaminating the field scattered from dielectric objects. In order to reduce those effects on the reconstructed permittivity profiles, some techniques such as regularization, iterative matrix inversion, and multiple incidence are applied to this problem.

• Journal of the Korean Data and Information Science Society
• /
• v.10 no.2
• /
• pp.415-428
• /
• 1999

In this paper, we study the saddlepoint approximations for the ratio of two independent sequences of random variables. In Section 2, we review the saddlepoint approximation to the probability density function. In section 3, we derive an saddlepoint approximation formular for the tail probability by following Daniels'(1987) method. In Section 4, we represent a numerical example which shows that the errors are small even for small sample size.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.51-56
• /
• 2008

The waveform inversion for isotropic media has ever been studied since the 1980s, but there has been few studies for anisotropic media. We present a seismic waveform inversion algorithm for 2-D heterogeneous transversely isotropic structures. A cell-based finite difference algorithm for anisotropic media in time domain is adopted. The steepest descent during the non-linear iterative inversion approach is obtained by backpropagating residual errors using a reverse time migration technique. For scaling the gradient of a misfit function, we use the pseudo Hessian matrix which is assumed to neglect the zero-lag auto-correlation terms of impulse responses in the approximate Hessian matrix of the Gauss-Newton method. We demonstrate the use of these waveform inversion algorithm by applying them to a two layer model and the anisotropic Marmousi model data. With numerical examples, we show that it's difficult to converge to the true model when we assumed that anisotropic media are isotropic. Therefore, it is expected that our waveform inversion algorithm for anisotropic media is adequate to interpret real seismic exploration data.

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

The large scale shaker can be employed to measure linear and nonlinear shear moduli of soils in situ as a function of shear strain. The method involves applying dynamic loads on a surface foundation measuring the dynamic response of the soil mass beneath the foundation with embedded instrumentation. This paper focuses on the development of a framework of the inverse analysis for the interpretation of test data to estimate linear and nonlinear shear moduli of soils along with the necessity of the inverse analysis. The suggested framework is based on the nonlinear least squares but it uses two iterative loops to account for the nonlinear behavior of soil that sensors are not located. The validity of the suggested inversion framework is tested through a series of numerical parametric studies. An example use of the suggested inversion framework is also shown. Because the field condition may affect the accuracy of suggested method, it is important to conduct a preliminary inverse analysis to quantify the discrepancy between the estimated modulus and the baseline.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.221-230
• /
• 2021

This paper proposes full waveform inversion (FWI) for estimating the physical properties of a layered half-space. An FWI solution is obtained using a genetic algorithm (GA), which is a well-known global optimization approach. The dynamic responses of a layered half-space subjected to a harmonic vertical disk load are measured and compared with those calculated using the estimated physical properties. The responses are calculated using the thin-layer method, which is accurate and efficient for layered media. Subsequently, a numerical model is constructed for a layered half-space using mid-point integrated finite elements and perfectly matched discrete layers. An objective function of the global optimization problem is defined as the L2-norm of the difference between the observed and estimated responses. A GA is used to minimize the objective function and obtain a solution for the FWI. The accuracy of the proposed approach is applied to various problems involving layered half-spaces. The results verify that the proposed FWI based on a GA is suitable for estimating the material properties of a layered half-space, even when the measured responses include measurement noise.

• Geomechanics and Engineering
• /
• v.36 no.6
• /
• pp.575-587
• /
• 2024

The contact pressure between the surrounding rock and the support is an important indicator of the surrounding rock pressure. There has been a bottleneck in the prediction of contact loads between surrounding rock and primary support in deep-buried mountain tunnels. The main reason is that a reliable method wasn't existed to quantify the contact loads. This study had been taken into account the flexible support role of the primary support, and the fitting curve of surrounding rock deformation for dynamic tunnel construction was proposed. New formulas for the calculation of contact loads between surrounding rock and primary support were obtained by inversion. Comparative analysis of the calculation results with numerical simulation verified the reliability of the calculation method in this study. It can be seen from the analyses that the contact load between surrounding rock and primary support increases, remains unchanged and decreases during acceleration, uniform velocity and deceleration, respectively, and the deformation of the surrounding rock in the acceleration and deceleration stages cannot completely converted into contact loads. The contact loads between surrounding rock and primary support of medium-strength and weak surrounding rock tunnels are generally within 150 kPa and 1 MPa, respectively. For tunnels with weak surrounding rock, advanced support can be installed to reduce the unique release coefficient λ₀ and the value of the constant D, with the purpose of reducing the contact loads between surrounding rock and primary support. Changes in support parameters have a small effect on the contact loads between surrounding rock and primary support, but increase or decrease the safety factor, resulting in a waste of resources or a situation that threatens the safety of the support. The results of this research provide guidance for the prediction of contact loads between surrounding rock and primary support for dynamic tunnel construction.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.2
• /
• pp.262-267
• /
• 2003

Error of the geometric series expansion method for the structural sensitivity analysis is estimated. Although the semi-analytic method has several advantages, accuracy of the method prevents it from practical application. One of the promising remedies is the use of geometric series formula for the matrix inversion. Its result of the sensitivity analysis converges that of the global difference method which is known as reliable one. To reduce computational efforts and to obtain reliable results, it is important to know how many terms need to expand. In this paper, the error formula is presented and Its usefulness is illustrated through numerical experiments.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.10
• /
• pp.3313-3321
• /
• 1996

A fundamental study in cooling and solidification process focused on ice storage was performed, including the interesting phenomena of density inversion, supercooling and dendritic ice. A numerical study was performed for natural convection and ice formation in the cooling and freezing processes with supercooling in a space between double cylinders. When water was cooled under the freezing point by a cooling wall in a cavity, solidification was not started at once, but a subcooled region was formed near the wall. Especially, when the cooling rate was low, subcooled region extended to a wide area. However, after a few minutes, supercooling is released by some triggers. Dendritic ice is suddenly formed within a subcooled region, and a dense ice layer begins to be developed from the cooling wall. Due to the difficulties, most previous studies on solidification process with numerical methods had not treated the supercooling phenomena, i.e. the case considering only the growth of dense ice. In this study, natural convection and ice formation considering existence of supercooling and dendritic ice were analyzed numerically with using finite difference method and boundary fixing method. The results of numerical analysis were well compared with the experimental results.

• Journal of the Korean earth science society
• /
• v.30 no.5
• /
• pp.550-564
• /
• 2009

When magnetotelluric (MT) data are obtained in the vicinity of the coast, the surrounding seas make it difficult to interpret subsurface structure, especially the deep part of the subsurface. We introduce an iterative method to correct the sea effect, based on the previous topographic correction method that removes the distortion due to topographic changes in seafloor MT data. The method first corrects the sea effect in observed MT impedance, and then inverts corrected response in a model space without the sea. Due to mutual coupling between the sea and the subsurface structure, the correction and inversion steps are iterated until the changes in each result become negligible. The method is tested for 1- and 2-D structures using synthetic MT data produced by 3-D forward modeling including surrounding seas. In all cases, the method closely recovers the true structure assumed to generate synthetic responses after a few iterations.