On the finite difference solution of two‐dimensional induction problems

Abstract

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. Copyright © 1976, Wiley Blackwell. All rights reserved