• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.13-18
• /
• 2000

We describe the concept of the singularity in the integral equation of electromagnetic (EM) modeling in comparison with that in the integral representation of electric fields in EM theory, which would clarify the singular integral problems of the Green's tensor. We have also derived and classified the singular integrals of the Green's tensors in 3-D, 2.5-D and 2-D as well as in the thin sheet integral equations of the EM scattering problem, which have the most important effect on the accuracy of the numerical solution of the problems.

• Geophysics and Geophysical Exploration
• /
• v.26 no.1
• /
• pp.1-7
• /
• 2023

This study derives the expressions of vector gravity and gravity gradient tensor due to an elliptical cylinder. The vector gravity for an arbitrary three-dimensional (3D) body is obtained by differentiating the gravitational potential, including the triple integral, according to the shape of the body in each axis direction. The vector gravity of the 3D body with axial symmetry is integrated along the axial direction and reduced to a double integral. The complex Green's theorem using complex conjugates subsequently converts the double integral into a one-dimensional (1D) closed-line integral. Finally, the vector gravity due to the elliptical cylinder is derived using 1D numerical integration by parameterizing a boundary of the elliptical cross-section as a closed line. Similarly, the gravity gradient tensor due to the elliptical cylinder is second-order differentiated from the gravitational potential, including the triple integral, and integrated along the vertical axis direction reducing it to a double integral. Consequently, all the components of the gravity gradient tensor due to an elliptical cylinder are derived using complex Green's theorem as used in the case of vector gravity.

• Geophysics and Geophysical Exploration
• /
• v.27 no.1
• /
• pp.51-56
• /
• 2024

In this paper, the vector gravity and gravity gradient tensor of an elliptical disk are derived. The vector gravity of an elliptical disk is defined by differentiating the gravitational potential due to the elliptical disk expressed by a double integral with respect to each axial direction. The vector gravity defined by the double integral is then transformed into a line integral of a closed curve along the elliptical disk boundary using the complex Green's theorem. Finally, vector gravity due to the elliptical disk is derived by 1D parametric numerical integration along the elliptical disk boundary. The xz, yz, zz components of the gravity gradient tensor due to the elliptical disk are obtained by differentiating the vector gravity with respect to vertical direction. The xx, yy, xy components are derived by differentiating the horizontal components of the vector gravity in the form of a double integral with respect to horizontal directions and then using the complex Green's theorem.

• Economic and Environmental Geology
• /
• v.27 no.1
• /
• pp.41-47
• /
• 1994

A fast Hankel transform (FHT) algorithm (Anderson, 1982) is applied to numerical evaluation of many Green's tensor integrals encountered in three-dimensional electromagnetic modeling using integral equations. Efficient computation of Hankel transforms is obtained by a combination of related and lagged convolutions which are available in the FHT. We express Green's tensor integrals for a layered half-space, and rewrite those to a form of related functions so that the FHT can be applied in an efficient manner. By use of the FHT, a complete or full matrix of the related Hankel transform can be rapidly and accurately calculated for about the same computation time as would be required for a single direct convolution. Computing time for a five-layer half-space shows that the FHT is about 117 and 4 times faster than conventional direct and multiple lagged convolution methods, respectively.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.1 no.3 s.3
• /
• pp.123-136
• /
• 2001

The purpose of application of moment tensor inversion method is to determine source parameters, such as, focal mechanism, seismic moment and source depth. This paper presents results of focal mechanism solutions of 14 recent events with magnitudes ranging from M3.3 to M4.8 by using moment tensor inversion method called TDMT_INV. The strike of faults is in the direction of NE-SW and NW-SE with the movement of strike-slip or strike-slip of minor reverse component. The compressional axis of the stress field is predominantly E-W or ENE-WSW except for some faults, which are distributed at Ryongnam Massif and Wonsan, they have a compressional axis of N-S or NNW-SSE.

• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.507-518
• /
• 2012

It is now widely accepted that the resulting displacement field within elastic, inhomogeneous, anisotropic solids subjected to equipartitioned, uniform illumination from uncorrelated sources, has intensities that follow diffusion-like equations. Typically, coda waves are invoked to illustrate this concept. These waves arrive later as a consequence of multiple scattering and appear at "the tail" (coda, in Latin) of seismograms and are usually considered an example of diffuse field. It has been demonstrated that the average correlations of motions within a diffuse field, in frequency domain, is proportional to the imaginary part of Green's function tensor. If only one station is available, the average autocorrelation is equal to the average squared amplitudes or the average power spectrum and this gives the Green's function at the source itself. Several works address this point from theoretical and experimental point of view. However, a complete and explicit analytical description is lacking. In this work we study analytically some properties of the Green's function, specifically the imaginary part of Green's function for 2D antiplane problems. This choice is guided by the fact that these scalar problems have a closed analytical solution (Kausel 2006). We assume the diffusiveness of the field and explore its analytical consequences.

• Geophysics and Geophysical Exploration
• /
• v.5 no.1
• /
• pp.64-67
• /
• 2002

• Geophysics and Geophysical Exploration
• /
• v.2 no.2
• /
• pp.86-95
• /
• 1999

This paper presents an efficient three-dimensional (3-D) modeling algorithm using the extended approximation to an electric field integral equation. Numerical evaluations of Green's tensor integral are performed in the spatial wavenumber domain. This approach makes it possible to reduce computing time, to handle smoothly varying conductivity model and to remove singularity problems encountered in the integration of Green's tensor at a source point. The responses obtained by 3-D modeling algorithm developed in this study are compared with those by the full integral equation for a thin-sheet EM scattering. The extensive analyses on the performance of modeling algorithm are made with the conductivity contrasts and source frequencies. These results show that the modeling algorithm are accurate up to the conductivity contrast of 1:16 and the frequency range of 100 Hz-100 kHz. The extended Born approximation, however, may produce inaccurate results for some source and model configurations in which the electric field is discontinuous across the conductivity boundary. We performed the modeling of a composite model of which conductivity varies continuously and this shows the modeling algorithm developed in this study is efficient for 3-D EM modeling. For a cross-hole source-receiver configuration a composite model of which conductivity varies continuously can be successfully simulated using this algorithm.

• The Journal of Engineering Geology
• /
• v.17 no.4
• /
• pp.673-680
• /
• 2007

The seismic source parameters of the Odaesan earthquake on 20 January 2007, including focal depth, focal mechanism, magnitude, and source characteristics, are analysed using seismic moment tensor inversion. The Green's function for different 3 crust models representing the southern Korean Peninsula are used. Final results show that the event, considering 6 seismic moment tensor elements, is caused by the typical strike slip fault with nearly NNE strike. The focal depth is estimated to be about 11km and 6 seismic moment tensor elements with 7.2% CLVD value shows typical double couple seismic source. The consistent characteristics of the strike and epicenter of the event with Odaesan fault imply that Odaesan earthquake is mainly caused by movement of the Odaesan fault.

• Advances in aircraft and spacecraft science
• /
• v.3 no.4
• /
• pp.379-397
• /
• 2016

In this manuscript, free vibrations of a unidirectional composite orthotropic Timoshenko beam based on finite strain have been studied. Using Green-Lagrange strain tensor and comprising all of the nonlinear terms of the tensor and also applying Hamilton's principle, equations of motion and boundary conditions of the beam are obtained. Using separation method in single-harmonic state, time and locative variables are separated from each other and finally, the equations of motion and boundary conditions are gained according to locative variable. To solve the equations, generalized differential quadrature method (GDQM) is applied and then, deflection and cross-section rotation of the beam in linear and nonlinear states are drawn and compared with each other. Also, frequencies of carbon/epoxy and glass/epoxy composite beams for different boundary conditions on the basis of the finite strain are calculated. The calculated frequencies of the nonlinear free vibration of the beam utilizing finite strain assumption for various geometries have been compared to von Karman one.