2018传热传质学年会集

:Pendry[1]教授料[9,11,13,19]或有效介质理论[5,2021]简化斗篷的参数分布，精心设计后实际制作以上提及的实现隐身斗篷所的是由于坐标变换法的本质导致的。本文 ,流隐身斗篷。3Kb。图1(a)二维热斗篷模型 图1(b)二维斗篷的截面示意图图1热流隐身斗篷示意图(K(x,y)T)(K(x,y)T)0 T(x,y)T(x,y)

在边界1 在边界4 -KTb

在边界2,3 Kb，可以求解得到原始背景温度场。J[K(T)]J[K(x,y)]T(x,y)Y(x,y)2 i1其中Ti与Yi分别为rR3边界上的计算温度值和目标温度值。将迭代过的导果。脚标“i”表示rR3边界上的计算网格点。“m”表示rR3边界上的计算网格点rR3边界上的温度分布能够恢复到原始背景温度，则rR3整个背景区域温度均可恢复。因此，为了减少计算量，仅需将rR3边界上的J[Kc(x,y)]越小，说明背景温度恢复的越好，则热斗篷功能越趋于完善。 Kn1(x,y)Kn(x,y)nPn(x,y)其中n0,1,2,... 其中，n代表迭代步数，βnnn+1步的迭代步长，Pn(xy)为下降方向（搜索方向。Pn(x,y)定义为：KPn(x,y)Jn(x,y)nPn1(x,y)with00 KKKJn(x,y)2x i n 。 Jn1(x,y)2xi

i KK(K(x,y)T)(K(x,y)T)(K(x,y)T)

(K(x,y)T)

，

T(x,y)

在边界1,4 -

T

在边界2,3 其中KxyPxyn+1J[Kc(x,y)]，可以通过改变函数（2）

m cciJ[Kn1(x,y)]T(KnnPn)cci

i1 cc m cc iJ[Kn1(x,y)]T(Kn)nTn(Pncc i

i1 TYT(P mn 。 mi[Tn(Pnii(K(x,y))

(K(x,y))

2(TY)(xx)(yy)

i

(x,y)0在边界1,4 -

b

J(x,y)(x,y)T(x,y)(x,y)T(x,y) m (TY)2在r=R 3iKc(x,y)。T(x,y)。Kλ(x,y)。第五步：通过方程（11）Jn(xy。K第九步：通过方程（9）βn。cc（1210cm×10cmR11cm，R21.3cm，R32cm，T1273K，T2300K。2（a）2（b）2区，计2区的存在，使外界温度场发生了，在背景中的观察者可以探测到2区的存在。这2（c2（a2（d10.1%。 图2（c）加入3区后的温度分布 图2设计结果rR3边界上的温度分布。用圆形符号表示的曲线代表放置斗篷后的背景温度在rR3边界上的温度分布。根据上文分析，当rR3边界上放置斗篷后的背景4（a（b0.25W/(m·K)。离散后的斗篷和两种材料近似后的斗篷对于温度场的影4（c0.62%,6°。

Pendry,J.B,Schurig,D,andSmith,D.R,2006,ControllingElectromagneticFields,J.Science,312,PP.Leonhardt,U,2006,OpticalConformalMap,J.Science,312(23),PP.1777-medium,J.Science,345(6195),PP.427-429.XuC,StiubianuG,andGorodetskyA,2018,Adaptiveinfrared-reflectingsystemsinspiredbycephalopods,Science,359(6383),PP.1495-1500.Zhang,S,Xia,C,andFang,N,2011,Broadbandacousticcloakforultrasoundwaves,J.Phys.Rev.Popa,B,Zigoneanu,L,andCummer,S.A,2011,Experimentalacousticgroundcloakinair,J.Phys.Rev.Lett.,106(25),PP.253901.Farhat,M,Guenneau,S,andEnoch,S,2009,Ultrabroadbandelasticcloakinginthintes,J.Phys.Rev.Lett.,103(2),PP.024301.Brun,M,Guenneau,S,andMovchan,A.B,2009,Achievingcontrolofin-neelasticwaves,J.Appl.Phys.Lett.,94(6),PP.061903.heatflux,J.Opt.Express,20(7),PP.8207.Fan,C.Z,Gao,Y,andHuang,J.P,2008,Shapedgradedmaterialswithanapparentnegativethermalconductivity,J.Appl.Phys.Lett.,92(25),PP.251907.heatlocalization,J.Sci.Rep.,3,PP.1593.Chen,T,Weng,C,andChen,J,2008,Cloakforcurvilinearlyanisotropicmediainconduction,J.Appl.Phys.Lett.,93(11),PP.114103.Guenneau,S,andPuvirajesinghe,T.M,2013,Fick'ssecondlawtransformed:onepathtocloakinginmassdiffusion,J.J.R.Soc.Interface,10(83),PP. Guenneau,S,Petiteau,D,Zerrad,M,Amra,C,andPuvirajesinghe,T,2015,Transformedfourierandfickequationsforthecontrolofheatandmassdifussion,J.AIPAdvances,5(5),PP.053404.Li,J.Y,Gao,Y,andHuang,J.P,2010,Abifunctionalcloakusingtransformationmedia,J.J.Appl.Han,T,Bai,X,Gao,D,Thong,J.T.L,Li,B,andQiu,C,2014,Experimentaldemonstrationofabilayerthermalcloak,J.Phys.Rev.Lett.,112(5),PP.054302.ShenX,LiY,JiangC,NiY,andHuangJ,2016,Thermalcloak-concentrator,J.Appl.Phys.Lett.,impedancetomography,J.InverseProblems,24(1),PP.015016.Narayana,S,andSato,Y,2012,Heatfluxmanipulationwithengineeredthermalmaterials,J.Phys.Rev.Lett.,108(21),PP.214303.Schittny,R,Kadic,M,Guenneau,S,andWegener,M,2013,Experimentsontransformationthermodynamics:moldingtheflowofheat,J.Phys.Rev.Lett.,110(19),PP.195901.Schurig,D,Mock,J.J,Justice,B.J,Cummer,S.A,Pendry,J.B,Starr,A.F,andSmith,D.R,2006,MetamaterialElectromagneticCloakatMicrowaveFrequencies,J.Science,314,PP.977-980.Ni,X,Wong,Z,Mrejen,M,Wang,Y,andZhang,X,2015,Anultrathininvisibilityskincloakforvisiblelight,J.Science,349(6254),PP.1310.Li,J,andPendry,J.B,2008,Hidingunderthecarpet:anewstrategyforcloaking,J.Phys.Rev.Liu,R,Ji,C,Mock,J.J,Chin,J.Y,Cui,T.J,andSmith,D.R,2009,BroadbandGround- neCloak,J.Science,323,PP.366-369.Huang,C,andWu,J,1994,Two‐dimensionalinverseprobleminestimatingheatfluxesofanenclosurewithunknowninternalheatsources,J.J.Appl.Phys,76(1),PP.133.Cheng,C,andChang,M,2005,PredictionsofinternaltemperaturedistributionofPEMFCbynon-destructiveinversemethod,J.J.PowerSources,139,PP.115-125.Huang,C,andYan,J,1995,Aninverseprobleminsimultaneouslymeasuringtemperature-dependentthermalconductivityandheatcapacity,J.Int.J.HeatMassTransfer,38(18),PP.3433-3441.Alifanov,O.M,andArtyukhin,E.A,1975,Regularizednumericalsolutionofnonlinearinverseheat-conductionproblem,J.J.Eng.Thermophys,29(1),PP.934-938.Chen,W,andYang,Y,2010,Inverseproblemofestimatingtheheatfluxattheroller/workpieceinterfaceduringarollingprocess,J.Appl.ThermalEng,30(10),PP.1247-1254.Borukhov,V.T,andTimoshpol'skii,V.I,2005,FunctionalIdentificationoftheNonlinearThermal-ConductivityCoefficientbyGradientMethods.I.ConjugateOp