Vol. 32

Latest Volume

All Volumes

PIER 176 [2023]
PIER 175 [2022]
PIER 174 [2022]
PIER 173 [2022]
PIER 172 [2021]
PIER 171 [2021]
PIER 170 [2021]
PIER 169 [2020]
PIER 168 [2020]
PIER 167 [2020]
PIER 166 [2019]
PIER 165 [2019]
PIER 164 [2019]
PIER 163 [2018]
PIER 162 [2018]
PIER 161 [2018]
PIER 160 [2017]
PIER 159 [2017]
PIER 158 [2017]
PIER 157 [2016]
PIER 156 [2016]
PIER 155 [2016]
PIER 154 [2015]
PIER 153 [2015]
PIER 152 [2015]
PIER 151 [2015]
PIER 150 [2015]
PIER 149 [2014]
PIER 148 [2014]
PIER 147 [2014]
PIER 146 [2014]
PIER 145 [2014]
PIER 144 [2014]
PIER 143 [2013]
PIER 142 [2013]
PIER 141 [2013]
PIER 140 [2013]
PIER 139 [2013]
PIER 138 [2013]
PIER 137 [2013]
PIER 136 [2013]
PIER 135 [2013]
PIER 134 [2013]
PIER 133 [2013]
PIER 132 [2012]
PIER 131 [2012]
PIER 130 [2012]
PIER 129 [2012]
PIER 128 [2012]
PIER 127 [2012]
PIER 126 [2012]
PIER 125 [2012]
PIER 124 [2012]
PIER 123 [2012]
PIER 122 [2012]
PIER 121 [2011]
PIER 120 [2011]
PIER 119 [2011]
PIER 118 [2011]
PIER 117 [2011]
PIER 116 [2011]
PIER 115 [2011]
PIER 114 [2011]
PIER 113 [2011]
PIER 112 [2011]
PIER 111 [2011]
PIER 110 [2010]
PIER 109 [2010]
PIER 108 [2010]
PIER 107 [2010]
PIER 106 [2010]
PIER 105 [2010]
PIER 104 [2010]
PIER 103 [2010]
PIER 102 [2010]
PIER 101 [2010]
PIER 100 [2010]
PIER 99 [2009]
PIER 98 [2009]
PIER 97 [2009]
PIER 96 [2009]
PIER 95 [2009]
PIER 94 [2009]
PIER 93 [2009]
PIER 92 [2009]
PIER 91 [2009]
PIER 90 [2009]
PIER 89 [2009]
PIER 88 [2008]
PIER 87 [2008]
PIER 86 [2008]
PIER 85 [2008]
PIER 84 [2008]
PIER 83 [2008]
PIER 82 [2008]
PIER 81 [2008]
PIER 80 [2008]
PIER 79 [2008]
PIER 78 [2008]
PIER 77 [2007]
PIER 76 [2007]
PIER 75 [2007]
PIER 74 [2007]
PIER 73 [2007]
PIER 72 [2007]
PIER 71 [2007]
PIER 70 [2007]
PIER 69 [2007]
PIER 68 [2007]
PIER 67 [2007]
PIER 66 [2006]
PIER 65 [2006]
PIER 64 [2006]
PIER 63 [2006]
PIER 62 [2006]
PIER 61 [2006]
PIER 60 [2006]
PIER 59 [2006]
PIER 58 [2006]
PIER 57 [2006]
PIER 56 [2006]
PIER 55 [2005]
PIER 54 [2005]
PIER 53 [2005]
PIER 52 [2005]
PIER 51 [2005]
PIER 50 [2005]
PIER 49 [2004]
PIER 48 [2004]
PIER 47 [2004]
PIER 46 [2004]
PIER 45 [2004]
PIER 44 [2004]
PIER 43 [2003]
PIER 42 [2003]
PIER 41 [2003]
PIER 40 [2003]
PIER 39 [2003]
PIER 38 [2002]
PIER 37 [2002]
PIER 36 [2002]
PIER 35 [2002]
PIER 34 [2001]
PIER 33 [2001]
PIER 32 [2001]
PIER 31 [2001]
PIER 30 [2001]
PIER 29 [2000]
PIER 28 [2000]
PIER 27 [2000]
PIER 26 [2000]
PIER 25 [2000]
PIER 24 [1999]
PIER 23 [1999]
PIER 22 [1999]
PIER 21 [1999]
PIER 20 [1998]
PIER 19 [1998]
PIER 18 [1998]
PIER 17 [1997]
PIER 16 [1997]
PIER 15 [1997]
PIER 14 [1996]
PIER 13 [1996]
PIER 12 [1996]
PIER 11 [1995]
PIER 10 [1995]
PIER 09 [1994]
PIER 08 [1994]
PIER 07 [1993]
PIER 06 [1992]
PIER 05 [1991]
PIER 04 [1991]
PIER 03 [1990]
PIER 02 [1990]
PIER 01 [1989]

All Issues

Covariant Isotropic Constitutive Relations in Clifford's Geometric Algebra

Constitutive relations for isotropic material media are formulated in a manifestly covariant manner. Clifford's geometric algebra is used throughout. Polarisable,c hiral and Tellegen medium are investigated. The investigation leads to the discovery of an underlying algebraic structure that completely classifies isotropic media. Variational properties are reviewed,sp ecial attention is paid to the imposed constraints on material parameters. Covariant reciprocity condition is given. Finally,dualit y transformations and their relevance to constitutive relations are investigated. Duality is shown to characterise ‘well-behavedness' of medium which has an interesting metric tensor related implication.

Implementing the Perfectly Matched Layer Absorbing Boundary Condition with Mimetic Differencing Schemes

This paper concerns the implementation of the perfectly matched layer (PML) absorbing boundary condition in the framework of a mimetic differencing scheme for Maxwell's Equations. We use mimetic versions of the discrete curl operator on irregular logically rectangular grids to implement anisotropic tensor formulation of the PML. The form of the tensor we use is fixed with respect to the grid and is known to be perfectly matched in the continuous limit for orthogonal coordinate systems in which the metric is constant, i.e. Cartesian coordinates, and quasi-perfectly matched (quasi-PML) for curvilinear coordinates. Examples illustrating the effectiveness and long-term stability of the methods are presented. These examples demonstrate that the grid-based coordinate implementation of the PML is effective on Cartesian grids, but generates systematic reflections on grids which are orthogonal but non-Cartesian (quasi-PML). On non-orthogonal grids progressively worse performance of the PML is demonstrated. The paper begins with a summary derivation of the anisotropic formulation of the perfectly matched layer and mimetic differencing schemes for irregular logically rectangular grids.

Frequency Domain Analysis of Waveguides and Resonators with Fit on Non-Orthogonal Triangular Grids

The focus of this paper is on the solution of Maxwell's equations for time-harmonic fields on triangular, possibly nonorthogonal meshes. The method is based on the well-known Finite Integration Technique (FIT) [33, 35] which is a proven consistent discretization method for the computation of electromagnetic fields. FIT on triangular grids was first introduced in [29, 31] for eigenvalue problems arising in the design of accelerator components and dielectric loaded waveguides. For many technical applications the 2D simulation on a triangular grid combines the advantages of FIT, as e.g. the consistency of the method or the numerical advantage of banded system matrices, with the geometrical flexibility of non-coordinate grids. The FIT-discretization on non-orthogonal 2D grids has close relations [26] to the N´ed´elec elements [14, 15] or edge elements in the Finite Element Method.

Computational Aspects of the Cell Method in Electrodynamics

A desire to unify the mathematical description of many physical theories, such as electrodynamics, mechanics, thermal conduction, has led us to understand that global (=integral) physical variables of eachth eory can be associated to spatial geometrical elements suchas points, lines, surfaces, volumes and temporal geometrical elements suchas instants and intervals. This association has led us to build a space-time classification diagram of variables and equations for each theory. Moreover, the possibility to express physical laws directly in a finite rather than differential form has led to the development of a computational methodology called Cell Method [12]. This paper discusses some practical aspects of this method and how it may overcome some of the main limitations of FDTD method in electrodynamics. Moreover, we will provide some numerical examples to compare the two methodologies.

Conservation of Discrete Energy and Related Laws in the Finite Integration Technique

We report some properties of the Finite Integration Technique (FIT), which are related to the definition of a discrete energy quantity. Starting with the well-known identities for the operator matrices of the FIT, not only the conservation of discrete energy in time and frequency domain simulations is derived, but also some important orthogonality properties for eigenmodes in cavities and waveguides. Algebraic proofs are presented, which follow the vectoranalytical proofs of the related theorems of the classical (continuous) theory. Thus, the discretization approach of the FIT can be considered as the framework for a consistent discrete electromagnetic field theory.

Higher Order Whitney Forms

The calculus of differential forms can be used to devise a unified description of discrete differential forms of any order and polynomial degree on simplicial meshes in any spatial dimension. A general formula for suitable degrees of freedom is also available. Fundamental properties of nodal interpolation can be established easily. It turns out that higher order spaces, including variants with locally varying polynomial order, emerge from the usual Whitneyforms by local augmentation. This paves the way for an adaptive pversion approach to discrete differential forms.

Discrete Hodge-Operators: an Algebraic Perspective

Discrete differential forms should be used to deal with the discretization of boundary value problems that can be stated in the calculus of differential forms. This approach preserves the topological features of the equations. Yet, the discrete counterparts of the metricdependent constitutive laws remain elusive. I introduce a few purely algebraic constraints that matrices associated with discrete material laws have to satisfy. It turns out that most finite element and finite volume schemes comply with these requirements. Thus convergence analysis can be conducted in a unified setting. This discloses basic sufficient conditions that discrete material laws have to meet in order to ensure convergence in the relevant energy norms.

Finite Element-Based Algorithms to Make Cuts for Magnetic Scalar Potentials: Topological Constraints and Computational Complexity

This paper outlines a generic algorithm to generate cuts for magnetic scalar potentials in 3-dimensional multiply-connected finite element meshes. The algorithm is based on the algebraic structures of (co)homology theory with differential forms and developed in the context of the finite element method and finite element data structures. The paper also studies the computational complexity of the algorithm and examines how the topology of the region can create an obstruction to finding cuts in O(m20) time and O(m0) storage, where m0 is the number of vertices in the finite element mesh. We argue that in a problem where there is no a priori data about the topology, the algorithm complexity is O(m20) in time and O(m4/30 ) in storage. We indicate how this complexity can be achieved in implementation and optimized in the context of adaptive mesh refinement.

Topological Approach to Computational Electromagnetism

Software systems designed to solve Maxwell's equations need abstractions that accurately explain what different kinds of electromagnetic problems really do have in common. Computational electromagnetics calls for higher level abstractions than what is typically needed in ordinary engineering problems. In this paper Maxwell's equations are described by exploiting basic concepts of set theory. Although our approach unavoidably increases the level of abstraction,it also simplifies the overall view making it easier to recognize a topological problem behind all boundary value problems modeling the electromagnetic phenomena. This enables us also to construct an algorithm which tackles the topological problem with basic tools of linear algebra.

Geometric Aspects of the Simplicial Discretization of Maxwell's Equations

Aspects of the geometric discretization of electromagnetic fields on simplicial lattices are considered. First, the convenience of the use of exterior differential forms to represent the field quantities through their natural role as duals (cochains) of the geometric constituents of the lattice (chains = nodes, edges, faces, volumes) is briefly reviewed. Then, the use of the barycentric subdivision to decompose the (ordinary) simplicial primal lattice together with the (twisted) non-simplicial barycentric dual lattice into simplicial elements is considered. Finally, the construction of lattice Hodge operators by using Whitney maps on the first barycentric subdivision is described. The objective is to arrive at a discrete formulation of electromagnetic fields on general lattices which better adheres to the underlying physics.

Data Structures for Geometric and Topological Aspects of Finite Element Algorithms

This paper uses simplicial complexes and simplicial (co)homology theory to expose a foundation for data structures for tetrahedral finite element meshes. Identifying tetrahedral meshes with simplicial complexes leads, by means of Whitney forms, to the connection between simplicial cochains and fields in the region modeled by the mesh. Furthermore, lumped field parameters are tied to matrices associated with simplicial (co)homology groups. The data structures described here are sparse, and the computational complexity of constructing them is O(n) where n is the number of vertices in the finite element mesh. Non-tetrahedral meshes can be handled by an equivalent theory. These considerations lead to a discrete form of Poincar´e duality which is a powerful tool for developing algorithms for topological computations on finite element meshes. This duality emerges naturally in the data structures. We indicate some practical applications of both data structures and underlying theory.

The Geometry of Time-Stepping

The space-time geometric structure of Maxwell's equations is examined and a subset of them is found to define a pair of exact discrete time-stepping relations. The desirability of adopting an approach to the discretization of electromagnetic problems which exploits this fact is advocated, and the name topological time-stepping for numerical schemes complying with it is suggested. The analysis of the equations leading to this kind of time-stepping reveals that these equations are naturally written in terms of integrated field quantities associated with space-time domains. It is therefore suggested that these quantities be adopted as state variables within numerical methods. A list of supplementary prescriptions for a discretization of electromagnetic problems suiting this philosophy is given, with particular emphasis on the necessity to adopt a space-time approach in each discretization step. It is shown that some existing methods already comply with these tenets, but that this fact is not explicitly recognized and exploited. The role of the constitutive equations in this discretization philosophy is briefly analyzed. The extension of this approach to more general kinds of space-time meshes, to other sets of basic time-stepping equations and to other field theories is finally considered.

Mimetic Finite Difference Methods for Maxwell's Equations and the Equations of Magnetic Diffusion

We have constructed mimetic finite difference methods for both the TE and TM modes for 2-D Maxwell's curl equations and equations of magnetic diffusion with discontinuous coefficients on nonorthogonal, nonsmooth grids. The discrete operators were derived using the discrete vector and tensor analysis to satisfy discrete analogs of the main theorems of vector analysis. Because the finite difference methods satisfy these theorems, they do not have spurious solutions and the "divergence-free" conditions for Maxwell's equations are automatically satisfied. The tangential components of the electric field and the normal components of magnetic flux used in the FDM are continuous even across discontinuities. This choice guarantees that problems with strongly discontinuous coefficients are treated properly. Furthermore on rectangular grids the method reduces to the analytically correct averaging for discontinuous coefficients. We verify that the convergence rate was between first and second order on the arbitrary quadrilateral grids and demonstrate robustness of the method in numerical examples.

Discrete Electromagnetism with the Finite Integration Technique

The Finite Integration Technique (FIT) is a consistent discretization scheme for Maxwell's equations in their integral form. The resulting matrix equations of the discretized fields can be used for efficient numerical simulations on modern computers. In addition, the basic algebraic properties of this discrete electromagnetic field theory allow to analytically and algebraically prove conservation properties with respect to energy and charge of the discrete formulation and gives an explanation of the stability properties of numerical formulations in the time domain.

'Generalized Finite Differences' in Computational Electromagnetics

The geometrical approach to Maxwell's equations promotes a way to discretize them that can be dubbed "Generalized Finite Differences", which has been realized independently in several computing codes. The main features of this method are the use of two grids in duality, the "metric-free" formulation of the main equations (Amp`ere and Faraday), and the concentration of metric information in the discrete representation of the Hodge operator. The specific role that finite elements have to play in such an approach is emphasized, and a rationale for Whitney forms is proposed, showing why they are the finite elements of choice.

Finite Formulation of the Electromagnetic Field

The objective of this paper is to present an approach to electromagnetic field simulation based on the systematic use of the global (i.e. integral) quantities. In this approach, the equations of electromagnetism are obtained directly in a finite form starting from experimental laws without resortingto the differential formulation. This finite formulation is the natural extension of the network theory to electromagnetic field and it is suitable for computational electromagnetics.