• 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.

• Economic and Environmental Geology
• /
• v.27 no.2
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.

• The Korean Journal of Rheology
• /
• v.11 no.1
• /
• pp.18-27
• /
• 1999

Predicting the deformation behaviors of sheets in thermoforming processes has been a daunting challenge due to the strong nonlinearities arising from very large deformations, mold-polymer contact condition and hyperelasticity constitutive equations. Nonlinear numerical analysis is always required to face this challenge especially for realistic processing conditions. In this study a 3-D algorithm and the membrane approximation are developed for thermoforming processes. The constitutive equation is expressed in terms of the 2nd Piola-Kirchhoff stress tensor and the Cauchy-Green deformation tensor. The 2-term Mooney-Rivlin model is used for the material model equation. The algorithm is established by the finite element formulation employing the total Lagrangian coordinate. The deformation behavior and the stress distribution results of 3-D algorithm with various point boundary conditions are compared to those of the membrane approximation algorithm. Also, the slip boundary condition and the no-slip boundary condition are applied for the systems that have molds. Finally, the effect of sheet temperatures on the final thickness distribution is investigated for the ABS material.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.491-498
• /
• 1994

The integral equation method is used for magnetotelluric (MT) modeling of a finite inhomogeneity in a two-layered earth. An integral equation relates the incident plane-wave field and the scattering currents in the three-dimensional (3-D) inhomogeneity through the electric tensor Green's function appropriate to a layered earth. This paper describes an effect of overburden and basement on MT responses of 3-D body. The effect of overburden is to reduce the detectability of target, and the reduction of detectability is more apparent for conductive overburden than for resistive one. The effect of basement, on the other hand, may enhance the anomaly due to 3-D body in the upper layer. In case of the resistive basement current perturbations about the body tend to be confined to the more conductive upper layer.