Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Thermal properties of field-assisted-sintered SiCN–Y 2 O 3 composites

Polymer-derived amorphous SiCN has excellent high-temperature stability and properties. To reduce the shrinkage during pyrolysis and to improve the high-temperature oxidation resistance, Y 2 O 3 was added as a filler. In this study, polymer-derived SiCN–Y 2 O 3 composites were fabricated by mixing a polymeric precursor of SiCN with Y 2 O 3 submicron powders in different ratios. The mixtures were cross-linked and pyrolyzed in argon. SiCN–Y 2 O 3 composites were processed using field-assisted sintering technology at 1350°C for 5 min under vacuum. Dense SiCN–Y 2 O 3 composite pellets were successfully made with relative density higher than 98% and homogeneous microstructure. Due to low temperature and short time of the heat-treatment, the grain growth of Y 2 O 3 was substantially inhibited. The Y 2 O 3 grain size was ~1 μm after sintering. The composites’ heat capacity, thermal diffusivity, and thermal expansion coefficients were characterized as a function of temperature. The thermal conductivity of the composites ceramics decreased as the amount of amorphous SiCN increased and the coefficient of thermal expansion (CTE) of the composites increased with Y 2 O 3 content. However, the thermal conductivity and CTE did not follow the rule of mixture. This is likely due to the partial oxidation of SiCN and the resultant impurity phases such as Y 2 SiO 5 , Y 2 Si 2 O 7 , and Y 4.67 (SiO 4 ) 3 O.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

Novel Ordered Crown Ether-Containing Polyimides for Ion Conduction

We report the synthesis and characterization of thermally-stable polyimides for use as battery and fuel cell electrolyte membranes. Dianhydrides used were 1,4,5,8- naphthalenetetracarboxylic dianhydride and/or 4,4'-(hexafluoroisopropylidene)diphthalic anhydride. Diamines used were anti-4,4-diaminodibenzo-l8-crown-6, 4,4'- diaminodibenzo-24-crown-8, 2,2-bis(4-aminophenyl)hexafluoropropane, and/or 2,5- diaminobenzenesulfonic acid. The polymers were characterized using electrochemical impedance spectroscopy (EIS), thermal analysis and X-ray diffraction. Polymers containing the hexafluoroisopropylidene (HFIP) group were soluble in common organic solvents, while polymers without the HFIP group were very poorly soluble. Sulfonation yields polymers that are sparingly soluble in aqueous base and/or methanol. Degree of sulfonation, determined by titration, was between one and three sulfonate groups per repeat unit. Proton conductivity was determined as a function of water content, with a maximum conductivity of l x 10(exp -2) per centimeter when fully hydrated. Crown ether-containing polymers exhibit a high degree of order that may be indicative of crown ether channel formation, which may facilitate Li(+) transport for use in battery membranes.

Irvin, Jennifer A.↗
DOE OSTI2022

Approximate bounds and temperature dependence of adiabatic connection integrands for the uniform electron gas

Thermal density functional theory is commonly used in simulations of warm dense matter, a highly energetic phase characterized by substantial thermal effects and by correlated electrons demanding quantum mechanical treatment. Methods that account for temperature dependence, such as Mermin–Kohn–Sham finite-temperature density functional theory and free energy density functional theory, are now employed with more regularity and available in many standard code packages. However, approximations from zero-temperature density functional theory are still often used in temperature-dependent simulations using thermally weighted electronic densities as an input to exchange–correlation functional approximations, a practice known to miss temperature-dependent effects in the exchange–correlation free energy of these systems. Here, in this work, the temperature-dependent adiabatic connection is demonstrated and analyzed using a well-known parameterization of the uniform electron gas free energy. Useful tools based on this formalism for analyzing and constraining approximations of the exchange–correlation at zero temperature are leveraged for the finite-temperature case. Inspired by the Lieb–Oxford inequality, which provides a lower bound for the ground-state exchange–correlation energy, bounds for the exchange–correlation at finite temperatures are approximated for various degrees of electronic correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Thermal conductance of pressed contacts at liquid helium temperatures

It is pointed out that the optimum design of cryogenic instruments requires accurate thermal models. The present models are limited by a lack of knowledge of the low temperature thermal conductance of the bolted joints which are typically used in the instrument-to-system interface. In connection with studies of pressed contacts, it has been found that the thermal conductance does not obey the Wiedemann-Franz law. The present investigation is concerned with the characterization of the thermal conductance of pressed contacts at liquid helium-4 temperatures, taking into account the dependence of thermal contact conductance on applied force and temperature. It is shown that for the 0.4 micron OFHC copper pressed contact pair, the thermal conductance varies roughly as the second power of the temperature, and increases with increasing applied force.

Salerno, L. J.↗
NASA NTRS2026

The Relevance of the VIPER Mission to NASA’s Artemis Human Exploration of the Moon

NASA has ambitious plans to send astronaut crews to the south polar region of the Moon and explore this lunar terrain with humans for the first time through a series of Artemis missions. Artemis will provide economic benefits, drive technological advancement, and inspire the next generation of explorers. Simultaneously, Artemis presents a unique opportunity for humans to conduct high-priority planetary science in situ in the south polar region of the Moon. Prior to human exploration, landed robotic precursor investigations can provide valuable information to reduce risk and maximize discovery and productivity of subsequent crewed missions. NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) is a robotic mission designed to explore and characterize the lunar south polar region prior to crewed Artemis landed missions. VIPER is designed to explore multiple thermal regimes to characterize the lunar polar environment and regolith properties, and to explore for surface and subsurface ices. Each of these objectives are key for both scientific investigation and to detect and assess volatile deposits to support in situ resource utilization (ISRU) and a sustained human presence on the Moon. VIPER would also operate with a real-time mission operations architecture with relevance and feed-forward to Artemis crewed operations on the lunar surface. VIPER would provide key information regarding the lunar environment, scientific exploration, technologies, and real-time operations to optimize the valuable surface exploration time of Artemis crew members while simultaneously reducing risk to crew and increasing health and human safety on the lunar surface.

Jennifer L Heldmann↗
NASA NTRS2024

Large Interferometer For Exoplanets (LIFE) XIII. The value of Combining Thermal Emission and Reflected Light for the Characterization of Earth Twins

Context. Following the recommendations to NASA (in the Astro2020 Decadal survey) and ESA (through the Voyage2050 process), the search for life on exoplanets will be a priority in the next decades. Two concepts for direct imaging space missions are being developed for this purpose: the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). These two concepts operate in different spectral regimes: HWO is focused on reflected light spectra in the ultraviolet, visible, and near-infrared (UV/VIS/NIR), while LIFE will operate in the mid-infrared (MIR) to capture the thermal emission of temperate exoplanets. Aims. In this study, we aim to assess the potential of HWO and LIFE to characterize a cloud-free Earth twin orbiting a Sun-like star at a distance of 10 pc, both as separate missions and in synergy with each other. We aim to quantify the increase in information that can be gathered by joint atmospheric retrievals on a habitable planet. Methods. We performed Bayesian retrievals on simulated data obtained by an HWO-like mission and a LIFE-like one separately, then jointly. We considered the baseline spectral resolutions currently assumed for these concepts and used two increasingly complex noise simulations, obtained using state-of-the-art noise simulators. Results. An HWO-like concept would allow one to strongly constrain H 2 O, O 2 , and O 3 in the atmosphere of a cloud-free Earth twin, while the atmospheric temperature profile is not well constrained (with an average uncertainty ≈100 K). LIFE-like observations would strongly constrain CO 2 , H 2 O, and O 3 and provide stronger constraints on the thermal atmospheric structure and surface temperature (down to ≈10 K uncertainty). For all the investigated scenarios, both missions would provide an upper limit on CH 4 . A joint retrieval on HWO and LIFE data would accurately define the atmospheric thermal profile and planetary parameters. It would decisively constrain CO 2 , H 2 O, O 2 , and O 3 and find weak constraints on CO and CH 4 . The significance of the detection is in all cases greater than or equal to the single-instrument retrievals. Conclusions. Both missions provide specific information that is relevant for the characterization of a terrestrial habitable exoplanet, but the scientific yield can be maximized by considering synergistic studies of UV/VIS/NIR+MIR observations. The use of HWO and LIFE together will provide stronger constraints on biosignatures and life indicators, with the potential to be transformative for the search for life in the Universe.

E Alei↗
DOE OSTI2026

A Physics-Informed Reinforcement Learning Framework for Economic-Thermal Co-Optimization of Crypto Mining Data Centers: Preprint

The rapid expansion of cryptocurrency mining has created a new class of high-density data centers characterized by extreme thermal flux and high sensitivity to volatile economic markets. Traditional thermal management strategies, typically reliant on rule-based control, maintain static setpoints that fail to account for fluctuating electricity prices and cryptocurrency values - factors critical to mining profitability. To address this, we present a physics-informed reinforcement learning (PIRL) framework for economic-thermal co-optimization in crypto mining data centers. This framework consists of a proximal policy optimization (PPO) agent, a virtual testbed powered by high-fidelity physics-based models, and an interactive frontend dashboard. The PPO agent is trained using the virtual testbed and strict hardware safety limits. This physics-informed approach allows the agent to learn a stochastic policy that dynamically balances mining revenue against operational costs by co-optimizing HVAC cooling setpoints and IT computational hashrate. The simulation results demonstrate that the integrated framework achieved an 8.62% increase in net operational profit compared to traditional baseline strategies while strictly adhering to safety-critical temperature constraints (coolant supply temperature < 32 degrees C). This work provides a scalable template for the deployment of reinforcement learning in mission critical facilities where economic volatility and physical safety must be managed simultaneously.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

In Situ Characterization of Ferroelectric HfO 2 During Rapid Thermal Annealing

Thermal annealing is critical in governing the phase distribution, and ultimately the electrical properties, of ferroelectric hafnium oxide films. In particular, rapid thermal annealing (RTA) has been shown to favor the formation of ferroelectric crystal phases, but the dynamic behavior of the film over the course of the anneal is not well understood. Herein, synchrotron X-ray diffraction is used to characterize the phase distribution of HfO 2 films deposited by atomic layer deposition (ALD) in situ during RTA, revealing complex phase transformations occurring on the scale of seconds. All samples investigated here transform into a nonmonoclinic phase, which is required for ferroelectric films. However, this phase often converts into the more stable monoclinic phase as annealing proceeds. The kinetic barrier to the transformation to the monoclinic phase during heating is higher than 1 eV f.u. -1 In this work, the initial crystallization into the ferroelectric phase and relatively high barrier to transformation to the nonferroelectric phase suggests that careful control of the thermal annealing profile can greatly increase the ferroelectric fraction of the film. By using a fast ramp rate and a short annealing time above the ferroelectric crystallization temperature, the remanent polarization of a pure HfO 2 film can be increased more than twofold.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Theory of local thermal instability in spherical systems

The gasdynamical properties of local thermal instability in optically thin astrophysical plasmas as it occurs in spherical accretion and winds is investigated. In a medium characterized by both thermal and hydrostatic equilibrium, if the cooling function is not an explicit function of position and does not display isentropic thermal instability, then isobaric thermal instability by the Field criterion is present if and only if convective instability is present by the Schwarzschild criterion. In this case, thermal overstability cannot occur. Convective instability by the Schwarzschild criterion will also occur in accretion flows locally dominated by external heating or in marginally unbound, radiatively cooling outflows. A very general Lagrangian equation for the development of nonradial thermal instability in flows with spherical symmetry is derived and is solved analytically in certain regimes. The results are applied to cluster X-ray cooling flows.

Balbus, Steven A.↗
NASA NTRS2001

Remote Sensing of Atlanta's Urban Sprawl and the Distribution of Land Cover and Surface Temperature

Between 1973 and 1992, an average of 20 ha of forest was lost each day to urban expansion of Atlanta, Georgia. Urban surfaces have very different thermal properties than natural surfaces-storing solar energy throughout the day and continuing to release it as sensible heat well after sunset. The resulting heat island effect serves as catalysts for chemical reactions from vehicular exhaust and industrialization leading to a deterioration in air quality. In this study, high spatial resolution multispectral remote sensing data has been used to characterize the type, thermal properties, and distribution of land surface materials throughout the Atlanta metropolitan area. Ten-meter data were acquired with the Advanced Thermal and Land Applications Sensor (ATLAS) on May 11 and 12, 1997. ATLAS is a 15-channel multispectral scanner that incorporates the Landsat TM bands with additional bands in the middle reflective infrared and thermal infrared range. The high spatial resolution permitted discrimination of discrete surface types (e.g., concrete, asphalt), individual structures (e.g., buildings, houses) and their associated thermal characteristics. There is a strong temperature contrast between vegetation and anthropomorphic features. Vegetation has a modal temperature at about 20 C, whereas asphalt shingles, pavement, and buildings have a modal temperature of about 39 C. Broad-leaf vegetation classes are indistinguishable on a thermal basis alone. There is slightly more variability (+/-5 C) among the urban surfaces. Grasses, mixed vegetation and mixed urban surfaces are intermediate in temperature and are characterized by broader temperature distributions with modes of about 29 C. Thermal maps serve as a basis for understanding the distribution of "hotspots", i.e., how landscape features and urban fabric contribute the most heat to the lower atmosphere.

Laymon, Charles A.↗
NASA NTRS2001

Remote Sensing of Atlanta's Urban Sprawl and the Distribution of Land Cover and Surface Temperatures

Between 1973 and 1992, an average of 20 ha of forest was lost each day to urban expansion of Atlanta, Georgia. Urban surfaces have very different thermal properties than natural surfaces-storing solar energy throughout the day and continuing to release it as sensible heat well after sunset. The resulting heat island effect serves as catalysts for chemical reactions from vehicular exhaust and industrialization leading to a deterioration in air quality. In this study, high spatial resolution multispectral remote sensing data has been used to characterize the type, thermal properties, and distribution of land surface materials throughout the Atlanta metropolitan area. Ten-meter data were acquired with the Advanced Thermal and Land Applications Sensor (ATLAS) on May 11 and 12, 1997. ATLAS is a 15-channel multispectral scanner that incorporates the Landsat TM bands with additional bands in the middle reflective infrared and thermal infrared range. The high spatial resolution permitted discrimination of discrete surface types (e.g., concrete, asphalt), individual structures (e.g., buildings, houses) and their associated thermal characteristics. There is a strong temperature contrast between vegetation and anthropomorphic features. Vegetation has a modal temperature at about 20 C, whereas asphalt shingles, pavement, and buildings have a modal temperature of about 39 C. Broad-leaf vegetation classes are indistinguishable on a thermal basis alone. There is slightly more variability (plus or minus 5 C) among the urban surfaces. Grasses, mixed vegetation and mixed urban surfaces are intermediate in temperature and are characterized by broader temperature distributions with modes of about 29 C. Thermal maps serve as a basis for understanding the distribution of "hotspots", i.e., how landscape features and urban fabric contribute the most heat to the lower atmosphere.

Laymon, Charles A.↗
NASA NTRS1987

Shuttle payload bay thermal environments: Summary and conclusion report for STS Flights 1-5

The thermal data for the payload bay of the first five shuttle flights is summarized and the engineering evaluation of that data is presented. After a general discussion on mission profiles and vehicle configurations, the thermal design and flight instrumentation systems of the payload bay are described. The thermal flight data sources and a categorization of the data are then presented. A thermal flight data summarization section provides temperature data for the five phases of a typical mission profile. These are: prelaunch, ascent, on-orbit, entry and postlanding. The thermal flight data characterization section encompasses this flight data for flight to flight variations, payload effects, temperature ranges, and other variations. Discussion of the thermal environment prediction models in use by industry and various NASA Centers, and the results predicted by these models, is followed by an evaluation of the correlation between the actual flight data and the results predicted by the models. Finally, the available thermal data are evaluated from the viewpoint of the user concerned with establishing the thermal environment in the payload bay. The data deficiencies are discussed and recommendations for their elimination are presented.

Fu, J. H.↗
NASA NTRS2002

Dynamic Characterization of an Inflatable Concentrator for Solar Thermal Propulsion

Solar-thermal propulsion is a concept for producing thrust sufficient for orbital transfers and requires innovative, lightweight structures. This note presents a description of an inflatable concentrator that consists of a torus, lens simulator, and three tapered struts. Modal testing was discussed for characterization and verification of the solar concentrator assembly. Finite element shell models of the concentrator were developed using a two-step nonlinear approach, and results were compared to test data. Reasonable model-to-test agreement was achieved for the torus, and results for the concentrator assembly were comparable to the test for several modes.

Leigh, Larry M.↗
NASA NTRSDate not supplied

Thermal Protection System P50 Cork Char Response Characterization to Ascent Flight Shear and Heating Environments

NASA’s Space Launch System (SLS) experienced deflagration of the acreage of the Core Stage base heat shield during the Green Run hot-fire test campaigns and Artemis I flight. The burning products from the base heat shield during Green Run hot-fire test operation led to concerns in aerothermal environments and thermal protection system (TPS) performance prior to Artemis I flight. Updated design thermal environments accounting for cork combustion were developed using Green Run measurements and assessed prior to Artemis I, resulting in increased TPS thickness for the base heat shield. This work documents a highly controlled ground test within the NASA Marshall Space Flight Center’s Propulsion Test Branch Hot Gas Facility (HGF) and investigates the P50 cork TPS response due to sensitivity changes in radiative heat rates and surface shear stress. Prior to testing at the HGF, most of the P50 cork panels had a thick layer of char due to the cork combustion process replicated in a thermal vacuum chamber. SLS ascent shear stress and heating rate profile were simulated within the HGF by running the combustion-driven wind tunnel at off-nominal conditions and Mach numbers below 1. This test data provided insight into aerothermal measurements from Artemis I flight as well as informing updated base heating design environments for Artemis II. An in-depth understanding of the TPS P50 cork char response and observations is presented for launch vehicle flight ascent environments.

aerothermodynamics↗
DOE OSTI2024

Exchange–correlation entropy from the generalized thermal adiabatic connection

Warm dense matter is a highly energetic phase characterized by strong correlations, thermal effects, and quantum mechanical electrons. Thermal density functional theory is commonly used in simulations of this challenging phase, driving the development of temperature-dependent approximations to the exchange–correlation free energy. Approaches using the adiabatic connection formula are well known at zero temperature and have been recently leveraged at non-zero temperatures as well. In this work, a generalized thermal adiabatic connection (GTAC) formula is proposed, introducing a fictitious temperature parameter. This allows extraction of the exchange–correlation entropy S XC using simulated interaction strength scaling. Here, this procedure uses a Hellmann–Feynman approach to express the exchange–correlation entropy in terms of a temperature- and interaction strength-dependent exchange–correlation potential energy. In addition, analysis of S XC as a function of interaction strength suggests new forms for approximations, and GTAC itself offers a new framework for exploring both the exact and approximate interplay of temperature, density, and interaction strength across a wide range of conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Comparative experimental study of heat transfer processes in accumulation energy recovery exchangers

In this study, an experimental comparison between three different energy accumulating ceramic heat exchangers for energy recovery in ventilation systems was presented. The units were selected to represent three different approaches on the energy recovery process: honeycomb structure with more accumulation mass (more energy can be stored in one unit)- Unit 1.1, honeycomb structure with lower accumulation mass- Unit 1.2, and a rectangular structure with expanded heat transfer surface- Unit 1.3. The achieved results are useful for the future application of such units in ventilation systems. It was established that all evaluated units demonstrated an acceptable effectiveness of thermal energy recovery from the exhaust air. Their average energy recovery efficiency ranged between 70 % and 80 %, aligning with expected performance benchmarks for regenerative heat exchangers employed in contemporary mechanical ventilation systems. It was also established that the factor which has the highest impact on thermal effectiveness is the heat transfer surface available in the tested heat exchangers. Units with the highest number of channels (i.e., with the highest amount of single channel walls, which can be used for heat transfer) achieve the highest thermal effectiveness. However, higher thermal effectiveness can negatively affect the ventilation potential of the units. Unit 1.3, which was characterized by the highest thermal efficiency, was also characterized by the lowest achievable flow rate. It was also found that the important technical aspect that should be taken into account when analyzing accumulation energy recovery units is also the fluctuation of the supply air temperature. The ability to ensure minimal temperature fluctuations is a significant operational advantage, as it ensures a higher level of safety at lower outside temperatures.

Energy recovery↗
NASA NTRS1990

Thermal Conductivity and Temperature Limits of Electron-Beam Physical-Vapor-Deposited ZrO2-7wt%Y2O3 Thermal Barrier Coating

The electron-beam physical-vapor-deposited (EB-PVD) ZrO2-7wt%Y2O3 thermal barrier coating has been widely used to protect engine hot section components in modem aircraft engines. Thermal conductivity and the conductivity increases due to sintering and phase changes are important coating design parameters. There is a need to characterize the coating thermal conductivity behavior and temperature limits, in order to potentially take full advantage of the coating capability. In addition, since the ZrO2-7wt%Y2O3 coating is often used as a baseline coating for the development of advanced lower thermal conductivity and higher temperature capability coatings, a thorough evaluation of the coating conductivity behavior at future higher engine operating temperatures will be useful for more accurately assessing the benefit gained from the new coating systems. In this study, thermal conductivity behavior of EB-PVD ZrO2-7wt%Y2O3 coating has been systematically investigated as a function of temperature (up to 2600 F) and time under steady-state and cyclic test conditions using a laser heat-flux technique. Thermal conductivity change kinetics were determined under realistic engine high temperature thermal gradient conditions, and correlated to the coating microstructural and phase changes. The thermal conductivity prediction models have been established in terms of heat flux, time, and testing temperatures.

Zhu, Dongming↗