arxiv: v1 [math-ph] 8 Jun 2011

Transcription

1 Electromagnetism on Anisotropic Fractals arxiv:06.49v [math-ph] 8 Jun 20 Martin Ostoja-Starzewski Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign Urbana, IL 680, U.S.A. October 25, 208 Abstract e derive basic equations of electromagnetic fields in fractal media which are specified by three indepedent fractal dimensions α i in the respective directions x i i =,2,3 of the Cartesian space in which the fractal is embedded. To grasp the generally anisotropic structure of a fractal, we employ the product measure, so that the global forms of governing equations may be cast in forms involving conventional integer-order integrals, while the local forms are expressed through partial differential equations with derivatives of integer order but containing coefficients involving the α i s. First, a formulation based on product measures is shown to satisfy the four basic identities of vector calculus. This allows a generalization of the Green-Gauss and Stokes theorems as well as the charge conservation equation on anisotropic fractals. Then, pursuing the conceptual approach, we derive the Faraday and Ampère laws for such fractal media, which, along with two auxiliary null-divergence conditions, effectively give the modified Maxwell equations. Proceeding on a separate track, we employ a variational principle for electromagnetic fields, appropriately adapted to fractal media, to independently derive the same forms of these two laws. It is next found that the parabolic for a conducting medium and the hyperbolic for a dielectric medium equations involve modified gradient operators, while the Poynting vector has the same form as in the non-fractal case. Finally, Maxwell s electromagnetic stress tensor is reformulated for fractal systems. In all the cases, the derived equations for fractal media depend explicitly on fractal dimensions and reduce to conventional forms for continuous media with Euclidean geometries upon setting the dimensions to integers. Introduction In a recent article Li & Ostoja-Starzewski, 2009 we formulated a continuum mechanics of anisotropic fractal elastic media using two independent approaches: a mechanical Newtonian one and a variational Lagrangian-Hamiltonian one. Proceeding separately on both paths, we derived the identical equations of

2 wave motion in one-, two- and three-dimensional media, and this provided a verification of the consistency of the product measure capable of grasping the anisotropic fractal media as well as the verification of the final equations. The question arises whether one can use the same strategy to derive the field equations of electromagnetism in fractal media, be they isotropic or anisotropic. Various studies generalizing the classical, quantum or relativistic field theories to fractal spaces Stillinger, 977; Svozil, 987; Palmer & Stavrinou, 2004; Nottale, 200 as well as the electromagnetism to fractal media appeared over the past two decades Tarasov, 2006, 200. They have been formulated either through a conceptual path in the case of electromagnetism, by developing the Gauss, Faraday and Ampère laws, so as to obtain the Maxwell equations or through variational principles. However, a derivation obtained on one path say, conceptual has never been verified through another, separate derivation respectively, variational. In the following we proceed under the premise that the conceptual derivation of Maxwell s equations for anisotropic fractal media should yield the same results as the derivation based on a variational principle, both derivations being based on the same dimensional regularization. To grasp the generally anisotropic structure of a fractal, we employ the recently introduced product measure, so that the global forms of governing equations may be cast in forms involving conventional integer-order integrals, while the local forms are expressed through partial differential equations with derivatives of integer order but containing coefficients involving fractal dimensions α i in three respective directions x i, i =,2,3. e first show that the formulation based on product measures satisfies the four basic identities of vector calculus. This readily leads to a generalization of the Green-Gauss and Stokes theorems involving the fractal gradient and curl operators, and hence to the charge conservation equation on anisotropic fractals. Then, pursuing the conceptual approach, we derive the Faraday and Ampère laws for such fractal media, which, along with two auxiliary null-divergence conditions, effectively give the modified Maxwell equations. Proceeding on a separate track, we employ a variational principle for electromagnetic fields, appropriately adapted to fractal media, to independently derive the same forms of these two laws. e next examine whether the parabolic for a conducting medium and the hyperbolic for a dielectric medium equations involve fractal gradient operators, and also whether the Poynting vector has the same form as in the non-fractal case. Finally, Maxwell s electromagnetic stress tensor is reformulated for fractal systems. In all the cases, the derived equations for fractal media depend explicitly on fractal dimensions and reduce to conventional forms for continuous media with Euclidean geometries upon setting the dimensions to integers. In order to make the presentation clear, in some places we give the tensorial relations in the symbolic as well as the index notations. 2

3 2 Background 2. Anisotropic Fractal Distributions The distribution of charges in an electrically and magnetically conducting isotropic fractal medium embedded in the Euclidea space E 3 obeys a power law QR R D, D < 3, where D is the fractal dimension and R is the spatial resolution. However, in an anisotropic fractal medium generally denoted by we write Qx,x 2,x 3 x α x 2 α 2 x 3 α 3, 2 whereeachα k playstheroleofafractaldimensionalongeachaxisx k i =,2,3 of the Cartesian space in which the fractal is embedded. In other directions, the fractal dimension is not necessarily the sum of projected fractal dimensions, however, as noted by Falconer 2003, Many fractals encountered in practice are not actually products, but are product-like. Hence, we expect the equality between the fractal dimension D of the total charge and α +α 2 +α 3 to hold for fractals encountered in practice. hen D 3, with each α i, the conventional concept of charge distribution is recovered. By analogy to the formulation of continuum mechanics on anisotropic fractals Li & Ostoja-Starzewski, 2009, we now express 2. through a product measure Q = ρx,x 2,x 3 dl α x dl α2 x 2 dl α3 x 3, 3 where ρ is the density of distribution. In 2.3, the length measure in each coordinate is given in terms of the transformation coefficients c k by dl αk x k = c k α k,x k dx k, k =,2,3 no sum, 4 where l k is the total length integral interval along x k and l k0 is the characteristic length in the given direction, like the mean pore size. The relation 2.4 impliesthat the infinitesimalfractalvolumeelement, dv D, isrelatedtothe cubic volume element, dv 3 = dx dx 2 dx 3 by dv D = dl α x dl α2 x 2 dl α3 x 3 = c c2 c3 dx dx 2 dx 3 = c 3 dv 3, 5 with c 3 = c c2 c3. Similarly, the infinitesimal fractal surface element, ds d, is related to the planar surface element, ds k d = dx i dx j with the normal vector along x k, according to ds k d = c c ds 2, with c k 2 = c i cj = c 3 /c k, i j, i,j k. 6 3

4 Thesumd k = α i +α j,i j,i,j k, isthefractaldimensionofthesurfaces k d alongthediagonals x i = x j ins k. Again, thisequalityisnotnecessarilytrue d elsewhere, but is expected to hold for fractals encountered in practice Falconer, The transformation coefficients c k showing up in 2.4 can be represented in terms of the modified Riemann-Liouville fractional integral of Jumarie 2005, 2009, thus αk c k lk x k = α k, k =,2,3, no sum, 7 l k0 although further developments will not depend on this explicit form. 2.2 Vector calculus on anisotropic fractals Following Li & Ostoja-Starzewski, 2009, we will extensively employ the fractal derivative fractal gradient operator grad D φ = e k D k φ or D k φ = c k φ x k no sum on k. 8 where e k are base vectors. Hence, the fractal divergence of a vector field divf = D f or D k f k = c k Note that this leads to a fractal curl operator of a vector field curlf = D f or D k f i = e jki f k x k. 9 c k f i x k. 0 e observe that the four fundamental identities of the conventional vector calculus carry over to these new operators: i The divergence of the curl of a vector field f: div curlf = D k f = [ f x j e jki c k f i x k ii The curl of the gradient of a scalar field φ: [ curl gradφ = D j D k φ = e ijk x j ] = e jki ] φ x k c k f i x j x k = 0. f i = e c k jki = 0. c k x j x k 2 In both cases above we can pull /c k in front of the gradient because the coefficient c k is independent of x j. 4

5 iii The divergence of the gradient of a scalar field φ is written in terms of the fractal gradient as [ ] [ ] div gradφ = D j D k φ = φ = φ, j, x j x j j 3 which gives an explicit form of the fractal Laplacian. iv The curl of the curl operating on a vector field f: [ ] curl curlf = D r D r f = e f i prj e c r jki x r c k x k [ ] [ ] f i = e prj e jki = δ c r x r c k kp δ ir δ kr δ ip f i x k c r x r c k x k [ ] [ ] = f r f p = D c r x r c p x p c r x r c r p D x r f D p D r f r 4 As further background, we now extend two integral theorems to generally anisotropic fractals. 2.3 Stokes and Green-Gauss theorems for anisotropic fractals Stokes theorem. e begin with A D f nds d and proceed in the index notation: n k D j f i ds d n k e kji f i, j ds d = n k e kji f i, j c k 2 ds 2 = n k e kji f i, j c i ds 2 = n k e kji f i c i, j ds 2 5 In the second equality above we employed the dimensional regularization, followed by taking note of 2.6, including the fact that c i is independent of x j. Next, by the Stokes theorem in the Euclidean space E 3, we have n k e kji f i c i, j ds 2 = f i c i dl i = A f i dl id 6 where, again, we used 2.4. Thus, we arrive at the Stokes theorem for fractals n curlf ds d = f dl α 7 A Formally, the above agrees with Tarasov s 2006 result in the case of isotropy, although our curl operator also grasps anisotropy. l 5

6 Green-Gauss theorem. e begin with A f nds d and proceed in the index notation: = f i c k 2 f i n k ds d =, k c 3 dv D = f i n k c k 2 ds 2 = f i, k c k 2 c 3 dv D = f i c k 2, k dv 2 f i, k c k dv D, 8 whereby in the second equality we employed the dimensional regularization, followed by taking note of 2.5, then followed by employing the Green-Gauss theorem in E 3, and in turn followed by noting that c k 2 is independent of x k. Thus, we arrive at the Green-Gauss theorem for fractals as f n ds d = D f dv D 9 The same comment as that following 2.7 applies here. 2.4 Charge conservation on anisotropic fractals Let J be the current density, n be the direction normal to the surface, and η be the charge density on the fractal. Then we have J n ds d = η dv D 20 On account of 2.8, we obtain global form D J dv D = d η dv D dt local form D J = 2 η The fact that D automatically appears in the above will be exploited below. Focusing on a linear electromagnetic response, with σ being the electric conductivity tensor, Ohm s law for an anisotropic medium reads J = σ E or J i = σ ij E j 22 the isotropic form following for σ ij = δ ij σ. The fact that the constitutive form above carries over directly from the non-fractal case is motivated by the analogy to elastic media, where Hooke s law is unchanged when going from nonfractal to fractal media Li & Ostoja-Starzewski, 2009; that result ensured the consistency of the Newtonian and Hamiltonian approaches to the derivation of governing equations - recall the beginning of the Introduction. 6

7 3 Formulation of Maxwell s equations via Stokes and Green-Gauss theorems 3. Faraday s law First, note that the Faraday law holds for a fractal embedded in E 3 d B nda d = E dl D dt A l On account of the Stokes theorem for fractals 2.7, we obtain d k n k da d = dt AB d kji dt Ae D j E i n k da d d e kji E i, j n k da d 2 dt A which, by localization leads, to 0 = B k +e kji D j E i B k +e kji or E i, j 3 0 = B+ D E 3.2 Ampère s law First, note that the Ampère law holds for a fractal embedded in E 3 C nda d = H dl D = n D HdA d 4 Noting we obtain A A l C = J+ D which, again by focusing on the integrand, yields l 5 J+ D D H da d = 0 6 J+ D D H = 0 7 On account ofthe speed oflight being related to the dielectric ǫ 0 and magnetic µ 0 constants by the well known relation c = / ǫ 0 µ 0, the above is equivalent to J+ E ǫ 0 c2 D B = 0 8 7

8 The equations 3.7 and 3.8 are subject to the constraints Gauss law for magnetism D B = 0 Gauss law D E = 0 9 where the presence of the fractal rather than the classical divergence operator D has clearly been suggested by Derivation from variational principle The variational principle for the Maxwell equations in conventionalnon-fractal setting E 3 is Seliger & hitham, 968 E = 0, B = 0, B + E = 0, E c2 B+ ǫ 0 J = 0, δ LdV 3 dt = 0, with [ L = ǫ 0 2 c2 B i B i 2 E Ei ie i +A i c2 e ijk B k, j + ] J i +χe i, i. ǫ 0 2 Here the vector potential A and the scalar potential χ are the Lagrange multipliers, while the integrand of 4.2 is slightly extended to include the current density J. By varying 4.2 with respect to E, B, A, and χ, gives, respectively, δe : E = A/ χ δb : B = A δa : E/ c 2 B+ ǫ 0 J = 0 3 δχ : E = 0 It is clear that 4.3,2 satisfy the first pair of Maxwell s equations 4. identically. Henceforth, our goal is to obtain the Maxwell equations for an anisotropic fractal, by using the same type of a variational approach. Forafractal embedded ine 3, thevariationalprincipleforelectromagnetic fields is written as t2 δ LdV D dt = 0 4 t 8

9 with L given by Varying 4.4 with respect to E, we obtain E = A/ D χc or E i = A i / c i χc i, i D i A i c i 5 while varying it with respect to B, we obtain the Gauss law for magnetism on fractals B = D Ac or B k = e ijk A i, j D j A i In both cases, the fractal gradient operator, D, shows up due to the presence of dv D in 4.4 ; recall 2.5. The next step is to verify whether substituting 4.5 into 3.3 would give a zero vector. Thus, 0 = B k +e kji E i, j = B k /+e kji c i [ A i c i χc i, i ], j [ = B k /+e kji ] A c i i c i, j χc i, ij [ ] = B k e kji A c i i c i, j e kji χc i c i, ij [ ] B k D j A i c i D j χc i, i 7 The first term in the square brackets in the last line of the above agrees with 4.6, while the last term vanishes, providing the fractal is isotropic c i =, i j, in which case 4.7 is satisfied identically. Thus, for an anisotropic fractal, there appears a source/disturbance an observation consistent with Tarasov s formulation, though his explicit form differs from ours. Also, substituting 4.6 into 3.9 yields 0 = e ijk c k A i, jk D k D j A i which vanishes providing the fractal is isotropic c i =, i j. Varying 4.4 with respect to A, we obtain E/ c 2 D B+ ǫ 0 J = 0 or E i / c 2 D j B k which agrees with 3.8 perfectly. Varying 4.4 with respect to χ, we obtain the Gauss law D E = 0 or D i E i = 9 c i ǫ 0 J i = 0 9 E i, i = 0 0

10 which perfectly agrees with Thus, we have derived by an entirely independent route the same set of Maxwell s equations modified to fractals as those in Section 2. In analogy to Li & Ostoja-Starzewski, 2009 which treated the continuum mechanics of fractal elastic media independently by Newtonian and Lagrangian approaches, this shows that the formulation of continuum physics equations based on the product measures is consistent. In the case of isotropy, our equations do not reduce to those of Tarasov Second order differential equations of electromagnetism In Gaussian units Eqns 3.9, 3.3 and 3.8 are D E = 0, c D B = 0, c B+ D E = 0, E D B+ 4π c J = 0, Taking derivatives gives the second-order Maxwell s equations for a fractal c p D D E D D E = c 2 2 E 2 t + 4π c 2 J [ ] E k, p /c p, p [ c k or E k, l /c l ], p = c 2 2 E k 2 + 4π c Now, consider two special cases: Conductor. ith Ohm s law for anisotropic media 2.22, the electrodynamics equations simplify to J k 2 D D = 0, D B = 0, c H+ D E = 0, σ c E D H = 0, 3 0

11 Upon using the fractal curl operation, we obtain a parabolic type equation D D H = σ c 2 H c p [ H k, p /c p or ], p = σ km H m c 2 Dielectric. From 5.3 we find a hyperbolic type equation 4 D D E = σ c 2 2 E 2 t c p or [ ] E k, p /c p, p = σ km 2 E m c Electromagnetic Energy and Stress Pontying vector. Now, we return to 3.3 and 3.7. Multiplying the first of these equations by H, and the second one by E, leads to D E H = B H Subtracting 6.2 from 6., on account of we obtain D H E = J E+ D E 2 H D E+E D H = D E H, 3 D E H+J E = D E B H or 4 e ijk E k H i, j +J i E i = E i D i H i B i Integrating 6.4 over the fractal s volume and applying the Green-Gauss theorem, yields [ G n ds d + J E dv D = D E ] B H dv D 5 S

12 where G = E H 6 is identified as the Poynting vector, which is seen to have the same form as in non-fractal media. Now, the electric and magnetic force densities are f E = ρe+p D E f M = J M B+M B D = c JE B+M B D 7 This has to be accompanied by an expression of the Ampère law in electrostatic units D H = 4π c J+ cḋ = 4π J+ 8 c cḋ 7 Conclusion e determine the equations governing electromagnetic fields in generally anisotropic fractal media using the dimensional regularization together with the recently introduced product measure allowing an independent characterization of anisotropy in three Cartesian directions of the Euclidean space in which an anisotropic fractal is embedded. This measure also gives rise to fractal gradient, divergence and curl operators, which are shown to satisfy the four fundamental identities of the vector calculus. Next, the conservation of the electric charge on a fractal shows that the fractal divergence has to be used. ith this Ansatz, in the first place, we obtain Maxwell equations modified to generally anisotropic fractals using two independent approaches: a conceptual one involving generalized Faraday and Ampère laws, and the one directly based on a variational principle for electromagnetic fields. In both cases the resulting equations are the same, thereby providing a self-consistent verification of our derivations. Just as in the previous works of Tarasov, we find that the presence of anisotropy in the fractal structure leads to a source/disturbance as a result of generally unequal fractal dimensions in various directions. However, our modified Mawell equations do not coincide with those of Tarasov 2006, 200. The modifications due to fractal geometry carry over to the parabolic for conductor and hyperbolic for dielectric equations, while the Poynting vector is found to have the same form as in the non-fractal case. Finally, the electromagnetic stress tensor is reformulated for fractal systems. Overall, all the derived relationsdepend explicitlyonthreefractaldimensionsα i inthe respectivecartesian directionsx i, i =,2,3, as wellas the spatial resolution; upon setting all α i =, these relations reduce to conventional forms of Maxwell equations. References Falconer, K.: Fractal Geometry: Mathematical Foundations and Applications, J. iley

13 Jumarie, G.:. On the representation of fractional Brownian motion as an integral with respect to dt a. Appl. Math. Lett. 8, Jumarie, G.: Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functions. Appl. Math. Lett. 223, Nottale, L. 200, Scale relativity and fractal space-time: theory and applications, Found. Sci. 5, Li, J. & Ostoja-Starzewski, M.: Fractal solids, product measures and fractional wave equations, Proc. R. Soc. A 465, doi:0.098/rspa ; Errata 200 doi:0.098/rspa Palmer, C. & Stavrinou, P.N. 2004, Equations of motion in a non-integerdimensional space, J. Phys. A: Math. Gen. 37, Seliger, R.L. & hitham, G.B. 968, Variational principles in continuum mechanics. Proc. R. Soc. A 305, -25. Stillinger, F.H. 977, Axiomatic basis for spaces with noninteger dimension, J. Math. Phys. 86, Svozil, K. 987, Quantum field theory on fractal spacetime: a new regularisation method, J. Phys. A: Math. Gen. 20, Tarasov, V.E Electromagnetic fields on fractals, Mod. Phys. Lett. A 220, Tarasov, V.E. 200 Fractional Dynamics, Springer. 3