04.08.2013 Views

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

assuming that the TM magnetic source field j.s due to a uniforrn<br />

sheet current at height z = -h, h > 0. T1lj.s assun~ptj-on, however,<br />

is immaterial for the following.<br />

In the sequel all field quantities are split into a normal and<br />

anoinalous part, denoted by the subscripts "nu and "a" , respectivel-y.<br />

The normal part refers to a one-dimensional conductivity structure.<br />

Let<br />

u(y,zl = un(z) + U,(Y,Z) (3.7)<br />

H(y,z) = Hn(z) + Ha(y,z)<br />

E and H are defined as solutions of the equations<br />

n n<br />

(3.3b)<br />

e<br />

AEn = k2 E + j (3.10a)<br />

n 11<br />

d l d<br />

-(- - 1-1 ) = H z > Or H (0) =<br />

dz ,?2 dz n n ' - n<br />

vanishing for z + m.<br />

In virtue of (3.5a,b), (3.9arb), and (3.103,b) Ea and H satisfy<br />

a<br />

1 d 1. 1 dIln<br />

div (- gradH ) = H a - - - , z > 0<br />

k<br />

a<br />

k * k2 dz -<br />

n<br />

If the anomalous domain is of fin2te extent, E has to vanish uni-<br />

a<br />

forinly at infinity. Under the same condition Ha has to vanish uni-<br />

formly in the lower half-space. At z=o H is zero.<br />

a<br />

If the anomalous domain is of infi.nite extent in horizontal direc-<br />

tion, we can demand only that Ear H +O for z + m.<br />

a<br />

For a numerical solution of (3.11 a) the following three clioices of<br />

a basic dornain are possible (boundaries hatched).<br />

In approach A, (3.11 a) is solved by Finite differences subject to<br />

the boundary condition Ea=O or better subject to an inpedance<br />

boundary collclition (below). In approach B (3.11a) is solved by<br />

finite differences only in the anomal.ous slab. At the ho~izontal<br />

boundaries boundary conditions involicity the normal structure<br />

above and below the slab are applied. I approach C (3.11~~) is re-<br />

duced to an integral equation over the anon~alous domain. These<br />

approaches will now be discussed in details.

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

Saved successfully!

Ooh no, something went wrong!