arxiv:physics/ v1 [physics.optics] 3 Jul 2003

Transcription

1 arxiv:physics/ v1 [physics.optics] 3 Jul 2003 Note on transformation to general curvilinear coordinates for Maxwell s curl equations D M Shyroki Department of Theoretical Physics, Belarusian State University, Fr. Skaryna Avenue 4, Mensk, Belarus shyroki@tut.by Abstract. Maxwell s curl equations are formulated in curvilinear non-orthogonal coordinates in three dimensions in a manner that enables to preserve conventional Cartesian mesh for the finite-difference schemes by means of mere redefinition of the permittivity and permeability coefficients. It is highlighted that in the previous works on this subject (Ward and Pendry 1996 J. Modern Opt ; 1998 Phys. Rev. B ) inaccurate transformation to curvilinear coordinates, rooted in questionable transformation laws assigned to electric and magnetic field vectors, led to physical dubiety of the final formulae. In this note the amended equations are presented, capable of being easily adopted for the accurate finite-difference time-domain (FDTD) modelling of electromagnetic propagation and scattering problems in complicated geometries within the existing computational codes. PACS numbers: De, Bf

2 Note on transformation to general curvilinear coordinates for Maxwell s Introduction Maxwell s equations span a long life of nearly one and a half century, yet it often remains a real challenge to the theorist to solve, say, the problem of light scattering from a complex surface or propagation in state-of-the-art photonic band-gap material. To tackle all those numerous problems, a vast array of numerical methods have been developed, in particular the finite-difference time-domain (FDTD) technique [1, 2] called otherwise the order-n method [3] due to its linear scaling with the size of the system. In its framework, the problem is formulated as an initial value problem leaned on discretized Maxwell s curl equations which can be integrated numerically. The FDTD algorithm was originally given by Yee [4] in rectangular Cartesian coordinates, i.e., on the cube mesh of points constituting the computational space. Meanwhile, many objects of interest have curved surfaces which are purely represented in a rectilinear mesh. Furthermore, severe problems associated with a uniform mesh arise when the refractive index varies strongly within the structure, for example in metallodielectric components. To circumvent these shortcomings, Holland [5] attempted to formulate the FDTD algorithm in general curvilinear coordinates, but extreme complexity of curvilinear mesh that exacts a toll in computer time [6] and makes erecting a system of discretized equations a programmer s nightmare [7] restrained the practical usage of Holland s result. Even in a two-dimensional case [6, 8], the awkwardness of the formulae leaves one beyond any hope to qualitatively analyse the equations. Fortunately, a few years ago Ward and Pendry [7, 9] noticed that considerable simplification of programming efforts can be achieved since transformation to curvilinear coordinates can be regarded as redefining the permittivity and permeability coefficients, ǫ and µ, so that the well-developed algorithm of computing on a rectangular mesh can be translated almost verbatim to the case of a non-uniform mesh, provided the mesh structure encoded in the effective ǫ and µ. This idea, being no doubt a really fruitful and natural one, was implemented afterwards in a computer program for calculating photonic band structures, Green s functions and transmission reflection coefficients [10]. An agreement with alternative numerical computations on orthogonal mesh and with plane-wave expansion method was believed confidently good if not perfect. However, one relevant inaccuracy in the cited works has been left out of sight: in Maxwellian electrodynamics, the electric and magnetic fields and displacements E, B, D, and H can be treated as vectors only with respect to pure rotations in threedimensional (3-D) Euclidean space, but care should be taken when transforming to 3-D Riemannian space referenced to arbitrary curvilinear coordinate system. Meanwhile, Ward and Pendry in a manner of early Holland s misconstruction grounded upon the highly questionable representation of affine-invariant Maxwell s equations and in fact came to a physically dubious formulation of final analytic results. The scope of this note is to clarify the origin of that misconception and present the amended formalism aimed to facilitate the non-orthogonal curvilinear FDTD technique.

3 Note on transformation to general curvilinear coordinates for Maxwell s Maxwell s equations in non-orthogonal curvilinear coordinates The sound theoretical basis for any problem within the scope of classical electrodynamics is constituted by Maxwell s equations. They describe electromagnetic field as a unified quantity invariant with respect to Lorentz group, i.e., rotations in Minkowski space, and hence their utmost consistent formulation involves four-dimensional electromagnetic field and flux tensors. However, if the medium in the laboratory frame is at rest, it is often found convenient to decompose Maxwell s equations into fully vectorial system for three-dimensional quantities E, B, D, and H, written in cgs units as E = 1 B, B = 0, (1) c t H = 1 D c t + 1 4πρv, D = 4πρ, (2) c where ρ is the charge density, v is the charge velocity, c is the vacuum speed of light. The accompanying constitutive relations are D = ǫe, B = µh. (3) Naturally the above equations are expected to be coordinate-free, that is preserving the formal structure in any coordinate system, under any spatial coordinate transformation. A principal point to be understood here is that, scholarly speaking, E, B, D, and H are vectors with respect to rotations in three-dimensional Euclidean space referenced to rectangular Cartesian system only. Even the point reflection is known to brake the seeming symmetry between, say, E and B vectors in free-space counterparts of equations (1), (2), since under r r operation we have E E for electric polar vector versus B B for magnetic pseudovector. Another example on the surface is r ar conversion: in contrast with true vectors, E should transform like E a 1 E to keep the work of electromagnetic force δa = ee dr invariant. The vague ideas of pseudoquantities etc. find an elegant and comprehensible mathematical conceptualization in the framework of conventional tensor calculus. To be specific, for more than half a century the affine-invariant form of Maxwell s equations (1), (2) is known in both research [11] and textbook [12] literature: B αβ 2 [α E β] = 1, [α B βγ] = 0, (4) c t α Ĥ αβ = 1 ˆD β 1 c t c 4πˆρvβ, ˆDα α = 4πˆρ. (5) Here E α is covariant vector, B αβ = B [αβ] is covariant bivector, Ĥ αβ = Ĥ[αβ] is contravariant bivector-density, and ˆD α is contravariant vector-density; the square brackets denote alternation in a usual manner: [αβ] = 1 2! ( αβ βα ), etc. Note that despite of someone s contradictory claims [13], the above equations are actually metricfree; all the metric information is hidden in constitutive relations (see [12, ch 6]) ˆD α = ǫ g 1 2g αβ E β, B αβ = µ g 1 2g αγ g βδ Ĥ γδ, (6)

4 Note on transformation to general curvilinear coordinates for Maxwell s... 4 where g αβ is the metric tensor of Riemannian space referenced to the right-handed coordinate system x α ; α, β, γ, δ = 1, 2, 3; g = det(g αβ ). In FDTD method, only the curl equations are needed, provided the divergence equations satisfied initially. Eliminating the flux quantities via (6) yields 2 [α E β] = µ c g 1 2g αγ g Ĥγδ βδ, α Ĥ αβ = ǫ t c g 2g 1 βγ E γ (7) t for source-free region. Allowing for invariant correspondence between bivector-density Ĥ αβ and pseudovector H γ [12], one can rewrite (7) as E β x α eαβγ = µγδ H δ c t, H β x α eαβγ = ǫγδ E δ c t, (8) where e ijk is the fully skew-symmetric Levi-Civita symbol, and in a manner of [7, 9] I introduce the effective permittivity and permeability quantities ǫ γδ = ǫg 1 2 g γδ, µ γδ = µg 1 2 g γδ. (9) Thus obtained equations are isomorphous to Maxwell s curl equations on a conventional Cartesian mesh. To restore fields in Euclidean space referenced to Cartesian system x i (designated by Roman indices in contrast with Greek ones reserved for curvilinear coordinates), the transformation matrix i x α should be used: E i E i = xα x E α, H i H i i H i = xα x H α. (10) i The circle above the equality sign ( ) means that this equality is not invariant but holds in the given coordinate system. I also assume here that coordinate transformation preserves the handedness of the system, which is likely to comprise all the practicable cases of interest. 3. Comparison with Ward and Pendry Formulae (8), (9), and (10) constitute the core of generalized non-orthogonal FDTD method. In [7, 9] Ward and Pendry came to some principally similar results, but for the effective permittivity and permeability they obtained (in current notation and units) ǫ γδ = ǫg 1 2 g γδ, µ γδ = µg 1 2 g γδ. (11) Note that the powers of g are different in (9), namely, 1 in the formula for ǫ vs. 2 1 in 2 that for µ, while in (11) they are 1 both. This discrepancy is because Ward and Pendry 2 assumed with no argument the symmetry between E and H fields [7, p 777] which is actually a questionable point when general coordinate transformations are considered. The difference between the powers of g in (9) definitely seems to be more reasonable and physically-intuitive due to the distinct difference between E and H fields as for their transformation characteristics. Altogether I should remark that the discrepancy between (9) and (11) vanishes when we consider the group of volume-preserving transformations, e.g., introducing the skew lattice instead of orthogonal one, but is essential for the general coordinate transformations with g 1, especially in the case of magnetic media with permeability µ comparable to or even exceeding the permittivity ǫ.

5 Note on transformation to general curvilinear coordinates for Maxwell s Conclusion The formulae by Ward and Pendry [7, 9] for the Maxwell s curl equations in curvilinear non-orthogonal coordinates are amended to allow for the marked difference between electric and magnetic fields with respect to their transformation laws. The formulae are considerably simple and can be easily adopted for an accurate non-orthogonal FDTD modelling of electromagnetic propagation and scattering problems in complicated geometries within the existing computational codes, since a mere cosmetic revision is actually needed. The hitherto demonstrated agreement between numerical results computed by Ward and Pendry and others may be accounted for two reasons: (i) special character of coordinate systems used, namely those with affine volume g = 1, and (ii) prevailing impact of effective ǫ over µ in dielectric media with high ǫ and unit µ. Acknowledgments I d like to thank Dr A V Lavrinenko who acquainted me with Ward s and Pendry s work, and Dr M N Polozov for a qualified check to the key ideas of this note. References [1] Taflove A 1995 Computational Electrodynamics: The Finite-Difference Time-Domain Method (Boston: Artech House) [2] Kunz K S, Luebbers R J 1993 The Finite Difference Time Domain Method for Electro-Magnetics (Boca Raton, FL: CRC Press) [3] Chan C T, Yu Q L, Ho K M 1995 Order-N spectral method for electromagnetic waves Phys. Rev. B [4] Yee K S 1966 Numerical solution of initial boundary value problems involving Maxwell s equations in isotropic media IEEE Trans. Antennas Propagat [5] Holland R 1983 Finite difference solutions of Maxwell s equations in generalized nonorthogonal coordinates IEEE Trans. Nuc. Sci [6] Fusco M 1990 FDTD algorithm in curvilinear coordinates IEEE Trans. Antennas Propagat [7] Ward A J, Pendry J B 1996 Refraction and geometry in Maxwell s equations J. Modern Opt [8] Qiu M, He S 2000 A nonorthogonal finite-difference time-domain method for computing the band structure of a two-dimensional photonic crystal with dielectric and metallic inclusions J. Appl. Phys [9] Ward A J, Pendry J B 1998 Calculating photonic Green s functions using a nonorthogonal finitedifference time-domain method Phys. Rev. B [10] Ward A J, Pendry J B 2000 A program for calculating photonic band structures, Green s functions and transmission/reflection coefficients using non-orthogonal FDTD method Comput. Phys. Commun [11] Van Dantzig D 1934 Electromagnetism, independent of metrical geometry Proc. Kon. Ned. Akad. v. Wet [12] Schouten J A 1951 Tensor Analysis for Physicists (Oxford: Clarendon) [13] Teixeira F L 2001 Geometric aspects of the simplical discretization of Maxwell s equations Progress In Electromagnetic Research