CN110826271B - 一种同时适用于传导、对流和辐射的热隐身斗篷装置 - Google Patents
一种同时适用于传导、对流和辐射的热隐身斗篷装置 Download PDFInfo
- Publication number
- CN110826271B CN110826271B CN201911026554.8A CN201911026554A CN110826271B CN 110826271 B CN110826271 B CN 110826271B CN 201911026554 A CN201911026554 A CN 201911026554A CN 110826271 B CN110826271 B CN 110826271B
- Authority
- CN
- China
- Prior art keywords
- thermal
- heat
- transformation
- cloak
- stealth
- 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
- 230000005855 radiation Effects 0.000 title claims abstract description 19
- 230000009466 transformation Effects 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 13
- 230000001052 transient effect Effects 0.000 claims abstract description 10
- 239000000463 material Substances 0.000 claims abstract description 9
- 238000009792 diffusion process Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000008033 biological extinction Effects 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 5
- 238000009795 derivation Methods 0.000 abstract description 3
- 230000035699 permeability Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 238000012795 verification Methods 0.000 abstract 1
- 230000004907 flux Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明属于红外隐身技术领域,具体为一种同时适用于传导、对流和辐射的热隐身斗篷装置。本发明的热隐身斗篷装置是通过坐标变换由圆形区域从圆心向外压缩而形成的圆环区域内部,并确定热隐身斗篷的材料参数,使任何物体在圆环内即可隐身,即放置在热隐身斗篷内的任何物体都不会被红外探测所发现。该装置不仅适用于稳态还适用于瞬态。具体实现的方法是:通过坐标变换将一个圆形区域压缩成一个环形区域,并将空间的变换转为材料的变换,从而得到热隐身斗篷的材料参数。热传导和热辐射分别是由Fourier定律和Rosseland扩散近似处理;热对流可以通过直接改变速度或改变多孔介质的渗透率控制。本发明的理论推导得到了有限元模拟的验证。
Description
技术领域
本发明属于红外隐身技术领域,具体涉及热隐身斗篷装置。
背景技术
热能传递有三种方式,即:传导、对流和辐射。然而,现有的理论没有一个能同时处理传导、对流和辐射问题,尽管这三种方式在自然界中普遍存在。这在很大程度上限制了实际应用,因为在大多数情况下,这三种方式总是同时存在。为了解决这一问题,我们提出了变换全热学理论,以此来同时控制传导、对流和辐射。热传导和热辐射分别是由Fourier定律和Rosseland扩散近似处理。热对流可以通过直接改变速度或改变多孔介质的渗透率控制。在理论的基础上,本发明进一步设计了不仅适用于稳态还适用于瞬态的热隐身斗篷器件作为模型应用,并通过瞬态有限元模拟进行了验证。这项工作实现了在理论框架下统一热传导、热对流和热辐射,提出了变换全热学,在热管理方面具有潜在应用。
发明内容
本发明的目的在于提出一种能够同时适用于传导、对流和辐射的热隐身斗篷装置。
本发明提供的能够同时适用于传导、对流和辐射的热隐身斗篷装置,是通过坐标变换由一个圆形区域从圆心压缩成的一个环形区域,并确定热隐身斗篷的材料参数,使任何物体在圆环内即可实现隐身,即:放置在热隐身斗篷内的任何物体都不会被红外探测所发现。
本发明不仅适用于稳态情形,也适用于瞬态情形。
本发明中,所述通过坐标变换将一个圆形区域从圆心向外压缩形成圆环区域的方式如下:
考虑二维情况,从虚拟空间(r,θ)到物理空间(r′,θ′)的坐标变化公式为:
其中,R1和R2分别是热隐身斗篷的内径和外径。公式(1)表示将一个圆形区域从圆心压缩成一个环形区域,参见图1。该坐标变换的Jacobian变换矩阵A可表示为:
本发明中,热隐身斗篷的材料参数确定方式如下:
本发明提供热隐身斗篷装置中,热传导通过Fourier定律处理,热辐射通过Rosseland扩散近似处理,热对流可以通过直接改变速度或改变多孔介质的渗透率控制。
考虑热传导、热对流和热辐射的纯流体,瞬态热输运过程的热力学演化过程可由以下公式主导:
其中,ρf,Cf,κf和vf分别表示流体材料的密度、热容、热导率和速度。为热辐射系数,βf为Rosseland平均消光系数,nf为相对折射率,σ为Stefan-Boltzmann常数(其值等于5.67×10-8Wm-2K-4)。T代表温度,t代表时间,为拉普拉斯算子。公式(3)和公式(4)分别表示热流守恒和质量守恒。
传导热流J1可以用傅里叶定律描述:
对流热流J2可以用公式(4)描述:
J2=ρfCfvfT (6)
辐射热流J3由Rosseland扩散近似描述:
总热流JT为:
JT=J1+J2+J3 (8)
公式(1)和公式(2)满足坐标变换形式不变性,由雅可比矩阵A描述,于是以上参数的变换规则为:
即,变换后的折射率不发生改变nf′=nf,变换后的Rosseland消光系数为变换后的热导率为变换后的密度和热容为变换后的速度为v′f=Avf;detA为Jacobian变换矩阵的行列式,Aτ为Jacobian变换矩阵的转置。
本发明的优点:
(1)本发明提出的方法具有普适性,可以同时调控对流、传导和辐射;
(2)本发明提出的方法利用坐标变换来进行调控,精准灵活;
(3)本发明提出的方法可适用于稳态和瞬态。
附图说明
图1是热隐身斗篷的二维示意图,其中,内外径对应的环形区域即为热隐身斗篷,内径以内区域可放置任意物体。
图2是热隐身斗篷的二维瞬态模拟图。颜色变化代表了温度分布的不同,白线为等温线。(a)-(d)的温度区间为300~360K,速度方向vf沿+x方向。(e)-(h)的温度区间为300~1200K,速度方向vf沿+x方向。(i)-(l)的温度区间为300~1200K,速度方向vf沿-x方向。模拟尺寸大小为10×10cm2,R1=2.4cm,R2=3.6cm。背景参数为ρfCf=106Jm-3K-1,nf=1,vf=10- 5m/s,βf=100m-1,κf=1Wm-1K–1。热隐身斗篷的参数设置是根据公式(7)来设计,其中Jacobian矩阵由公式(9)决定。
图3是背景热通量沿x轴的分布。其中,(a)、(b)、(c)为图2中第一、第二、第三列相同边界条件下的稳态热流分布。
具体实施方式
下面结合具体实施例和附图进一步详细说明本发明,但本发明并不仅限于此。
图1展示了热隐身斗篷的二维示意图,内径与外径之间的环形区域即为热隐身斗篷,内径内的区域可放置任意被隐身的物体。在热传导、热对流和热辐射三种传热方式共同作用下,本发明设计的热隐身斗篷器件可以实现完美的热隐身效果。
在进行理论推导后,可以得到坐标变换后的参数设置及雅可比矩阵。为了验证理论推导的正确性,本发明利用有限元模拟软件COMSOL Multiphysics进行模拟。图2为二维瞬态的模拟结果。在三种情况(a-d),(e-h),(i-l)达到瞬态的过程中,无论初始温度和速度方向如何设置,背景的等温线都没有被扭曲且斗篷内部没有热流进入,达到了热学隐身的效果。边界条件设置为:左右边界分别为高温热源和低温冷源,上下边界绝热,中间白色区域的边界同样设置为绝热。
为了定量分析,还绘制了背景中热传导、热对流和热辐射的总热流通量图,见图3。其中,(a),(b)和(c)分别对应图2(a-d),(e-h),(i-l)所示的三种情况。在三种情况下,我们都可以观察到总通量(呈现水平的直线)是沿x轴的常数,因此总热通量是如公式(1)所述的保守量。同样可以验证我们的理论分析。
Claims (1)
1.一种同时适用于传导、对流和辐射的热隐身斗篷装置,其特征在于,为一个圆环区域内部,是通过坐标变换将一个圆形区域从圆心向外压缩而形成,并确定热隐身斗篷的材料参数,使任何物体在圆环内即可隐身;
所述通过坐标变换将一个圆形区域从圆心向外压缩形成圆环区域的方式如下:
考虑二维情况,对于热隐身斗篷,从虚拟空间(r,θ)到物理空间(r′,θ′)的坐标变化公式为:
其中,R1和R2分别是热隐身斗篷的内径和外径;公式(1)表示将一个圆形区域从圆心压缩成一个环形区域,该坐标变换的Jacobian变换矩阵A可表示为:
热隐身斗篷的材料参数确定方式如下:
考虑热传导、热对流和热辐射的纯流体,瞬态热输运过程的热力学演化过程的公式为:
其中,ρf,Cf,κf和vf分别表示流体材料的密度、热容、热导率和速度;为热辐射系数,βf为Rosseland平均消光系数,nf为相对折射率,σ为Stefan-Boltzmann常数;T代表温度,t代表时间,为拉普拉斯算子;公式(1)和公式(2)分别表示热流守恒和质量守恒;
传导热流J1用傅里叶定律描述:
对流热流J2用公式(4)描述:
J2=ρfCfvfT (6)
辐射热流J3由Rosseland扩散近似描述:
总热流JT为:
JT=J1+J2+J3 (8)
公式(1)和公式(2)满足坐标变换形式不变性,由雅可比矩阵A描述,于是以上参数的变换规则为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911026554.8A CN110826271B (zh) | 2019-10-26 | 2019-10-26 | 一种同时适用于传导、对流和辐射的热隐身斗篷装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911026554.8A CN110826271B (zh) | 2019-10-26 | 2019-10-26 | 一种同时适用于传导、对流和辐射的热隐身斗篷装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110826271A CN110826271A (zh) | 2020-02-21 |
CN110826271B true CN110826271B (zh) | 2023-05-30 |
Family
ID=69550645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911026554.8A Active CN110826271B (zh) | 2019-10-26 | 2019-10-26 | 一种同时适用于传导、对流和辐射的热隐身斗篷装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110826271B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111457236B (zh) * | 2020-03-21 | 2022-05-20 | 复旦大学 | 一种呈现红外热幻像且可见光下隐身的全热学超表面 |
CN111877073B (zh) * | 2020-07-26 | 2021-09-28 | 复旦大学 | 一种用于冻土路基保护热防护结构及方法 |
CN113656993B (zh) * | 2021-07-01 | 2023-10-03 | 复旦大学 | 基于温度依赖变换热电学设计的热电隐身斗篷 |
CN113688495B (zh) * | 2021-07-01 | 2024-04-26 | 复旦大学 | 基于温度依赖变换热电场理论的热电转换器及其设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106959053A (zh) * | 2017-01-31 | 2017-07-18 | 大连理工大学 | 一种基于多层石墨烯环层的可调控二维热隐身斗篷 |
CN108959712A (zh) * | 2018-05-31 | 2018-12-07 | 西安交通大学 | 一种被隐区域有内热源的二维热隐身斗篷设计方法 |
CN109033495A (zh) * | 2018-05-31 | 2018-12-18 | 西安交通大学 | 一种基于主动调控热源的二维对流热隐身斗篷设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10895421B2 (en) * | 2016-12-14 | 2021-01-19 | University Of South Florida | Thermoelectric flow cloaking via metamaterials |
-
2019
- 2019-10-26 CN CN201911026554.8A patent/CN110826271B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106959053A (zh) * | 2017-01-31 | 2017-07-18 | 大连理工大学 | 一种基于多层石墨烯环层的可调控二维热隐身斗篷 |
CN108959712A (zh) * | 2018-05-31 | 2018-12-07 | 西安交通大学 | 一种被隐区域有内热源的二维热隐身斗篷设计方法 |
CN109033495A (zh) * | 2018-05-31 | 2018-12-18 | 西安交通大学 | 一种基于主动调控热源的二维对流热隐身斗篷设计方法 |
Non-Patent Citations (1)
Title |
---|
变换热学:热超构材料及其应用;沈翔瀛;黄吉平;《物理学报》;第65卷(第17期);106-132 * |
Also Published As
Publication number | Publication date |
---|---|
CN110826271A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826271B (zh) | 一种同时适用于传导、对流和辐射的热隐身斗篷装置 | |
CN110826272B (zh) | 一种同时适用于传导、对流和辐射的热旋转装置 | |
CN110826265B (zh) | 基于变换热辐射和热传导理论设计的热隐身斗篷 | |
CN113656993B (zh) | 基于温度依赖变换热电学设计的热电隐身斗篷 | |
CN109492341A (zh) | 表面等离激元波导的光热效应仿真方法 | |
CN110600087A (zh) | 呈现热学变色龙现象的各向同性双壳层结构及其实现方法 | |
CN110807280B (zh) | 一种同时适用于传导、对流和辐射的热聚集装置 | |
Zhou et al. | Study on PID parameters tuning based on particle swarm optimization | |
Ghaddar et al. | Steady thermal comfort by radiant heat transfer: The impact of the heater position | |
CN111062160B (zh) | 能够同时调控热传导、热对流和热辐射的热透明装置 | |
CN111177957B (zh) | 能够同时调控热传导、热对流和热辐射的热隐身斗篷 | |
CN113656990B (zh) | 基于温度依赖坐标变换理论的热电旋转器及其设计方法 | |
CN110826266B (zh) | 基于变换热辐射和热传导理论设计的热旋转器装置 | |
Wen et al. | Development and validation of online parameter estimation for HVAC systems | |
Zheng et al. | Numerical Simulation and Experimental Study of Comfort Air Conditioning Influenced by Bottom-supply and Stratum Ventilation Modes | |
Xu et al. | A Fast Computational Fluid Dynamics Model for the Flow and Heat Transfer Characteristics Analysis of Indoor Substation Rooms | |
CN113656992B (zh) | 基于温度依赖变换热电学设计的热电聚集器 | |
Aleshin et al. | Verification of CFD modeling for rotor type of recuperation energy | |
Dobrego et al. | Simulation of cooling tower and influence of aerodynamic elements on its work under conditions of wind | |
Zhao et al. | Influence of covering layer on surface temperature of floor radiant heating system | |
Zhang | Technical-Environmental-Economical Evaluation of the Implementation of a Highly Efficient District Heating System in China | |
Back et al. | Combining CFD and GIS software capabilities to enhance rapid fine-scale urban micro-and bioclimatic modelling | |
Xia et al. | Numerical Simulation on Heat Pipe Heat Exchanger: Effects of Different Ambient Temperature | |
Khrustalyov et al. | Modeling of convective flows in pneumobased objects. Part 1 | |
CN111460604A (zh) | 蓄热类电采暖设备的建模处理方法及装置 |
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 |