3D finite-difference frequency-domain code for electromagnetic induction tomography

The effect of shrapnel on target chamber components and experiments at large lasers such as the National Ignition Facility at LLNL and the Megajoule Laser at CESTA in France is an important issue in fielding targets and exposure samples. Modeling calculations are likely to be an important component of this effort. Some work in this area has been performed by French workers, who are collaborating with the LLNL on many issues relating to target chamber, experiment-component, and diagnostics survival. Experiments have been performed at the PhCbus laser in France to measure shrapnel produced by laser-driven targets; among these shots were experiments that accelerated spheres of a size characteristic of some of the more damaging shrapnel. These spheres were stopped in polyethylene witness plates. The penetration depth is characteristic of the velocity of the shrapnel. Experimental calibration of steel sphere penetration into polyethylene was performed at the CESTA facility. The penetration depth has been reported (ref. 1) and comparisons with modeling calculations have been made (ref. 2). There was interest in a comparison study of the modeling of these experiments to provide independent checks of the calculations. This work has been approved both by DOE headquarters and by the French Atomic Energy Commission (CEA); it is task number 99-3.2 of the 1999 ICF agreement between the DOE and the CEA. Daniel Gogny of the CEA who is on a long-term assignment to LLNL catalyzed this collaboration. This report contains the initial results of our modeling effort.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  Douglas LaBrecque,et al.  Monitoring an underground steam injection process using electrical resistance tomography , 1993 .

[3]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[4]  James G. Berryman,et al.  Electromagnetic induction tomography field experiment at Lost Hills, CA , 1998 .

[5]  Michael S. Zhdanov,et al.  Methods for modelling electromagnetic fields Results from COMMEMI—the international project on the comparison of modelling methods for electromagnetic induction , 1997 .

[6]  Jin-Fa Lee,et al.  A perfectly matched anisotropic absorber for use as an absorbing boundary condition , 1995 .

[7]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[8]  T. Habashy,et al.  Rapid 2.5‐dimensional forward modeling and inversion via a new nonlinear scattering approximation , 1994 .

[9]  Jon F. Claerbout,et al.  Fundamentals of Geophysical Data Processing: With Applications to Petroleum Prospecting , 1985 .

[10]  W. Heinrich,et al.  Improved finite-difference formulation in frequency domain for three-dimensional scattering problems , 1992 .

[11]  Alex Becker,et al.  A borehole-to-surface electromagnetic survey , 1998 .

[12]  M. Zhdanov,et al.  Quasi‐linear approximation in 3-D electromagnetic modeling , 1996 .

[13]  Jin-Fa Lee,et al.  A comparison of anisotropic PML to Berenger's PML and its application to the finite-element method for EM scattering , 1997 .

[14]  Michael S. Zhdanov,et al.  Underground imaging by frequency-domain electromagnetic migration , 1996 .