Open AccessOn the finite difference solution of two‐dimensional induction problems
Author(s) -
BrewittTaylor C. R.,
Weaver J. T.
Publication year - 1976
Publication title -
geophysical journal of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0016-8009
DOI - 10.1111/j.1365-246x.1976.tb01280.x
Subject(s) - finite difference , generalization , finite difference method , mathematics , simple (philosophy) , grid , mathematical analysis , polarization (electrochemistry) , position (finance) , numerical analysis , electromagnetic induction , geometry , physics , philosophy , chemistry , epistemology , finance , quantum mechanics , economics , electromagnetic coil
Summary The numerical solution by finite differences of two‐dimensional problems in electromagnetic induction is reexamined with a view to generalizing the method to three‐dimensional models. Previously published work, in which fictitious values were used to derive the finite difference equations, is discussed and some errors in the theory which appear to have gone undetected so far, are pointed out. It is shown that the previously published B‐polarization formulas are incorrect at points where regions of different conductivity meet, and that the E ‐polarization formulas are inaccurate when the step sizes of the numerical grid around the point are uneven. An appropriately‐modified version of the two‐dimensional theory is developed on the assumption that the Earth's conductivity is a smoothly‐varying function of position, a method which naturally lends itself to three‐dimensional generalization. All the required finite‐difference formulas are derived in detail, and presented in a form which is suitable for programming. A simple numerical calculation is given to illustrate the application of the method and the results are compared with those obtained from previous work.