CN110826271B - 一种同时适用于传导、对流和辐射的热隐身斗篷装置 - Google Patents

一种同时适用于传导、对流和辐射的热隐身斗篷装置 Download PDF

Info

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
Application number
CN201911026554.8A
Other languages
English (en)
Other versions
CN110826271A (zh
Inventor
黄吉平
须留钧
杨帅
戴高乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fudan University
Original Assignee
Fudan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fudan University filed Critical Fudan University
Priority to CN201911026554.8A priority Critical patent/CN110826271B/zh
Publication of CN110826271A publication Critical patent/CN110826271A/zh
Application granted granted Critical
Publication of CN110826271B publication Critical patent/CN110826271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational 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′,θ′)的坐标变化公式为:
Figure BDA0002248794290000011
其中,R1和R2分别是热隐身斗篷的内径和外径。公式(1)表示将一个圆形区域从圆心压缩成一个环形区域,参见图1。该坐标变换的Jacobian变换矩阵A可表示为:
Figure BDA0002248794290000012
本发明中,热隐身斗篷的材料参数确定方式如下:
本发明提供热隐身斗篷装置中,热传导通过Fourier定律处理,热辐射通过Rosseland扩散近似处理,热对流可以通过直接改变速度或改变多孔介质的渗透率控制。
考虑热传导、热对流和热辐射的纯流体,瞬态热输运过程的热力学演化过程可由以下公式主导:
Figure BDA0002248794290000021
Figure BDA0002248794290000022
其中,ρf,Cff和vf分别表示流体材料的密度、热容、热导率和速度。
Figure BDA0002248794290000023
为热辐射系数,βf为Rosseland平均消光系数,nf为相对折射率,σ为Stefan-Boltzmann常数(其值等于5.67×10-8Wm-2K-4)。T代表温度,t代表时间,
Figure BDA00022487942900000210
为拉普拉斯算子。公式(3)和公式(4)分别表示热流守恒和质量守恒。
传导热流J1可以用傅里叶定律描述:
Figure BDA0002248794290000024
对流热流J2可以用公式(4)描述:
J2=ρfCfvfT (6)
辐射热流J3由Rosseland扩散近似描述:
Figure BDA0002248794290000025
总热流JT为:
JT=J1+J2+J3 (8)
公式(1)和公式(2)满足坐标变换形式不变性,由雅可比矩阵A描述,于是以上参数的变换规则为:
Figure BDA0002248794290000026
即,变换后的折射率不发生改变nf′=nf,变换后的Rosseland消光系数为
Figure BDA0002248794290000027
变换后的热导率为
Figure BDA0002248794290000028
变换后的密度和热容为
Figure BDA0002248794290000029
变换后的速度为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-1f=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′,θ′)的坐标变化公式为:
Figure FDA0004151929400000011
其中,R1和R2分别是热隐身斗篷的内径和外径;公式(1)表示将一个圆形区域从圆心压缩成一个环形区域,该坐标变换的Jacobian变换矩阵A可表示为:
Figure FDA0004151929400000012
热隐身斗篷的材料参数确定方式如下:
考虑热传导、热对流和热辐射的纯流体,瞬态热输运过程的热力学演化过程的公式为:
Figure FDA0004151929400000013
Figure FDA0004151929400000014
其中,ρf,Cff和vf分别表示流体材料的密度、热容、热导率和速度;
Figure FDA0004151929400000015
为热辐射系数,βf为Rosseland平均消光系数,nf为相对折射率,σ为Stefan-Boltzmann常数;T代表温度,t代表时间,
Figure FDA0004151929400000016
为拉普拉斯算子;公式(1)和公式(2)分别表示热流守恒和质量守恒;
传导热流J1用傅里叶定律描述:
Figure FDA0004151929400000017
对流热流J2用公式(4)描述:
J2=ρfCfvfT (6)
辐射热流J3由Rosseland扩散近似描述:
Figure FDA0004151929400000018
总热流JT为:
JT=J1+J2+J3 (8)
公式(1)和公式(2)满足坐标变换形式不变性,由雅可比矩阵A描述,于是以上参数的变换规则为:
Figure FDA0004151929400000021
即,变换后的折射率不发生改变nf′=nf,变换后的Rosseland消光系数为
Figure FDA0004151929400000022
变换后的热导率为
Figure FDA0004151929400000023
变换后的密度和热容为
Figure FDA0004151929400000024
变换后的速度为v′f=Avf;detA为Jacobian变换矩阵的行列式,Aτ为Jacobian变换矩阵的转置。
CN201911026554.8A 2019-10-26 2019-10-26 一种同时适用于传导、对流和辐射的热隐身斗篷装置 Active CN110826271B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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