Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal modeling techniques”

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

Towards inverse microstructure-centered materials design using generative phase-field modeling and deep variational autoencoders

The field of Integrated Computational Materials Engineering (ICME) combines a broad range of methods to study materials’ responses over a spectrum of length scales. A relatively unexplored aspect of microstructure-sensitive materials design is uncertainty propagation and quantification (UP/UQ) of materials’ microstructure, as well as establishing process-structure–property (PSP) relationships for inverse material design. In this study, an efficient UP technique built on the idea of changing probability measures and a deep generative unsupervised representative machine learning method for microstructure-based design of thermal conductivity of materials is proposed. Probability measures are used to represent microstructure space, and Wasserstein metrics are used to test the efficiency of the UP method. By using deep Variational AutoEncoder (VAE), we identify the correlations between the material/process parameters and the thermal conductivity of heterogeneous dual-phase microstructures. Through high-throughput screening, UP, and the deep-generative VAE method, PSP relationships that are too complex can be revealed by exploiting the materials’ design space with an emphasis on microstructures. As a last point, we demonstrate generative machine learning serves as a useful tool for inverse microstructure-centered materials design, and we demonstrate this by examining the inverse design of thermal conductivity in nano-structured materials. Here, the results reveal the effects of morphology, volume fraction, characteristic length scale, and the individual thermal diffusivity of phases on the thermal conductivity of dual-phase alloys. Our findings emphasize the advantages of high-throughput phase-field modeling and generative deep learning for linking PSP and inverse microstructure-centered materials design.

36 MATERIALS SCIENCE↗

Arctic Critical Infrastructure: Assessing and Predicting the Risk to Critical Permafrost Infrastructure from Climate Change: A New Thermomechanical Approach

This study presents the development of a computational framework designed to predict the interaction between permafrost and infrastructure, addressing potential failure modes and mitigation strategies in the context of climate change. The framework, rooted in advanced modeling and simulation (mod/sim) techniques, integrates thermomechanical coupling to account for the complex interplay between heat flow, ice content, and mechanical behavior in permafrost. Existing models fail to fully capture these dynamics, particularly as they relate to the effects of ice saturation on structural integrity. Our innovative Arctic Coastal Erosion (ACE) framework fills this gap by coupling thermal and mechanical models to accurately simulate subsidence and deformation in permafrost environments. We applied the ACE framework to a representative runway, demonstrating its capability to predict settlement due to rising temperatures and subsequent permafrost thaw. This proof-of-concept showcases the potential of the framework to evaluate risks to Arctic infrastructure, which supports over four million people and 70% of existing permafrost-based structures. By simulating various infrastructure types and environmental conditions, our research offers insights into failure mechanisms and evaluates structural solutions to mitigate risk. The anticipated deliverables, including a prototype runway exemplar, position this project as a critical advancement in permafrost infrastructure modeling, with applications in national security and resilience planning.

54 ENVIRONMENTAL SCIENCES↗

Subvoxel Control of Fiber Orientation via Multidirectional Shearing in 3D Printing

Anisotropy, the characteristic of materials exhibiting different properties based on their direction, is widespread in nature. Conventional manufacturing techniques often fall short in recreating such complex anisotropy. While 3D printing allows for precise fiber deposition in anisotropic composites, previous studies have only achieved bulk reorientation of fibers at the voxel level in two dimensions. This limits the replication of localized 3D anisotropies found in natural materials. To address this, a novel 3D printing technique is presented that enables subvoxel control of fiber orientation in all three directions using multidirectional shearing via nozzle rotation and inclination. The fiber orientation control in these experimental tests is driven by a numerical model, enabling a fully digital approach to program microstructures within a strand. This programmability is demonstrated through mechanical and thermal tests, illustrating a localized and controllable response to external stimuli. Achieving such complex anisotropy holds potential across several fields, including wearables and biomedical implants, lightweight composite structures, and energy storage technologies such as batteries and supercapacitors.

3D printing↗

Low‐Density 2D Superlattices Assembled via Directional DNA Bonding

Abstract It is critical to assemble nanoparticles (NPs) into superlattices with controlled symmetries and spacings on substrates for metamaterials applications, where such structural parameters dictate their properties. Here, we use DNA to assemble anisotropic NPs of three shapes—cubes, octahedra, and rhombic dodecahedra—on substrates and investigate their thermally induced reorganization into two‐dimensional (2D) crystalline films. We report two new low‐density 2D structures, including a honeycomb lattice based on octahedral NPs. The low‐density lattices favored here are not usually seen when particles are crystallized via other bottom‐up assembly techniques. Furthermore, we show that, consistent with the complementary contact model, a primary driving force for crystallization is the formation of directional, face‐to‐face DNA bonds between neighboring NPs and between NPs and the substrate. Our results can be used to deliberately prepare crystalline NP films with novel morphologies.

Miao, Ziyi↗

Thermophysical Properties of NaCl–UCl 3 –PuCl 3 Molten Salts: A Combined Computational and Experimental Study

Actinide-bearing molten salts for use as fuels are an essential part of next generation molten salt reactors. Yet, numerous multicomponent salt mixtures are underdeveloped or have not been investigated. Here, this study, based on a combination of experimental and modeling techniques, is dedicated to determining and understanding a variety of properties of the ternary system of NaCl–UCl 3 –PuCl 3 , which represents a scenario for burnup of NaCl–UCl 3 fuel, at two compositions (∼10 and 5 mol % PuCl 3 in eutectic NaCl–UCl 3 pseudobinary) and a range of temperatures. Evaluation of the heat flow and mass loss data showed the 0.61NaCl–0.30UCl 3 –0.09PuCl 3 salt had a melting temperature of 551 ± 5 °C. Two additional thermal effects were observed occurring at approximately 410 and 494 °C. The transition occurring at 410 °C may be due to the presence of oxide in the salt. Extrapolation of thermodynamic data indicates the transition occurring at 494 °C is due to the formation of a liquid phase. Experimental testing determined the density of this system is a linear function of temperature and can be represented by the equation ρ = 4.014–0.0010T(°C), R 2 = 0.992. Additionally, by using atomistic modeling, we found that increasing the PuCl 3 content from 5 to 10 mol % led to the formation of larger Pu 3+ clusters and slower transport of ions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantifying thermal barrier coating delamination through luminescence modeling

Thermal barrier coatings (TBCs) are widely used in the hot sections of gas turbine systems as they are remarkably efficient in insulating the underlying superalloys, leading to higher operating temperatures and therefore enhancing the performance of the engines. However, this benefit is only possible if the integrity of the TBC under aggressive thermo-mechanical environments is ensured. Delamination is a common but hard to detect failure mode. We present, in this work, supported by experimentation, the implementation of a modeling approach applying the Kubelka-Munk theory to provide numerical quantification of luminescence contrast and intensity due to top coat delamination in TBCs. The method relies on the drastic change in reflectivity when a delamination forms, exploiting it for high-contrast luminescence mapping. Two distinct TBC configurations containing an erbium-doped yttria-stabilized zirconia (YSZ:Er 3+ ) layer for delamination sensing were used to validate this model. A delamination zone induced by Rockwell indentation was successfully tracked by measuring an increase of the intensity of the erbium emission line at 562 nm. Luminescence-based methods for delamination detection can provide a revolutionary non-invasive technique, with potential for both off-line and on-line engine monitoring.

36 MATERIALS SCIENCE↗

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗

Experimental and Computational Investigation of Integrated Internal and Film Cooling Designs Incorporating a Thermal Barrier Coating

Abstract Few studies in the open literature have studied the effect of thermal barrier coatings (TBC) when used in combination with shaped hole film cooling and enhanced internal cooling techniques. The current study presents Reynolds-averaged Navier–Stokes (RANS) conjugate heat transfer simulations that identify trends in cooling design performance as well as experimental measurements of overall effectiveness using a flat-plate matched-Biot number model with a simulated TBC layer of 0.42D thickness, where D is the film cooling hole diameter. Coolant is fed to the film cooling holes in a co-flow configuration, and the results of both the smooth and rib-turbulated channels are compared. At a constant coolant flow rate, enhanced internal cooling was found to provide a 44% increase in spatially-averaged overall effectiveness, $\bar{\bar{\phi}}$, without a TBC. The results show that the addition of a TBC can raise $\bar{\bar{\phi}}$ on a film-cooled component surface by 47%. The optimum velocity ratio was found to decrease with the addition of enhanced cooling techniques and a TBC as the film provided minimal benefit at the expense of reduced internal cooling. While the computational results closely identified trends in overall system performance without a TBC, the model over-predicted effectiveness on the metal–TBC interface. The results of this study will inform turbine component design as material science advances increase the reliability of the TBC.

Engineering↗

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗

Coordinated Thermal Safety Attack and Defense on EV Battery Management Systems

Battery temperature sensor and battery current sensor data which are key sensing inputs to the Battery Management Controllers in electric vehicles, are vulnerable to possible cyber/ physical manipulation due to known vulnerabilities inherited from CAN bus technology that is used for in-vehicle communications between electronic control units that transfer sensing and control data. In this paper, we first create a simulation that enables us to evaluate impact of cyber physical attacks on electric vehicle battery management system in a controlled environment that violates thermal safety. Specifically, we emulate a Level 3 - DC fast charging system with SAE J1772/CCS, integrated with standard charging controls and thermal safety controls on EVs, and various sensing data flows. Second, we propose a coordinated current and battery temperature attack that has crippling economic, and safety impacts. Third, we quantify the usability, economic and safety impacts of such attacks as a function of the extent of data manipulation. Finally, we propose a physics model driven detection technique to detect presence of such attacks.

25 ENERGY STORAGE↗

Thermal Conductivity Degradation in High Burnup U-Pu-Zr Fuel

Recent advancements in the characterization of irradiated U-Pu-Zr fuels have revealed complexities that challenge existing understanding of constituent redistribution. Traditionally, models have proposed three concentric regions within the fuel, each characterized by distinctive phases and porosity. However, through detailed analysis of high burnup U-Pu-Zr, we discovered the presence of four distinct constituent redistribution regions. Particularly novel is the observation of significant Pu redistribution, a previously unreported phenomenon that necessitates a reevaluation of current models. This work aims to delve deeper into these findings, seeking to correlate mesoscale measurements of thermal diffusivity and respective thermal conductivity with the phases present in each redistribution region. To achieve this objective, we employed mesoscale thermoreflectance methods using the unique, Idaho National Laboratory (INL) developed, Thermal Conductivity Microscope (TCM) at INL’s Irradiated Materials Characterization Laboratory. The TCM employs two tightly focused lasers: one for heating to generate periodic thermal waves in the substrate, and another spatially separated probe laser to detect changes in the optical reflectivity of the gold-coated substrate resulting from thermal wave diffusion. We conducted several thermal diffusivity measurements within each region of constituent redistribution of a U-19Pu-10Zr fuel pin cross section irradiated to 11 at. % burnup. The TCM measurement positions strategically aligned with transmission electron microscopy (TEM) lift-out locations previously collected from the fuel sample. Complementary microstructural analysis techniques such as optical and scanning electron microscopy (OM/SEM), electron probe microanalysis for chemical compositions, and TEM-based selective area electron diffraction (SAED) analysis for crystallographic insights into each phase were also utilized. This comprehensive approach allowed us to correlate local thermal diffusivity data with microstructural characteristics, enabling the computation of local thermal conductivity at each position. The significance of this contribution lies in its pioneering use of the TCM for ternary fuel mesoscale examination, shedding light on the previously overlooked effects of Pu redistribution on local thermal conductivity. By informing current models capturing constituent redistribution and heat transfer, our findings pave the way for more accurate predictions of metallic fuel performance. Moreover, this work sets the stage for future comparisons with similar TCM examinations on U-19Pu-10Zr fuels at ultra-low burnup, facilitating a comprehensive understanding of thermal property changes across different burnup levels. Ultimately, our study not only enriches our understanding of the thermophysical properties of individual redistribution regions within U-Pu-Zr fuel but also offers valuable insights for the design and operational parameters of proposed next-generation fast reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

New Materials and an Efficient Processing Approach for Materials for Harsh Environments – Continuous Electric Field Assisted Sintering

Electric-field assisted sintering (EFAS) is an advanced manufacturing method for the consolidation and processing of ceramic and metal materials. It is a particularly excellent candidate for efficient materials processing as it has greater than 90% energy cost savings and improved CO2 emissions compared to traditional sintering in a furnace. This is due to direct rapid heating of the materials and molds leveraging an electrification phenomenon called Joule heating where electrical current is passed through the sample/mold assemblies causing direct heating. The technique does have limitations though as it is traditionally a batch-process and suffers from scale up difficulties. To overcome this, researchers at Idaho National Laboratory changed the form factor of traditional EFAS instruments to develop and construct a continuous rolling electric-field assisted sintering (CEFAS) device. Its unique design creates highly localized electrical, thermal, and mechanical process parameters for continuous samples of limitless length. Advanced core shell powder materials were also created using fluidized powder bed atomic layer deposition of core/shell solid oxide electrochemical electrode/electrolyte BaCe0.7Zr0.1Y0.1Yb0.1O3 (BCZYYb) and yttrium doped BaZrO3 (BZY) to be a feedstock for the new CEFAS instrument. The innovative core/shell powder is hypothesized to have high protonic conductivity with enhanced environmental resistance. The completed project covers innovations in process control and design enabled by control theory, finite element modeling and mechanical design, three dimensionally printed carbon-carbon composite manufacturing and utilization, and green materials processing techniques.

36 MATERIALS SCIENCE↗

Estimating Soil Thermal Inertia Profiles From the Passive Equilibration of a Temperature Probe

Knowledge of the distribution of soil thermal properties is important for understanding subsurface hydrological and biogeochemical processes. This study describes and evaluates quick thermal profiling (QTP), a new measurement technique aimed at providing rapid, depth-resolved measurements of soil thermal inertia at numerous locations across the landscape. A cylindrical probe with temperature sensors at multiple depths is quickly inserted into the ground, and soil thermal inertia is estimated from how quickly the probe temperature equilibrates with the soil. To this end, a finite volume heat transfer model is used to generate temperature equilibration time series across combinations of controlling factors, and a gridded search inversion approach is applied to infer soil thermal inertia. Field tests in the Arctic indicate that QTP measurements have a minimum uncertainty of 0.14 J m −2 K −1 s −1/2 and covary with dual-probe heat pulse thermal analyzer measurements (concordance correlation coefficient = 0.56) with a root-mean-square error of 0.40 J m −2 K −1 s −1/2 . Besides demonstrating the value of QTP for estimating thermal inertia, this study identifies various sources of measurement uncertainty, particularly probe-soil contact resistance and frictional heating. Further, analysis of soil samples indicates that thermal inertia can be used to estimate thermal conductivity and dry bulk density in the studied area, although such inferences are highly site-specific. Overall, the QTP method holds promise to generate thermal inertia data products and to complement other characterization approaches for advancing understanding of soil properties across far more locations than is currently possible.

Lamb, J. R. [Lawrence Berkeley National Laboratory↗

Characterizing the nonequilibrium response of FeRh thin films using time-domain thermoreflectance

Time-domain thermoreflectance (TDTR) characterization of FeRh throughout its first-order antiferromagnetic (AF) to ferromagnetic (FM) transition shows that the transient reflectance ΔR⁡(t)/R strongly depends on the magnetic order of the sample. Using TDTR, which uses optical pulses to induce small temperature excursions, we have found that Δ⁢R⁡(t)/R of the AF phase exhibits a large negative response, while the response of the FM phase is positive. This magnetic phase sensitivity has allowed us to study the transient response of both the AF and FM phases to the pump-pulse excitation and the mixed phase of the material. These results are significant since the ultrafast properties of antiferromagnetic materials and mixed antiferromagnetic and ferromagnetic materials are difficult to detect using other conventional techniques. We have found that the AF phase exhibits a strong subpicosecond decaying signal not observed in the FM phase. The magnetic phase dependence of the sign of Δ⁢R⁡(t)/R is qualitatively explained using the results of ab initio density functional theory calculations. Using the two-temperature model, we found that the change in the thermalization time across the transition is caused by differences in both the electronic heat capacity and the electron-phonon coupling factor of the AF and FM phases. The electron-phonon coupling constant in the AF phase is also determined using the two-temperature model conducted using the ntmpy code package. For the FM phase, we provide boundaries for the magnitude of the electron-phonon coupling factor for the FM phase. Finally, these results indicate that TDTR can be used to study the transient properties of magnetic materials that are otherwise challenging to probe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Abstract Low-level jets (LLJs), in which the wind speed attains a local maximum at low altitudes, have been found to occur in the U.S. mid-Atlantic offshore, a region of active wind energy deployment as of 2023. In contrast to widely studied regions such as the U.S. southern Great Plains and the California coastline, the mechanisms that underlie LLJs in the U.S. mid-Atlantic are poorly understood. This work analyzes floating lidar data from buoys deployed in the New York Bight to understand the characteristics and causes of coastal LLJs in the region. Application of the Hilbert–Huang transform, a frequency analysis technique, to LLJ case studies reveals that mid-Atlantic jets frequently occur during times of adjustment in synoptic-scale motions, such as large-scale temperature and pressure gradients or frontal passages, and that they do not coincide with motions at the native inertial oscillation frequency. Subsequent analysis with theoretical models of inertial oscillation and thermal winds further reveals that these jets can form in the stationary geostrophic wind profile from horizontal temperature gradients alone—in contrast to canonical LLJs, which arise from low-level inertial motions. Here, inertial oscillation can further modulate the intensity and altitude of the wind speed maximum. Statistical evidence indicates that these oscillations arise from stable stratification and the associated frictional decoupling due to warmer air flowing over a cold sea surface during the springtime land–sea breeze. These results improve our conceptual understanding of mid-Atlantic jets and may be used to better predict low-level wind speed maxima. Significance Statement The purpose of this work is to identify and characterize the atmospheric mechanisms that result in an occasional low-level maximum in the wind speed off the U.S. mid-Atlantic coastline. Our findings show that these low-level jets form due to horizontal temperature gradients arising from fronts and synoptic systems, as well as from the land–sea breeze that forces warmer air over the cold ocean surface. This work aids predictability of such jets, improves our understanding of this coastal environment, and has implications for future deployment of offshore wind energy in this region.

17 WIND ENERGY↗

Helium Bubble Growth and Retention Characteristics in Aging Palladium Tritide

The helium-related characteristics of aged palladium tritide, including tritium loss, the formation of He bubbles from tritium decay-induced He atoms, and He thermal desorption behavior are computed using a continuum model of bubble evolution that incorporates the log-normal bubble spacing distribution deduced from 3 He nuclear magnetic resonance (NMR) measurements. This spacing distribution produces significant differences between mean (expectation) values of the bubble size and pressure and the average (median) values of these aging characteristics obtained by a previous calculation. The new calculations find much of the He is retained in a large quantity of smaller bubbles with higher bubble pressures, improving the overall stability of the bubbles and extending the age of He retention. By contrast, the integrated tritide material characteristics of swelling and tritium pressure are found relatively insensitive to details of the bubble spacing distribution. Inter-bubble fracture is predicted to begin with small, closely-spaced bubbles and progress to include larger bubbles near the onset of rapid He release. The critical concentration for this onset increases almost linearly with material tensile strength and decreases with shear modulus. Release of He from the aging solid is modeled as a growing network of fractured inter-bubble ligaments. For small particle or small grain material, the physical size of the fracture cluster, relative to the material dimensions, becomes important. In contrast with normal aging, inter-bubble fracture occurring during thermal desorption appears to begin with large, widely-spaced bubbles, indicating elevated temperature techniques may be of limited use for evaluating development of the bubble fracture network.

74 ATOMIC AND MOLECULAR PHYSICS↗

Confronting Dark Matter with the Multiwavelength Sky (Final Technical Report)

The existence of dark matter is a key clue to new fundamental physics, yet its nature and non-gravitational interactions remain mysterious. In the course of this five-year project (spanning July 2015 - July 2020), we have developed novel analysis techniques to separate dark matter signals from astrophysical backgrounds in the inner Galaxy, and measured hitherto-unknown properties of the Fermi Bubbles and the Galactic Center excess in gamma-rays from the inner Milky Way. We have improved predictions for the spatial distribution and energy spectra of photon signals from dark matter, and set new limits from astrophysical and cosmological datasets on dark matter inhabiting complex dark sectors, multi-particle dark matter annihilation, and tiny primordial black holes. We have successfully adapted effective field theory techniques from collider physics to make high-precision predictions for signals from heavy weakly-interacting DM, allowing us to essentially exclude the possibility that DM is a thermal wino, and have discovered novel theoretical and observational implications for dark matter bound states and self-interactions. We have developed a comprehensive public code package for modeling the effects of exotic energy injections, such as from dark matter annihilation and decay, on the ionization and temperature history of the early cosmos, and have used these tools to set robust, powerful and broadly-applicable constraints on models of light dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low-Density 2D Superlattices Assembled via Directional DNA Bonding

It is critical to assemble nanoparticles (NPs) into superlattices with controlled symmetries and spacings on substrates for metamaterials applications, where such structural parameters dictate their properties. Here, we use DNA to assemble anisotropic NPs of three shapes-cubes, octahedra, and rhombic dodecahedra-on substrates and investigate their thermally induced reorganization into two-dimensional (2D) crystalline films. We report two new low-density 2D structures, including a honeycomb lattice based on octahedral NPs. The low-density lattices favored here are not usually seen when particles are crystallized via other bottom-up assembly techniques. Furthermore, we show that, consistent with the complementary contact model, a primary driving force for crystallization is the formation of directional, face-to-face DNA bonds between neighboring NPs and between NPs and the substrate. Finally, our results can be used to deliberately prepare crystalline NP films with novel morphologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗