Computation of Maxwell s equations on Manifold using DEC
|
|
- Jennifer Shields
- 5 years ago
- Views:
Transcription
1 arxiv: v6 [math.na] 9 Dec 009 Computation of Maxwell s equations on Manifold using DEC Zheng Xie Yujie Ma. Center of Mathematical Sciences, Zhejiang University (3007),China. Key Laboratory of Mathematics Mechanization, Chinese Academy of Sciences, (00090), China Abstract In this paper, the method of discrete exterior calculus for numerically solving Maxwell s equations in space manifold and the time is discussed, which is a kind of lattice gauge theory. The analysis of its stable condition and error is also accomplished. This algorithm has been implemented on C++ plateform for simulating TE/M waves in vacuum. Keywords: Discrete exterior calculus, Discrete variation, Maxwell s equations, Lattice gauge theory. PACS(00): 4.0.Jb, 0.30.Jr, 0.40.Sf, 0.60.Cb. Introduction Computational electromagnetism is concerned with the numerical study of Maxwell s equations. The Yee scheme is known as finite difference time domain and is one of the most successful numerical methods, particularly in the area of microwave problems []. It preserves important structural lenozhengxie@yahoo.com.cn yjma@mmrc.iss.ac.cn This work is partially supported by CPSFFP (No ), NKBRPC (No. 004CB38000), and NNSFC (No )
2 features of Maxwell s equations [ 6]. Bossavit et al present the Yee-like scheme and extend Yee scheme to unstructured grids. This scheme combines the best attributes of the finite element method (unstructured grids) and Yee scheme (preserving geometric structure) [7, 8]. Stern et al [9] generalize the Yee scheme to unstructured grids not just in space, but in 4-dimensional spacetime by discrete exterior calculus(dec) [0 0]. This relaxes the need to take uniform time steps. In this paper, we generalize the Yee scheme to the discrete space manifold and the time. The spacetime manifold used here is split as a product of D time and D or 3D space manifold. The space manifold can be approximated by triangular and tetrahedrons depending on dimension, and the time by segments. So the spacetime manifold is approximated by prism lattice, on which the discrete Lorentz metric can be defined.. With the technique of discrete exterior calculus, the R value discrete connection, curvature and Bianchi identity are defined on prim lattice. With discrete variation of an inner product of discrete forms and their dual forms, the discrete source equation and continuity equation are derived.. Those equations compose the discrete Maxwell s equations in vacuum case, which just need the local information of triangulated manifold such as length and area. The discrete Maxwell s equations here can be re-grouped into two sets of explicit iterative schemes for TE and TM waves, respectively. Those schemes can directly use acute triangular, rectangular, regular polygon and their combination, which has been implemented on C++ plateform to simulate the electromagnetic waves propagation and interference on manifold. DEC for Maxwell s equations Maxwell s equations can be simply expressed once the language of exterior differential forms is used. The electric and magnetic fields are jointly described by a curvature form F in a 4-D spacetime manifold. The Maxwell s equations reduce to the Bianchi identity and the source equation df = 0 d F = J ()
3 where d denotes the exterior differential operator, denotes the Hodge star operator, and -form J is called the electric current form satisfying the continuity equation d J = 0. As the exterior derivative is defined on any manifold, the differential form version of the Bianchi identity makes sense for any 3D or 4D spacetime manifold, whereas the source equation is defined if the manifold is oriented and has a Lorentz metric. Now, we introduce the discrete counterpart of those differential geometric objects to derive the numerical computational schemes for Maxwell s equations. Discrete Lorentz metric The spacetime manifold used here is split as a product of D time and D or 3D space manifold. The D or 3D space manifold can be approximated by triangular or tetrahedrons, and the time by segments. The length of edge and area of triangular and volume of tetrahedrons gives the discrete Riemann metric on space grids. The metric on time grid is the minus of length square. The spacetime manifold is approximated by prism lattice, on which the discrete Lorentz metric can be defined as the product of discrete metric on space and time. Discrete exterior calculus A discrete differential k-form, k Z, is the evaluation of the differential k- form on all k-simplices. Dual forms, i.e., forms that we evaluate on the dual cell. Suppose each simplex contains its circumcenter. The circumcentric dual cell D(σ 0 ) of simplex σ 0 is D(σ 0 ) := Int(c(σ 0 )c(σ ) c(σ r )), σ 0 σ σ r where σ i is all the simplices which contains σ 0,..., σ i, and c(σ i ) is the circumcenter of σ i. The two operators in Eqs.() can be discretized as follows:. Discrete exterior differential operator d, this operator is the transpose of the incidence matrix of k-cells on k +-cells. 3
4 . Discrete Hodge Star, the operator scales the cells by the volumes of the corresponding dual and primal cells. Discrete connection and curvature Discrete connection form or gauge field A assigns to each element in the set of edges E an element of the gauge group R: A : E R. Discrete curvature form is the discrete exterior derivative of the discrete connection form F = da : P R. The value of F on each element in the set of triangular P is the coefficient of Holonomy group of this face. The form F automatically satisfies the discrete Bianchi identity df = 0. () Note that since the gauge group R used here is abelian, we need not pick a starting vertex for the loop. We may traverse the edges in any order, so long as taking orientations into account. Discrete Maxwell s equations For source case, we need discrete current form J. Let A = E Lagrangian functional be A i and the L(A,J) = da,da + A,J, where da,da := (A) E (d) E F ( ) F F (d) T F E (A)T E A,J := (A) E ( ) E E (J T ) E Supposing that there is a variation of A i, vanishing on the boundary, we 4
5 have Ai L(A,J) = Ai ( da,da + A,J ) = Ai ( (A) E (d) E F ( ) F F (d) T F E (A)T E +(A) E ( ) E E (J T ) E ) = (0,..., }{{},...,0) E (d) E F ( ) F F (d) T F E (A)T E i (A) E (d) E F ( ) F F (d) T F E (0,..., }{{},...,0) T E i +(0,..., }{{},...,0) E ( ) E E (J T ) E i = (0,..., }{{},...,0) E (d) E F ( ) F F (d) T F E (A)T E i +(0,..., }{{},...,0) E ( ) E E (J T ) E i The Hamilton s principle of stationary action states that this variation must equal zero for any such vary of A i, implying the Euler-Lagrange equations (d) E F ( ) F F (d) T F E (A)T E +( ) E E (J T ) E = 0, which is the discrete source equation δf = J, (3) where δ = d T. Since (d T ) = 0, the discrete continuity equation can express as: d T J = 0. (4) The equations of discrete Bianchi identity (), source equation (3), and continuity equation (4) are called discrete Maxwell s equations. Discrete Gauge transformations Discrete gauge transformations are maps A A+df for any 0 form or scalar function f on vertex. Since the discrete exterior derivative maps F F +d f = F, 5
6 the discrete Maxwell s equations(-4) are invariant under discrete gauge transformations. Since the discrete continuity equation (4) ensures so we have df,j = (f) V (d T ) V E ( ) E E (J T ) E = 0, L(A+df,J) = d(a+df),d(a+df) + A+df,J = L(A,J). That is to say the Lagrangian function is also invariant under discrete gauge transformations. 3 Explicit schemes Schemes for TE wave The discrete current form, discrete curvature form and its dual can be written as J = ( ρ e dt,j e ) F n+ = E n+ dt+b n F n = H n dt D n, where n and n+ denote the coordinate of the time, E = E i e i (electric E field) is discrete form on space, B = B i P i (magnetic field) is discrete P form on space, H = H i P i (magnetizing field) is the dual of B on P space, D = D i e i (electric displacement field) is the dual of E on space, E ρ e dt (charge density) is the discrete form on time, J e = J ei e i (electric E current density) is the discrete form on space. The discrete Maxwell s equations can be rewritten as d s B n = 0 d s E n+ dt = d t B n d T sd n = (ρ e dt) n d T sh n dt = d T t D n + Je, n where d s, d T s are the restriction of d, dt on space, and d t B n := Bn+ B n dt d T t D n := Dn+ D n 6 dt. (5)
7 If allowing for the possibility of magnetic charges and current discrete 3 form J = (ρ m, J m dt), the symmetric scheme can be written as where d s B n = ρ n m d s E n+ dt = d t B n J n+ m dt d T sd n = (ρ e dt) n d T sh n dt = d T t D n + Je, n ρ m = Tetρ mi T i (magnetic charges) is discrete 3-form on space, (6) J m = P J mi P i (current) is discrete form on space. The compact form of Eqs.(6) can be written as df = J d T F = J, with discrete continuity equations or integrability conditions d J = 0 d T J = 0. Proposition 3. If the initial condition satisfies the first and third equations in Eqs.(6), the solution of the second and fourth equations in Eqs.(6) automatically satisfy Eqs.(6). Proof. Because the dimension of spacetime is 3 +, therefore d T s (ρ e dt) = 0 d T t J e = 0 d s ρ m = 0 d t J m dt = 0, and the continuity equations can be reduced to So we have d T t ( ρ edt n ) d T s J n e = 0 d sj n+ m dt+d t ρ n m = 0. d T t dt s Dn d T t (ρ e dt) n = d T t (ρ e dt) n d T s (dt s Hn dt J n e ) = 0 d t d s B n d t ρ n m = d t ρ n m +d s(d s E n+ dt+j n+ m = 0. 7 dt)
8 Nowwe show thescheme (5) ontheproduct ofddiscrete spacemanifold and time. The second and fourth equations in Eqs.(5) based on Fig. are D n+ D n Bn+ B n +J n e = Hn Hn e = En+ e +E n+ e +E n+ 3 e 3. P where denotes the measure of forms and dual. The summation on the right is orient, that is to say, inverse the orientation of e i, then multiply with Ēi. Eqs.(7) can be implemented on D discrete manifold directly(see Fig.). Eq.(7) on rectangular gird is just the Yee scheme. (7) Figure : edge and face with direction In the absence of magnetic or dielectric materials, the relations are simple: D i = ε 0 E i B i = µ 0 H i, (8) where ε 0 and µ 0 are two universal constants, called the permittivity of free space and permeability of free space, respectively. With relations(8), Eqs.(7) can be rewritten into an explicit iterative scheme for TE wave. E n+ E n ǫ 0 µ 0 H n+ H n +J n e = Hn Hn e = En+ e +E n+ e +E n+ 3 e 3 P TE (9) 8
9 The symmetric TE wave scheme induced from Eqs.(6) can be written as follows. E n+ E n ǫ 0 +Je n = Hn Hn e µ 0 H n+ H n m = En+ e +E n+ e +E n+ 3 e 3 P +J n+ Schemes for TM wave If writing F n+ = H n+ dt D n F n = E n dt B n J = ( ρ e,j e dt) J = ( ρ m dt,j m ), TE (0) where H = H i e i is the discrete form on space, D = H i P i is the E P discrete form on space, E = D i P i is the dual of D on space, B = P B i e i is the dual of H on space, ρ e = ρ ei T i is the discrete 3 form E Tet on space, J e = J ei P i is the discrete form on space, ρ m dt is the discrete P form on time, J m = J mi E i is the discrete form on space, the discrete E Maxwell s equations can be rewritten as d s D n = ρ n e d s H n+ dt = d t D n +J n+ e dt d T s Bn = (ρ m dt) n d T s En dt = d T t Bn Jm n. () Proposition 3. If the initial condition satisfies the first and third equations in Eqs.(), the solution of the second and fourth equations in Eqs.() automatically satisfy Eqs.(). Proof.Because the dimension of spacetime is 3+ or +, therefore d T s (ρ m dt) = 0 d T t J m = 0 d s ρ e = 0 d t J e dt = 0, 9
10 and the continuity equations can be reduced to So we have d T t (ρ m dt) n d T s ( J n m ) = 0 d sj n+ e dt d t ρ n e = 0. d t d s D n d t ρ n e = d t ρ n e d s(d s H n+ dt J n+ e dt) = 0 d T t d T sb n d T t (ρ m dt) n = d T t (ρ m dt) n +d T s(d T se n+ dt+ Jm) n = 0. Now we show the scheme () on the product of D discrete space manifold and time. The second and fourth equations in Eqs.() based on Fig. are B n+ B n D n+ D n +J n m = En En e +J n+ e = Hn+ e +H n+ e +H n+ 3 e 3. P () With relations (8), Eqs.() can be rewritten into an explicit iterative scheme for TM wave. E n+ E n ǫ 0 +J n+ e = Hn+ e +H n+ e +H n+ 3 e 3 P H n+ H n µ 0 General schemes +J n m = En En e TM (3) For real world materials, the constitutive relations are not simple proportionalities, except approximately. The relations can usually still be written: D = εe B = µh, but ε and µ are not, in general, simple constants, but rather functions. With Ohm s law E = σ J, J m = σ m H, 0
11 where σ is the electrical conductivity and σ m is magnetic conductivity. The DEC schemes can be written as ǫ En+ E n +σ En+ +E n = Hn Hn e µ Hn+ H n ǫ En+ E n µ Hn+ H n +σ m H n+ +H n +σ En+ +E n H n+ +H n +σ m = En+ e +E n+ e +E n+ 3 e 3, P = Hn+ e +H n+ e +H n+ 3 e 3 P = En En e 4 Stability, convergence and accuracy Stability. TE TM The Courant-Friedrichs-Lewy condition is a necessary condition for convergence while solving certain partial differential equations numerically. Now, we find this condition for scheme (0). Condition for scheme (3) can be induced in the same way. First, we decompose DEC algorithm into temporal and spacial eigenvalue problems. The temporal eigenvalue problem: H n 0 t = ΛH n 0 It can approximated by difference equation Supposing H n+ 0 H n 0 +Hn 0 () = ΛH n 0. (4) H n+ 0 = H n 0 cos(n ) H n 0 = H n 0 cos(n ) and H n+ 0 = H n 0 sin(n ) H n 0 = H n 0 sin(n ),
12 and substituting those into Eq.(4), we obtain therefore cos(n )+cos(n ) () = Λ, sin(n )+sin(n ) () = Λ, 4 () Λ 0. This is the stabile condition for the temporal eigenvalue problem. The spacial eigenvalue problem: c H = ΛH It can be approximated by difference equation (5) based on Fig.. P 3 c ΛH 0 = l 3 l A0 (H A H 0 )+ l 34 l B0 (H B H 0 )+ l 45 l C0 (H C H 0 ) (5) Let Figure : Face and dual face H i = H 0 cos(cl 0i ) or H i = H 0 sin(cl 0i ), and substitute into Eq.(5) to obtain P 3 c Λ = l 3 l A0 (cos(cl 0A ) )+ l 34 l B0 (cos(cl 0B ) )+ l 45 l C0 (cos(cl 0C ) ) P 3 c Λ = l 3 l A0 (sin(cl 0A ) )+ l 34 l B0 (sin(cl 0B ) )+ l 45 l C0 (sin(cl 0C ) ).
13 So we have ( c l3 + l 34 + l ) 45 Λ 0. P 3 l A0 l B0 l C0 In order to keep the stability of scheme (3), we need ( () c l3 + l 34 + l ) 45, (6) P 3 l A0 l B0 l C0 or Convergence P Mim P3 P ( 3 c l3 + l 34 + l ) 45. l A0 l B0 l C0 By the definition of truncation error, the exact solution of Maxwell s equations satisfy the same relation as DEC scheme except for an additional term O(() + e ). This expresses the consistency, and so convergence for DEC scheme by Lax equivalence theorem (consistency + stability = convergence). Accuracy The derivative of Maxwell s equations is approximated by first order difference in schemes (0) and (3). Equivalently, H and E are approximated by linear interpolation functions. Consulting the definition about accuracy of finite volume method, we can also say that schemes (0) and (3) have first order spacial and temporal accuracy, and have second order spacial and temporal accuracy on rectangular grid with equivalent space and time steps. 3
14 5 Implementation The DEC algorithm of Maxwell s equations was implemented in C++ platform. The Fig.3 shows the flowchart of DEC schemes for Maxwell s equations. Figure 3 In the common practice, not every simulation step needs to be visualized, especially when the time step size is too small. Fig.4 exhibit Gaussian pluses 4
15 waveforms simulated by DEC. Figure 4: Simulation of Gaussian pluse on Stanford bunny by DEC Fig.5 exhibits two sources Gaussian pluses interference simulated by DEC. Our algorithm can simulate more complex situation on surface and 3D-space manifold. Figure 5: Simulation of Gaussian pluse interference on sphere by DEC References [] K.S. Yee, Numerical solution of inital boundary value problems involving Maxwell s equations in isotropic media. IEEE Trans. Ant. Prop. 4(3), (966). [] A. Bondeson, T. Rylander, P. Ingelstrom, Computational electromagnetics, Texts in Applied Mathematics, vol. 5. Springer, New York (005). [3] M. Clemens, T. Weiland, Magnetic field simulation using conformal FIT formulations. IEEE Trans. Magn. 38(), (00). 5
16 [4] P.W. Gross, P.R. Kotiuga, Electromagnetic theory and computation: a topological approach, Mathematical Sciences Research Institute Publications, vol. 48. Cambridge University Press, Cambridge, (004). [5] E. Hairer, C. Lubich, G. Wanner, Geometric numerical integration, Springer Series in Computational Mathematics, vol. 3. Springer-Verlag, Berlin (00). [6] A. Nicolet, F. Zolla, Y.O. Agha and S. Guenneau, Geometrical transformations and equivalent materials in computational electromagnetism, International Journal for Computation and Mathematics in Electrical and Electronic Engineering COMPEL, vol. 7, , (008). [7] A. Bossavit, L. Kettunen, : Yee-like schemes on a tetrahedral mesh, with diagonal lumping. Int. J. Numer. Modell. (-), 9 4 (999). [8] A. Bossavit, L. Kettunen, Yee-like schemes on staggered cellular grids: A synthesis between FIT and FEM approaches. IEEE Trans. Magn. 36(4), (000). [9] A. Stern, Computational Electromagnetism with Variational Integrators and Discrete Differential Forms. arxiv: [0] H. Whitney, Geometric integration theory. Princeton University Press, Princeton, (957). [] D.N. Arnold, R.S. Falk, R. Winther, Finite element exterior calculus, homological techniques, and applications. Acta Numer. 5, -55, (006). [] S. P. Novikov, Discrete connection on the triangulated manifold and difference linear equations arxiv: math-ph/ [3] B. Auchmann, S. Kurz, A geometrically defined discrete hodge operator on simplicial cells. IEEE Trans. Magn. 4(4), , (006). [4] A. Dimakis, F. Muller-Hoissen, Discrete Riemannian geometry. J Math Phys, 999, 40(3): [5] M. Desbrun, A.N. Hirani, M. Leok, J. E. Marsden, Discrete exterior calculus arxiv: math.dg/
17 [6] M. Meyer, M. Desbrun, P. Schröder, A.H. Barr, Discrete differential geometry operators for triangulated -manifolds. In InternationalWorkshop on Visualization and Mathematics, VisMath, (00). [7] J. M. Hyman, M. Shashkov, Natural discretizations for the divergence, gradient, and curl on logically rectangular grids. Comput. Math. Appl., 33(4):8-04, (997). [8] R. Hiptmair, Discrete Hodge operators, Numer. Math., 90():65-89, (00). [9] D.K. Wise, p-form electromagnetism on discrete spacetimes. Classical Quantum Gravity 3(7), (006). [0] M. Leok, Foundations of computational geometric mechanics. Ph.D. thesis, California Institute of Technology (004). 7
Variational Integrators for Maxwell s Equations with Sources
PIERS ONLINE, VOL. 4, NO. 7, 2008 711 Variational Integrators for Maxwell s Equations with Sources A. Stern 1, Y. Tong 1, 2, M. Desbrun 1, and J. E. Marsden 1 1 California Institute of Technology, USA
More informationarxiv: v3 [math.na] 27 May 2009
GEOMETRIC COMPUTATIONAL ELECTRODYNAMICS WITH VARIATIONAL INTEGRATORS AND DISCRETE DIFFERENTIAL FORMS arxiv:0707.4470v3 [math.na] 7 May 009 ARI STERN, YIYING TONG, MATHIEU DESBRUN, AND JERROLD E. MARSDEN
More informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
More informationNew Model Stability Criteria for Mixed Finite Elements
New Model Stability Criteria for Mixed Finite Elements Andrew Gillette Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin http://www.math.utexas.edu/users/agillette
More informationSpace-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories
21st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 2014. Space-Time Finite-Element Exterior Calculus and Variational Discretizations of Gauge Field Theories Joe Salamon
More informationDiscrete Exterior Calculus
Discrete Exterior Calculus and the Averaged Template Matching Equations Anil N. Hirani Jet Propulsion Laboratory California Institute of Technology hirani at caltech dot edu http://www.cs.caltech.edu/
More informationElectromagnetic Theory with Discrete Exterior Calculus
Progress In Electromagnetics Research, Vol. 159, 59 78, 2017 Electromagnetic Theory with Discrete Exterior Calculus Shu C. Chen 1 and Weng C. Chew 2, * Abstract A self-contained electromagnetic theory
More informationIs a dual mesh really necessary?
Is a dual mesh really necessary? Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICIAM, Vancouver, 2011. Joint work with Ari Stern, UCSD. 19 July 2011
More informationarxiv:physics/ v2 [physics.comp-ph] 14 Nov 2005
Geometric Finite Element Discretization of Maxwell Equations in Primal and Dual Spaces arxiv:physics/0503013v2 [physics.comp-ph] 14 Nov 2005 Bo He and F. L. Teixeira ElectroScience Laboratory and Department
More informationFrom point cloud data to the continuum model of geometry
From point cloud data to the continuum model of geometry J. Harrison University of California, Berkeley July 22, 2007 Continuum model of geometry By the continuum model of geometry we refer to smooth manifolds
More informationA simplicial gauge theory on spacetime
A simplicial gauge theory on spacetime Tore G. Halvorsen, NTNU, Trondheim Norway. Joint work with Snorre H. Christiansen, UiO. 1 Abstract In this talk I will introduce the variational form of the SU(N)
More informationGeometry of Electromagnetism and its Implications in Field and Wave Analysis
Geometry of Electromagnetism and its Implications in Field and Wave Analysis L. Kettunen, T. Tarhasaari, Tampere University of Technology, Institute of Electromagnetics, FIN-33101 Tampere, Finland Abstract
More informationGeometric finite element discretization of Maxwell equations in primal and dual spaces
Physics Letters A 349 (2006) 1 14 www.elsevier.com/locate/pla Geometric finite element discretization of Maxwell equations in primal and dual spaces Bo He,F.L.Teixeira Electroscience Laboratory and Department
More informationGlobal formulation of 3D magnetostatics using flux and gauged potentials
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2000; 00:1 6 [Version: 2000/01/19 v2.0] Global formulation of 3D magnetostatics using flux and gauged potentials M.
More informationDifferential Geometry: Discrete Exterior Calculus
Differential Geometry: Discrete Exterior Calculus [Discrete Exterior Calculus (Thesis). Hirani, 23] [Discrete Differential Forms for Computational Modeling. Desbrun et al., 25] [Build Your Own DEC at Home.
More informationDifferential Form Approach to the Analysis of Electromagnetic Cloaking and Masking
MOTL DRAFT 1 Differential Form Approach to the Analysis of Electromagnetic Cloaking and Masking F. L. Teixeira Abstract We bring attention to the relationship between (1) electromagnetic masking or cloaking
More informationNotes on Discrete Exterior Calculus
Notes on Discrete Exterior Calculus Andrew Gillette August 27, 2009 These notes are aimed at mathematics graduate students familiar with introductory algebraic and differential topology. The notes were
More informationFinite Difference Solution of Maxwell s Equations
Chapter 1 Finite Difference Solution of Maxwell s Equations 1.1 Maxwell s Equations The principles of electromagnetism have been deduced from experimental observations. These principles are Faraday s law,
More informationHigh Order Differential Form-Based Elements for the Computation of Electromagnetic Field
1472 IEEE TRANSACTIONS ON MAGNETICS, VOL 36, NO 4, JULY 2000 High Order Differential Form-Based Elements for the Computation of Electromagnetic Field Z Ren, Senior Member, IEEE, and N Ida, Senior Member,
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationHOMOGENEOUS AND INHOMOGENEOUS MAXWELL S EQUATIONS IN TERMS OF HODGE STAR OPERATOR
GANIT J. Bangladesh Math. Soc. (ISSN 166-3694) 37 (217) 15-27 HOMOGENEOUS AND INHOMOGENEOUS MAXWELL S EQUATIONS IN TERMS OF HODGE STAR OPERATOR Zakir Hossine 1,* and Md. Showkat Ali 2 1 Department of Mathematics,
More informationDISCRETE EXTERIOR CALCULUS
DISCRETE EXTERIOR CALCULUS MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN Abstract. We present a theory and applications of discrete exterior calculus on simplicial complexes of arbitrary
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationProjective geometry and spacetime structure. David Delphenich Bethany College Lindsborg, KS USA
Projective geometry and spacetime structure David Delphenich Bethany College Lindsborg, KS USA delphenichd@bethanylb.edu Affine geometry In affine geometry the basic objects are points in a space A n on
More informationIn collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract
73 Chapter 3 Discrete Exterior Calculus In collaboration with Mathieu Desbrun, Anil N. Hirani, and Jerrold E. Marsden. Abstract We present a theory and applications of discrete exterior calculus on simplicial
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationGeometric Discretization of Lagrangian Mechanics and Field Theories
Geometric Discretization of Lagrangian Mechanics and Field Theories Thesis by Ari Stern In Partial Fulfullment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology
More informationGlobal formulation for 3D magneto-static using flux and gauged potential approaches
Global formulation for 3D magneto-static using flux and gauged potential approaches Maurizio Repetto, Francesco Trevisan 1 Dep. of Ingegneria Elettrica, Industriale, Politecnico di Torino C. Duca degli
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationExplicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems. Abstract
Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems Jianyuan Xiao, 1, 2 Hong Qin, 1, 3, Jian Liu, 1, 2 Yang He, 1, 2 Ruili Zhang, 1, 2 and Yajuan Sun 4 arxiv:1510.06972v2
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationSPECIAL RELATIVITY AND ELECTROMAGNETISM
SPECIAL RELATIVITY AND ELECTROMAGNETISM MATH 460, SECTION 500 The following problems (composed by Professor P.B. Yasskin) will lead you through the construction of the theory of electromagnetism in special
More informationA Discrete Vector Calculus in Tensor Grids
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS Vol. (00,, pp. 44 c 00 Institute of Mathematics, National Academy of Sciences A Discrete Vector Calculus in Tensor Grids Nicolas Robidoux Stanly Steinberg Abstract
More informationA Generalization for Stable Mixed Finite Elements
A Generalization for Stable Mixed Finite Elements Andrew Gillette joint work with Chandrajit Bajaj Department of Mathematics Institute of Computational Engineering and Sciences University of Texas at Austin,
More informationLattice electromagnetic theory from a topological viewpoint a
JOURNAL OF MATHEMATICAL PHYSICS VOLUME 40, NUMBER 1 JANUARY 1999 Lattice electromagnetic theory from a topological viewpoint a F. L. Teixeira b) and W. C. Chew Center for Computational Electromagnetics,
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationOverview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title A doubled discretization of Abelian Chern-Simons theory Author(s) Adams, David H. Citation Adams, D.
More informationNumerical Analysis of Electromagnetic Fields in Multiscale Model
Commun. Theor. Phys. 63 (205) 505 509 Vol. 63, No. 4, April, 205 Numerical Analysis of Electromagnetic Fields in Multiscale Model MA Ji ( ), FANG Guang-You (ྠ), and JI Yi-Cai (Π) Key Laboratory of Electromagnetic
More informationr r 1 r r 1 2 = q 1 p = qd and it points from the negative charge to the positive charge.
MP204, Important Equations page 1 Below is a list of important equations that we meet in our study of Electromagnetism in the MP204 module. For your exam, you are expected to understand all of these, and
More informationClassical Physics. SpaceTime Algebra
Classical Physics with SpaceTime Algebra David Hestenes Arizona State University x x(τ ) x 0 Santalo 2016 Objectives of this talk To introduce SpaceTime Algebra (STA) as a unified, coordinate-free mathematical
More informationCosmology on Simplicial Complexes
Gravitation and Regge Calculus Astro Coee, Frankfurt, April 2015 Outline Gravitation and Regge Calculus 1 Gravitation and Regge Calculus Foundations of General Relativity Geometric Structure of Regge Calculus
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationMath. 460, Sec. 500 Fall, Special Relativity and Electromagnetism
Math. 460, Sec. 500 Fall, 2011 Special Relativity and Electromagnetism The following problems (composed by Professor P. B. Yasskin) will lead you through the construction of the theory of electromagnetism
More informationChapter Three: Propagation of light waves
Chapter Three Propagation of Light Waves CHAPTER OUTLINE 3.1 Maxwell s Equations 3.2 Physical Significance of Maxwell s Equations 3.3 Properties of Electromagnetic Waves 3.4 Constitutive Relations 3.5
More informationarxiv: v1 [physics.class-ph] 8 Apr 2019
Representation Independent Boundary Conditions for a Piecewise-Homogeneous Linear Magneto-dielectric Medium arxiv:1904.04679v1 [physics.class-ph] 8 Apr 019 Michael E. Crenshaw 1 Charles M. Bowden Research
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationPHYSICAL REVIEW B 71,
Coupling of electromagnetic waves and superlattice vibrations in a piezomagnetic superlattice: Creation of a polariton through the piezomagnetic effect H. Liu, S. N. Zhu, Z. G. Dong, Y. Y. Zhu, Y. F. Chen,
More informationInvariant differential operators on the sphere
Third CoE Talk at the University of Tokyo p. 1/21 Invariant differential operators on the sphere Michael Eastwood Australian National University Third CoE Talk at the University of Tokyo p. 2/21 The round
More informationA posteriori error estimates in FEEC for the de Rham complex
A posteriori error estimates in FEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationHow to recognise a conformally Einstein metric?
How to recognise a conformally Einstein metric? Maciej Dunajski Department of Applied Mathematics and Theoretical Physics University of Cambridge MD, Paul Tod arxiv:1304.7772., Comm. Math. Phys. (2014).
More informationJustin Solomon MIT, Spring 2017
Fisher et al. Design of Tangent Vector Fields (SIGGRAPH 2007) Justin Solomon MIT, Spring 2017 Original version from Stanford CS 468, spring 2013 (Butscher & Solomon) Divergence Theorem Green s Theorem
More informationE. YABLONOVITCH photonic crystals by using level set methods
Appl. Phys. B 81, 235 244 (2005) DOI: 10.1007/s00340-005-1877-3 Applied Physics B Lasers and Optics C.Y. KAO 1, Maximizing band gaps in two-dimensional S. OSHER 2 E. YABLONOVITCH photonic crystals by using
More informationComparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 1014 1020 Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses Jian-Bo Chen, Yan Shen,
More informationRICCI SOLITONS ON COMPACT KAHLER SURFACES. Thomas Ivey
RICCI SOLITONS ON COMPACT KAHLER SURFACES Thomas Ivey Abstract. We classify the Kähler metrics on compact manifolds of complex dimension two that are solitons for the constant-volume Ricci flow, assuming
More informationGravitation: Tensor Calculus
An Introduction to General Relativity Center for Relativistic Astrophysics School of Physics Georgia Institute of Technology Notes based on textbook: Spacetime and Geometry by S.M. Carroll Spring 2013
More informationModern Geometric Structures and Fields
Modern Geometric Structures and Fields S. P. Novikov I.A.TaJmanov Translated by Dmitry Chibisov Graduate Studies in Mathematics Volume 71 American Mathematical Society Providence, Rhode Island Preface
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationInside-out electromagnetic cloaking
Inside-out electromagnetic cloaking Nina A. Zharova 1,2, Ilya V. Shadrivov 1, and Yuri S. Kivshar 1 1 Nonlinear Physics Center, Research School of Physical Sciences and Engineering, Australian National
More informationarxiv: v1 [physics.plasm-ph] 15 Sep 2013
Symplectic Integration of Magnetic Systems arxiv:309.373v [physics.plasm-ph] 5 Sep 03 Abstract Stephen D. Webb Tech-X Corporation, 56 Arapahoe Ave., Boulder, CO 80303 Dependable numerical results from
More informationarxiv: v1 [physics.optics] 19 Jun 2008
arxiv:0806.3231v1 [physics.optics] 19 Jun 2008 Necessary and sufficient conditions for reflectionless transformation media in an isotropic and homogenous background Wei Yan, Min Yan, Min Qiu Laboratory
More informationFinite Formulation of Electromagnetic Field
Finite Formulation of Electromagnetic Field Enzo ON Dept.Civil Engin., Univ. of rieste, Piazzale Europa 1, 34127 rieste, talia. e-mail: tonti@univ.trieste.it October 16, 2000 Abstract he paper shows that
More informationarxiv: v1 [physics.optics] 30 Mar 2010
Analytical vectorial structure of non-paraxial four-petal Gaussian beams in the far field Xuewen Long a,b, Keqing Lu a, Yuhong Zhang a,b, Jianbang Guo a,b, and Kehao Li a,b a State Key Laboratory of Transient
More informationOn Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material Surfaces
Vol. 114 2008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology On Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More informationarxiv:hep-th/ v3 21 Jul 1997
CERN-TH/96-366 hep-th/9612191 Classical Duality from Dimensional Reduction of Self Dual 4-form Maxwell Theory in 10 dimensions arxiv:hep-th/9612191v3 21 Jul 1997 David Berman Theory Division, CERN, CH
More informationt, H = 0, E = H E = 4πρ, H df = 0, δf = 4πJ.
Lecture 3 Cohomologies, curvatures Maxwell equations The Maxwell equations for electromagnetic fields are expressed as E = H t, H = 0, E = 4πρ, H E t = 4π j. These equations can be simplified if we use
More informationarxiv:math-ph/ v1 10 Oct 2005
DISCRETE MODELS O THE SEL-DUAL AND ANTI-SEL-DUAL EQUATIONS arxiv:math-ph/0510041v1 10 Oct 2005 Volodymyr Sushch Lviv, Uraine; Koszalin, Poland Abstract. In the case of a gauge-invariant discrete model
More informationA simple and compact approach to hydrodynamic using geometric algebra. Abstract
A simple and compact approach to hydrodynamic using geometric algebra Xiong Wang (a) Center for Chaos and Complex Networks (b) Department of Electronic Engineering, City University of Hong Kong, Hong Kong
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationLagrangian Description for Particle Interpretations of Quantum Mechanics Single-Particle Case
Lagrangian Description for Particle Interpretations of Quantum Mechanics Single-Particle Case Roderick I. Sutherland Centre for Time, University of Sydney, NSW 26 Australia rod.sutherland@sydney.edu.au
More informationCLASS OF BI-QUADRATIC (BQ) ELECTROMAGNETIC MEDIA
Progress In Electromagnetics Research B, Vol. 7, 281 297, 2008 CLASS OF BI-QUADRATIC (BQ) ELECTROMAGNETIC MEDIA I. V. Lindell Electromagnetics Group Department of Radio Science and Engineering Helsinki
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationarxiv: v2 [math.na] 8 Sep 2015
OMPLEXES OF DISRETE DISTRIBUTIONAL DIFFERENTIAL FORMS AND THEIR HOMOLOGY THEORY MARTIN WERNER LIHT arxiv:1410.6354v2 [math.na] 8 Sep 2015 Abstract. omplexes of discrete distributional differential forms
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationarxiv: v1 [math.dg] 10 Dec 2015
Preprint Curvature corrected estimates for geodesic arc-length. arxiv:1512.03461v1 [math.dg] 10 Dec 2015 Leo Brewin School of Mathematical Sciences Monash University, 3800 Australia 19-Oct-2015 Abstract
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationEffective theory of quadratic degeneracies
Effective theory of quadratic degeneracies Y. D. Chong,* Xiao-Gang Wen, and Marin Soljačić Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA Received 28
More informationNew Concept Conformal Antennas Utilizing Metamaterial and Transformation Optics
New Concept Conformal Antennas Utilizing Metamaterial and Transformation Optics The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation
More informationThe Coulomb And Ampère-Maxwell Laws In The Schwarzschild Metric, A Classical Calculation Of The Eddington Effect From The Evans Unified Field Theory
Chapter 6 The Coulomb And Ampère-Maxwell Laws In The Schwarzschild Metric, A Classical Calculation Of The Eddington Effect From The Evans Unified Field Theory by Myron W. Evans, Alpha Foundation s Institutute
More information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationFOUNDATIONS OF COMPUTATIONAL GEOMETRIC MECHANICS
FOUNDATIONS OF COMPUTATIONAL GEOMETRIC MECHANICS MELVIN LEOK Received the SIAM Student Paper Prize and the Leslie Fox Prize in Numerical Analysis (second prize). Abstract. Geometric mechanics involves
More informationLecture 13 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 13 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Covariant Geometry - We would like to develop a mathematical framework
More informationOverview in Images. 5 nm
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) S. Lin et al, Nature, vol. 394, p. 51-3,
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationHigh-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials
Progress In Electromagnetics Research Letters, Vol. 77, 103 107, 01 High-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials Wei-Jun Chen
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationA Symmetric and Low-Frequency Stable Potential Formulation for the Finite-Element Simulation of Electromagnetic Fields
A Symmetric and Low-Frequency Stable Potential Formulation for the Finite-Element Simulation of Electromagnetic Fields Martin Jochum, Ortwin Farle, and Romanus Dyczij-Edlinger Abstract A low-frequency
More informationTowards Computational Electromagnetics in Spacetime
Tampereen teknillinen yliopisto. Julkaisu 953 Tampere University of Technology. Publication 953 Janne Keränen Towards Computational Electromagnetics in Spacetime Thesis for the degree of Doctor of Science
More informationAposteriorierrorestimatesinFEEC for the de Rham complex
AposteriorierrorestimatesinFEEC for the de Rham complex Alan Demlow Texas A&M University joint work with Anil Hirani University of Illinois Urbana-Champaign Partially supported by NSF DMS-1016094 and a
More informationGRAPH QUANTUM MECHANICS
GRAPH QUANTUM MECHANICS PAVEL MNEV Abstract. We discuss the problem of counting paths going along the edges of a graph as a toy model for Feynman s path integral in quantum mechanics. Let Γ be a graph.
More informationAn abstract Hodge Dirac operator and its stable discretization
An abstract Hodge Dirac operator and its stable discretization Paul Leopardi Mathematical Sciences Institute, Australian National University. For presentation at ICCA, Tartu,. Joint work with Ari Stern,
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More informationSupersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College, Cambridge. October 1989
Supersymmetric Quantum Mechanics and Geometry by Nicholas Mee of Trinity College Cambridge October 1989 Preface This dissertation is the result of my own individual effort except where reference is explicitly
More informationarxiv:physics/ v1 19 Feb 2004
Finite Difference Time Domain (FDTD) Simulations of Electromagnetic Wave Propagation Using a Spreadsheet David W. Ward and Keith A. Nelson Department of Chemistry Massachusetts Institute of Technology,
More informationDiscrete Differential Geometry. Peter Schröder with help from Eitan Grinspun, Mathieu Desbrun and the rest of the DDG crew.
Discrete Differential Geometry Peter Schröder with help from Eitan Grinspun, Mathieu Desbrun and the rest of the DDG crew ps@cs.caltech.edu A Bit of History Geometry is the key! studied for centuries Hermann
More informationWilson and Polyakov loops of gravitational gauge fields in Rindler space
Wilson and Polyakov loops of gravitational gauge fields in Rindler space E.Koorambas 8A Chatzikosta, 5 Ampelokipi, Athens, Greece E-mail:elias.koor@gmail.com (August 4, 04) Abstract: We will study the
More informationA Review of Basic Electromagnetic Theories
A Review of Basic Electromagnetic Theories Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820)
More information