Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

mtnet.dias.ie
from mtnet.dias.ie More from this publisher
04.08.2013 Views

assuming that the TM magnetic source field j.s due to a uniforrn sheet current at height z = -h, h > 0. T1lj.s assun~ptj-on, however, is immaterial for the following. In the sequel all field quantities are split into a normal and anoinalous part, denoted by the subscripts "nu and "a" , respectivel-y. The normal part refers to a one-dimensional conductivity structure. Let u(y,zl = un(z) + U,(Y,Z) (3.7) H(y,z) = Hn(z) + Ha(y,z) E and H are defined as solutions of the equations n n (3.3b) e AEn = k2 E + j (3.10a) n 11 d l d -(- - 1-1 ) = H z > Or H (0) = dz ,?2 dz n n ' - n vanishing for z + m. In virtue of (3.5a,b), (3.9arb), and (3.103,b) Ea and H satisfy a 1 d 1. 1 dIln div (- gradH ) = H a - - - , z > 0 k a k * k2 dz - n If the anomalous domain is of fin2te extent, E has to vanish uni- a forinly at infinity. Under the same condition Ha has to vanish uni- formly in the lower half-space. At z=o H is zero. a If the anomalous domain is of infi.nite extent in horizontal direc- tion, we can demand only that Ear H +O for z + m. a For a numerical solution of (3.11 a) the following three clioices of a basic dornain are possible (boundaries hatched). In approach A, (3.11 a) is solved by Finite differences subject to the boundary condition Ea=O or better subject to an inpedance boundary collclition (below). In approach B (3.11a) is solved by finite differences only in the anomal.ous slab. At the ho~izontal boundaries boundary conditions involicity the normal structure above and below the slab are applied. I approach C (3.11~~) is re- duced to an integral equation over the anon~alous domain. These approaches will now be discussed in details.

i I Earth Anomalous domain 3.2. Air half-space and conductor. as basic domain (Finite difference method) For the TE- and TM-mode we have to solve the differential equa- tions (3.Sa,b), i.e. with the boundary condltion that the differences E =E-E and a n Ha=H-HI> vanish' at infinity. E has to be computed for any given n two-dimensional external source field al.ong the lines of Sec.2.3. The H -field belongs to a uniform external magnetic field. n In the finite difference method, the differential operators in (3.12a,b) are reduced to finite differences. For simplicity a d square grid.with grid ~11th h is assumed. Consider the following configuration of a nodal point 0 and its four neighbours:

i<br />

I Earth<br />

Anomalous<br />

domain<br />

3.2. Air half-space and conductor. as basic domain<br />

(Finite difference method)<br />

For the TE- and TM-mode we have to solve the differential equa-<br />

tions (3.Sa,b), i.e.<br />

with the boundary condltion that the differences E =E-E and<br />

a n<br />

Ha=H-HI> vanish' at infinity. E has to be computed for any given<br />

n<br />

two-dimensional external source field al.ong the lines of Sec.2.3.<br />

The H -field belongs to a uniform external magnetic field.<br />

n<br />

In the finite difference method, the differential operators in<br />

(3.12a,b) are reduced to finite differences. For simplicity a<br />

d<br />

square grid.with grid ~11th h is assumed. Consider the following<br />

configuration of a nodal point 0 and its four neighbours:

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!