• Journal of Industrial Technology
• /
• v.24 no.B
• /
• pp.163-170
• /
• 2004

In this study, First, the results of travel-time inversion (first arrival inversion using the travel-time of the first arrival) were compared with those of full-wave inversion for numerical data. Numerical experiments to find key parameters other than initial velocity model showed that the frequency of source has a great effect on the result of full-wave inversion. Finally, this research presented the corrected full-wave inversion applying the correction term to the final result of full-wave inversion. The corrected full-wave inversion depicted cavities inside concretes even when the inversion started with 20% error in an initial velocity model for cavities. However, full-wave inversion did not reveal cavities.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1201-1206
• /
• 2004

This paper deals with the numerical analysis of natural convection flow induced by the density inversion effect of water inside horizontal pipe. The numerical method is based on SIMPLE/PWIM in general coordinate for its wide applicabilities. The numerical tool was validated through the comparison with the previous results concerning the density inversion effect of water It is shown that the developed numerical tool could predict the flow pattern and the heat transfer phenomena qualitatively And it is also found that the density inversion effect of water has significant effects on the flow pattern.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• v.3B no.2
• /
• pp.67-73
• /
• 2003

In this paper, the impedance of a piezoelectric transducer is calculated using the three-dimensional finite element method. The validity of numerical routine is confirmed experimentally. Using this numerical routine, the effects of material coefficients on piezoelectric actuators characteristics are analyzed. The material constants, which make significant effects, are selected and the relations between material constants are studied. Using these processes, three variables of material constants for a piezoelectric transducer are selected and the design sensitivity method is adopted as an inversion scheme. The validity of the inversion scheme for a piezoelectric transducer is confirmed by applying the proposed method to the sample piezoelectric transducer.

• Kyungpook Mathematical Journal
• /
• v.52 no.4
• /
• pp.383-395
• /
• 2012

In this paper, we use a computational algorithm for the inversion of the one and two-dimensional $\mathcal{L}_2$-transform based on the Bromwich's integral and the Fourier series. The new inversion formula can evaluate the inverse of the $\mathcal{L}_2$-transform with considerable accuracy over a wide range of values of the independent variable and can be devised for the functions which are not Laplace transformable and have damping motion in small interval near origin.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed, when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• Geophysics and Geophysical Exploration
• /
• v.6 no.4
• /
• pp.199-206
• /
• 2003

Characteristics of the static shift are discussed by comparing the three-dimensional MT inversion with/without static shift parameterization. The galvanic distortion by small-scale shallow feature often leads severe distortion in inverted resistivity structures. The new inversion algorithm is applied to four numerical data sets contaminated by different amount of static shift. In real field data interpretations, we generally do not have any a-priori information about how much the data contains the static shift. In this study, we developed an algorithm for finding both Lagrangian multiplier for smoothness and the trade-off parameter for static shift, simultaneously in 3-D MT inversion. Applications of this inversion routine for the numerical data sets showed quite reasonable estimation of static shift parameters without any a-priori information. The inversion scheme is successfully applied to all the four data sets, even when the static shift does not obey the Gaussian distribution. Allowing the static shift parameters have non-zero degree of freedom to the inversion, we could get more accurate block resistivities as well as static shifts in the data. When inversion does not consider the static shift as inversion parameters (conventional MT inversion), the block resistivities on the surface are modified considerably to match possible static shift. The inhomogeneous blocks on the surface can generate the static shift at low frequencies. By those mechanisms, the conventional 3-D MT inversion can reconstruct the resistivity structures to some extent in the deeper parts even when moderate static shifts are in the data. As frequency increased, however, the galvanic distortion is not frequency independent any more, and thus the conventional inversion failed to fit the apparent resistivity and phase, especially when strong static shift is added. Even in such case, however, reasonable estimation of block resistivity as well as static shift parameters were obtained by 3-D MT inversion with static shift parameterization.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.4
• /
• pp.155-165
• /
• 2005

Love wave and Rayleigh wave are the major elastic waves belonging to the category of the surface wave. Those waves are used to determine the ground stiffness profile using their dispersion characteristics. The fact that Love wave is not contaminated by P-wave makes Love wave superior to Rayleigh wave and other body waves. Therefore, the information that Love wave carries is more distinct and clearer than that of others. Based on theoretical research, the joint inversion analysis that uses the dispersion information of both Love and Rayleigh wave was proposed. Numerical analysis, theoretical model test, and field test were performed to verify the joint inversion analysis. Results from 2D, 3D finite element analysis were compared with those from the transfer matrix method in the numerical analysis. On the other hand, the difference of results from each inversion analysis was investigated in the theoretical model analysis. Finally, practical applicability of the joint inversion analysis was verified by performing field test. As a result, it is confirmed that considering dispersion information of each wave simultaneously prevents excessive divergence and improves accuracy.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.1
• /
• pp.41-53
• /
• 2008

Polynomials are one of most important and widely used numerical tools in dealing with a smooth function on a bounded domain and trigonometric functions work for smooth periodic functions. However, they are not the best choice if a function has a bounded support in space and in frequency domain. The Prolate Spheroidal wave function (PSWF) of order zero has been known as a best candidate as a basis for band-limited functions. In this paper, we review some basic properties of PSWFs defined as eigenfunctions of bounded Fourier transformation. We also propose numerical inversion schemes based on PSWF and present some numerical examples to show their feasibilities as signal processing tools.

• Geomechanics and Engineering
• /
• v.35 no.2
• /
• pp.109-119
• /
• 2023

Ground subsidence in urban areas due to excessive development and degraded underground facilities is a serious problem. Geophysical surveys have been conducted to estimate the distribution and scale of cavities and subsidence. In this study, electrical resistivity tomography (ERT) was performed near an area of road subsidence in an urban area. The subsidence arose due to groundwater leakage that carried soil into a neighboring excavation site. The ERT survey line was located between the main subsidence area and an excavation site. Because ERT data are affected by rapid topographic changes and surrounding structures, the influence of the excavation site on the data was analyzed through field-scale numerical modeling. The effect of an excavation should be considered when interpreting ERT data because it can lead to wrong anomalous results. A method for performing 2D inversion after correcting resistivity data for the effect of the excavation site was proposed. This method was initially tested using a field-scale numerical model that included the excavation site and subsurface anomaly, which was a loosened zone, and was then applied to field data. In addition, ERT data were interpreted using an existing in-house 3D algorithm, which considered the effect of excavation sites. The inversion results demonstrated that conductive anomalies in the loosened zone were greater compared to the inversion that did not consider the effects of excavation.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.129-135
• /
• 2012

The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.