Topology optimization for negative permeability metamaterials using level-set algorithm

This paper aims to develop a level-set-based topology optimization approach for the design of negative permeability electromagnetic metamaterials, where the topological configuration of the base cell is represented by the zero-level contour of a higher-dimensional level-set function. Such an implicit expression enables us to create a distinct interface between the free space and conducting phase (metal). By seeking for an optimality of a Lagrangian functional in terms of the objective function and the governing wave equation, we derived an adjoint system. The normal velocity (sensitivity) of the level-set model is determined by making the Eulerian derivative of the Lagrangian functional non-positive. Both the governing and adjoint systems are solved by a powerful finite-difference time-domain algorithm. The solution to the adjoint system is separated into two parts, namely the self-adjoint part, which is linearly proportional to the solution of the governing equation; and the non-self-adjoint part, which is obtained by swapping the locations of the incident wave and the receiving planes in the simulation model. From the demonstrative examples, we found that the well-known U-shaped metamaterials might not be the best in terms of the minimal value of the imaginary part of the effective permeability. Following the present topology optimization procedure, some novel structures with desired negative permeability at the specified frequency are obtained.

[1]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[2]  Willie J Padilla,et al.  Composite medium with simultaneously negative permeability and permittivity , 2000, Physical review letters.

[3]  Chun Lu,et al.  Level set simulation of dislocation dynamics in thin films , 2006 .

[4]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[5]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[6]  Dominique Lesselier,et al.  Shape reconstruction of buried obstacles by controlled evolution of a level set: from a min-max formulation to numerical experimentation , 2001 .

[7]  P Y Chen,et al.  Synthesis design of artificial magnetic metamaterials using a genetic algorithm. , 2008, Optics express.

[8]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[9]  J. Pendry,et al.  Magnetism from conductors and enhanced nonlinear phenomena , 1999 .

[10]  Ole Sigmund,et al.  A topology optimization method for design of negative permeability metamaterials , 2010 .

[11]  Stanley Osher,et al.  A survey on level set methods for inverse problems and optimal design , 2005, European Journal of Applied Mathematics.

[12]  Qing Li,et al.  A level-set procedure for the design of electromagnetic metamaterials. , 2010, Optics express.

[13]  Jiangtao Huangfu,et al.  Left-handed materials composed of only S-shaped resonators. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  Lin He,et al.  Reconstruction of shapes and impedance functions using few far-field measurements , 2009, J. Comput. Phys..

[15]  J. Pendry,et al.  Negative refraction makes a perfect lens , 2000, Physical review letters.

[16]  David R. Smith,et al.  Electromagnetic parameter retrieval from inhomogeneous metamaterials. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Willie J. Padilla,et al.  Electrically resonant terahertz metamaterials: Theoretical and experimental investigations , 2007 .

[18]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[19]  Rolf Schuhmann,et al.  Extraction of effective metamaterial parameters by parameter fitting of dispersive models , 2007 .

[20]  Roger F. Harrington,et al.  Field computation by moment methods , 1968 .

[21]  V. Veselago.  The Electrodynamics of Substances with Simultaneously Negative Values of ∊ and μ , 1968 .

[22]  James K. Guest,et al.  Level set topology optimization of fluids in Stokes flow , 2009 .

[23]  K. Svanberg.  The method of moving asymptotes—a new method for structural optimization , 1987 .

[24]  Nader Engheta,et al.  A reciprocal phase shifter using novel pseudochiral or ω medium , 1992 .

[25]  J. Ming,et al.  Finite element method in electromagnetics , 2002 .

[26]  M. Bendsøe.  Optimal shape design as a material distribution problem , 1989 .

[27]  F. Santosa.  A Level-set Approach Inverse Problems Involving Obstacles , 1995 .

[28]  David R. Smith,et al.  Controlling Electromagnetic Fields , 2006, Science.

[29]  M. Zhou,et al.  Generalized shape optimization without homogenization , 1992 .

[30]  Yang Xiang,et al.  A level set method for dislocation dynamics , 2003 .

[31]  H. Giessen,et al.  Three-dimensional metamaterials at optical frequencies , 2008, 2008 Conference on Lasers and Electro-Optics and 2008 Conference on Quantum Electronics and Laser Science.

[32]  Sergei A. Tretyakov,et al.  Resonance Properties of Bi-Helix Media at Microwaves , 1997 .

[33]  Wei Li,et al.  Level-set based topology optimization for electromagnetic dipole antenna design , 2010, J. Comput. Phys..

[34]  O. Dorn,et al.  Level set methods for inverse scattering , 2006 .

[35]  G. Allaire,et al.  Structural optimization using sensitivity analysis and a level-set method , 2004 .

[36]  Xiaoming Wang,et al.  A level set method for structural topology optimization , 2003 .

[37]  Qing Li,et al.  A variational level set method for the topology optimization of steady-state Navier-Stokes flow , 2008, J. Comput. Phys..

[38]  M. Cuer,et al.  Control of singular problem via differentiation of a min-max , 1988 .

[39]  Salvatore Torquato,et al.  A variational level set approach for surface area minimization of triply-periodic surfaces , 2007, J. Comput. Phys..