Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat map”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Boron coordination change in barium borate melts and glasses and its contribution to configurational heat capacity, entropy, and fragility

High-energy x-ray diffraction from molten and glassy BaB 2 O 4 and BaB 4 O 7 has been performed using aerodynamic levitation and laser heating over a wide range of temperatures. Remarkably, even in the presence of a heavy metal modifier dominating x-ray scattering, it was possible to extract accurate values for the tetrahedral, sp 3 , boron fraction, N 4 , which declines with increasing temperature, using bond valence-based mapping from the measured mean B–O bond lengths while accounting for vibrational thermal expansion. These are used within a boron-coordination-change model to extract enthalpies, ΔH, and entropies, ΔS, of isomerization between sp 2 and sp 3 boron. The results for BaB 4 O 7 , ΔH = 22(3) kJ mol -1 boron, ΔS = 19(2) J mol -1 boron K -1 , agree quantitatively with those found previously for Na 2 B 4 O 7 . Analytical expressions for N 4 (J, T) and associated configurational heat capacity, $C_{P}^{conf}$(J, T), and entropy, S conf (J, T), contributions are extended to cover a wide composition range 0 ≤ J = BaO/B 2 O 3 ≤ 3 using a model for ΔH(J) and ΔS(J) derived empirically for lithium borates. Maxima in the $C_{P}^{conf}$(J, T g ) and fragility index contributions are thereby predicted for J ≲ 1, higher than the maximum observed and predicted in N 4 (J, T g ) at J ≃ 0.6. We discuss the utility of the boron-coordination-change isomerization model in the context of borate liquids containing other modifiers and the prospect of neutron diffraction to aid in empirical determinations of modifier-dependent effects, illustrated by new neutron diffraction data on Ba 11 B 4 O 7 glass, its well-known α-polymorph, and lesser-known δ-phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Perovskite grain wrapping by converting interfaces and grain boundaries into robust and water-insoluble low-dimensional perovskites

Stabilizing perovskite solar cells requires consideration of all defective sites in the devices. Substantial efforts have been devoted to interfaces, while stabilization of grain boundaries received less attention. Here, we report on a molecule tributyl(methyl)phosphonium iodide (TPI), which can convert perovskite into a wide bandgap one-dimensional (1D) perovskite that is mechanically robust and water insoluble. Mixing TPI with perovskite precursor results in a wrapping of perovskite grains with both grain surfaces and grain boundaries converted into several nanometer-thick 1D perovskites during the grain formation process as observed by direct mapping. The grain wrapping passivates the grain boundaries, enhances their resistance to moisture, and reduces the iodine released during light soaking. The perovskite films with wrapped grains are more stable under heat and light. The best device with wrapped grains maintained 92.2% of its highest efficiency after light soaking under 1-sun illumination for 1900 hours at 55°C open-circuit condition.

14 SOLAR ENERGY↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

XtalCAMP: a comprehensive program for the analysis and visualization of scanning Laue X-ray micro-/nanodiffraction data

XtalCAMP is a software package based on the MATLAB platform, which is suitable for, but not limited to, the analysis and visualization of scanning Laue X-ray micro-/nanodiffraction data. The main objective of the software is to provide complementary functionalities to the Laue indexing software packages used at several synchrotron beamlines. Here, the graphical user interfaces allow the easy analysis of characteristic microstructure features, including real-time intensity mapping for a quick examination of phase, grain and defect distribution, 2D color-coded mapping of microstructural properties from the output of other Laue indexing software, crystal orientation visualization, grain boundary characterization based on orientation/misorientation calculation, principal strain/stress analysis, and strain ellipsoid representation, as well as a series of additional toolkits. As an example, XtalCAMP is applied to the microstructural investigation of a solution-heat-treated Ni-based superalloy manufactured using a laser 3D-printing technique, and a deformed natural quartzite from Val Bregaglia in the Central Alps.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial stabilization for inverted perovskite solar cells with long-term stability

Perovskite solar cells (PSCs) commonly exhibit significant performance degradation due to ion migration through the top charge transport layer and ultimately metal electrode corrosion. Here, we demonstrate an interfacial management strategy using a boron chloride subphthalocyanine (Cl 6 SubPc)/fullerene electron-transport layer, which not only passivates the interfacial defects in the perovskite, but also suppresses halide diffusion as evidenced by multiple techniques, including visual element mapping by electron energy loss spectroscopy. As a result, we obtain inverted PSCs with an efficiency of 22.0% (21.3% certified), shelf life of 7000 h, Τ 80 of 816 h under damp heat stress (compared to less than 20 h without Cl 6 SubPc), and initial performance retention of 98% after 2000 h at 80 °C in inert environment, 90% after 2034 h of illumination and maximum power point tracking in ambient for encapsulated devices and 95% after 1272 h outdoor testing ISOS-O-1. Our strategy and results pave a new way to move PSCs forward to their potential commercialization solidly.

14 SOLAR ENERGY↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-throughput investigation of structural evolution upon solid-state in Cu–Cr–Co combinatorial multilayer thin-film

Cu-Cr-Co combinatorial multilayer thin-films were prepared by a high-throughput ion beam sputtering system. Based on the thickness ratio among the individual nanoscale monolayers (Cu, Cr, Co), the resulting stoichiometry covered the entire phase diagram. The chemical composition and structure of Cu-Cr-Co combinatorial chip upon solid-state reaction were studied by lab-based micro-X-ray fluorescence (mu-XRF) and high-throughput synchrotron X-ray diffraction (XRD), respectively. A composition-structure map for Cu-Cr-Co combinatorial chip was developed through automated data analysis employing hierarchical clustering techniques. The structural evolution of Cu-Cr-Co combinatorial chip as a function of heat-treatment temperature, time, and modulation period were studied systematically. Furthermore, the effect of the elemental distribution in-depth direction was investigated to gain more insights regarding phase transformation. This work provides an efficient method and new perspectives for the design and optimization of the composition and structure of high-performance thin-films.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Boiling Transitions During Droplet Contact on Superheated Nano/Micro-Structured Surfaces

Manipulating surface topography is one of the most promising strategies for increasing the efficiency of numerous industrial processes involving droplet contact with superheated surfaces. In such scenarios, the droplets may immediately boil upon contact, splash and boil, or could levitate on their own vapor in the Leidenfrost state. In this work we report the outcomes of water droplets coming in gentle contact with designed nano/micro-textured surfaces at a wide range of temperatures as observed using high-speed optical and X-ray imaging. Further, we report a paradoxical increase in the Leidenfrost temperature ($T_{\text{LFP}}$) as the texture spacing is reduced below a critical value (~10 μm) that represents a minima in $T_{\text{LFP}}$. Although droplets on such textured solids appear to boil upon contact, our studies suggest that their behavior is dominated by hydrodynamic instabilities implying that the increase in $T_{\text{LFP}}$ may not necessarily lead to enhanced heat transfer. On such surfaces, the droplets display a new regime characterized by splashing accompanied by a vapor jet penetrating through the droplets before they transition to the Leidenfrost state. We provide a comprehensive map of boiling behavior of droplets over a wide range of texture spacings that may have significant implications towards applications such as electronics cooling, spray cooling, nuclear reactor safety and containment of fire calamities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Heat Wave Impacts on Western US Electricity System

Visualization of hourly generation, load, unserved load, and locational marginal energy prices in the western United States during a July 22-28 heatwave event in 2018 and 2058. In addition to hourly time series data of each of the parameters, choropleth maps showing the hourly value for each balancing authority are provided as well as a county-level choropleth map of temperature.

Mongird, Kendall↗
DOE OSTI2024

Neutron Diffraction Residual Stress Study of a AA2219-T87 Self-Reacting Friction Stir Weld

2219-T87 is a precipitation hardenable aluminum-copper alloy which sees wide use in structural aerospace components. Thick panels of this alloy are joined via self-reacting friction welding (SRFSW); however, this thermomechanical process can result in significant loss of the alloy’s strengthening precipitates and large residual stresses which are detrimental to mechanical behavior and performance. High-resolution maps of the residual strain and stress states along the normal, transverse, and longitudinal directions of 2219-T87 SRFSW were obtained using neutron diffraction measurements. Residual stress had the highest tensile value in the heat-affected zone (HAZ) and the largest compressive stress in the base metal region of the normal and transverse directions. Line profile residual stress distributions displayed an “M”-shaped distribution in each direction with the longitudinal direction being the most pronounced while 2D residual stresses displayed an hourglass-shaped pattern. Further, average maximal longitudinal residual stress values ranged from 60.6 to 85.6% of yield. Microhardness testing across the transverse weld section produced a “V”-shaped curve with slight hardness recovery in the stir zone. Maximal and minimal microhardness values were observed in the base metal (154 HV) and thermomechanically affected zone (82 HV), respectively. Microstructural evolution was recorded using optical microscopy and showed decreasing grain size from the HAZ to the stir zone.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploratory analysis of machine learning techniques in the Nevada geothermal play fairway analysis

Play fairway analysis (PFA) is commonly used to generate geothermal potential maps and guide exploration studies, with a particular focus on locating and characterizing blind geothermal systems. This study evaluates the application of machine learning techniques to PFA in the Great Basin region of Nevada. Following the evaluation of various techniques, we identified two approaches to PFA that produced promising results, 1) supervised Bayesian probabilistic neural networks to generate geothermal potential maps with confidence intervals, and 2) unsupervised principal component analysis paired with k-means clustering to generate both cluster maps to help identify spatial patterns, as well as new combined feature inputs. We applied these techniques to perform a comparative analysis between two principal sets of geological and geophysical features related to permeability and heat and a set of positive (known geothermal resources) and negative training sites (known drill sites with unsuitable geothermal conditions). We found that these methods constrain previously unrecognized feature controls on geothermal favorability, many of which are spatially organized within the extent of cluster groups and the major structural-hydrologic domains of the study area. Furthermore, we utilized exploratory unsupervised modeling to highlight spatial relationships between input data and predictive output results of our supervised modeling. As a result, we demonstrate how our models compare to the previous Nevada PFA and how the rapid insights these machine learning techniques offer may support future assessments of both known and undiscovered blind geothermal systems in the Great Basin region of Nevada and beyond.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Long-lived period-doubled edge modes of interacting and disorder-free Floquet spin chains

Floquet spin chains have been a venue for understanding topological states of matter that are qualitatively different from their static counterparts by, for example, hosting π edge modes that show stable period-doubled dynamics. However the stability of these edge modes to interactions has traditionally required the system to be many-body localized in order to suppress heating. In contrast, here we show that even in the absence of disorder, and in the presence of bulk heating, π edge modes are long lived. Their lifetime is extracted from exact diagonalization and is found to be non-perturbative in the interaction strength. A tunneling estimate for the lifetime is obtained by mapping the stroboscopic time-evolution to dynamics of a single particle in Krylov subspace. In this subspace, the π edge mode manifests as the quasi-stable edge mode of an inhomogeneous Su-Schrieffer-Heeger model whose dimerization vanishes in the bulk of the Krylov chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Multi-Source Machine Learning and Thermoplastics Enhanced Aerostructure Manufacturing (mTEAM)

RTX Technology Research Center (RTRC), together with Collins Aerospace (Collins) and Oak Ridge National Laboratory (ORNL) has developed an Artificial Intelligence (AI) / Machine Learning (ML) guided solution to advance the manufacturing and assembly of high performance and lightweight thermoplastic composite (TPC) aerospace products. The solution aims to lower risk, cost and lead time for induction heating based welding and consolidation processes for TPC structure. The cost and lead time of part and material specific process development for induction welding (IW) and induction consolidation will be reduced by replacing traditional empirical methods with optimization methods that merge AI/ML and physics-based process simulations and process experiments with sensing and controls. TPC-IW process development is empirical in nature, and uncertainties in material & process behavior exist near & far from the induction coil. Physics-based simulations can be leveraged directly for process optimization but can be too computationally expensive to run in high fidelity and real time to do robust process optimization. The key impact of successful TPC induction consolidation and welding is cost & lead time reduction for part & material specific consolidation and welding recipes. This is an enabler for more rapid deployment of TPC structures via joining assembly, which can reduce energy & cost intensive usage of autoclaves & ovens. The solution aimed to advance the U.S. Department of Energy’s interests in using thermoplastics and automation in composite manufacturing for improvement of products for existing markets via increased production speeds, reduced costs, and lowered use of energy. Welded TPC structures can offer significant weight & energy savings for high-value commercial aerospace & industrial applications compared to metal & thermoset composite structures assembled by mechanical fastening and/or adhesive bonding. The project was organized into two Budget Periods. Budget Period 1 (BP1) was 15 months and its goal was to perform ML process optimization framework development & deployment on lab-coupon aerostructure components. A Go/No-Go Review was performed at the end of BP1 to verify fulfilment of key tasks & milestones to justify a Go Decision to move into the next Budget Period. Budget Period 2 (BP2) was 12 months and its goal was the deployment of the ML framework for ML process optimization of pilot industrial scale aerostructure components. The overall project aim was to develop & demonstrate ML-enhanced modeling framework that learns process-property mapping from multiple data sources at different fidelities. During BP1, the team accomplished key tasks & milestones to demonstrate the concept of multi-source ML for TPC aerostructure consolidation and assembly. First, the team completed documentation of induction based TPC heating requirements including baseline metrics to compare measured results against. Next the team completed demonstration of data generation from physics-based simulations for ML surrogate model generation and demonstrated the integration of physics-based simulation data into multi-source AI/ML algorithms. In parallel, the team established the lab-coupon scale induction welding system and completed a process to label and reduce generated data from physics-based simulation and experiments for ML surrogate models to enable multi-source ML model training & testing. To complete BP1, the team integrated physics-based simulation data and experimental data into multi-source ML algorithms. This was based on the team completing ML deployment of the induction welding on a lab system at RTRC and AI/ML deployment on existing induction welding line at Collins. ORNL visited both Collins and RTRC sites to witness the TPC induction welding process. Then, ORNL designed and constructed a new version of their vision-based sensing system better adapted to acquire process signals of the TPC induction welding process for process anomaly and defect detection. In BP2, the team accomplished key tasks & milestones to scale up multi-source ML for TPC aerostructure consolidation and assembly from the lab-coupon scale to the pilot-industrial scale. In BP2, the team demonstrated real time anomaly & defect detection via experiments performed by ORNL & RTRC. The team completed ML-optimization heating trials for TPC induction consolidation at Collins, and the team confirmed pilot industrial scale experimental data from Collins was compatible with the developed ML pipeline from RTRC. The team completed sub-element scale ML process optimization demonstration at RTRC, where the team leveraged RTRC’s robotic TPC welding setup to de-risk the ML process optimization by performing ML analysis of recorded temperatures to account for complex part features. Then, the team applied its ML-derived control strategies and ML process optimization framework at Collins to the pilot-industrial scale on a demo skin-stiffener part representative of a nacelle aerostructure fan cowl section. The key innovation is the AI/ML framework enabling effective process development of high performance, lightweight, energy efficient TPCs for composite aircraft structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mapping of a Novel Zero-Liquid Discharge Desalination System Based on Humidification–Dehumidification onto the Field of Existing Desalination Technologies

It is well-established that increasing demands for fresh water are paving the way for desalination technologies. However, this correlates with an increase in brine production whose treatment and disposal can be complicated and expensive. This paper presents a thermodynamic model to bound the operation and development of a novel Humidification–Dehumidification-based system featuring Zero-Liquid Discharge and off-grid capabilities. The model employs conservation laws to find feasible state points to meet a baseline operation of 10 kg/h of product water separated from a hypersaline feed stream with 100 g/kg salt concentration. The system incurs in a 1039 kWh/m3 energy intensity that can be supplied completely by an electric source or in combination with heating steam. Follow-up sensitivity analysis highlights the robustness of the system in handling variations of 25% in product flowrate and 75% in feed salinity, practically without incurring any additional energy demands. The proposed system operating costs between 72 USD/m3 and 96 USD/m3 are comparable to those of existing brine disposal techniques. Furthermore, an operational map of existing desalination technologies suggests a niche characterized by high recovery rates and high feed salinities that are generally unfulfilled by conventional desalination methods. Overall, the proposed system shows potential for off-grid hypersaline brine treatment. This study sets the stage for future development of physics-based and data-driven predictive models as the proposed system iterates into a pilot plant deployment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Machine Learning for Predicting Multipactor Susceptibility in Planar RF Structures

Multipactor discharge is a persistent challenge in high-power microwave (HPM) and accelerator systems, where secondary electron avalanches can cause heating, vacuum degradation, and failure. This work presents the first supervised machine learning (ML) framework for multipactor prediction, trained on high-fidelity 3D Particle-in-Cell (PIC) simulation data in planar geometries. The model maps operational, geometric, and material-dependent secondary electron yield (SEY) parameters to the time-averaged electron growth rate, enabling rapid reconstruction of susceptibility charts. Among the models evaluated, tree-based ensemble methods such as Random Forest and Extra Trees demonstrate superior generalization to unseen materials compared to neural networks such as multilayer perceptron (MLP). Performance metrics, including Intersection over Union (IoU), Structural Similarity Index Measure (SSIM), and Pearson correlation, show close agreement with simulation benchmarks. Principal Component Analysis attributes generalization limits to material feature-space disjointedness.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Paleoclimatic implications of glacial fluctuations in the Sierra Nevada del Cocuy, northern Andes, Colombia, during the Lateglacial and Holocene

The reconstruction of former mountain glaciers from geomorphic mapping and cosmogenic-nuclide surface-exposure dating provides a unique opportunity to infer patterns of past terrestrial climate variability. Tropical mountain glaciers are particularly valuable as there are comparatively few terrestrial climate proxies at equatorial latitudes relative to higher latitudes. As the single largest climate zone on Earth, the tropics play an outsized role in mediating global climate via the ocean-atmosphere transfer of latent heat and water vapor. Nonetheless, there remains a persistent gap in our understanding of how the tropics influenced – or were influenced by – the high-magnitude climate shifts of the Late Pleistocene, and whether this high-energy region simply responded to extratropical forcing or was itself a driver of global climatic change. To help address this knowledge gap, we analyzed geologic evidence for past glacial fluctuations in three adjacent valleys in the Sierra Nevada del Cocuy, the highest subrange of the Eastern Cordillera in the Colombian Andes, to provide a terrestrial record of atmospheric temperature during the latter part of Termination 1. Coupled with geomorphic mapping and paleo-snowline reconstructions, our beryllium-10 glacial chronology indicates that glaciers in the humid inner tropics underwent pronounced growth and gradual decay during the Antarctic Cold Reversal (14.5–12.8 ka) and Younger Dryas (12.8–11.7 ka) periods, respectively, following a trend that, according to directly dated moraine records from throughout both polar hemispheres, appears to have been global. While the specific mechanism(s) behind this large-scale behavior remains to be corroborated, we revisit the hypothesis that ocean-atmosphere heat transfer and water vapor flux are key drivers of abrupt Lateglacial temperature fluctuations. Subsequent to the Lateglacial, deglaciation of the Sierra Nevada del Cocuy accelerated during the Early Holocene, a pattern also observed in other tropical glacier records. More recently, the magnitude of snowline rise and glacier retreat over the last two centuries supports the view that modern tropospheric warming is anomalously strong at least relative to the last ∼16,000 years.

Andes↗