A self-adaptive goal-oriented hp-finite element method with electromagnetic applications. Part II: Electrodynamics

Abstract We present the formulation, implementation, and applications of a self-adaptive, goal-oriented, hp -Finite Element (FE) Method for Electromagnetic (EM) problems. The algorithm delivers (without any user interaction) a sequence of optimal hp -grids. This sequence of grids minimizes the error in a prescribed quantity of interest with respect to the problem size, and it converges exponentially in terms of the relative error in a user-prescribed quantity of interest against the CPU time, including problems involving high material contrasts, boundary layers, and/or several singularities. The goal-oriented refinement strategy is an extension of a fully automatic, energy-norm based, hp -adaptive algorithm. We illustrate the efficiency of the method with 2D numerical simulations of Maxwell’s equations using both H 1 -conforming (continuous) elements and H(curl)-conforming (Nedelec edge) elements. Applications include alternate current (AC) resistivity logging instruments in a borehole environment with steel casing for the assessment of rock formation properties behind casing. Logging instruments, steel casing, and rock formation properties are assumed to exhibit axial symmetry around the axis of a vertical borehole. For the presented challenging class of problems, the self-adaptive goal-oriented hp -FEM delivers results with 5–7 digits of accuracy in the quantity-of-interest.

[1]  Leszek Demkowicz,et al.  H1, H(curl) and H(div)-conforming projection-based interpolation in three dimensionsQuasi-optimal p-interpolation estimates , 2005 .

[2]  Mark Ainsworth,et al.  Topics in Computational Wave Propagation , 2003 .

[3]  David Pardo.  Integration of hp-adaptivity with a two grid solver: applications to electromagnetics , 2004 .

[4]  David Pardo,et al.  A goal‐oriented hp‐adaptive finite element method with electromagnetic applications. Part I: electrostatics , 2006 .

[5]  Leszek Demkowicz,et al.  Goal-oriented hp-adaptivity for elliptic problems , 2004 .

[6]  Thomas J. R. Hughes,et al.  Encyclopedia of computational mechanics , 2004 .

[7]  L. Demkowicz,et al.  An hp‐adaptive finite element method for electromagnetics. Part 3: A three‐dimensional infinite element for Maxwell's equations , 2003 .

[8]  Serge Prudhomme,et al.  On goal-oriented error estimation for elliptic problems: application to the control of pointwise errors , 1999 .

[9]  Rolf Rannacher,et al.  A posteriori error control in finite element methods via duality techniques: Application to perfect plasticity , 1998 .

[10]  Leszek F. Demkowicz,et al.  A Fully Automatic hp-Adaptivity , 2002, J. Sci. Comput..

[11]  Franck Assous,et al.  Theoretical tools to solve the axisymmetric Maxwell equations , 2002 .

[12]  R. RannacherInstitut,et al.  Weighted a Posteriori Error Control in Fe Methods , 1995 .

[13]  Leszek Demkowicz,et al.  An hp‐adaptive finite element method for electromagnetics—part II: A 3D implementation , 2002 .

[14]  R. Garbacz,et al.  Introduction to electromagnetic fields , 1982, IEEE Antennas and Propagation Society Newsletter.

[15]  J. Bladel.  Singular electromagnetic fields and sources , 1996 .

[16]  Leszek Demkowicz,et al.  Integration of hp-adaptivity and a two-grid solver for elliptic problems , 2006 .

[17]  J. Oden,et al.  Goal-oriented error estimation and adaptivity for the finite element method , 2001 .

[18]  Rolf Rannacher,et al.  Duality-based adaptivity in the hp-finite element method , 2003, J. Num. Math..

[19]  J. R. Lovell.  Finite Element Methods in Resistivity Logging , 1993 .

[20]  Anthony T. Patera,et al.  A hierarchical duality approach to bounds for the outputs of partial differential equations , 1998 .

[21]  R. Harrington.  Time-Harmonic Electromagnetic Fields , 1961 .

[22]  Leszek Demkowicz,et al.  hp-Adaptive Finite Elements for Time-Harmonic Maxwell Equations , 2003 .

[23]  L. Demkowicz,et al.  An hp-adaptive finite element method for electromagnetics: Part 1: Data structure and constrained approximation , 2000 .

[24]  L. Demkowicz,et al.  Integration of hp-adaptivity and a two grid solver for electromagnetic problems , 2006 .

[25]  Leszek Demkowicz,et al.  Fully automatic hp-adaptivity in three dimensions , 2006 .