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

4. Model calculations. for tl1.r,ee--d~nei~sioi~a1. structures - - 4.1. Introduction - In the three-dimensional case the TE- and TM-mode become mixed and cannot longer be treated separately. Now the differential equation for a vector field instead of a scalar field is to be solved. In numerical solutions questions of storage and computer time become important. Assume as example that in approach A a basic domain with 20 cells in each direction is chosen. In this case, only the storage of the electric field vector would require 48000 locations. For an iterative improvement of one field component at least 0.0005 sec are needed for each cell. This yields 12 sec for a complete itera- tion, and 20 min for 100 iterations. This appears to be the ].east time required for this model. Hence methods for a reduction of com- puter time and.storage are particularly appreciated in this case. The equation to be solved is cur1213 - (r) - -t k2 (r) E (r) = - i ~ l(g) i ~ ~ ~ --- where k2.(r) = i~p~c(g). - 1 (r) is the source current density. -e - After the splitting where E is that solution of -11 curlZE (r) + k2(r)F, (r) = -impo& -11 - n--n- (4.3) which vanishes at infinity, we obtain for the anomalous field the two alternative equations c ~ r l . + ~ k2E ~ = - k2E -a n-a a Eq. (4.4a) is the starting point for the volume integral or integral equation approach, Eq. (4.4b) is the point or-' starting for the sur- face integral approach.

4.2. In3ral - --- equation method Let -1 G.(r --a [r), -- i = 1,2,3 be a solution of cur12c. .-I (r --0 lr)+k2 - n (r)G. - --I (r --0 Ir) ,. A - = zi6(g-%), (4.5) vanishing at infinity. Here, the x. are unit vectors along the A h -l h A A A cartesian coordinate axes: x =x, x = -1 - -2 1, z3 = - 2. 14u'iul&31y (4.5) by E (r) and (4.4a) by G. -a - (r I?:) and integrate the difference with re- -1 -0 - spect to - r over the whole space. Green's vector theorem where d-r i.s a volume element, dA a surface element and - 1 the out- ward normal vector, yields since and G vanish at infinity. After combining all three compo-- -i nents and introducing - E instead of E , the vector integral equation -a 1 -gtg)d~ (4.7) -- is obtained. Here is Green's tensor being defined as wj(rJ % --O 3 3 A h A t(r Ir)= C x.G. (r [r)= Z G.. (r .Ir)x: x. - -1-1 -0 - ir j=1 13 -0 -- -1 -3 i= I (using dyadic notation). The tensor elements G.. admit a simple = 7 physical interpretation: Gij (r _o [r) - is the j-th electric field coin- ponent of an oscillating electric di-pole of unit moment pointlng in x.-direction, placed in the normal conductivity structure at &; 1 -- the point of observation is - r. Note that the first subscript and argument refer to the source, the second subscript and argument to the receiver. Because of the fundamental reci.p~:ocity in electro- mag~~etl.sm, source and observer parameters axe interchangeable, i.e. For a proof replace in (4.5) - r by - r', write an analogous equation for G.(rlrl), multiply cross-wise by G and Gi, integrate the -3 - - - j - difference wit11 respect to - r' over the whole space, and obtain (4.9) on using (4 .G) . Due to (4.9) , the equation (4.7) is alternatj.vely written

4.2. In3ral - --- equation method<br />

Let -1 G.(r --a [r), -- i = 1,2,3 be a solution of<br />

cur12c. .-I (r --0 lr)+k2 - n (r)G. - --I (r --0 Ir)<br />

,.<br />

A<br />

- = zi6(g-%), (4.5)<br />

vanishing at infinity. Here, the x. are unit vectors along the<br />

A h -l h A A A<br />

cartesian coordinate axes: x =x, x =<br />

-1 - -2 1, z3 = - 2. 14u'iul&31y (4.5) by<br />

E (r) and (4.4a) by G.<br />

-a - (r I?:) and integrate the difference with re-<br />

-1 -0 -<br />

spect to - r over the whole space. Green's vector theorem<br />

where d-r i.s a volume element, dA a surface element and - 1 the out-<br />

ward normal vector, yields<br />

since and G vanish at infinity. After combining all three compo--<br />

-i<br />

nents and introducing - E instead of E , the vector integral equation<br />

-a<br />

1 -gtg)d~ (4.7)<br />

--<br />

is obtained. Here is Green's tensor being defined as<br />

wj(rJ<br />

% --O<br />

3 3 A h<br />

A<br />

t(r Ir)= C x.G. (r [r)= Z G.. (r .Ir)x: x.<br />

- -1-1 -0 - ir j=1<br />

13 -0 -- -1 -3<br />

i= I<br />

(using dyadic notation). The tensor elements G.. admit a simple<br />

= 7<br />

physical interpretation: Gij (r _o [r) - is the j-th electric field coin-<br />

ponent of an oscillating electric di-pole of unit moment pointlng<br />

in x.-direction, placed in the normal conductivity structure at &;<br />

1 --<br />

the point of observation is - r. Note that the first subscript and<br />

argument refer to the source, the second subscript and argument to<br />

the receiver. Because of the fundamental reci.p~:ocity in electro-<br />

mag~~etl.sm, source and observer parameters axe interchangeable, i.e.<br />

For a proof replace in (4.5) - r by - r', write an analogous equation<br />

for G.(rlrl), multiply cross-wise by G and Gi, integrate the<br />

-3 - - - j -<br />

difference wit11 respect to - r' over the whole space, and obtain (4.9)<br />

on using (4 .G) . Due to (4.9) , the equation (4.7) is alternatj.vely<br />

written

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

Saved successfully!

Ooh no, something went wrong!