劉冶，李竹影，俞翔，曹文康，趙亞（海軍工程大學(xué)，武漢 430033）

形狀任意的超材料電磁隱身波長(zhǎng)變換器的設(shè)計(jì)

劉冶，李竹影，俞翔，曹文康，趙亞
（海軍工程大學(xué)，武漢 430033）

目的 研究形狀任意的超材料電磁隱身波長(zhǎng)變換器的設(shè)計(jì)方法。方法 基于坐標(biāo)變換理論，通過(guò)對(duì)空間中的不同區(qū)域先后進(jìn)行壓縮和擴(kuò)張變換，設(shè)計(jì)出該裝置不同區(qū)域的電磁參數(shù)張量。結(jié)果 一方面，入射波經(jīng)過(guò)該裝置時(shí)，在其界面不產(chǎn)生散射，并且入射波穿過(guò)裝置進(jìn)入外部自由空間時(shí)又能恢復(fù)其初始傳播狀態(tài)；另一方面，該裝置在特定區(qū)域改變了電磁波的波長(zhǎng)。結(jié)論 有限元軟件的仿真結(jié)果顯示波長(zhǎng)變換器實(shí)現(xiàn)了預(yù)期的電磁功能，從而驗(yàn)證了該裝置電磁參數(shù)張量表達(dá)式的正確性。

超材料；坐標(biāo)變換；波長(zhǎng)變換器；隱身

人工超材料實(shí)現(xiàn)了自然材料并不具備的奇異物理性能，除了負(fù)折射、逆多普勒效應(yīng)、完美透鏡、逆Cerenkov輻射效應(yīng)等反常物理特性[1]外，而以坐標(biāo)變換理論[2—3]為基礎(chǔ)的超材料[4]在理論上可實(shí)現(xiàn)對(duì)電磁波的完美隱身。自Schurig等人在2006年完成了世界上第一個(gè)簡(jiǎn)化參數(shù)的圓柱形電磁隱身斗篷的參數(shù)設(shè)計(jì)和微波仿真驗(yàn)證[4]后，基于坐標(biāo)變換理論來(lái)設(shè)計(jì)電磁隱身斗篷及其拓展裝置的研究很快成為超材料領(lǐng)域的又一個(gè)重要研究方向。除了隱身斗篷[5—11]，人們應(yīng)用坐標(biāo)變換理論還設(shè)計(jì)了電磁聚焦器[12]、電磁彎折器[13]、規(guī)則形狀的波長(zhǎng)變換器[14]、波形轉(zhuǎn)換器[15]等多種功能型超材料電磁裝置。

在超材料裝置的研究中，某課題組先后推導(dǎo)出了圓柱形和橢圓柱形波長(zhǎng)變換器的電磁參數(shù)張量，發(fā)展了新型超材料波長(zhǎng)變換器的設(shè)計(jì)方法[14]。在此基礎(chǔ)上，文中將波長(zhǎng)變換器的構(gòu)型拓展到任意形狀，推導(dǎo)了任意形狀波長(zhǎng)變換器的電磁參數(shù)張量，使波長(zhǎng)變換器在構(gòu)型設(shè)計(jì)上的更靈活。

1　電磁參數(shù)張量的推導(dǎo)

文中僅討論如何用波長(zhǎng)變換器來(lái)放大電磁波的波長(zhǎng)，對(duì)波長(zhǎng)縮小的討論同理。在設(shè)計(jì)波長(zhǎng)變換器時(shí)需要進(jìn)行兩次空間變換，第一次變換是將原空間中的某一部分區(qū)域進(jìn)行壓縮，即將圖1a中的區(qū)域II，III壓縮成如圖1b所示的區(qū)域III，以實(shí)現(xiàn)入射波的無(wú)散射傳播；第二次變換是將原空間中剩下的區(qū)域I擴(kuò)張成如圖1c所示的區(qū)域II′，以使得裝置可以放大入射電磁波的波長(zhǎng)。

1.1空間壓縮

圖1為形狀任意的波長(zhǎng)變換器的二維構(gòu)型，該構(gòu)型由三條任意的共形曲線組成。設(shè)曲線的歸一化邊界方程為R（0θ），按照從內(nèi)到外的順序，圖1中三條封閉曲線的方程依次為aR（0θ），bR（0θ），cR（0θ），a，b，c為比例系數(shù)。則有區(qū)域I為0≤r（θ）≤aR0（θ），區(qū)域II為aR（0θ）＜r（θ）≤bR（0θ），區(qū)域III 為bR（0θ）＜r（θ）≤cR（0θ），區(qū)域II′為0＜r（θ）≤bR（0θ）。圖1b表示的壓縮過(guò)程，對(duì)應(yīng)的坐標(biāo)變換式為：

圖1　形狀任意的波長(zhǎng)變換器的空間變換過(guò)程Fig.1　Procedures of space transformations in designing of the invisible wavelength transformer with arbitrary shape

如果令c=1，cR（0θ）=R（0θ）是該波長(zhǎng)變換器的外邊界曲線方程，此時(shí)a，b轉(zhuǎn)化為R（0θ）所圍區(qū)域的壓縮系數(shù)，可令a=n，b=m，并設(shè)：

可以得到：

在直角坐標(biāo)系中：

由式（2）—（4）以及三角函數(shù)的定義可以得到：

由式（5），根據(jù)坐標(biāo)變換理論[2，6]，可以推導(dǎo)出波長(zhǎng)變換器中空間壓縮區(qū)域的電磁參數(shù)張量分別為：

1.2空間擴(kuò)張

在壓縮過(guò)程后，為滿足變換空間的連續(xù)性，需要將原空間中的區(qū)域I擴(kuò)張至區(qū)域II′，如圖1c所示。該擴(kuò)張過(guò)程對(duì)應(yīng)的坐標(biāo)變換式在直角坐標(biāo)系下可寫為：

根據(jù)坐標(biāo)變換理論，可以很快得到波長(zhǎng)變換器在該區(qū)域的電磁參數(shù)張量為：

2　建模與仿真

任意形狀波長(zhǎng)變換器的仿真模型如圖2所示。

圖2　任意形狀波長(zhǎng)變換器的仿真模型Fig.2　The sketch map of simulation model of the invisible wavelength transformer with arbitrary shape

在圖2中，取壓縮系數(shù)m=0.6，n=0.3，并設(shè)波長(zhǎng)變換器的歸一化曲線方程為R0（θ）=0.04·（cos 2θ+ sin 3θ+cosθ+3）。采用基于有限元的COMSOL Multiphysics商業(yè)仿真軟件進(jìn)行數(shù)值仿真計(jì)算。仿真條件中，波源為單位強(qiáng)度的TM平面波，工作頻率為10 GHz，矩形求解域?yàn)?.45 m×0.36 m，求解域的四周設(shè)置為吸收邊界，波長(zhǎng)變換器的各個(gè)邊界均為連續(xù)邊界，TM波從矩形域的左邊界垂直入射。最終的仿真結(jié)果如圖3所示。

圖3　任意形狀波長(zhǎng)變換器的磁場(chǎng)分布Fig.3　The magnetic field distribution of the invisible wavelength-transformer with arbitrary shape

從圖3中可知，當(dāng)TM波傳播至波長(zhǎng)變換器的壓縮區(qū)域III時(shí)，其等相位面開(kāi)始平滑地沿著裝置的擴(kuò)張區(qū)域II′的邊界彎曲，而當(dāng)TM波傳播至擴(kuò)張區(qū)域II′的內(nèi)部時(shí)，TM波的波長(zhǎng)明顯被放大，并且該入射波在擴(kuò)張區(qū)域傳播時(shí)并沒(méi)有出現(xiàn)散射和反射的現(xiàn)象。當(dāng)TM波從波長(zhǎng)變換器內(nèi)部透射至外部自由空間時(shí)，其等相位面又恢復(fù)到原來(lái)的入射波形態(tài)。如果從外部自由空間來(lái)觀測(cè)TM波傳播的情況，則會(huì)發(fā)現(xiàn)TM波好像一直處于初始狀態(tài)，說(shuō)明波長(zhǎng)變換器的存在不影響TM波在該裝置外部的傳播方式。從圖3中可以看到，形狀任意的波長(zhǎng)變換器同時(shí)具備對(duì)電磁波隱身和變換電磁波波長(zhǎng)的功能，該結(jié)果也驗(yàn)證了式（6）和（8）的正確性。

文中從數(shù)值仿真計(jì)算的結(jié)果中還提取了波長(zhǎng)變換器外部自由空間區(qū)域和內(nèi)部擴(kuò)張空間區(qū)域的場(chǎng)強(qiáng)變化圖，并比較了這兩區(qū)域中的電磁波波長(zhǎng)，如圖4所示。圖4a為自由空間區(qū)域，TM波的波長(zhǎng)為0.03 m，傳播過(guò)程中一直保持不變；圖4b為波長(zhǎng)變換器的擴(kuò)張空間區(qū)域，TM波的波長(zhǎng)為0.06 m，恰為自由空間中波長(zhǎng)的2倍，這個(gè)放大倍數(shù)正好就是圖1c中擴(kuò)張區(qū)域II′的擴(kuò)張倍數(shù)，也是波長(zhǎng)變換器在建模過(guò)程中比例系數(shù)b相對(duì)于a的倍數(shù)。因此，圖4的仿真計(jì)算結(jié)果再一次驗(yàn)證了任意形狀波長(zhǎng)變換器的理想隱身和變換入射波波長(zhǎng)的電磁功能。

圖4　求解域不同區(qū)域的場(chǎng)強(qiáng)變化Fig.4　Diagram of field density in different regions of solution domain

3　結(jié)語(yǔ)

應(yīng)用坐標(biāo)變換理論，設(shè)計(jì)出了同時(shí)具備對(duì)電磁波隱身和變換入射電磁波波長(zhǎng)等功能的形狀任意的超材料電磁隱身波長(zhǎng)變換器，推導(dǎo)出了波長(zhǎng)變換器不同區(qū)域的電磁參數(shù)張量表達(dá)式，并采用有限元仿真軟件對(duì)該裝置的電磁特性進(jìn)行數(shù)值仿真，仿真結(jié)果顯示所設(shè)計(jì)的波長(zhǎng)變換器完全實(shí)現(xiàn)了上述電磁功能，從而驗(yàn)證了該裝置電磁參數(shù)張量推導(dǎo)的正確性。關(guān)于該裝置的設(shè)計(jì)方法不僅增加了超材料電磁隱身波長(zhǎng)變換器在構(gòu)型設(shè)計(jì)上的靈活性，還拓展了坐標(biāo)變換理論在設(shè)計(jì)超材料隱身電磁裝置上的應(yīng)用。

[1]SHELBY R A，SMITH D R，SCHULTZ S.Experimental Verification of a Negative Index of Refraction[J].Science，2001，292（5514）：77—79.

[2]SCHURIG D，PENDRY J B，SMITH D R.Calculation of Material Properties and Ray Tracing in Transformation Media[J]. Optics Express，2006，14（21）：9794—9804.

[3]PENDRY J B，SCHURIG D，SMITH D R.Controlling Electromagnetic Fields[J].Science，2006，312（5781）：1780—1782.

[4]SCHURIG D，MOCK J J，JUSTICE B J，et al.Metamaterial Electromagnetic Cloak at Microwave Frequencies[J].Science，2006，314（5801）：977—980.

[5]WU Q，ZHANG K，MENG F Y，et al.Material Parameters Characterization for Arbitrary N-sided Regular Polygonal Invisible Cloak[J].J Phys D：Appl Phys，2009，42：035408.

[6]LI C，LI F.Two-dimensional Electromagnetic Cloaks with Arbitrary Geometries[J].Optics Express，2008，16（17）：13414—13420.

[7]MA H，QU S B，XU Z，et al.Material Parameter Equation for Elliptical Cylindrical Cloaks[J].Phys Rev A，2008，77（1）：013825.

[8]劉冶，李竹影，張旺洲，等.組合型電磁隱身斗篷的超材料設(shè)計(jì)與仿真[J].功能材料，2013，44（15）：2235—2238.

LIU Ye，LI Zhu-ying，ZHANG Wang-zhou，et al.Design and Emulation of Combined-shaped Electromagnetic Stealthy Cloak Made of Metamaterials[J].Journal of Functional Materials，2013，44（15）：2235—2238.

[9]劉冶，李竹影，趙林，等.一種無(wú)奇異參數(shù)橢圓形電磁隱身斗篷的設(shè)計(jì)基礎(chǔ)理論[J].裝備環(huán)境工程，2015，12（1）：6—10.

LIU Y，LI Zhu-ying，ZHAO Lin，et al.A Design Fundamental Theory of Elliptic Cylindrical Invisible Cloak with Non-singular Electromagnetic Tensors[J].Equipment Environmental Engineering，2015，12（1）：6—10.

[10]JIANG W X，CUI T J，YU G X，et al.Arbitrarily Elliptical-cylindrical Invisible Cloaking[J].J Phys D：Appl Phys，2008，41 （8）：085504.

[11]CHEN H S，ZHENG B，SHEN L，et al.Ray-optics Cloaking Devices for Large Object in Incoherent Natural Light[J].Nature Communications，2013（4）：2652.

[12]YU G X，JIANG W X，ZHOU X Y，et al.Non-rotationally Invariant Invisibility Cloaks and Concentrators of Em Wavesz [J].Eur Phys J Appl Phys，2008，44：181—185.

[13]JIANG W X，CUI T J，ZHOU X Y，et al.Arbitrary Bending of Electromagnetic Waves Using Realizable Inhomogeneous and Anisotropic Materials[J].Phys Rev E，2008，78（6）：066607.

[14]劉冶，李竹影，張旺洲，等.超材料電磁隱身波長(zhǎng)變換器的設(shè)計(jì)與仿真[J].功能材料，2014，45（11）：11027—11031.

LIU Ye，LI Zhu-ying，ZHANG Wang-zhou，et al.Design and Simulation of Metamaterial Electromagnetic Invisible Wavelength-transformer[J].Journal of Functional Materials，2014，45（11）：11027—11031.

[15]MA H，QU S B，XU Z，et al.General Method for Designing Wave Shape Transformers[J].Opt Express，2008，16（26）：22072—22082.

Design of Metamaterial Electromagnetic Invisible Wavelength Transformer with Arbitrary Shape

LIU Ye，LI Zhu-ying，YU Xiang，CAO Wen-kang，ZHAO Ya
（Naval University of Engineering，Wuhan 430033，China）

Objective To study the design method of metamaterial electromagnetic invisible wavelength transformer with arbitrary shape.Methods Based on coordinate transformation theory,by means of compressing certain region and then expanding another one in the space,the electromagnetic tensors of the device in different regions were deduced.Results On one hand,when the incoming waves propagated into the device,there was no scattering on the boundary of the device.And the waves returned to their original propagation when they propagated into the free space from the inner region of the device.On the other hand,the device could change the wavelength of electromagnetic wave in the particular region.Conclusion The simulation results of the software based on finite element method showed that the wavelength transformer realized the anticipated electromagnetic function,which validated the electromagnetic tensor expressions of the device.

metamaterial；coordinate transformation；wavelength-transformer；invisibility

2015-07-07；Revised：2015-08-09

10.7643/issn.1672-9242.2016.01.018

TJ99；TM25

A

1672－9242（2016）01－0098-04

2015-07-07；

2015-08-09

國(guó)防預(yù)研基金項(xiàng)目（9140A270304）

Fund：Supported by the National Defense Pre-research Foundation（9140A270304）

劉冶（1986—），男，湖北武漢人，博士，工程師，主要研究方向?yàn)槌牧想[身裝置的構(gòu)型設(shè)計(jì)與性能優(yōu)化。

Biography：LIU Ye（1986—），Male，from Wuhan，Hubei，Ph.D.，Engineer，Research focus：configuration design and property optimization of metamaterial invisible device.