Three-dimensional wideband electromagnetic modeling on massively parallel computers

被引:101
作者
Alumbaugh, DL
Newman, GA
Prevost, L
Shadid, JN
机构
[1] Sandia National Laboratories, Albuquerque, NM
[2] Sandia National Laboratory, M.S. 0750, Albuquerque, NM 87185
关键词
D O I
10.1029/95RS02815
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
A method is presented for modeling the wideband, frequency domain electromagnetic (EM) response of a three-dimensional (3-D) earth to dipole sources operating at frequencies where EM diffusion dominates the response (less than 100 kHz) up into the range where propagation dominates (greater than 10 MHz). The scheme employs the modified form of the vector Helmholtz equation for the scattered electric fields to model variations in electrical conductivity, dielectric permitivity and magnetic permeability. The use of the modified form of the Helmholtz equation allows for perfectly matched layer (PML) absorbing boundary conditions to be employed through the use of complex grid stretching. Applying the finite difference operator to the modified Helmholtz equation produces a linear system of equations for which the matrix is sparse and complex symmetrical. The solution is obtained using either the biconjugate gradient (BICG) or quasi-minimum residual (QMR) methods with preconditioning; in general we employ the QMR method with Jacobi scaling preconditioning due to stability. In order to simulate larger, more realistic models than has been previously possible, the scheme has been modified to run on massively parallel (MP) computer architectures. Execution on the 1840-processor Intel Paragon has indicated a maximum model size of 280 x 260 x 200 cells with a maximum hop rate of 14.7 Gflops. Three different geologic models are simulated to demonstrate the use of the code for frequencies ranging from 100 Hz to 30 MHz and for different source types and polarizations. The simulations show that the scheme is correctly able to model the air-earth interface and the jump in the electric and magnetic fields normal to discontinuities. For frequencies greater than 10 MHz, complex grid stretching must be employed to incorporate absorbing boundaries while below this normal (real) grid stretching can be employed.
引用
收藏
页码:1 / 23
页数:23
相关论文
共 26 条
[1]  
ALUMBAUGH DL, 1995, APR SAGEEP 95 S APPL
[2]  
ALUMBAUGH DL, 1994, 64TH ANN INT M SOC E
[3]  
Barrett Richard, 1994, Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods
[4]  
Berenger J., 1993, J COMPUT PHYS, V114, P185
[5]  
Chew W., 1994, Waves and Fields in Inhomogeneous Media
[6]   A 3D PERFECTLY MATCHED MEDIUM FROM MODIFIED MAXWELLS EQUATIONS WITH STRETCHED COORDINATES [J].
CHEW, WC ;
WEEDON, WH .
MICROWAVE AND OPTICAL TECHNOLOGY LETTERS, 1994, 7 (13) :599-604
[7]   SPECTRAL APPROACH TO SOLVING 3-DIMENSIONAL MAXWELL DIFFUSION-EQUATIONS IN THE TIME AND FREQUENCY DOMAINS [J].
DRUSKIN, V ;
KNIZHERMAN, L .
RADIO SCIENCE, 1994, 29 (04) :937-953
[8]  
Druskin V. L., 1988, Izvestiya Academy of Sciences USSR, Physics of the Solid Earth, V24, P641
[9]   CONJUGATE GRADIENT-TYPE METHODS FOR LINEAR-SYSTEMS WITH COMPLEX SYMMETRICAL COEFFICIENT MATRICES [J].
FREUND, RW .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1992, 13 (01) :425-448
[10]   BEYOND THE BORN AND RYTOV APPROXIMATIONS - A NONLINEAR APPROACH TO ELECTROMAGNETIC SCATTERING [J].
HABASHY, TM ;
GROOM, RW ;
SPIES, BR .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 1993, 98 (B2) :1759-1775