Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator

• J. Torquil Smith

  Geophysics Program, AK-50, University of Washington, Seattle, WA 98195

抄録

The preceding paper derives a staggered‐grid, finite‐difference approximation applicable to electromagnetic induction in the Earth. The staggered‐grid, finite‐difference approximation results in a linear system of equations [Formula: see text]x = b, where [Formula: see text] is symmetric but not Hermitian. This is solved using the biconjugate gradient method, preconditioned with a modified, partial Cholesky decomposition of [Formula: see text]. This method takes advantage of the sparsity of [Formula: see text], and converges much more quickly than methods used previously to solve the 3-D induction problem. When simulating a conductivity model at a number of frequencies, the rate of convergence slows as frequency approaches 0. The convergence rate at low frequencies can be improved by an order of magnitude, by alternating the incomplete Cholesky preconditioned biconjugate gradient method with a procedure designed to make the approximate solutions conserve current.

収録刊行物

• GEOPHYSICS

  GEOPHYSICS 61 (5), 1319-1324, 1996-09

  Society of Exploration Geophysicists

キーワード

• Geochemistry and Petrology
• Geophysics

詳細情報

• CRID: 1362262943954478976
• NII論文ID: 30022770741
• DOI: 10.1190/1.1444055
• ISSN: 19422156
• ISSN: 00168033
• Web Site: https://library.seg.org/doi/pdf/10.1190/1.1444055
• データソース種別:
  • Crossref
  • CiNii Articles