CN110701956B - 一种基于热偶极子的热隐身方法 - Google Patents
一种基于热偶极子的热隐身方法 Download PDFInfo
- Publication number
- CN110701956B CN110701956B CN201910946020.0A CN201910946020A CN110701956B CN 110701956 B CN110701956 B CN 110701956B CN 201910946020 A CN201910946020 A CN 201910946020A CN 110701956 B CN110701956 B CN 110701956B
- Authority
- CN
- China
- Prior art keywords
- thermal
- thermocouple
- field
- particles
- pole
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41H—ARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
- F41H3/00—Camouflage, i.e. means or methods for concealment or disguise
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
本发明属于热隐身技术领域,具体为一种基于热偶极子的热隐身方法。本发明方法包括:构建热偶极子结构,精确计算热偶极子产生的热场;调节热偶极子的热偶极矩使其产生的热场与颗粒在均匀热场中产生的效果相消,从而实现热隐身。本发明方法适用于二维情形,也适用三维情形;而且在稳态和瞬态条件下都可以实现热隐身效果。本发明方法简单易行制对材料没有严苛要求。本发明也为其他领域实现隐身、伪装等性能提供了新的方法。
Description
技术领域
本发明属于热隐身技术领域,具体涉及一种基于热偶极子的热隐身方法。
背景技术
随着人们对能源问题的日益关注,许多研究人员已将研究重点转向热能管理这一领域。近十年来新兴的热超构材料领域更是推动了热能管理这一领域的发展。最具代表性的例子是热隐身,它几乎贯穿了热超构材料的发展历程。热隐身的特征是基质(背景)的热场均匀。为了实现这一目标,人们提出了许多方案,最初的方法是基于变换热学实现热隐身。然而,变换热学存在几个问题会限制其实际应用。首先是各向异性,这需要不同的径向和切向分量的张量热导率。第二种是不均匀性,即为空间分布的热导率。第三个是奇点,热导率为零或无穷大。第四个是需要非常大的热导率,然而,自然界中材料的热导率的分布范围一般从0.026Wm-1K-1(空气)到430Wm-1K-1(银),其他范围的热导率都是不存在的。这几个问题虽然限制了材料的实际应用,但也促进了超材料的发展。幸运的是,各向异性和非均匀性问题很快得到了解决。然而,奇点和热导率不存在的问题仍然不能同时解决。这两个问题(奇点和不存在)在很大程度上限制了热超材料的发展。
为了完全解决这两个问题,本发明提出了一种基于热偶极子实现热隐身的理论,它可以同时消除奇点和热导率不存在的要求。事实上,本发明甚至不需要设计任何壳层,一个热偶极子就足够了,见图1。因为热偶极子热场的特殊性,通过设计热偶极子矩(M)就可以抵消粒子的影响。接下来,本发明建立了基于热偶极子的二维热学理论并成功的推广到三维。通过有限元模拟和实验室实验验证了该理论的正确性。除此之外,本发明还证明了热偶极子实现热隐身不仅在稳态情况下成立,在瞬态情况下也依然成立。
发明内容
本发明的目的在于提出一种在稳态和瞬态情况下均能呈现热隐身的基于热偶极子的热隐身方法。
本发明提供的基于热偶极子的热隐身方法,首先构建热偶极子结构,精确计算热偶极子产生的热场,通过调节热偶极子的热偶极矩使其产生的热场与颗粒在均匀热场中产生的效果相消,从而实现热隐身。本发明方法是基于热偶极子实现热隐身,是一种的热隐身新机制。
本发明提供的方法,其中,二维情形可以直接推广到三维。
本发明方法,无论是稳态还是瞬态条件下,都可以实现热隐身的效果。
本发明提供的方法,其主要原理为,通过计算颗粒在外部均匀热场中产生的效果和热偶极子产生的热场来确定热偶极子的热偶极矩。其中热传导过程由傅里叶定律描述。下面具体推导本发明方法中,热偶极子实现热隐身满足的条件:
设颗粒的热导率和半径为kp和rp,热偶极子的热偶极矩为M,基质(背景)的热导率为km。根据傅里叶定律,当颗粒嵌在一个基质(背景)中时,基质(背景)的热场Gme可表示为:
其中,M=Q·l,当r>>l时,公式(4)成立。根据矢量场的叠加原理,基质(背景)的热场 (由均匀外场和热偶极子产生)Gs可表示为:
热隐身的特点是基质(背景)的热场应均匀,所以公式(6)中的第二项应为零,即:
根据公式(7),可得热偶极子的偶极矩为:
M=(κm-κp)fG0(8)
热偶极子实现热隐身的二维理论可以直接拓展到三维。在三维情况下,设颗粒的热导率和半径为κ′p和r′p,热偶极子的热偶极矩为M′,基质(背景)的热导率为κ′m。类似于二维情形的推导,当热偶极子的热偶极矩M′满足:
M′=(κ′m-κ′p)f′G′0 (9)
其中,G′0为均匀热场,f′=4πr′p 3/3为颗粒的体积。
本发明的优点:
(1)本发明提出的方法是精确解;
(2)本发明提出的方法简单,利用自然材料即可实现;
(3)本发明提出的方法不仅适用于稳态还适用于瞬态。
附图说明
图1为二维结构示意图和模拟图。其中,(a)为颗粒在热场中的示意图,(d)为相应的模拟图。(b)为热偶极子在颗粒中的示意图,(e)为相应的模拟图。(c)为颗粒在热场和热偶极子激发的场中的示意图。(f)相应的模拟图。其他参数:整体尺寸大小:20×20cm2,颗粒半径rp=6cm,偶极子间距l=2cm,颗粒和基质(背景)的热导率κp与κm分别为200Wm-1K-1与400Wm-1K-1。热偶极矩M=452.4Wm产生的热量为Q=22620W,热偶极子半径rd=0.5cm,初始温度T0=303K。
图2为二维结构热偶极子实现热隐身的实验图。其中,(a)为实验示意图,(b)、(c)分别为颗粒在有无热偶极子的情况下在热场中的实验温度分布图。(d)、(e)分别为对应于(b)、(c) 的模拟图。其他参数:铜的热导率为400Wm-1K-1,空气的热导率为0.026Wm-1K-1,样品中共256个半径为0.22cm的空气洞。
图3展示了二维情况下利用热偶极子器件实现热隐身的效果。其中,(a)展示了当有热偶极子存在时的温度分布,(b)展示了均匀温度分布,(c)展示了(a)图和(b)图的基质(背景)区域的温度差。
图4(a)和(b)展示了热偶极子的半径rd和距离l对热隐身效果的影响。
图5为瞬态情况下二维结构模拟图。其中,(a-d)分别展示了颗粒在无热偶极子的情况下在均匀热场中10s,20s,40s,80s时的温度分布图。(e-h)分别展示了颗粒在有热偶极子的情况下在均匀热场中10s,20s,40s,80s时的温度分布图。
图6展示了三维情况下利用热偶极子器件实现热隐身的效果。其中,(a)展示了当有热偶极子存在时的温度分布,(b)展示了均匀温度分布,(c)展示了(a)图和(b)图的基质(背景)区域的温度差。其他参数:整体大小:20×20×20cm3,颗粒半径r′p=6cm,偶极子间距l′=2cm,颗粒和基质(背景)的热导率k′p与k′m分别为200Wm-1K-1与400Wm-1K-1。热偶极矩M′=36.2 Wm产生的热量为Q′=1810W,热偶极子半径r′d=0.5cm,初始温度T0=303K。
图7为瞬态情况下三维结构模拟图。其中,(a-d)分别展示了颗粒在无热偶极子的情况下在均匀热场中10s,20s,40s,80s时的温度分布。(e-h)分别展示了颗粒在有热偶极子的情况下在均匀热场中10s,20s,40s,80s时的温度分布。
具体实施方式
下面结合具体实例和附图来详细说明本发明,但本发明并不仅限于此。
首先,我们利用有限元模拟来验证我们的推论,利用商业软件COMSOL进行模拟。其中图 1(a)与(d)中左端温度323K,右端温度283K,上下边界绝热。如果一个热导率与基质(背景) 不同的颗粒处在均匀基质(背景)的中心,其示意图见图1(a),其温度分布图见图1(d)。扭曲的外场分布使颗粒可见,然后,我们将一对热偶极子放置在颗粒中间,见图1(b)与(e),这时不再添加外场,所有边界绝热,作为参考温度令θ=±π/2处的温度为303K,其示意图见图 1(b),其温度分布图见图1(e)。最后,将图1(a)图与(b)图合并,见图1(c),其温度分布图见图1(f),因为热偶极子的添加,扭曲的外场分布可以恢复。即利用热偶极子实现了热隐身。
然后我们做了实验室实验。其中,图2(a)为实验示意图,实验样品是由经过激光雕刻的铜板构成,为了使样特定区域的等效热导率为200Wm-1K-1,我们在热导率为400Wm-1K-1的铜板上雕刻出热导率为0.026Wm-1K-1的空气洞。并将样品置于均匀热场中,左端加热源,右端为冷源。由于特定区域的热导率与基质(背景)不同,所以基质(背景)的温度分布不均匀,图2(b)为用红外热像仪Flir E60拍摄的样品的温度分布图。然后在样品上放置一对热偶极子,热偶极子的功率分别为22620W和-22620W,由于热偶极子的加入,基质(背景)的温度分布重新变得均匀,图2(c)为我们用红外热像仪Flir E60拍摄的样品的温度分布图。图2(d)和(e) 分别是严格对应于实验情况(b)和(c)的模拟图。理论结果,实验结果和模拟结果都验证了利用热偶极子实现热隐身的有效性。
在计算过程中,本发明利用到了近似条件:r>>l。因此我们需要验证此近似条件的有效性。首先比较了有无热偶极子存在时的结果,见图3(a)和(b)。图3(c)为图3(a)和(b)基质(背景)的温度差分布,即ΔT=T1-T2。其中温度差最大值为ΔTmax=0.04K,相比于冷热源温差ΔT0=40K,相对误差仅为0.1%。可见相对误差可以评价热偶极子对于热隐身的效果。
接下来本发明继续模拟了瞬态情况下热偶极子实现热隐身,见图5。当初始温度被设置为冷热源的中间温度303K时,热隐身的效果很好。图5(a-d)是没有热偶极子情况下用于作为参考的模拟图,分别对应于10、20、40、80s时的温度分布,基质(背景)温度分布不均匀。图5(e-h)是添加了热偶极子情况下的模拟图,和图5(a-d)一一对应,并且基质(背景)温度均匀,实现了热隐身的效果。
以上均是基于二维情况下的讨论。为了使研究结果更贴近实用,本发明对三维情况下利用热偶极子实现热隐身的情况继续做了模拟,见图6。根据三维理论推导,可以利用公式(9)设置热偶极子的热偶极矩。图6(a)和(b)分别对应有无热偶极子存在时的情况,基质(背景)相同的均匀温度分布验证了热偶极子实现热隐身在三维情况下的有效性。类似于二维讨论,我们绘制了图6(a)和(b)基质(背景)区域的温度差ΔT′=T′1-T′2定量的进行讨论,见图6(c)。温度差最大值为0.07K,相比于冷热源温差40K,相对误差仅为0.18%,显示了热偶极子的有效性。
类似于二维瞬态模拟,我们还进行了三维情况下的瞬态模拟,见图7。当初始温度被设置为冷热源的中间温度303K时,热隐身的效果很好。图7(a-d)是没有热偶极子情况下用于作为参考的模拟图,分别对应于10、20、40、80s时的温度分布,基质(背景)温度分布不均匀。图7(e-h)是添加了热偶极子情况下的模拟图,和图7(a-d)一一对应,并且基质(背景)温度均匀,实现了热隐身的效果。
Claims (2)
1.一种基于热偶极子的热隐身方法,其特征在于,首先构建热偶极子结构,精确计算热偶极子产生的热场;调节热偶极子的热偶极矩使其产生的热场与颗粒在均匀热场中产生的效果相消,从而实现热隐身;具体为:
设颗粒的热导率和半径为κp和rp,热偶极子的热偶极矩为M,基质即背景的热导率为κm;根据傅里叶定律,当颗粒嵌在一个基质中时,基质的热场Gme表示为:
其中,M=Q·l,当r>>l时,公式(4)成立;根据矢量场的叠加原理,由均匀外场和热偶极子产生的基质热场Gs表示为:
热隐身的特点是基质的热场应均匀,公式(6)中的第二项应为零,即:
所述调节热偶极子的热偶极矩使其产生的热场与颗粒在均匀热场中产生的效果相消;根据公式(7),对于二维情形,实现的条件为热偶极子的热偶极矩M满足:
M=(km-kp)fG0
2.根据权利要求1所述基于热偶极子的热隐身方法,其特征在于,对于三维情形,实现的条件为热偶极子的热偶极矩M′满足:
M′=(k′m-k′p)f′G′0
其中,f′=4πr′p 3/3为颗粒的体积;κ′m为基质的热导率,κ′p和r′p为颗粒的热导率和半径,G′0为均匀热场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910946020.0A CN110701956B (zh) | 2019-10-01 | 2019-10-01 | 一种基于热偶极子的热隐身方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910946020.0A CN110701956B (zh) | 2019-10-01 | 2019-10-01 | 一种基于热偶极子的热隐身方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110701956A CN110701956A (zh) | 2020-01-17 |
CN110701956B true CN110701956B (zh) | 2021-09-17 |
Family
ID=69196506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910946020.0A Active CN110701956B (zh) | 2019-10-01 | 2019-10-01 | 一种基于热偶极子的热隐身方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110701956B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5955175A (en) * | 1996-09-20 | 1999-09-21 | W. L. Gore & Associates, Inc. | Infra-red reflective coverings |
DE19704070C2 (de) * | 1997-02-04 | 1998-12-24 | Buck Chem Tech Werke | Tarn- und/oder Täuscheinrichtung |
RU2178136C2 (ru) * | 1999-07-06 | 2002-01-10 | Войсковая Часть 25840 | Устройство маскировки подвижных объектов |
JP5449036B2 (ja) * | 2009-08-05 | 2014-03-19 | 日本アンテナ株式会社 | アンテナおよびアンテナ装置 |
-
2019
- 2019-10-01 CN CN201910946020.0A patent/CN110701956B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110701956A (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu et al. | Encrypted thermal printing with regionalization transformation | |
Papini et al. | Can Hall magnetohydrodynamics explain plasma turbulence at sub-ion scales? | |
Ji et al. | Designing thermal energy harvesting devices with natural materials through optimized microstructures | |
Guenneau et al. | Fick's second law transformed: one path to cloaking in mass diffusion | |
Hu et al. | Local heating realization by reverse thermal cloak | |
Guenneau et al. | Transformation thermodynamics: cloaking and concentrating heat flux | |
CN109684740A (zh) | 一种基于混合网格及时间步长的电磁学多尺度计算方法 | |
Meng et al. | Classical and semirelativistic magnetohydrodynamics with anisotropic ion pressure | |
CN110701956B (zh) | 一种基于热偶极子的热隐身方法 | |
Zhu et al. | Field-coupling topology design of general transformation multiphysics metamaterials with different functions and arbitrary shapes | |
Ji et al. | Deep learning based design of thermal metadevices | |
Zhu et al. | Arbitrary-shape transformation multiphysics cloak by topology optimization | |
Zhang et al. | A bilayer thermal-electric camouflage device suitable for a wide range of natural materials | |
Zohdi | Estimation of electrical heating load-shares for sintering of powder mixtures | |
Parashar et al. | Dissipation of kinetic Alfvénic turbulence as a function of ion and electron temperature ratios | |
Skugoreva et al. | Asymptotic cosmological regimes in scalar–torsion gravity with a perfect fluid | |
Ji et al. | Achieving thermal concentration based on fiber reinforced composite microstructures design | |
Li et al. | Realization and analysis of an Intelligent flux transfer regulator by allocating thermal and DC electric fields | |
Mackay | On the sensitivity of directions that support Voigt wave propagation in infiltrated biaxial dielectric materials | |
Dang et al. | Quasi-periodic array modeling using reduced basis from elemental array | |
Xu et al. | Investigating the thermodynamic performances of to-based metamaterial tunable cells with an entropy generation approach | |
Feng et al. | Design of an omnidirectional camouflage device with anisotropic confocal elliptic geometry in thermal-electric field | |
Xing et al. | Multifunctional thermal rotating cloak with nonconformal geometry | |
Zhang et al. | Three-dimensional thermal illusion devices with arbitrary shape | |
Bai et al. | Magnification device for Laplace equation by using homogeneous and isotropic media with positive values |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |