Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

3D electromagnetic modeling using staggered finite differences

Conference ·
OSTI ID:501494

The method of finite differences has been employed to solve a variety of 3D electromagnetic (EM) forward problems arising in geophysical applications. Specific sources considered include dipolar and magnetotelluric (MT) field excitation in the frequency domain. In the forward problem, the EM fields are simulated using a vector Helmholtz equation for the electric field, which are approximated using finite differences on a staggered grid. To obtain the fields, a complex-symmetric matrix system of equations is assembled and iteratively solved using the quasi minimum method (QMR) method. Perfectly matched layer (PML) absorbing boundary conditions are included in the solution and are necessary to accurately simulate fields in propagation regime (frequencies > 10 MHZ). For frequencies approaching the static limit (< 10 KHz), the solution also includes a static-divergence correction, which is necessary to accurately simulate MT source fields and can be used to accelerate convergence for the dipolar source problem.

Research Organization:
Sandia National Labs., Albuquerque, NM (United States)
Sponsoring Organization:
USDOE Office of Energy Research, Washington, DC (United States)
DOE Contract Number:
AC04-94AL85000
OSTI ID:
501494
Report Number(s):
SAND--97-1295C; CONF-970847--3; ON: DE97006366
Country of Publication:
United States
Language:
English