OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 5 — Mar. 2, 2009
  • pp: 3581–3586
« Show journal navigation

Invisible cloak design with controlled constitutive parameters and arbitrary shaped boundaries through Helmholtz’s equation

Xi Chen, Yunqi Fu, and Naichang Yuan  »View Author Affiliations


Optics Express, Vol. 17, Issue 5, pp. 3581-3586 (2009)
http://dx.doi.org/10.1364/OE.17.003581


View Full Text Article

Acrobat PDF (240 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

An approach to design an invisible cloak with controlled constitutive parameters and arbitrary shaped boundaries is presented. Helmholtz’s equation is adopted to establish a mapping between original and transformed coordinates inside the cloak. Then the constitutive parameters are obtained by the established mapping. The analytical solution of a regular cloak and the numerical solution of an irregular cloak both verify that that our method will guide electromagnetic wave efficiently and control the constitutive parameters of the cloak conveniently. It has great significance in realizing a cloak practically.

© 2009 Optical Society of America

1. Introduction

Since the theory of invisible cloak was firstly presented in 2006 [1–2

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

], this realm has got much attention. The basic principle of this theory is to keep Maxwell’s equations form-invariant, when the original coordinate system has altered to a new topology. A distorted space with a hollow region can be created by coordinate transformation. If the constitutive parameters have changed accordingly, the EM wave will propagate around the hollow region smoothly. This phenomenon has been verified experimentally by Schurig et al. [3

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial Electromagnetic Cloak at Microwave Frequencies,” Science 314, 977–980(2006). [CrossRef] [PubMed]

] in 2006 and the cloak they realized is cylindrical-form. Cloaks with more complicated shapes have been studied gradually. Axes symmetrical cloaks, such as polygonal [4–5

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transformations of Maxwell’s Equations,” Photon. Nanostruc. Fundam. Appl. 6, 87–95(2008). [CrossRef]

] and elliptical [6–8

H. Ma, S. B. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825(2008). [CrossRef]

] shapes have been designed by analytical functions. In recent months, various methods on achieving arbitrary shaped cloaks have been proposed. For instance, Yan et al [9

W. Yan, M. Yan, Z. Ruan, and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” New J. Phys. 10, 043040 (2008). [CrossRef]

] gives a method based on analyzing the wave vectors on cloak boundaries. Li et al [10

C. Li and F. Li, “Two-dimensional electromagnetic cloaks with arbitrary geometries,” Opt. Express 16, 13414–13420 (2008). [CrossRef] [PubMed]

] and Ma et al [11–12

H. Ma, S. Qu, Z. Xu, and J. Wang, “Numerical method for designing approximate cloaks with arbitrary shapes,” Phys. Rev. E 78, 036608 (2008). [CrossRef]

] distort polar axes into interior region of the cloak to realize a compressed space. Jiang et al [13

W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, “Analytical design of conformal invisible cloaks for arbitrarily shaped objects,” Phys. Rev. E 77, 066607 (2008). [CrossRef]

] and Nicolet et al. [14

A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Lett. 33, 1584–1586(2008). [CrossRef] [PubMed]

] respectively use different particular functions to describe the irregular boundaries. Hu et al. [15

J. Hu, X. Zhou, and G. Hu, “Design Method for Electromagnetic Cloak with Arbitrary Shapes Based on Laplace’s Equation,” Arxiv.0807.3837.

] solves Laplace’s equation to get the transformed coordinates inside cloak. Although those processes can guide EM wave efficiently, they are not flexible to realize a parameter controllable cloak with arbitrary shapes. Nonlinear transformation is valid to control the variation of coordinates of the cloak with regular shape [16

W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Design of optical cloak with nonlinear transformation,” Opt. Express 16, 5444–5452(2008). [CrossRef] [PubMed]

], but once utilized to the cloak with arbitrary shape, this approach will cause a few inconvenience in programming. The methods given by [9–15

W. Yan, M. Yan, Z. Ruan, and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” New J. Phys. 10, 043040 (2008). [CrossRef]

] are utilized to realize an arbitrary cloak, but due to the fixed boundary conditions, the transformed coordinates can not be altered, which means the material parameters are fixed or uncontrollable. It is inconvenient when those constitutive parameters are not desired or hard to implement practically.

Toward implement a material parameters controllable electromagnetic cloak with arbitrary shaped boundaries, this paper proposes utilizing Helmholtz’s equation to get the relationship between original and transformed coordinates inside the cloak. The constitutive parameters can be easily obtained by the Jacobian transformation matrix. Full wave simulation results indicate that by using this method, the electromagnetic wave will propagate around the obstacle shielded by the cloak with no scattering. The constitutive parameters of the cloak are verified being controlled by the propagation constant of Helmholtz’s equation.

2. Helmholtz’s equation method for parameter controllable cloak

According to the theory of coordinate transformation [1

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

], if the original coordinates (x 1, x 2, x 3) are distorted to the new coordinates (x 1′(x 1, x 2, x 3), x 2′(x 1, x 2, x 3), x 3′(x 1, x 2, x 3)), in order to keep Maxwell’s equations form-invariant, the relationship between original constitutive parameters and transformed constitutive parameters must satisfy the following expression

εrij= A ji A ji εr ij [det (A)] 1
(1)
μrij= A ji A ji μr ij [det (A)] 1
(2)

where εij r and μij r are the permittivity and permeability in original space respectively, εr ij and μr ij are those two in distorted space, Ai i ′ = ∂xi /∂xi ′. is the gradient coefficient of inverse mapping xi ′ → xi , A = [Ai i′ ] is the Jacobian transformation matrix. Rewrite Eqs. (1) and (2) as

εrij= ( A ii) 1 ( A ji) 1 εrij [det (A)]
(3)
μri'j'= ( A ii) 1 ( A ji) 1 μrij [det (A)]
(4)

where (Ai i′ )-1 is the element of inverse matrix (A)-1 . Figure 1 shows a section of irregular cloak. Region I and II respectively denote free space and a section of the cloak with arbitrary shaped boundary curves Γ outer and Γ inner . Region III is shielded by the cloak, which can be seen as the expansion of point P(x 1 p , x 2 p , x 3 p ).

Fig. 1. Geometrical graph of the section of irregular invisible cloak

To ensure electromagnetic wave propagating into cloak with no reflection, the outer boundary condition should be x 1Γ outer = xi ′. In region III, no EM field exists, thus the inner boundary condition is chosen as x 1 p = x 1Γ outer . In region II , we utilize homogeneous Helmholtz’s equation and specified boundary conditions to describe the original coordinates as

{ ( 2 x1 2+ 2 x2 2+ 2 x3 2) xi+ k2 xi =0 xi Γouter= xi, xi Γinner= xi P
(5)

where k is the propagation constant. From Eq. (5) we can obtain original coordinate xi straightly. If boundary curves Γ outer and Γ inner are both regular, xi will be given as an analytical function xi = fi (x 1′, x 2′, x 3′). If those boundaries are irregular, it’s difficult to obtain the analytical solution of Eq. (5), and the original coordinate xi will be solved numerically. If k equal to zero, Eq. (5) will change to a Laplace’s equation. This situation has been discussed by Hu et al [15

J. Hu, X. Zhou, and G. Hu, “Design Method for Electromagnetic Cloak with Arbitrary Shapes Based on Laplace’s Equation,” Arxiv.0807.3837.

]. Insert xi to Eqs. (3) and (4), we can get the permittivity and permeability in region II . For a regular spherical cloak, the outer and inner boundaries are both spherical surface, the distortion is just in the radial direction, and Eq. (5) will be changed to the following form

{ 1 r2 [ r ( r 2 r r)]+ k2r =0 r r = r1= 0,r r= r2= r2
(6)

where r 1, r 2 are the radius of the inner and outer boundary of the cloak respectively, the solution of Eq. (6) is given by

r= A· sin kr kr B· cos kr kr
(7)

with

A= k r22 cosk r1 sin [k ( r2 r1)],B= k r22 sink r1 sin [k ( r2 r1)]
(8)

where k(r 2-r 1) ≠ nπ, n is an arbitrary integer. Equation (7) indicates that the original coordinate r is not only relative to the transformed coordinate r′ but also to the propagation constant k. Besides, the changing of r is nonlinear. Thus although boundary conditions are fixed, r is still alterable with k. That means the constitutive parameters of cloak will be controllable by the propagation constant k. We can obtain gradient coefficient as

A rr= r r= A ( cos kr r sin kr kr2)+B ( sin kr r+ cos kr kr2)
(9)

Let θ = θ′ , φ = (φ′, the other transformation coefficient are

A θθ= A φφ= r r=A sin kr kr2B cos kr kr2
(10)

From Eq. (3), (4), (9) and (10), the constitutive parameters of spherical cloak are

εr= μr= ( A sin kr kr2B cos kr kr2)2 A ( cos kr r sin kr kr2)+B ( sin kr r+ cos kr kr2)
(11)
εθ= εφ= μθ= μφ= A ( cos kr r sin kr kr2)+B ( sin kr r+ cos kr kr2)
(12)

Where A, B are given by Eq. (8). The curve clusters of Eq. (11) and Eq. (12) for different k are shown respectively in Fig. 2, where r 1 = 0.2 m, r 2 = 0.25 m.

Fig. 2. Curves of constitute parameters with variable k inside the cloak

Because the value of μ is equal to ε, we deal with s alone in the next discussion. In Fig. 2(a), εr near r 2 varies obviously with different k. For k < 6π, εr increases steadily and the maximum of which is under 0.41. When k = 8π, the value of εr jumps intensely to 1 in the region between 0.24m and 0.25m. In Fig. 2(b), the varying range of εθ and εφ also rely on k. When k shifts from 2π to 8π, this range changes from [6.3, 3.8] to [8.2, 1]. Thus we can adopt different groups of constitutive parameters to form a cloak, which may give us a chance to choose appropriate parameters to realize a real cloak.

3. Designing a parameter controllable cloak with arbitrary boundaries

For the reason of simplifying, we choose a 2D-cloak shown in Fig. 1 as the model to make full-wave simulation by COMSOL Multi-physics Software, which is based on finite element method (FEM). In the simulation, we fill region III with perfect electric conductor (PEC) and apply scattering boundary conditions to the outside of calculation area. Perfectly matched layers (PML) are used to absorb TM plane wave with z-polarized electric fields generated by E-field boundary. The frequency of plane wave is 7.5 GHz and the calculation area is set to be 0.48m×0.48m. Similar to the processing Hu et al. [15

J. Hu, X. Zhou, and G. Hu, “Design Method for Electromagnetic Cloak with Arbitrary Shapes Based on Laplace’s Equation,” Arxiv.0807.3837.

] adopted, we use PDE model of Helmholtz’ s equation to calculate original coordinates in region II firstly. The constitutive parameters of the cloak are obtained secondly. Then utilize TE wave model to verify the property of cloak we designed, this whole processing will compute 1704399 unknowns for each k. Figure 3 shows the distribution of electric field near the cloak with different k. Figure 3 indicates that in any case of k, TM wave will penetrate the cloak with very low refraction and propagate around the PEC region smoothly. When k increases from 10π to 50π, in the area near inner boundary, the electric field become more homogeneous, while in the region near outer boundary, the field seems been compressed and relies on the shape of cloak.

Fig. 3. Distribution of electric-field near the cloak with different k. (a) k = 10π. (b) k = 30π. (c) k = 50π.

Figure 4 displays the distribution of εzz with different k. As can be seen, the region of εzz < 1 expands from the space near inner boundary to the most area inside the cloak, when k increases from 10π to 50π. In the region near outer boundary, εzz increase rapidly to the value above 4. For k choose 50π, this value even equal to 10. The other constitutive parameters μxx , μyy and μzz can be easily obtained and analyzed directly through the distribution figures similar to Fig. 4, which is not mentioned here. The above discussion has verified that the material parameters of the cloak are controlled by varying k. The electromagnetic cloak designed by Schurig et al [3

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial Electromagnetic Cloak at Microwave Frequencies,” Science 314, 977–980(2006). [CrossRef] [PubMed]

] is constituted by 10 layers of split-ring resonators (SRRs). The relative permeability μr of this cloak increases from 0.003 to 0.279 and the maximum change to the SRRs structure is less than 0.09 mm. As the shift of μr relies on tiny alteration to the structure, the precision in fabricating SRRs should be strict. Utilizing our approach, the constitutive parameters will be chosen in a wide range for a certain k. It can effectively reduce the restriction of fabricating precision, especially in the case that the structure is very small. Thus it is convenient and flexible in accomplishing a real invisible cloak.

Fig. 4. Distribution of εzz inside the cloak with different k. (a) k = 10π. (b) k = 30π. (c) k = 50π.

4. Conclusion

A method using Helmholtz’s equation to fabricate a parameter controllable cloak is demonstrated. The constitutive parameters of a regular spherical cloak have been obtained and verified to be controlled by the propagation constant of Helmholtz’s equation. An arbitrary shaped 2D-cloak has been designed by our method and simulated by FEM. The results coincide with what we expect. Although the cloak simulated in this paper is given as 2D, our method is still usable in designing a 3D irregular shaped cloak. In conclusion, this approach can guide wave efficiently and control the constitutive parameters conveniently, which will make the designing of arbitrary shaped cloak more flexible. It has great significance in engineering and should be given more attention.

Acknowledgment

This work was partially supported by the National Science Foundation of China (Grant No.60871069) and Aeronautical Science Foundation of China (Grant No.20070188002).

References and links

1.

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling Electromagnetic Fields,” Science 312, 1780–1782 (2006). [CrossRef] [PubMed]

2.

U. Leonhardt, “Optical conformal mapping,” Science 312, 1777–1780(2006). [CrossRef] [PubMed]

3.

D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, “Metamaterial Electromagnetic Cloak at Microwave Frequencies,” Science 314, 977–980(2006). [CrossRef] [PubMed]

4.

M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, “Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transformations of Maxwell’s Equations,” Photon. Nanostruc. Fundam. Appl. 6, 87–95(2008). [CrossRef]

5.

Q. Wu, K. Zhang, F. Y. Meng, and L. W. Li, “Arbitrary N-side Regular Polygonal Invisible Cloak,” in Proceedings of the 2008 International Workshop on Metamaterials, (Academic, Nanjing, China, 2008), 60–63.

6.

H. Ma, S. B. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, “Material parameter equation for elliptical cylindrical cloaks,” Phys. Rev. A 77, 013825(2008). [CrossRef]

7.

D. Kwon and D. H. Werner, “Two-dimensional eccentric elliptic electromagnetic cloaks,” Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]

8.

W. X. Jiang, T, J. Cui, G. X. Yu, X. Q. Lin, Q. Cheng, and J. Y. Chin, “Arbitrarily elliptical-cylindrical invisible cloaking,” J. Phys. D: Appl. Phys. 41, 085504 (2008). [CrossRef]

9.

W. Yan, M. Yan, Z. Ruan, and M. Qiu, “Coordinate transformations make perfect invisibility cloaks with arbitrary shape,” New J. Phys. 10, 043040 (2008). [CrossRef]

10.

C. Li and F. Li, “Two-dimensional electromagnetic cloaks with arbitrary geometries,” Opt. Express 16, 13414–13420 (2008). [CrossRef] [PubMed]

11.

H. Ma, S. Qu, Z. Xu, and J. Wang, “Numerical method for designing approximate cloaks with arbitrary shapes,” Phys. Rev. E 78, 036608 (2008). [CrossRef]

12.

H. Ma, S. Qu, Z. Xu, and J. Wang, “Approximation approach of designing practical cloaks with arbitrary shapes,” Opt. Express 16, 15449–15454 (2008). [CrossRef] [PubMed]

13.

W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, “Analytical design of conformal invisible cloaks for arbitrarily shaped objects,” Phys. Rev. E 77, 066607 (2008). [CrossRef]

14.

A. Nicolet, F. Zolla, and S. Guenneau, “Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section,” Opt. Lett. 33, 1584–1586(2008). [CrossRef] [PubMed]

15.

J. Hu, X. Zhou, and G. Hu, “Design Method for Electromagnetic Cloak with Arbitrary Shapes Based on Laplace’s Equation,” Arxiv.0807.3837.

16.

W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Design of optical cloak with nonlinear transformation,” Opt. Express 16, 5444–5452(2008). [CrossRef] [PubMed]

OCIS Codes
(160.1190) Materials : Anisotropic optical materials
(230.0230) Optical devices : Optical devices
(230.3205) Optical devices : Invisibility cloaks

ToC Category:
Physical Optics

History
Original Manuscript: January 12, 2009
Revised Manuscript: February 6, 2009
Manuscript Accepted: February 6, 2009
Published: February 23, 2009

Citation
Xi Chen, Yunqi Fu, and Naichang Yuan, "Invisible cloak design with controlled constitutive parameters and arbitrary shaped boundaries through Helmholtz’s equation," Opt. Express 17, 3581-3586 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-5-3581


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. B. Pendry, D. Schurig, and D. R. Smith, "Controlling Electromagnetic Fields," Science 312, 1780-1782 (2006). [CrossRef] [PubMed]
  2. U. Leonhardt, "Optical conformal mapping," Science 312, 1777-1780(2006). [CrossRef] [PubMed]
  3. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, "Metamaterial Electromagnetic Cloak at Microwave Frequencies," Science 314, 977-980 (2006). [CrossRef] [PubMed]
  4. M. Rahm, D. Schurig, D. A. Roberts, S. A. Cummer, and D. R. Smith, "Design of Electromagnetic Cloaks and Concentrators Using Form-Invariant Coordinate Transformations of Maxwell’s Equations," Photon. Nanostruct. Fundam. Appl. 6, 87-95 (2008). [CrossRef]
  5. Q. Wu, K. Zhang, F. Y. Meng, and L. W. Li, "Arbitrary N-side Regular Polygonal Invisible Cloak," in Proceedings of the 2008 International Workshop on Metamaterials, (Academic, Nanjing, China, 2008), 60-63.
  6. H. Ma, S. B. Qu, Z. Xu, J. Q. Zhang, B. W. Chen, and J. F. Wang, "Material parameter equation for elliptical cylindrical cloaks," Phys. Rev. A 77, 013825 (2008). [CrossRef]
  7. D. Kwon and D. H. Werner, "Two-dimensional eccentric elliptic electromagnetic cloaks," Appl. Phys. Lett. 92, 013505 (2008). [CrossRef]
  8. W. X. Jiang, T, J. Cui, G. X. Yu, X. Q. Lin, Q. Cheng, and J. Y. Chin, "Arbitrarily elliptical-cylindrical invisible cloaking," J. Phys. D: Appl. Phys. 41, 085504 (2008). [CrossRef]
  9. W. Yan, M. Yan, Z. Ruan, and M. Qiu, "Coordinate transformations make perfect invisibility cloaks with arbitrary shape," New J. Phys. 10, 043040 (2008). [CrossRef]
  10. C. Li and F. Li, "Two-dimensional electromagnetic cloaks with arbitrary geometries," Opt. Express 16, 13414-13420 (2008). [CrossRef] [PubMed]
  11. H. Ma, S. Qu, Z. Xu, and J. Wang, "Numerical method for designing approximate cloaks with arbitrary shapes," Phys. Rev. E 78, 036608 (2008). [CrossRef]
  12. H. Ma, S. Qu, Z. Xu, and J. Wang, "Approximation approach of designing practical cloaks with arbitrary shapes," Opt. Express 16, 15449-15454 (2008). [CrossRef] [PubMed]
  13. W. X. Jiang, J. Y. Chin, Z. Li, Q. Cheng, R. Liu, and T. J. Cui, "Analytical design of conformal invisible cloaks for arbitrarily shaped objects," Phys. Rev. E 77, 066607 (2008). [CrossRef]
  14. A. Nicolet, F. Zolla, and S. Guenneau, "Electromagnetic analysis of cylindrical cloaks of an arbitrary cross section," Opt. Lett. 33, 1584-1586 (2008). [CrossRef] [PubMed]
  15. J. Hu, X. Zhou, and G. Hu, "Design Method for Electromagnetic Cloak with Arbitrary Shapes Based on Laplace’s Equation," Arxiv.0807.3837.
  16. W. Cai, U. K. Chettiar, A. V. Kildishev and V. M. Shalaev, "Design of optical cloak with nonlinear transformation," Opt. Express 16, 5444-5452 (2008). [CrossRef] [PubMed]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Figures

Fig. 1. Fig. 2. Fig. 3.
 
Fig. 4.
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited