Abstract:Quantitative interpretation of the data from controlled‐source electromagnetic methods, whether via forward modelling or inversion, requires solving a considerable number of forward problems, and multigrid methods are often employed to accelerate the solving process. In this study, a new extrapolation cascadic multigrid method is employed to solve the large sparse complex linear system arising from the finite element approximation of Maxwell's equations using secondary potentials. The equations using secondary potentials are discretized by the classic nodal finite element method on nonuniform rectilinear grids. The resulting linear systems are solved by the extrapolation cascadic multigrid method with a new prolongation operator and preconditioned Stabilized bi‐conjugate gradient method smoother. High‐order interpolation and global extrapolation formulas are utilized to construct the multigrid prolongation operator. The extrapolation cascadic multigrid method with the new prolongation operator is easier to implement and more flexible in application than the original one. Finally, several synthetic examples including layered models, models with anisotropic anomalous bodies or layers, are used to validate the accuracy and efficiency of the proposed method. Numerical results show that the extrapolation cascadic multigrid method improves the efficiency of 3D controlled‐source electromagnetic forward modelling a lot, compared with traditional iterative solvers and some state‐of‐the‐art methods or software (e.g., preconditioned flexible generalized minimal residual method, emg3d) in the considered models and grid settings. The efficiency benefit is more evident as the number of unknowns increases, and the proposed method is more efficient at low frequencies. The extrapolation cascadic multigrid method can also be used to solve systems of equations arising from related applications, such as induction logging, airborne electromagnetic, etc.

Inverse Problem of the Multislice Method in Retrieving Projected Complex Potentials from the Exit-Wave Function.

Exit-wave Phase Retrieval from a Single High-Resolution Transmission Electron Microscopy Image of a Weak-Phase Object

An Improved Method Assigning Three-Dimensional Atomic Potentials to Multiple Slices in Exit-Wave Simulations of Transmission Electron Microscopy.

A Novel Pseudospectral Time Domain Method for Complex EM Simulations

Far-field Characterization of Material Permittivity Based on Nonlinear Inversions

Retrieving Material Permittivity Under Oblique Incidence

Experimental Investigation on Through-Wall Imaging Based on Non-Linear Inversions

Reconstruction of Complex Single-Particle Images Using Charge-Flipping Algorithm

Exact inversion of partially coherent dynamical electron scattering for picometric structure retrieval

Linearized Inverse Schrödinger Potential Problem at a Large Wavenumber

Theory and Application of Wavelet Based Bi-Orthonormal Decomposition Method in the Solution of Linear Inverse Problems in Electromagnetics

Contracting Electromagnetic Full-Wave Inversion of 2-D Inhomogeneous Objects With Irregular Shapes Based on the Hybrid SESI Forward Solver

An Iterative Numerical Algorithm for Electromagnetic Imaging

Image reconstruction from electron and X-ray diffraction patterns using iterative algorithms: experiment and simulation

3D Controlled‐source electromagnetic modelling in anisotropic media using secondary potentials and a cascadic multigrid solver

Robust Determination Of The Chemical Potential In The Pole Expansion And Selected Inversion Method For Solving Kohn-Sham Density Functional Theory

The Modeling of 3-D DC Resistivity Based on Integral Equation of Potential

Variational projector-augmented wave method: A full-potential approach for electronic structure calculations in solid-state physics

Simple Algorithm for Partial Wave Expansion of Plasmonic and Evanescent Fields.

A New Scheme for the Exit-Wave Reconstruction from A Small Set of Focus Series of Hrem Images

A Synergistic GPR Approach of Back Projection Algorithm and Full-Wave Inversion for Reconstruction in Unknown Multilayered Environments