Abstract:With geophysical surveys evolving from traditional 2D to 3D models, the large volume of data adds challenges to inversion, especially when aiming to resolve complex 3D structures. An iterative forward solver for the controlled-source electromagnetic (CSEM) method requires less memory than that for a direct solver; however, it is not easy to iteratively solve an ill-conditioned linear system of equations arising from finite-element discretization of Maxwell’s equations. To solve this problem, we have developed efficient and robust iterative solvers for frequency- and time-domain CSEM modeling problems. For the frequency-domain problem, we first transform the linear system into its equivalent real-number format, and then introduce an optimal block-diagonal preconditioner. Because the condition number of the preconditioned linear equation system has an upper bound of [Formula: see text], we can achieve fast solution convergence when applying a flexible generalized minimum residual solver. Applying the block preconditioner further results in solving two smaller linear systems with the same coefficient matrix. For the time-domain problem, we first discretize the partial differential equation for the electric field in time using an unconditionally stable backward Euler scheme. We then solve the resulting linear equation system iteratively at each time step. After the spatial discretization in the frequency domain, or space-time discretization in the time domain, we exploit the conjugate-gradient solver with auxiliary-space preconditioners derived from the Hiptmair-Xu decomposition to solve these real linear systems. Finally, we check the efficiency and effectiveness of our iterative methods by simulating complex CSEM models. The most significant advantage of our approach is that the iterative solvers we adopt have almost the same accuracy and robustness as direct solvers but require much less memory, rendering them more suitable for large-scale 3D CSEM forward modeling and inversion.

The SCBiCG class of algorithms for complex symmetric linear systems with applications in several electromagnetic model problems

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

An Improved Quasi-Static Finite-Difference Scheme for Induced Field Evaluation in Mri Based on the Biconjugate Gradient Method

An Efficient BiCGstab Solved Impedance Method for Induced Field Evaluation with a Hyperthermia Applicator.

An efficient impedance method for induced field evaluation based on a stabilized Bi-conjugate gradient algorithm.

A Novel Pseudospectral Time Domain Method for Complex EM Simulations

Lanczos-type Variants of the COCR Method for Complex Nonsymmetric Linear Systems

A Scalable HPC-Based Domain Decomposition Method for Multiphysics Modeling of RF Devices

Electromagnetic-Circuital-Thermal Co-Simulation for Integrated Circuits and Microwave Applications

Analyzing Electromagnetic Scattering Using Characteristic Basis Function Method With Compressed Sensing

A semi-conjugate gradient method for solving unsymmetric positive definite linear systems

Solution of large-scale 3D controlled-source electromagnetic modeling problem using efficient iterative solvers

A global variant of the COCR method for the complex symmetric Sylvester matrix equation AX + XB = C

Application of the SSOR preconditioned CG algorithm to the vector FEM for 3D full-wave analysis of electromagnetic-field boundary-value problems

Advanced Parallel Algorithm for System-Level EMC Modeling of High-Speed Electronic Package

Modified CRI iteration methods for complex symmetric indefinite linear systems

Comparison the Iterative Solvers for Large Sparse Matrix in 3D Electromagnetic Forward Modelling

Scalable solvers for complex electromagnetics problems

Generalized Krylov Recycling Methods For Solution Of Multiple Related Linear Equation Systems In Electromagnetic Analysis

Iterative Methods by Space Decomposition and Subspace Correction

Simulation of multiscale problems using equivalence principle algorithm