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.

Complex Matrix Equation Solver Based on Computational Metasurface

A Scattering-Matrix-Based Field Stitching Technique for Accelerating the Simulation of Large Area Metasurfaces

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

Leveraging Parallel Speedup for Intelligent Design of Metasurface

Multiphysics Computing of Challenging Antenna Arrays under a Supercomputer Framework

Massively Parallel Simulation Of Large-Scale Electromagnetic Problems Using One High-Performance Computing Scheme And Domain Decomposition Method

Electromagnetic-Thermal Co-simulation of Large Antenna Array on Platform Using Enhanced Finite Element Solver with Massively Parallel Computing Capability

A Novel Matrix-Free Finite Element Method for Time-Harmonic Maxwell’s Equations

Massively Parallel Multiphysics Solver for Simulating Large-Scale Antenna Array and Other 3-D Structures

Solving Multivariable Equations with Tandem Metamaterial Kernels

Reconfigurable metamaterial processing units that solve arbitrary linear calculus equations

A Nested Iterative Solver for Sub-structure CMA of Composite Metallic-Dielectric Objects

Fullwave design of cm-scale cylindrical metasurfaces via fast direct solvers

A Novel Parallel Direct Solver Using Adaptive Cross Approximation for Analysis of Electromagnetic Radiation Problems With Complex Structures

Reconfigurable application-specific photonic integrated circuit for solving partial differential equations

Solving partial differential equations with waveguide-based metatronic networks

H-Matrix Accelerated Direct Matrix Solver for Maxwell's Equations using the Chebyshev-based Nyström Boundary Integral Equation Method

Arbitrary Jones matrix on-demand design in metasurface using multiple meta-atoms

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

Complex-Amplitude Programmable Versatile Metasurface Platform Driven by Guided Wave

Inverse design of dielectric metasurface by spatial coupled mode theory