scholarly journals Electromagnetic modeling of three‐dimensional bodies in layered earths using integral equations

Geophysics ◽  
1984 ◽  
Vol 49 (1) ◽  
pp. 60-74 ◽  
Author(s):  
Philip E. Wannamaker ◽  
Gerald W. Hohmann ◽  
William A. SanFilipo
Keyword(s):  
Integral Equations ◽  
Green's Functions ◽  
Green’S Functions ◽  
Three Dimensional ◽  
Computation Time ◽  
Vector Current ◽  
Hankel Transform ◽  
Transverse Magnetic ◽  
Electromagnetic Modeling ◽  
Method Of Integral Equations

We have developed an algorithm based on the method of integral equations to simulate the electromagnetic responses of three‐dimensional bodies in layered earths. The inhomogeneities are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric tensor Green’s function appropriate to a layered earth, and it is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic tensor Green’s functions over the scattering currents. Efficient evaluation of the tensor Green’s functions is a major consideration in reducing computation time. We find that tabulation and interpolation of the six electric and five magnetic Hankel transforms defining the secondary Green’s functions is preferable to any direct Hankel transform calculation using linear filters. A comparison of responses over elongate three‐dimensional (3-D) bodies with responses over two‐dimensional (2-D) bodies of identical cross‐section using plane wave incident fields is the only check available on our solution. Agreement is excellent; however, the length that a 3-D body must have before departures between 2-D transverse electric and corresponding 3-D signatures are insignificant depends strongly on the layering. The 2-D transverse magnetic and corresponding 3-D calculations agree closely regardless of the layered host.

THREE‐DIMENSIONAL INDUCED POLARIZATION AND ELECTROMAGNETIC MODELING

Geophysics ◽  
1975 ◽  
Vol 40 (2) ◽  
pp. 309-324 ◽  
Author(s):  
Gerald W. Hohmann
Keyword(s):  
Integral Equation ◽  
Green's Functions ◽  
Green’S Functions ◽  
Three Dimensional ◽  
Induced Polarization ◽  
The Body ◽  
Scale Model ◽  
Two Dimensional ◽  
Electromagnetic Modeling ◽  
The Earth

The induced polarization (IP) and electromagnetic (EM) responses of a three‐dimensional body in the earth can be calculated using an integral equation solution. The problem is formulated by replacing the body by a volume of polarization or scattering current. The integral equation is reduced to a matrix equation, which is solved numerically for the electric field in the body. Then the electric and magnetic fields outside the inhomogeneity can be found by integrating the appropriate dyadic Green’s functions over the scattering current. Because half‐space Green’s functions are used, it is only necessary to solve for scattering currents in the body—not throughout the earth. Numerical results for a number of practical cases show, for example, that for moderate conductivity contrasts the dipole‐dipole IP response of a body five units in strike length approximates that of a two‐dimensional body. Moving an IP line off the center of a body produces an effect similar to that of increasing the depth. IP response varies significantly with conductivity contrast; the peak response occurs at higher contrasts for two‐dimensional bodies than for bodies of limited length. Very conductive bodies can produce negative IP response due to EM induction. An electrically polarizable body produces a small magnetic field, so that it is possible to measure IP with a sensitive magnetometer. Calculations show that horizontal loop EM response is enhanced when the background resistivity in the earth is reduced, thus confirming scale model results.

Electromagnetic modeling of 3-D structures by the method of system iteration using integral equations

Geophysics ◽  
1992 ◽  
Vol 57 (12) ◽  
pp. 1556-1561 ◽  
Author(s):  
Zonghou Xiong
Keyword(s):  
Integral Equations ◽  
Three Dimensional ◽  
Computation Time ◽  
Iterative Solution ◽  
Matrix Inversion ◽  
Electromagnetic Modeling ◽  
Direct Inversion ◽  
Conductivity Structure ◽  
The Matrix ◽  
Earth Conductivity

A new approach for electromagnetic modeling of three‐dimensional (3-D) earth conductivity structures using integral equations is introduced. A conductivity structure is divided into many substructures and the integral equation governing the scattering currents within a substructure is solved by a direct matrix inversion. The influence of all other substructures are treated as external excitations and the solution for the whole structure is then found iteratively. This is mathematically equivalent to partitioning the scattering matrix into many block submatrices and solving the whole system by a block iterative method. This method reduces computer memory requirements since only one submatrix at a time needs to be stored. The diagonal submatrices that require direct inversion are defined by local scatterers only and thus are generally better conditioned than the matrix for the whole structure. The block iterative solution requires much less computation time than direct matrix inversion or conventional point iterative methods as the convergence depends on the number of the submatrices, not on the total number of unknowns in the solution. As the submatrices are independent of each other, this method is suitable for parallel processing.

Three-Dimensional Green’s Functions in Anisotropic Elastic Bimaterials With Imperfect Interfaces

2003 ◽  
Vol 70 (2) ◽  
pp. 180-190 ◽  
Author(s):  
E. Pan
Keyword(s):  
Green's Functions ◽  
Green’S Functions ◽  
Three Dimensional ◽  
Imperfect Interface ◽  
Boundary Integral ◽  
Stress Fields ◽  
Superposition Method ◽  
Interface Conditions ◽  
Complementary Part ◽  
Anisotropic Elastic

In this paper, three-dimensional Green’s functions in anisotropic elastic bimaterials with imperfect interface conditions are derived based on the extended Stroh formalism and the Mindlin’s superposition method. Four different interface models are considered: perfect-bond, smooth-bond, dislocation-like, and force-like. While the first one is for a perfect interface, other three models are for imperfect ones. By introducing certain modified eigenmatrices, it is shown that the bimaterial Green’s functions for the three imperfect interface conditions have mathematically similar concise expressions as those for the perfect-bond interface. That is, the physical-domain bimaterial Green’s functions can be obtained as a sum of a homogeneous full-space Green’s function in an explicit form and a complementary part in terms of simple line-integrals over [0,π] suitable for standard numerical integration. Furthermore, the corresponding two-dimensional bimaterial Green’s functions have been also derived analytically for the three imperfect interface conditions. Based on the bimaterial Green’s functions, the effects of different interface conditions on the displacement and stress fields are discussed. It is shown that only the complementary part of the solution contributes to the difference of the displacement and stress fields due to different interface conditions. Numerical examples are given for the Green’s functions in the bimaterials made of two anisotropic half-spaces. It is observed that different interface conditions can produce substantially different results for some Green’s stress components in the vicinity of the interface, which should be of great interest to the design of interface. Finally, we remark that these bimaterial Green’s functions can be implemented into the boundary integral formulation for the analysis of layered structures where imperfect bond may exist.

Digital Modeling of Electromagnetic Processes in Technological Induction Devices

2021 ◽  
Vol 7 (7) ◽  
pp. 26-32
Author(s):  
Viktor B. DEMIDOVICH ◽  
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Integral Equations ◽  
Transverse Magnetic ◽  
Combined Method ◽  
Method Of Integral Equations ◽  
Leading Role ◽  
Casting Aluminum ◽  
Point To Point ◽  
Element Method

Development of an electrical calculation method plays the leading role in simulating induction devices. In modeling electrical devices and complexes, it is often necessary to simultaneously solve both chain and field problems, i.e., to deal with both lumped and distributed parameters. The article considers the method of integral equations for induction systems with non-magnetic and ferromagnetic loading, which is based on the theory of long-range action. The method’s key statement is that the field at any point is determined as the sum of the fields produced by all sources, including primary and secondary ones. Another finite element method is based on the theory of short-range action, which describes the electromagnetic wave propagation from point to point, its refraction and reflection at the boundaries of media. The article substantiates the development of a combined method based on using the method of integral equations for calculating the input parameters of inductors (an external problem) and the finite element method for calculating the field distribution in the load (an internal problem). The combined method has well proven itself in modeling induction heating and melting of metals and oxides, heating a tape in a transverse magnetic field, induction plasmatrons, and casting aluminum into an electromagnetic crystallizer.

scholarly journals On the unique solvability of a class of modified boundary integral equations

1984 ◽  
Vol 27 (3) ◽  
pp. 303-311 ◽  
Author(s):  
R. E. Kleinman ◽  
G. F. Roach
Keyword(s):  
Helmholtz Equation ◽  
Integral Equations ◽  
Unique Solvability ◽  
Green's Functions ◽  
Green’S Functions ◽  
Boundary Integral Equations ◽  
Boundary Integral ◽  
Transmission Problem

In a recent paper the authors considered the transmission problem for the Helmholtz equation by using a reformulation of the problem in terms of a pair of coupled boundary integral equations with modified Green's functions as kernels. In this note we settle the question of the unique solvability of these modified boundary integral equations.

Sign in / Sign up

Export Citation Format

Share Document