Inherent thermal-responsive strategies for safe lithium batteries

Safe batteries are the basis for next-generation application scenarios such as portable energy storage devices and electric vehicles,which are crucial to achieving carbon neutralization.Electrolytes,separators,and electrodes as main components of lithium batteries strongly affect the occurrence of safety accidents.Responsive materials,which can respond to external stimuli or environmental change,have triggered extensive attentions recently,holding great promise in facilitating safe and smart batteries.This review thoroughly discusses recent advances regarding the construction of high-safety lithium batteries based on internal thermal-responsive strategies,together with the corresponding changes in electrochemical performance under external stimulus.Furthermore,the existing challenges and outlook for the design of safe batteries are presented,creating valuable insights and proposing directions for the practical implementation of safe lithium batteries.

Mechanism of internal thermal runaway propagation in blade batteries

Blade batteries are extensively used in electric vehicles,but unavoidable thermal runaway is an inherent threat to their safe use.This study experimentally investigated the mechanism underlying thermal runaway propagation within a blade battery by using a nail to trigger thermal runaway and thermocouples to track its propagation inside a cell.The results showed that the internal thermal runaway could propagate for up to 272 s,which is comparable to that of a traditional battery module.The velocity of the thermal runaway propagation fluctuated between 1 and 8 mm s^(-1),depending on both the electrolyte content and high-temperature gas diffusion.In the early stages of thermal runaway,the electrolyte participated in the reaction,which intensified the thermal runaway and accelerated its propagation.As the battery temperature increased,the electrolyte evaporated,which attenuated the acceleration effect.Gas diffusion affected thermal runaway propagation through both heat transfer and mass transfer.The experimental results indicated that gas diffusion accelerated the velocity of thermal runaway propagation by 36.84%.We used a 1D mathematical model and confirmed that convective heat transfer induced by gas diffusion increased the velocity of thermal runaway propagation by 5.46%-17.06%.Finally,the temperature rate curve was analyzed,and a three-stage mechanism for internal thermal runaway propagation was proposed.In Stage I,convective heat transfer from electrolyte evaporation locally increased the temperature to 100℃.In Stage II,solid heat transfer locally increases the temperature to trigger thermal runaway.In StageⅢ,thermal runaway sharply increases the local temperature.The proposed mechanism sheds light on the internal thermal runaway propagation of blade batteries and offers valuable insights into safety considerations for future design.

Personal Thermal Management by Radiative Cooling and Heating

Maintaining thermal comfort within the human body is crucial for optimal health and overall well-being.By merely broadening the setpoint of indoor temperatures,we could significantly slash energy usage in building heating,ventilation,and air-conditioning systems.In recent years,there has been a surge in advancements in personal thermal management(PTM),aiming to regulate heat and moisture transfer within our immediate surroundings,clothing,and skin.The advent of PTM is driven by the rapid development in nano/micro-materials and energy science and engineering.An emerging research area in PTM is personal radiative thermal management(PRTM),which demonstrates immense potential with its high radiative heat transfer efficiency and ease of regulation.However,it is less taken into account in traditional textiles,and there currently lies a gap in our knowledge and understanding of PRTM.In this review,we aim to present a thorough analysis of advanced textile materials and technologies for PRTM.Specifically,we will introduce and discuss the underlying radiation heat transfer mechanisms,fabrication methods of textiles,and various indoor/outdoor applications in light of their different regulation functionalities,including radiative cooling,radiative heating,and dual-mode thermoregulation.Furthermore,we will shine a light on the current hurdles,propose potential strategies,and delve into future technology trends for PRTM with an emphasis on functionalities and applications.

Heat transfer enhanced inorganic phase change material compositing carbon nanotubes for battery thermal management and thermal runaway propagation mitigation

Developing technologies that can be applied simultaneously in battery thermal management(BTM)and thermal runaway(TR)mitigation is significant to improving the safety of lithium-ion battery systems.Inorganic phase change material(PCM)with nonflammability has the potential to achieve this dual function.This study proposed an encapsulated inorganic phase change material(EPCM)with a heat transfer enhancement for battery systems,where Na_(2)HPO_(4)·12H_(2)O was used as the core PCM encapsulated by silica and the additive of carbon nanotube(CNT)was applied to enhance the thermal conductivity.The microstructure and thermal properties of the EPCM/CNT were analyzed by a series of characterization tests.Two different incorporating methods of CNT were compared and the proper CNT adding amount was also studied.After preparation,the battery thermal management performance and TR propagation mitigation effects of EPCM/CNT were further investigated on the battery modules.The experimental results of thermal management tests showed that EPCM/CNT not only slowed down the temperature rising of the module but also improved the temperature uniformity during normal operation.The peak battery temperature decreased from 76℃to 61.2℃at 2 C discharge rate and the temperature difference was controlled below 3℃.Moreover,the results of TR propagation tests demonstrated that nonflammable EPCM/CNT with good heat absorption could work as a TR barrier,which exhibited effective mitigation on TR and TR propagation.The trigger time of three cells was successfully delayed by 129,474 and 551 s,respectively and the propagation intervals were greatly extended as well.

Development status and prospect of underground thermal energy storage technology

Underground Thermal Energy Storage(UTES)store unstable and non-continuous energy underground,releasing stable heat energy on demand.This effectively improve energy utilization and optimize energy allocation.As UTES technology advances,accommodating greater depth,higher temperature and multi-energy complementarity,new research challenges emerge.This paper comprehensively provides a systematic summary of the current research status of UTES.It categorized different types of UTES systems,analyzes the applicability of key technologies of UTES,and evaluate their economic and environmental benefits.Moreover,this paper identifies existing issues with UTES,such as injection blockage,wellbore scaling and corrosion,seepage and heat transfer in cracks,etc.It suggests deepening the research on blockage formation mechanism and plugging prevention technology,improving the study of anticorrosive materials and water treatment technology,and enhancing the investigation of reservoir fracture network characterization technology and seepage heat transfer.These recommendations serve as valuable references for promoting the high-quality development of UTES.

Prediction of the thermal conductivity of Mg–Al–La alloys by CALPHAD method

Mg-Al alloys have excellent strength and ductility but relatively low thermal conductivity due to Al addition.The accurate prediction of thermal conductivity is a prerequisite for designing Mg-Al alloys with high thermal conductivity.Thus,databases for predicting temperature-and composition-dependent thermal conductivities must be established.In this study,Mg-Al-La alloys with different contents of Al2La,Al3La,and Al11La3phases and solid solubility of Al in the α-Mg phase were designed.The influence of the second phase(s) and Al solid solubility on thermal conductivity was investigated.Experimental results revealed a second phase transformation from Al_(2)La to Al_(3)La and further to Al_(11)La_(3)with the increasing Al content at a constant La amount.The degree of the negative effect of the second phase(s) on thermal diffusivity followed the sequence of Al2La>Al3La>Al_(11)La_(3).Compared with the second phase,an increase in the solid solubility of Al in α-Mg remarkably reduced the thermal conductivity.On the basis of the experimental data,a database of the reciprocal thermal diffusivity of the Mg-Al-La system was established by calculation of the phase diagram (CALPHAD)method.With a standard error of±1.2 W/(m·K),the predicted results were in good agreement with the experimental data.The established database can be used to design Mg-Al alloys with high thermal conductivity and provide valuable guidance for expanding their application prospects.

Pre‑hatch thermal manipulation of embryos and post‑hatch baicalein supplementation mitigated heat stress in broiler chickens

Background High environmental temperatures induce heat stress in broiler chickens,affecting their health and pro-duction performance.Several dietary,managerial,and genetics strategies have been tested with some success in mitigating heat stress(HS)in broilers.Developing novel HS mitigation strategies for sustaining broiler production is critically needed.This study investigated the effects of pre-hatch thermal manipulation(TM)and post-hatch baica-lein supplementation on growth performance and health parameters in heat-stressed broilers.Results Six hundred fertile Cobb 500 eggs were incubated for 21 d.After candling on embryonic day(ED)10,238 eggs were thermally manipulated at 38.5℃ with 55%relative humidity(RH)from ED 12 to 18,then transferred to the hatcher(ED 19 to 21,standard temperature)and 236 eggs were incubated at a controlled temperature(37.5℃)till hatch.After hatch,180-day-old chicks from both groups were raised in 36 pens(n=10 birds/pen,6 replicates per treatment).The treatments were:1)Control,2)TM,3)control heat stress(CHS),4)thermal manipulation heat stress(TMHS),5)control heat stress supplement(CHSS),and 6)thermal manipulation heat stress supplement(TMHSS).All birds were raised under the standard environment for 21 d,followed by chronic heat stress from d 22 to 35(32–33℃ for 8 h)in the CHS,TMHS,CHSS,and TMHSS groups.A thermoneutral(22–24℃)environment was maintained in the Control and TM groups.RH was constant(50%±5%)throughout the trial.All the data were analyzed using one-way ANOVA in R and GraphPad software at P<0.05 and are presented as mean±SEM.Heat stress significantly decreased(P<0.05)the final body weight and ADG in CHS and TMHS groups compared to the other groups.Embryonic TM significantly increased(P<0.05)the expression of heat shock protein-related genes(HSP70,HSP90,and HSPH1)and antioxidant-related genes(GPX1 and TXN).TMHS birds showed a significant increment(P<0.05)in total cecal volatile fatty acid(VFA)concentration compared to the CHS birds.The cecal microbial analysis showed significant enrichment(P<0.05)in alpha and beta diversity and Coprococcus in the TMHSS group.Conclusions Pre-hatch TM and post-hatch baicalein supplementation in heat-stressed birds mitigate the detrimental effects of heat stress on chickens’growth performance,upregulate favorable gene expression,increase VFA produc-tion,and promote gut health by increasing beneficial microbial communities.

Ultrahigh performance passive radiative cooling by hybrid polar dielectric metasurface thermal emitters

Real-world passive radiative cooling requires highly emissive,selective,and omnidirectional thermal emitters to maintain the radiative cooler at a certain temperature below the ambient temperature while maximizing the net cooling power.Despite various selective thermal emitters have been demonstrated,it is still challenging to achieve these conditions sim-ultaneously because of the extreme difficulty in controlling thermal emission of photonic structures in multidimension.Here we demonstrated hybrid polar dielectric metasurface thermal emitters with machine learning inverse design,en-abling a high emissivity of~0.92 within the atmospheric transparency window 8-13μm,a large spectral selectivity of~1.8 and a wide emission angle up to 80 degrees,simultaneously.This selective and omnidirectional thermal emitter has led to a new record of temperature reduction as large as~15.4°C under strong solar irradiation of~800 W/m2,signific-antly surpassing the state-of-the-art results.The designed structures also show great potential in tackling the urban heat island effect,with modelling results suggesting a large energy saving and deployment area reduction.This research will make significant impact on passive radiative cooling,thermal energy photonics and tackling global climate change.

Thermal safety boundary of lithium-ion battery at different state of charge

Thermal runaway(TR)is a critical issue hindering the large-scale application of lithium-ion batteries(LIBs).Understanding the thermal safety behavior of LIBs at the cell and module level under different state of charges(SOCs)has significant implications for reinforcing the thermal safety design of the lithium-ion battery module.This study first investigates the thermal safety boundary(TSB)correspondence at the cells and modules level under the guidance of a newly proposed concept,safe electric quantity boundary(SEQB).A reasonable thermal runaway propagation(TRP)judgment indicator,peak heat transfer power(PHTP),is proposed to predict whether TRP occurs.Moreover,a validated 3D model is used to quantitatively clarify the TSB at different SOCs from the perspective of PHTP,TR trigger temperature,SOC,and the full cycle life.Besides,three different TRP transfer modes are discovered.The interconversion relationship of three different TRP modes is investigated from the perspective of PHTP.This paper explores the TSB of LIBs under different SOCs at both cell and module levels for the first time,which has great significance in guiding the thermal safety design of battery systems.

Efficient Electromagnetic Wave Absorption and Thermal Infrared Stealth in PVTMS@MWCNT Nano‑Aerogel via Abundant Nano‑Sized Cavities and Attenuation Interfaces

Pre-polymerized vinyl trimethoxy silane(PVTMS)@MWCNT nano-aerogel system was constructed via radical polymerization,sol-gel transition and supercritical CO_(2)drying.The fabricated organic-inorganic hybrid PVTMS@MWCNT aerogel structure shows nano-pore size(30-40 nm),high specific surface area(559 m^(2)g^(−1)),high void fraction(91.7%)and enhanced mechanical property:(1)the nano-pore size is beneficial for efficiently blocking thermal conduction and thermal convection via Knudsen effect(beneficial for infrared(IR)stealth);(2)the heterogeneous interface was beneficial for IR reflection(beneficial for IR stealth)and MWCNT polarization loss(beneficial for electromagnetic wave(EMW)attenuation);(3)the high void fraction was beneficial for enhancing thermal insulation(beneficial for IR stealth)and EMW impedance match(beneficial for EMW attenuation).Guided by the above theoretical design strategy,PVTMS@MWCNT nano-aerogel shows superior EMW absorption property(cover all Ku-band)and thermal IR stealth property(ΔT reached 60.7℃).Followed by a facial combination of the above nano-aerogel with graphene film of high electrical conductivity,an extremely high electromagnetic interference shielding material(66.5 dB,2.06 mm thickness)with superior absorption performance of an average absorption-to-reflection(A/R)coefficient ratio of 25.4 and a low reflection bandwidth of 4.1 GHz(A/R ratio more than 10)was experimentally obtained in this work.

Self‑Assembly of Binderless MXene Aerogel for Multiple‑Scenario and Responsive Phase Change Composites with Ultrahigh Thermal Energy Storage Density and Exceptional Electromagnetic Interference Shielding

The severe dependence of traditional phase change materials(PCMs)on the temperature-response and lattice deficiencies in versatility cannot satisfy demand for using such materials in complex application scenarios.Here,we introduced metal ions to induce the self-assembly of MXene nanosheets and achieve their ordered arrangement by combining suction filtration and rapid freezing.Subsequently,a series of MXene/K^(+)/paraffin wax(PW)phase change composites(PCCs)were obtained via vacuum impregnation in molten PW.The prepared MXene-based PCCs showed versatile applications from macroscale technologies,successfully transforming solar,electric,and magnetic energy into thermal energy stored as latent heat in the PCCs.Moreover,due to the absence of binder in the MXene-based aerogel,MK3@PW exhibits a prime solar-thermal conversion efficiency(98.4%).Notably,MK3@PW can further convert the collected heat energy into electric energy through thermoelectric equipment and realize favorable solar-thermal-electric conversion(producing 206 mV of voltage with light radiation intensity of 200 mw cm^(−2)).An excellent Joule heat performance(reaching 105℃with an input voltage of 2.5 V)and responsive magnetic-thermal conversion behavior(a charging time of 11.8 s can achieve a thermal insulation effect of 285 s)for contactless thermotherapy were also demonstrated by the MK3@PW.Specifically,as a result of the ordered arrangement of MXene nanosheet self-assembly induced by potassium ions,MK3@PW PCC exhibits a higher electromagnetic shielding efficiency value(57.7 dB)than pure MXene aerogel/PW PCC(29.8 dB)with the same MXene mass.This work presents an opportunity for the multi-scene response and practical application of PCMs that satisfy demand of next-generation multifunctional PCCs.

Early warning method for thermal runaway of lithium-ion batteries under thermal abuse condition based on online electrochemical impedance monitoring

Early warning of thermal runaway(TR)of lithium-ion batteries(LIBs)is a significant challenge in current application scenarios.Timely and effective TR early warning technology is urgently required considering the current fire safety situation of LIBs.In this work,we report an early warning method of TR with online electrochemical impedance spectroscopy(EIS)monitoring,which overcomes the shortcomings of warning methods based on traditional signals such as temperature,gas,and pressure with obvious delay and high cost.With in-situ data acquisition through accelerating rate calorimeter(ARC)-EIS test,the crucial features of TR were extracted using the RReliefF algorithm.TR mechanisms corresponding to the features at specific frequencies were analyzed.Finally,a three-level warning strategy for single battery,series module,and parallel module was formulated,which can successfully send out an early warning signal ahead of the self-heating temperature of battery under thermal abuse condition.The technology can provide a reliable basis for the timely intervention of battery thermal management and fire protection systems and is expected to be applied to electric vehicles and energy storage devices to realize early warning and improve battery safety.

Microstructure and thermal properties of dissimilar M300–CuCr1Zr alloys by multi-material laser-based powder bed fusion

Multi-material laser-based powder bed fusion (PBF-LB) allows manufacturing of parts with 3-dimensional gradient and additional functionality in a single step. This research focuses on the combination of thermally-conductive CuCr1Zr with hard M300 tool steel.Two interface configurations of M300 on CuCr1Zr and CuCr1Zr on M300 were investigated. Ultra-fine grains form at the interface due to the low mutual solubility of Cu and steel. The material mixing zone size is dependent on the configurations and tunable in the range of0.1–0.3 mm by introducing a separate set of parameters for the interface layers. Microcracks and pores mainly occur in the transition zone.Regardless of these defects, the thermal diffusivity of bimetallic parts with 50vol% of CuCr1Zr significantly increases by 70%–150%compared to pure M300. The thermal diffusivity of CuCr1Zr and the hardness of M300 steel can be enhanced simultaneously by applying the aging heat treatment.

Flexible and Robust Functionalized Boron Nitride/Poly(p‑Phenylene Benzobisoxazole)Nanocomposite Paper with High Thermal Conductivity and Outstanding Electrical Insulation

With the rapid development of 5G information technology,thermal conductivity/dissipation problems of highly integrated electronic devices and electrical equipment are becoming prominent.In this work,“high-temperature solid-phase&diazonium salt decomposition”method is carried out to prepare benzidine-functionalized boron nitride(m-BN).Subsequently,m-BN/poly(pphenylene benzobisoxazole)nanofiber(PNF)nanocomposite paper with nacremimetic layered structures is prepared via sol–gel film transformation approach.The obtained m-BN/PNF nanocomposite paper with 50 wt%m-BN presents excellent thermal conductivity,incredible electrical insulation,outstanding mechanical properties and thermal stability,due to the construction of extensive hydrogen bonds andπ–πinteractions between m-BN and PNF,and stable nacre-mimetic layered structures.Itsλ∥andλ_(⊥)are 9.68 and 0.84 W m^(-1)K^(-1),and the volume resistivity and breakdown strength are as high as 2.3×10^(15)Ωcm and 324.2 kV mm^(-1),respectively.Besides,it also presents extremely high tensile strength of 193.6 MPa and thermal decomposition temperature of 640°C,showing a broad application prospect in high-end thermal management fields such as electronic devices and electrical equipment.

Thermo-hydro-mechanical (THM) coupled simulation of the land subsidence due to aquifer thermal energy storage (ATES) system in soft soils

Aquifer thermal energy storage(ATES)system has received attention for heating or cooling buildings.However,it is well known that land subsidence becomes a major environmental concern for ATES projects.Yet,the effect of temperature on land subsidence has received practically no attention in the past.This paper presents a thermo-hydro-mechanical(THM)coupled numerical study on an ATES system in Shanghai,China.Four water wells were installed for seasonal heating and cooling of an agriculture greenhouse.The target aquifer at a depth of 74e104.5 m consisted of alternating layers of sand and silty sand and was covered with clay.Groundwater level,temperature,and land subsidence data from 2015 to 2017 were collected using field monitoring instruments.Constrained by data,we constructed a field scale three-dimensional(3D)model using TOUGH(Transport of Unsaturated Groundwater and Heat)and FLAC3D(Fast Lagrangian Analysis of Continua)equipped with a thermo-elastoplastic constitutive model.The effectiveness of the numerical model was validated by field data.The model was used to reproduce groundwater flow,heat transfer,and mechanical responses in porous media over three years and capture the thermo-and pressure-induced land subsidence.The results show that the maximum thermoinduced land subsidence accounts for about 60%of the total subsidence.The thermo-induced subsidence is slightly greater in winter than that in summer,and more pronounced near the cold well area than the hot well area.This study provides some valuable guidelines for controlling land subsidence caused by ATES systems installed in soft soils.

Miscibility of light oil and flue gas under thermal action

The miscibility of flue gas and different types of light oils is investigated through slender-tube miscible displacement experiment at high temperature and high pressure.Under the conditions of high temperature and high pressure,the miscible displacement of flue gas and light oil is possible.At the same temperature,there is a linear relationship between oil displacement efficiency and pressure.At the same pressure,the oil displacement efficiency increases gently and then rapidly to more than 90% to achieve miscible displacement with the increase of temperature.The rapid increase of oil displacement efficiency is closely related to the process that the light components of oil transit in phase state due to distillation with the rise of temperature.Moreover,at the same pressure,the lighter the oil,the lower the minimum miscibility temperature between flue gas and oil,which allows easier miscibility and ultimately better performance of thermal miscible flooding by air injection.The miscibility between flue gas and light oil at high temperature and high pressure is more typically characterized by phase transition at high temperature in supercritical state,and it is different from the contact extraction miscibility of CO_(2) under conventional high pressure conditions.

Improvement strategy on thermophysical properties of A_(2)B_(2)O_(7)-type rare earth zirconates for thermal barrier coatings applications:A review

The A_(2)B_(2)O_(7)-type rare earth zirconate compounds have been considered as promising candidates for thermal barrier coating(TBC) materials because of their low sintering rate,improved phase stability,and reduced thermal conductivity in contrast with the currently used yttria-partially stabilized zirconia (YSZ) in high operating temperature environments.This review summarizes the recent progress on rare earth zirconates for TBCs that insulate high-temperature gas from hot-section components in gas turbines.Based on the first principles,molecular dynamics,and new data-driven calculation approaches,doping and high-entropy strategies have now been adopted in advanced TBC materials design.In this paper,the solid-state heat transfer mechanism of TBCs is explained from two aspects,including heat conduction over the full operating temperature range and thermal radiation at medium and high temperature.This paper also provides new insights into design considerations of adaptive TBC materials,and the challenges and potential breakthroughs are further highlighted for extreme environmental applications.Strategies for improving thermophysical performance are proposed in two approaches:defect engineering and material compositing.

Coupling effect on the thermal hazard assessment of hazardous chemical materials via calorimetric technologies and simulation approaches

The coupling effect of heat absorption and release exists in the thermal decomposition of a few chemical materials.However,the impact of the above coupling on thermal hazard assessment is not considered in the literature studies.In this work,nitroguanidine(NQ)and 1,3,5-trinitro-1,3,5-triazine(RDX)are selected as representative materials to explore the influence of the coupling effect on the thermal hazard assessment of chemical materials.The linear heating experiments of NQ and RDX are carried out by a microcalorimeter and synchronous thermal analyser.The thermal decomposition curves are decoupled by advanced thermokinetics software.The thermal decomposition and kinetic parameters before and after decoupling are calculated.The results of TG experiment show that both NQ and RDX began to lose mass during the endothermic stage.The endothermic melting and exothermic decomposition of NQ and RDX are coupled within this stage.The coupling effect has different degrees of influence on its initial decomposition temperature and safety parameters.Compared with the parameters in the coupling state,the initial decomposition temperature and adiabatic induction period after decoupling decrease.The self-accelerating decomposition temperature increases,and internal thermal runaway time decreases.In the thermal hazard assessment of chemical materials with coupling effects,the calculated parameters after decoupling should be taken as an important safety index。

General three-dimensional thermal illusion metamaterials

Thermal illusion aims to create fake thermal signals or hide the thermal target from the background thermal field to mislead infrared observers,and illusion thermotics was proposed to regulate heat flux with artificially structured metamaterials for thermal illusion.Most theoretical and experimental works on illusion thermotics focus on two-dimensional materials,while heat transfer in real three-dimensional(3D)objects remains elusive,so the general 3D illusion thermotics is urgently demanded.In this study,we propose a general method to design 3D thermal illusion metamaterials with varying illusions at different sizes and positions.To validate the generality of the 3D method for thermal illusion metamaterials,we realize thermal functionalities of thermal shifting,splitting,trapping,amplifying and compressing.In addition,we propose a special way to simplify the design method under the condition that the size of illusion target is equal to the size of original heat source.The 3D thermal illusion metamaterial paves a general way for illusion thermotics and triggers the exploration of illusion metamaterials for more functionalities and applications.

A graph neural network approach to the inverse design for thermal transparency with periodic interparticle system

Recent years have witnessed significant advances in utilizing machine learning-based techniques for thermal metamaterial-based structures and devices to attain favorable thermal transport behaviors.Among the various thermal transport behaviors,achieving thermal transparency stands out as particularly desirable and intriguing.Our earlier work demonstrated the use of a thermal metamaterial-based periodic interparticle system as the underlying structure for manipulating thermal transport behavior and achieving thermal transparency.In this paper,we introduce an approach based on graph neural network to address the complex inverse design problem of determining the design parameters for a thermal metamaterial-based periodic interparticle system with the desired thermal transport behavior.Our work demonstrates that combining graph neural network modeling and inference is an effective approach for solving inverse design problems associated with attaining desirable thermal transport behaviors using thermal metamaterials.