Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal response functions”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗
DOE OSTI2023

Multi-Functional Applications of H-Glass Embedded with Stable Plasmonic Gold Nanoislands

Metal nanoparticles (MNPs) are synthesized using various techniques on diverse substrates that significantly impact their properties. However, among the substrate materials investigated, the major challenge is the stability of MNPs due to their poor adhesion to the substrate. Herein, it is demonstrated how a newly developed H-glass can concurrently stabilize plasmonic gold nanoislands (GNIs) and offer multifunctional applications. The GNIs on the H-glass are synthesized using a simple yet, robust thermal dewetting process. The H-glass embedded with GNIs demonstrates versatility in its applications, such as i) acting as a room temperature chemiresistive gas sensor (70% response for NO 2 gas); ii) serving as substrates for surface-enhanced Raman spectroscopy for the identifications of Nile blue (dye) and picric acid (explosive) analytes down to nanomolar concentrations with enhancement factors of 4.8 × 10 6 and 6.1 × 10 5 , respectively; and iii) functioning as a nonlinear optical saturable absorber with a saturation intensity of 18.36 × 10 15 W m -2 at 600 nm, and the performance characteristics are on par with those of materials reported in the existing literature. In conclusion, this work establishes a facile strategy to develop advanced materials by depositing metal nanoislands on glass for various functional applications.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Emerging satellite observations for diurnal cycling of ecosystem processes

Diurnal cycling of plant carbon uptake and water use, and their responses to water and heat stresses, provide direct insight into assessing ecosystem productivity, agricultural production and management practices, carbon and water cycles, and feedbacks to the climate. Temperature, light, atmospheric water demand, soil moisture and leaf water potential vary over the course of the day, leading to diurnal variations in stomatal conductance, photosynthesis and transpiration. Earth observations from polar-orbiting satellites are incapable of studying these diurnal variations. Here, we review the emerging satellite observations that have the potential for studying how plant functioning and ecosystem processes vary over the course of the diurnal cycle. The recently launched ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) and Orbiting Carbon Observatory-3 (OCO-3) provide land surface temperature, evapotranspiration (ET), gross primary production (GPP) and solar-induced chlorophyll fluorescence data at different times of day. New generation operational geostationary satellites such as Himawari-8 and the GOES-R series can provide continuous, high-frequency data of land surface temperature, solar radiation, GPP and ET. Future satellite missions such as GeoCarb, TEMPO and Sentinel-4 are also planned to have diurnal sampling capability of solar-induced chlorophyll fluorescence. We explore the unprecedented opportunities for characterizing and understanding how GPP, ET and water use efficiency vary over the course of the day in response to temperature and water stresses, and management practices. We also envision that these emerging observations will revolutionize studies of plant functioning and ecosystem processes in the context of climate change and that these observations and findings can inform agricultural and forest management and lead to improvements in Earth system models and climate projections.

Plant functioning↗
NASA NTRS1971

A cardiovascular system model for lower-body negative pressure response

Mathematical models used to study complex physiological control systems are discussed. Efforts were made to modify a model of the cardiovascular system for use in studying lower body negative pressure. A computer program was written which allows orderly, straightforward expansion to include exercise, metabolism (thermal stress), respiration, and other body functions.

Mitchell, B. A., Jr.↗
NASA NTRS2001

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗
NASA NTRS2003

Radial Clearance Found To Play a Key Role in the Performance of Compliant Foil Air Bearings

Compliant foil air bearings are at the forefront of the Oil-Free turbomachinery revolution, which supports gas turbine engines with hydrodynamic bearings that use air instead of oil as the working fluid. These types of bearings have been around for almost 50 years and have found a home in several commercial applications, such as in air cycle machines, turbocompressors, and microturbines, but are now being aggressively pursued for use in small and midrange aircraft gas turbine engines. Benefits include higher operating speeds and temperatures, lower maintenance costs, and greater reliability. The Oil-Free Turbomachinery team at the NASA Glenn Research Center is working to foster the transition of Oil-Free technology into gas turbine engines by performing in-house experiments on foil air bearings in order to gain a greater insight into their complex operating principles. A research program recently undertaken at Glenn focused on the concept of radial clearance and its influence on bearing performance. The tests were conducted on foil bearings with different radial clearances. As defined for a foil bearing, radial clearance is a measure of the small amount of shaft radial motion that is present from play that exists in the elastic support structure, such as between the top and bump foils and the bump foils and bearing shell (see the drawing). With an insufficient amount of radial clearance, the bearing imparts a high preload on the shaft, which when excessive, can reduce the loadcarrying capability of the bearing. On the other hand, systems using foil bearings with excessive radial clearance may experience rotordynamic instabilities because of low bearing preload. Therefore, without a more thorough understanding of radial clearance, it is difficult to accurately predict the performance of a given bearing design. The test program demonstrated that there is a direct correlation between radial clearance and the performance of foil air bearings. As shown in the graph, an optimum radial clearance exists that will maximize the amount of load that the bearing is capable of supporting. With respect to this optimum, two different performance regimes were observed that are a function of the amount of radial clearance. Tests showed that bearings with radial clearances below the optimum in regime I were susceptible to sudden seizure, a failure mode indicative of thermal runaway caused by high preload. The high preload is in response to an insufficient amount of radial clearance available to accommodate the thermal growth of the bearing and shaft. However, radial clearances greater than the optimum in regime II resulted in low bearing preloads that did not cause any heat-related problems, and the failure mode was due to fluid-film breakdown. In fact, bearings operating with radial clearances twice as much as the optimum suffered a decrease in the maximum load capacity of only about 20 percent. Therefore, special attention has to be given to the range of operating conditions expected in the bearing/shaft system since changes in temperature, centrifugal, and hydrodynamic effects can all affect radial clearance. This enhanced understanding of foil air bearing behavior will greatly aid our efforts to transition Oil-Free technology to future aircraft engines.

Radil, Kevin C.↗
DOE OSTI2025

Undercooling minimization in ultrasound coupled DTA measurements of molten salts

Molten salts are ionic liquids that are used for the electrolytic pyroprocessing of metals and as heat transfer fluids in very high temperature processes. Recently, halide-based molten salt reactors (MSR) have gained momentum for high density and environmentally responsible electricity generation. The function of these reactors and their fuel cycle depend on a knowledge of the halide salt’s thermodynamic properties [1]. Therefore, high accuracy phase equilibria of MSR relevant base halide and actinide containing salts are needed. Modern thermal analysis of the phase transitions of halide salts is usually done during heating at relatively high scan rates with commercial devices. Normally such measurements are adequate for pure or pseudo-binary salts. However, as the number of components in the mixture increases, accurately resolving the liquidus becomes increasingly difficult. Reversing the scanning mode greatly increases the sensitivity of phase transition measurements but can decrease accuracy due to undercooling [2] — a common occurrence in molten halide systems [3]. This work presents the development of a differential thermal analysis (DTA) cell intended for use in radiological gloveboxes. Non-contact ultrasonic agitation of the halide salt is implemented to limit kinetic limitations on crystallization during cooling to minimize undercooling [4]. Measurements on halide salts are also presented to elucidate the effect of non-contact mixing on undercooling. Reference [1] S. Boyd and C. Taylor, “3 - Chemical fundamentals and applications of molten salts,” in Molten Salt Reactors and Thorium Energy, T. J. Dolan, Ed., Woodhead Publishing, 2017, pp. 29–91. doi: 10.1016/B978-0-08-101126-3.00003-8. [2] K. Nitsch, A. Cihlár, and M. Rodová, “Molten state and supercooling of lead halides,” J. Cryst. Growth, vol. 264, no. 1, pp. 492–498, Mar. 2004, doi: 10.1016/j.jcrysgro.2004.01.011. [3] L. Rycerz, “Practical remarks concerning phase diagrams determination on the basis of differential scanning calorimetry measurements,” J. Therm. Anal. Calorim., vol. 113, no. 1, pp. 231–238, Jul. 2013, doi: 10.1007/s10973-013-3097-0. [4] Md. H. Zahir, S. A. Mohamed, R. Saidur, and F. A. Al-Sulaiman, “Supercooling of phase-change materials and the techniques used to mitigate the phenomenon,” Appl. Energy, vol. 240, pp. 793–817, Apr. 2019, doi: 10.1016/j.apenergy.2019.02.045.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Laser Doppler vibrometry for piezoelectric coefficient ($d_{33}$) measurements in irradiated aluminum nitride

Sensors used for experiments in advanced reactors must survive in harsh environments. Few material systems can be used to construct sensors viable for extreme conditions. Aluminum nitride (AlN) is one such material because it has high thermal stability and radiation resistance and sustains good piezoelectric and dielectric properties at high temperatures. In this work, the piezoelectric coefficient $d_{33}$ of the AlN single-crystal with thermal and irradiation damage was investigated using an indirect method with a laser Doppler vibrometer (LDV). Surface electrodes were deposited on the AlN samples, and the vibration response of the samples to an applied voltage was monitored using the LDV as a function of the excitation frequency. The $d_{33}$ estimation was based on the excitation voltage and the thickness-mode displacement extracted from the LDV measurements. Further, six AlN substrates were irradiated with 8 MeV Al 2+ at three fluences (10 15 , 10 16 , and 10 17 ions/cm 2 ) and two temperatures (300 °C and 500 °C). The $d_{33}$ for the six irradiated samples and one pristine sample were measured, and the measurement uncertainty was estimated based on five repeated tests. All samples were also measured by a commercial piezometer for comparison. The experimental results demonstrate that the piezoelectric coefficients obtained by LDV were about 0.8–1.16 pm/V lower than those obtained by the piezometer. With the compensation of the clamping effect, the corrected LDV values are similar to the piezometer results. Both show similar trends in all samples, which validates the feasibility of the proposed method for $d_{33}$ measurement. Based on the LDV results, the irradiated samples show a 12%–22% decrease in $d_{33}$ compared with the pristine samples. The samples irradiated under the same fluence at a higher temperature (500 °C) demonstrated a lower $d_{33}$ than those at 300 °C. The effect of the retro-reflective tape, sample temperature, and sample size on the $d_{33}$ measurement were also studied.

36 MATERIALS SCIENCE↗
NASA NTRS1985

On the nonlinear characteristics of the axisymmetric flow regime: Cylindrical and spherical systems

The physical relationship between steady axisymmetric flows that might be observed in the atmosphere and in laboratory vessels is investigated theoretically. This is accomplished by comparing both the nonlinear structure and the thermal forcing mechanisms in two truncated spectral models of flow in the atmosphere and the rotating laboratory cylinder, respectively. Under statically stable conditions, the response of the internally forced spherical model (which is developed here from a set of new orthonormal basis functions) exhibits steady behavior different from that in the externally forced cylindrical model. Two regions of multiple steady solutions occur in the cylindrical model, under stable conditions, that are not found in the spherical one. The possible physical relevance of these multiple solutions is investigated by determining their location in parameter space with respect to the classical Hadley-Rossby transition curve. The results suggest that the wave flow regime, in an annulus, might develop catastrophically when an upper symmetric flow ceases to exist.

Higgins, R. W.↗
DOE OSTI2025

Electrically‐Driven Metal‐Insulator Transitions Emerging from Localizing Current Density and Temperature

Negative differential resistance (NDR) is a key electronic response enabling two‐terminal artificial neurons that can be achieved through different physical phenomena, including phase‐homogeneous current density and temperature (electro‐thermal) localizations and spatially‐localized metal‐insulator phase transitions (MITs). These two effects have been observed to occur sequentially in select electrically‐biased transition metal oxides. However, it is unknown why and under what conditions localizing behaviors precede MITs, particularly as a function of device length scale. To this end, the interplay between phase‐homogeneous electro‐thermal localizations and MITs is investigated in a 3D multiphysics simulation of a lateral thin film device, using the material properties of the prototype MIT material VO 2 . These findings demonstrate that the MIT is nucleated through dynamically localizing current density and temperature. A critical device width (≈0.7 µm in this study) is identified, below which both the electrically‐induced electro‐thermal and phase inhomogeneities cease to appear. It is demonstrated that the formation of spatial inhomogeneities directly relates to device dimensions, and demonstrate the decoupling of NDR from the MIT through device scaling relationships. These results provide insight into the material phenomena underlying the material's electrical responses, clarifying conditions under which spatial inhomogeneities form in electrically‐biased MIT materials.

artificial neuron↗
DOE OSTI2020

Frequency-dependent sensitivity of AC elastocaloric effect measurements explored through analytical and numerical models

We present a comprehensive study of the frequency-dependent sensitivity for measurements of the AC elastocaloric effect by applying both exactly soluble models and numerical methods to the oscillating heat flow problem. These models reproduce the finer details of the thermal transfer functions observed in experiments, considering here representative data for single-crystal Ba(Fe1−xCox)2As2. Based on our results, we propose a set of practical guidelines for experimentalists using this technique. This work establishes a baseline against which the frequency response of the AC elastocaloric technique can be compared and provides intuitive explanations of the detailed structure observed in experiments.

Straquadine, J. A. W. (ORCID:0000000165309541)↗
NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗
NASA NTRS2021

Ambient Performance Testing of the CubeSat Infrared Atmospheric Sounder (CIRAS)

Hyperspectral infrared measurements of Earth’s atmosphere from space have proven their value for weather forecasting, climate science and atmospheric composition. The CubeSat Infrared Atmospheric Sounder (CIRAS) instrument will demonstrate a fully functional infrared temperature, water vapor and carbon monoxide sounder in a CubeSat sized volume for at least an order of magnitude lower cost than legacy systems. Design for a CubeSat significantly reduces cost of access to space and enables flight in a constellation to reduce revisit time and enable new measurements including 3D winds. A technology demonstration of CIRAS is currently under development at JPL. The effort has completed integration and ambient testing of a high fidelity brassboard, complete with the flight configured optics assembly developed by Ball Aerospace with a JPL Immersion Grating and Black Silicon Entrance Slit. The brassboard includes a flight-configured High Operating Temperature Barrier Infrared Detector (HOT-BIRD) mounted in an Integrated Dewar Cryocooler Assembly (IDCA), enabling testing in the ambient environment. Ambient testing included radiometric testing of the system to characterize the instrument operability and NEdT. Spatial testing was performed to characterize the system line spread function (LSF) in two axes and report FWHM of the LSF. Spectral testing involved an air path test to characterize the spectral/spatial transformation matrix, and an etalon was used to measure the Spectral Response Functions (SRFs). Results of the testing show the CIRAS performs exceptionally well and meets the key performance required of the system. The end result of testing is the CIRAS instrument now meets TRL 4 with confidence in a brassboard configuration ready for thermal vacuum (TVac) testing necessary to achieve TRL 5 for the system.

Wilson, Robert C.↗
NASA NTRS2024

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles With Interrupted Relaxation Periods of TIMETAL 21S at 650 oC

The predictive performance of a multimechanism GVIPS model with saturating hardening function is assessed against a set of high temperature cyclic experiments under both strain-control and stress-control at varying rates of loading as well as several complex interrupted cyclic responses. The material specimens were composed of the titanium alloy, TIMETAL 21S, and tested at 650oC, see Lissenden et al., 2007. The model did a very good job of predicting the variety of tests, particularly given the fact that only monotonic loading tests (tensile, creep, relaxation) and a single fully reversed cycle was used for characterization. It was particularly interesting that the model was capable of reasonably capturing the rate of accumulated strain during ratchetting under tensile mean stress. Given complex interrupted cyclic response, the model was able to both qualitatively and even quantitatively predict reasonably well both the cyclic and associated relaxation periods (in all quadrants of the stress-strain space) thus confirming the validity of the chosen functional forms for both hardening and thermal recovery. The quantitative inaccuracy is associated with the fact that the material specimens tested by Lissenden et al., 2007 was significantly “softer” than those tested by Castelli in Saleeb et al., 1994 which were used for characterization of the GVIPS model parameters.

Viscoplaticity↗
NASA NTRS2026

GVIPS Prediction of Cyclic Ratchetting Behavior and Cycles with Interrupted Relaxation Periods of TIMETAL 21S at 650 °C

The predictive performance of a multimechanism GVIPS model with saturating hardening function is assessed against a set of high temperature cyclic experiments under both strain-control and stress-control at varying rates of loading as well as several complex interrupted cyclic responses. The material specimens were composed of the titanium alloy, TIMETAL 21S, and tested at 650 °C, see Lissenden et al., 2007. The model did a very good job of predicting the variety of tests, particularly given the fact that only monotonic loading tests (tensile, creep, relaxation) and a single fully reversed cycle was used for characterization. It was particularly interesting that the model was capable of reasonably capturing the rate of accumulated strain during ratchetting under tensile mean stress. Given complex interrupted cyclic response, the model was able to both qualitatively and even quantitatively predict reasonably well both the cyclic and associated relaxation periods (in all quadrants of the stress-strain space) thus confirming the validity of the chosen functional forms for both hardening and thermal recovery. The quantitative inaccuracy is associated with the fact that the material specimens tested by Lissenden et al., 2007 was significantly “softer” than those tested by Castelli in Saleeb et al., 1994 which were used for characterization of the GVIPS model parameters.

Viscoplaticity↗
DOE OSTI2022

High Performance Computing to Quantify the Evolution of Microscopic Concentration Gradients During Flash Processing

During the Flash process, the cross section of a plain-carbon or a low-alloy steel is austenitized through rapid heating and transformed on rapid cooling to a predominantly martensite + bainite structure with small amounts of retained austenite. Unlike conventional heat treating, homogeneity is intentionally avoided during Flash processing of steels. The Flash process assembly consists of a pair of rolls that transfer the steel sheets through the heating and cooling stage of the thermal cycle. The initial microstructure of the steel consists of ferrite (body-centered cubic iron) + carbide ((Fe,X)mCn) mixture. The heating rate through the peak temperature is a function of temperature and reaches a peak of about 300-400°C/s and the cooling rate has a maximum value of 3,000-4,000°C/s. The on-heating phase transformations include carbide dissolution, austenite (face-centered cubic iron) nucleation and growth, and diffusion of carbon and other substitutional elements in the steel. The on-cooling phase transformations include formation of martensite (body-centered tetragonal phase containing supersaturated solute) and bainite (ferrite plates with or without fine carbides). In this project, the focus is on Fe-C-Cr steels that are currently Flash processed for armor applications. The modeling effort proposed here will help optimize the Flash thermal cycle for these low alloy steels to achieve the target performance, which is an ongoing effort at SFP Works. A significant feature of Flash processed Fe-C-Cr steels is the presence of scatter in the through-thickness in the sheet. The variability in hardness results from a variability in the bainite + martensite microstructure that is sensitive to the local chemical concentration of C and Cr. Such a chemical inhomogeneity is intentionally obtained in the Flash process. Although such a microstructural gradient is presumably responsible for the exceptional properties of the Flash processed steel, it is very important to quantify the gradients as a function of Flash variabilities in processing parameters and the input microstructure. Understanding the mechanistic pathway that leads to microstructural gradients could be ground-breaking and instrumental for achieving better process control and optimized microstructural state to meet application-specific strength-ductility requirements. Since the final microstructure depends on setting up precise solute concentration gradients through a rapid heating process, and transforming these regions into various phases, it is important to understand how small changes in steel chemistry, input microstructure (carbide size and distribution), and process variables (Flash thermal cycle) will impact the solute concentration gradients.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Efficient Reformulation of HOTFGM: Heat Conduction with Variable Thermal Conductivity

Functionally graded materials (FGMs) have become one of the major research topics in the mechanics of materials community during the past fifteen years. FGMs are heterogeneous materials, characterized by spatially variable microstructure, and thus spatially variable macroscopic properties, introduced to enhance material or structural performance. The spatially variable material properties make FGMs challenging to analyze. The review of the various techniques employed to analyze the thermodynamical response of FGMs reveals two distinct and fundamentally different computational strategies, called uncoupled macromechanical and coupled micromechanical approaches by some investigators. The uncoupled macromechanical approaches ignore the effect of microstructural gradation by employing specific spatial variations of material properties, which are either assumed or obtained by local homogenization, thereby resulting in erroneous results under certain circumstances. In contrast, the coupled approaches explicitly account for the micro-macrostructural interaction, albeit at a significantly higher computational cost. The higher-order theory for functionally graded materials (HOTFGM) developed by Aboudi et al. is representative of the coupled approach. However, despite its demonstrated utility in applications where micro-macrostructural coupling effects are important, the theory's full potential is yet to be realized because the original formulation of HOTFGM is computationally intensive. This, in turn, limits the size of problems that can be solved due to the large number of equations required to mimic realistic material microstructures. Therefore, a basis for an efficient reformulation of HOTFGM, referred to as user-friendly formulation, is developed herein, and subsequently employed in the construction of the efficient reformulation using the local/global conductivity matrix approach. In order to extend HOTFGM's range of applicability, spatially variable thermal conductivity capability at the local level is incorporated into the efficient reformulation. Analytical solutions to validate both the user-friendly and efficient reformulations am also developed. Volume discretization sensitivity and validation studies, as well as a practical application of the developed efficient reformulation are subsequently carried out. The presented results illustrate the accuracy and implementability of both the user-friendly formulation and the efficient reformulation of HOTFGM.

Zhong, Yi↗
DOE OSTI2021

Studying drought-induced forest mortality using high spatiotemporal resolution evapotranspiration data from thermal satellite imaging

Drought can have pervasive and wide-spread impacts to forest health, as evidenced in several severe events occurring over the recent decades. Extensive forest die-off due to drought can impair the ecological functioning of forests, impacting habitat, water yield and quality from forested lands, and altering forest fire dynamics and intensity. Satellite remote sensing provides an effective means for detecting and monitoring spatial patterns of forest mortality over large areas, exploiting free and open long-term image archives available at a range in spatial and temporal resolutions. While remotely sensed surface reflectances and vegetation indices have been widely used to study optical response of forest canopies to drought events, retrievals of evapotranspiration (ET) derived from thermal satellite imagery – particularly at resolutions approaching crown scale - can provide insights into cumulative tree stresses that can incite disease and trigger mortality. In this study, we applied a multi-sensor satellite data fusion approach to estimate daily 30-m resolution ET and an associated Evaporative Stress Index (ESI) to study drought-induced mortality in a temperate forest at the Missouri Ozark AmeriFlux (MOFLUX) site, located in the central United States. The study covered the period from 2010 to 2014, including an exceptional drought year of 2012. Modeled ET agreed well with eddy flux measurements from the MOFLUX tower, with average monthly relative errors of 15%. Plot-scale ESI, describing temporal anomalies in the ratio of actual-to-reference ET, was used as an index of relative forest health to investigate relationships between forest mortality and drought severity. ESI showed good agreement with observed predawn leaf water potential, especially during the drought year. Furthermore, plot-scale ESI was also correlated with the subsequent year's tree mortality, suggesting the importance of considering the forest health condition prior to drought when studying drought-induced forest impacts. This study demonstrates the utility of multi-year ET remote sensing data at the stand or plot scale as an indicator of forest health and as a predictor of future mortality due to drought.

54 ENVIRONMENTAL SCIENCES↗