Search NASA⌕ Search

SEARCH · Search NASA

Results for “hard technologies”

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 595 records · Page 33

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗

CLARION2-TRINITY: A Compton-suppressed HPGe and GAGG:Ce-Si-Si array for absolute cross-section measurements with heavy ions

The design and performance of a new Compton-suppressed HPGe and charged-particle array, CLARION2-TRINITY, are described. The TRINITY charged-particle array is comprised of 64 Cerium-doped Gadolinium Aluminum Gallium Garnet (GAGG:Ce) crystals configured into five rings spanning 7–54 degrees, and two annular silicon detectors that can shadow or extend the angular coverage to backward angles with minimal $\gamma$-ray attenuation. GAGG:Ce is a non-hygroscopic, bright, and relatively fast scintillator with a light distribution well matched to SiPMs. Count rates up to 40 kHz per crystal are sustainable. Fundamental characteristics of GAGG:Ce are measured and presented, including light- and heavy-ion particle identification (PID) capability, pulse-height defects, radiation hardness, and emission spectra. The CLARION2 array consists of up to 16 Compton-suppressed HPGe Clover detectors (efficiency at 1 MeV) configured into four rings (eight HPGe crystal rings) using a non-Archimedean geometry that suppresses back-to-back coincident 511-keV gamma rays. The entire array is instrumented with 100- and 500-MHz (14 bit) waveform digitizers which enable triggerless operation, pulse-shape discrimination, fast timing, and pileup correction. Lastly, two examples of experimental data taken during the commissioning of the CLARION2-TRINITY system are given: a PID spectrum from 16 O + 18 O fusion-evaporation, and PID and Doppler-corrected -ray spectra from 48 Ti + 12 C Coulomb excitation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Coexistence of ordered and disordered vacancies in tungsten diboride with a broken hexagonal symmetry

As an important member of transition-metal borides, tungsten diboride (i.e., WB 2+𝑥 ) contains complex atomic vacancies with the originally assigned hexagonal symmetries, exhibiting many fascinating properties such as superconductivity and superior hardness. However, due to the difficulties in exploring atomic vacancies of transition-metal borides, the actual structure and composition of WB 2+𝑥 have been long-standing unsettled issues, which have impeded in-depth understanding of its structural stability and origins of such extraordinary properties. Here, we report a systematic investigation of the crystal structure of high-pressure synthesized WB 2+𝑥 samples by combination of state-of-the-art diffraction techniques and microstructural observations, leading to the discovery of an unusual coexistence of both the ordered and disordered atomic vacancies in WB 2+𝑥 with a broken hexagonal symmetry that can be well described by a dual-phase model involving hP12-W 0.70 ⁢B 1.73 and P1-W 0.75⁢ B 1.64 . Superposition of thin sample layers with disordered vacancies along the [001] or [110] direction is revealed to produce ordered vacancies, rationalizing the observed size-dependent symmetry breaking. In conclusion, these findings not only provide solid foundations for studying the phase stability and properties of this boride but also give powerful insights into how the intricate atomic vacancies can influence crystal structures of transition-metal borides.

high pressure↗

EXPERIMENTAL INVESTIGATION OF IMPACT AND ABRASION EROSIONS AT LOW VELOCITIES AND COMPARISON TO CFD SIMULATIONS

Wear and erosion damage of materials due to interaction of the surface with solid particles is a major issue in various industries. This phenomenon is common in mining, and oil and gas production, causing damages worth millions of dollars annually. Erosion is becoming a critical issue in renewable technologies as well such as concentrated solar thermal power (CSP) systems. In some CSP systems, solid particles can be used to absorb solar energy and enhance the heat storage. However, these particles may cause significant amount of wear to system components while falling down the power tower at high temperatures. This damage can be costly, and therefore, investigating solid particle erosion in CSP systems is of utmost importance. In collaboration with Sandia National Laboratories, the erosion durability of particle-based Generation 3 (GEN3) CSP system are being investigated both experimentally and numerically. In order to develop models to predict the life time of such systems, erosion under low temperature and low particle impact velocity conditions are investigated. Impact erosion and abrasion erosion experiments are conducted for low velocity and low temperature conditions for various candidate particles and containment materials. The impact erosion testing is conducted for three different containment material including Stainless Steel 316, special refractory material, and Inconel 740H with CARBOBEAD HSP 40/70 ceramic particles and UNIWEST 430 silica quartz erodent particles. All this data will enable us to better model the erosion behavior of these materials and provides a data base to investigate erosion at high temperatures experienced in CSP systems. At low temperatures and low velocities the softer refractory material exhibited greater erosion by an order-of-magnitude when compared to the harder Inconel 740H and Stainless Steel 316 materials. The overall erosion ratio for all the candidate materials was of the order of 10 -7 kg/kg. We also investigated wear in containment materials from sliding of particles over the surfaces at very low angles, also referred to as abrasion. For the velocities expected in a typical 1MW CSP plant, visible surface wear was noted, but the overall mass loss per unit surface area was found to be dependent on the hardness of the containment material tested. In this manuscript, we also present Computational Fluid Dynamics (CFD) simulations and erosion prediction of low velocity impact erosion. These results are compared with the experimental data to evaluate the ability of available erosion models for predicting erosion under these conditions.

14 SOLAR ENERGY↗

Tensile behavior of dual-phase titanium alloys under high-intensity proton beam exposure: radiation-induced omega phase transformation in Ti-6Al-4V

A high-intensity proton beam exposure with 181 MeV energy has been conducted at Brookhaven Linac Isotope Producer facility on various material specimens for accelerator targetry applications, including titanium alloys as a beam window material. The radiation damage level of the analyzed capsule was 0.25 dpa at the beam center region with an irradiation temperature around 120 oC. Tensile tests showed increased hardness and a large decrease in ductility for the dual a+ß-phase Ti-6Al-4V Grade-5 and Grade-23 extra low interstitial alloys, with the near a-phase Ti-3Al-2.5V Grade-9 alloy still exhibiting uniform elongation of a few % after irradiation. Transmission Electron Microscope analyses on Ti-6Al-4V indicated clear evidence of a high-density of defect clusters with size less than 2 nm in each a-phase grain. The ß-phase grains did not contain any visible defects such as loops or black dots, while the diffraction patterns clearly indicated ?-phase precipitation in an advanced formation stage. The radiation-induced ?-phase transformation in the ß-phase could lead to greater loss of ductility in Ti-6Al-4V alloys in comparison with Ti-3Al-2.5V alloy with less ß-phase.google scho

Tensile Test, titanium, acclerator, proton beam↗

Estimates of the wavenumber wavelet power spectrum of magnetic fluctuations during magnetic reconnection

Fluctuation analyses of experimental observations generally lack high temporal resolution and are in frequency-space f, contrary to theoretical efforts in wavenumber-space k. This is due to the inherent limits of the Fourier transform, though it is prominent due to the ease of diagnostic implementation. Advances in wavelet-based analysis have provided relief due to its temporal resolution, but in its common use, is still hard to compare to theoretical models. By using the two-point correlation technique in conjunction with large data sets, a wavelet power spectrum in wavenumber-space can be created. Dubbed the wavenumber wavelet power spectrum, this spectrum relates wavenumber to power in time. Further, this analysis technique more closely connects characterizations of experimentally observed fluctuations with other system parameters and theoretical predictions. In this article, we develop the wavenumber wavelet power spectrum using magnetic fluctuations caused by tearing instability driven magnetic reconnection in reproducible, high temperature laboratory plasmas. These dynamic magnetic fluctuations generated in reversed field pinch plasmas are broadband, ranging from the low frequency, 10's of kHz, up to the ion gyroradii frequencies, 100's of kHz. The dominant fluctuations have poloidal and toroidal mode numbers (m,n)=(1,6−10) and can grow to 2%–3% of the mean magnetic field. During these reconnection events, ions, and electrons are energized, magnetic fluctuation amplitudes increase, plasma flow is halted, and the toroidal magnetic flux increases, all on a semi-periodic basis. The newly developed spectrum provides better temporal resolution of spectrum characteristics to correlate with these particle energization phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗

Transient Photocurrent From High-Voltage Vertical GaN Diodes Irradiated With Electrons: Experiments and Simulations

Radiation-hard high-voltage vertical GaN p-n diodes are being developed for use in power electronics subjected to ionizing radiation. Here, we present a comparison of the measured and simulated photocurrent response of diodes exposed to ionizing irradiation with 70 keV and 20 MeV electrons at dose rates in the range of 1.4x10 7 - 5.0x10 8 rad(GaN)/s. The simulations correctly predict the trend in measured steady-state photocurrent and agree with experimental results within a factor of 2. Furthermore, simulations of the transient photocurrent response to dose rates with uniform and non-uniform ionization depth profiles uncover the physical processes involved that cannot be otherwise experimentally observed due to orders of magnitude larger RC time constant of the test circuit. The simulations were performed using an Exploratory Physics Development code developed at Sandia National Laboratories. The code offers the capability to include defect physics under more general conditions, not included in commercially available software packages, extending the applicability of the simulations to different types of radiation environments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Low Level Calibration of Liquid Scintillation Counting

Liquid Scintillation Counting (LSC) is an analytical laboratory method that quantifies the concentrations of alpha and beta emitting isotopes. The LSC method can be utilized for the detection of tritium, both environmentally and at nuclear facilities. Most LSC samples are aqueous-based and prepared as a radioactive sample dissolved or suspended in a high efficiency scintillation cocktail that converts radioactive energy into light pulses. When placed in counting vials, the liquid scintillation cocktail allows for optical coupling in which energy is released from the samples radioactive decay, and transferred to the scintillator. That energy is then absorbed to produce excited electrons that fall back to their ground-state and emit a pulse of light that is detected by the liquid scintillation analyzers photomultiplier tubes. Quench is a signal attenuation due to physical, chemical, or color interferences. To compensate for quench, quench standards of a known concentration are implemented at increasing increments of quench per vial to create a quench curve. Quench standards can be purchased pre-made or customized from the manufacturer, however, there are several issues that arise from purchasing the pre-made quench standards. The main issues with ordering pre-made quench standards is that they use higher count values, between 30,000 and 250,000 disintegrations per minute (dpm) of tritium per vial and only come in glass vials. Using higher counting standards in glass vials makes it hard to see lower concentrations in plastic vials. Additionally, when using higher energy standards the accuracy of the calibration readings tend to be lower. On the other hand, making a set of quench standards in-house allows for the customization of the dpm concentration to roughly 1,500-3,000 and preparation in plastic vials. This allows for lower activity readings that results in closer precision of the samples' true concentration. Technical Objectives: Create a quench curve and calibrate the LSA for tritium. Verify the calibration of the quench curve prior to cross sample runs. Check calibration against Quality Control (QC) reference value. Conclusion: a tritium quench curve calibration was performed using a scintillation cocktail with ultima gold (AB and AB/F), acetone (only with AB), and water (only with AB/F). By performing the tritium quench curve calibration, the counting efficiency was correlated to various levels of quench. Calibration checks with quality control standards were also performed to verify that the calibration was correct. The observed results provided acceptable values that matched the quality control standard reference values with an acceptable uncertainty. The next step will be fine tuning and cross comparison checks with legacy instrumentation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Adhesion of Self-Complementary, Sinusoidal Surfaces Fabricated Using Two-Photon Polymerization

Microscale, pick-and-place assembly is a non-lithographic assembly method poised to impact diverse fields including flexible electronics, microfluidics and robotics. However, a major technological challenge is the need to deterministically control adhesion between parts. Here, switchable adhesion involving 3D-printed, self-complementary surfaces is demonstrated. Mechanical properties of metasurfaces pressed against flat, rigid substrates are modeled using finite element methods. A series of flat slabs and metastructured slabs with 2D sinusoidal surfaces are printed using two-photon polymerization (2PP) of a shape-memory resin. The surface frequency of featured slabs was varied between $3.\bar3$ mm −1 and $26.\bar6$ mm −1 with similar amplitudes. Adhesion between printed metasurfaces and glass and between printed, self-complementary metasurfaces is studied above and below the cured resin’s glass transition temperature (∼45 °C). Simple heating of adhering surfaces to above 60 °C lowers adhesion, and compression of surfaces while above the glass transition temperature followed by cooling to room temperature elevates adhesion. The nominal adhesive strength between printed, self-complementary surfaces, as determined by the maximum observable pull-off stress, exceeds 3 MPa. Further tailoring complementary surfaces for adhesion control may facilitate microscale disassembly for recovery of components or precious metals.

amorphous materials↗

Salt-Induced Polymorphs Observed in Colloidal Single Crystals

Polymorphs are solid materials with the same chemical composition but different crystallographic structures. A unique aspect of polymorphs is that they exhibit different physical properties, such as solubility, melting point, density, color, hardness, and bioavailability. Here, we synthesized polymorphs of colloidal crystals engineered with DNA by slow-cooling gold nanoparticle-core programmable atom equivalents (PAEs, particles with DNA sequences that control their bonding characteristics) under salt concentrations ranging from 0.5 to 4 M NaCl. This approach yielded a diverse set of single-crystalline phases with cubic, tetragonal, and hexagonal lattice symmetries. The structural transitions observed here arise solely from the modulation of interparticle repulsion via ionic strength and thermal processing. Notably, in certain cases, we observed diffusionless phase transformations, wherein the superlattices evolve from cubic to lower-symmetry tetragonal lattices. By tuning the thermal stability and salt concentration, we captured intermediate, metastable body-centered tetragonal structures during the slow-cool process, indicating that subtle changes in free energy can direct crystallization to low-symmetry phases. In conclusion, this study demonstrates that thermal and ionic parameters can be tuned to access and stabilize colloidal crystal polymorphs with emergent structures and interesting functional properties.

Colloidal crystallization↗

Proton-Radiation Tolerant All-Perovskite Multijunction Solar Cells

Radiation-resistant but cost-efficient, flexible, and ultralight solar sheets with high specific power (W g -1 ) are the "holy grail" of the new space revolution, powering private space exploration, low-cost missions, and future habitats on Moon and Mars. Herein, this study investigates an all-perovskite tandem photovoltaic (PV) technology that uses an ultrathin active layer (1.56 micrometers) but offers high power conversion efficiency, and discusses its potential for high-specific-power applications. This study demonstrates that all-perovskite tandems possess a high tolerance to the harsh radiation environment in space. The tests under 68 MeV proton irradiation show negligible degradation (<6%) at a dose of 10 13 p + cm -2 where even commercially available radiation-hardened space PV degrade >22%. Using high spatial resolution photoluminescence (PL) microscopy, it is revealed that defect clusters in GaAs are responsible for the degradation of current space-PV. By contrast, negligible reduction in PL of the individual perovskite subcells even after the highest dose studied is observed. Studying the intensity-dependent PL of bare low-gap and high-gap perovskite absorbers, it is shown that the V OC , fill factor, and efficiency potentials remain identically high after irradiation. Radiation damage of all-perovskite tandems thus has a fundamentally different origin to traditional space PV.

14 SOLAR ENERGY↗

First Measurement of Z Opacity Sample Evolution near Solar Interior Conditions Using Time-Resolved Spectroscopy

Opacity model differences with Fe opacity measurements at high temperature (T>1⁢8⁢0 eV ) and high electron density (𝑛 𝑒 >3×10 22 cm −3 ) at nearly solar interior conditions have remained unresolved [Bailey et al., Nature 517, 56 (2015) and Nagayama et al., Phys. Rev. Lett. 122, 235001 (2019)]. Systematic errors from temporal gradients are one hypothesis for the discrepancy. Past data recorded on x-ray film provided spectral measurements over a time determined by the 2.8-ns backlighter duration. Here, we present the first measurements of opacity sample temporal evolution using novel hCMOS ultrafast x-ray camera technology. The measured conditions, measured backlighter time history, and modeled opacities are used to show that temporal gradients do not resolve the model-data discrepancy. The methods demonstrated provide potential advantages, including opacities at more extreme conditions, spectral line shift measurements, and improved measurements at other facilities.

atomic spectra↗

Molten Chloride Thermophysical Properties, Chemical Optimization, and Purification

Next-generation concentrating solar power (CSP) technology requires high-temperature heat-transfer fluids and thermal energy storage that can operate in the temperature range of 500°–720°C. This high temperature demands a new chemistry for the heat-transfer fluids that must have high-temperature thermal and chemical stabilities, high energy density, and low corrosion rate on metallic materials used in the heat exchangers, piping, and thermal energy storage tanks. A ternary MgCl 2 -KCl-NaCl chloride-salt system has been proposed to achieve such goals, given its high-temperature thermal stability. The chloride-salt-based systems have been used by the metallurgy industry as the electrolyte to electrolytically extract metallic magnesium. They have also attracted attentions from the nuclear industry for developing the molten salt fast reactors that can operate at close to atmospheric pressure and at higher temperature than traditional light-water reactors. However, the corrosion properties of the ternary salt are not well known. The corrosion mechanism is believed to be correlated to the presence of MgOHCl impurity at the operation temperature. MgOHCl is fundamentally caused by the hygroscopic nature of the MgCl 2 component. Therefore, the first key objective of this project is to understand the nature of the corrosiveness of the ternary chloride-salt system and its prevention by designing new lab-scale purification processes. Because of its importance, a working group called the Chloride Salt Collective across multiple U.S. Department of Energy national laboratories is established to gain experimental agreement on the corrosion behaviors. In addition, water generation during salt dehydration could lead to highly corrosive HCl/water mixture via a hydrolysis process between water and anhydrous MgCl 2 or MgCl 2 hydrates. MgOHCl will eventually be converted to MgO particulates regardless of the chemical path chosen (i.e., by thermal decomposition, chemical purification with an active metal, or electrochemical purification utilizing the reduction power of electron flows). Large MgO particles with high hardness can wear out bearing and bushing and even clog the delicate channels in the primary heat exchanger. Because of this, a second objective is to understand the chemistry behind salt purification and transfer the scientific understanding to a set of engineering know-how for large-scale salt melter design. Because there is limited information on the thermophysical properties and the sensitivity of these properties around the targeted salt composition, a third objective is to measure thermophysical properties relevant to system and component design and support Topic 1 Liquid Pathway development. For example, composition shift due to salt reaction with water can be a potential issue because it can lead to changes of thermophysical properties, such as (1) melting point (or freeze point), which can jeopardize safe operation of the supercritical CO 2 power cycle proposed for next-generation CSP, and (2) heat capacity and density, which can decrease the overall power-generation capability of the plants. The value of the thermal conductivity of the selected salt composition is also a critical criterion for the solar receiver design.

14 SOLAR ENERGY↗

Effect of Heat Treatment on Microstructure and Mechanical Property of 316L Stainless Steel Produced by Laser Powder Bed Fusion

The advanced non-light water reactor designs (Gen IV reactors), including molten salt/ very high temperature/ sodium-cooled and lead-cooled fast reactors, typically operate at higher temperatures and more extreme radiation conditions than light water reactors. An intrinsic part of the deployment and progress of Gen IV reactor designs is selecting the most suitable structural material for a specific application. Additive manufacturing (AM), a fairly new process of making physical, three-dimensional objects from a computer design file, is going to completely change the way of design, build and certify nuclear systems. It offers a range of opportunities to produce complex geometries from existing materials, offers new routes for processing of previously difficult to process materials, allows for design of new high-performance materials, and finally facilitates hybridization of dissimilar materials. This emerging technology has successfully produced cars, wind turbine blade molds and even live cells. It could also open up big opportunities for the nuclear industry to quickly deploy technologies at a fraction of the cost. So far, AM techniques have been preliminarily applied in the field of nuclear reactors, including the classical parts such as the pressure vessel of a small reactor with 508-III steel, the bottom nozzle of a fuel assembly with 304L steel, the fuel cladding with zirconium alloy and the integrated impeller of a pump and the multi-channel valve body with 316L steel [6,7]. The AM applications for operating nuclear reactors started in auxiliary plant components and have slowly migrated to metallic reactors and core components, but many of these are not safety critical components. Although many parts used for nuclear reactors have been fabricated by AM techniques, practical applications in engineering are still a long way off due to the uncertainty factors focused on the processing, material properties, analysis methods and application standards, which feeds the safety and life-cycle of the nuclear reactor. Due to rapid, repeated heating and cooling during production, a high dislocation density was present in the AM material. This microstructure feature is unstable at elevated temperature while high temperature is one of the typical operation environments for nuclear reactors. Thus, it is important to understand the thermal effect on the microstructure of AM material. The objectives of this study are to investigate the effect of heat treatment on the microstructure and mechanical properties of 316L stainless steel produced by laser powder bed fusion additive manufacturing, and to determine an appropriate heat treatment practice that will be applied to the lightweight AM lattice-structured material with the same chemistry. The heat treatment study consisted of annealing the samples at a temperature range of 800 to 1200 oC with a 50 oC increment for different times (1-24 hours), followed by vacuum or air cooling. Microstructural characterization was carried out by Scanning Electron Microscope (SEM). Grain size and crystallographic orientation were investigated by Electron Backscatter Diffraction (EBSD). Vickers hardness tests with a 0.5 kg load were employed to determine the hardness of samples after different heat treatments. After heat treatment, the random crystallographic orientation was preserved, and the volume fraction of high-angle grain boundaries (grain boundary misorientation =15 oC) remained the same. The dislocation density decreased with annealing temperature due to recovery. The fine subgrain structures in the as-printed specimen were quite stable up to 1200 oC. Minimal recrystallization was observed up to 1200 oC. Recrystallization initiated only after 8.5 hours at 1200 oC. The SEM images did not show obvious dependence of microstructure on cooling rate. The hardness of the specimens decreased with increasing annealing temperature as a result of the decrease in dislocation density. It is interesting to note that the AM material showed very similar hardness to the wrought material when annealing at similar temperature, although the microstructures are very different. Annealing at 1050 oC for 1 hour followed by air cooling was selected as the heat treatment procedure for the lattice designed lightweight AM 316L material.

36 MATERIALS SCIENCE↗

Parallel quantum annealing

Quantum annealers of D-Wave Systems, Inc., offer an efficient way to compute high quality solutions of NP-hard problems. This is done by mapping a problem onto the physical qubits of the quantum chip, from which a solution is obtained after quantum annealing. However, since the connectivity of the physical qubits on the chip is limited, a minor embedding of the problem structure onto the chip is required. In this process, and especially for smaller problems, many qubits will stay unused. We propose a novel method, called parallel quantum annealing, to make better use of available qubits, wherein either the same or several independent problems are solved in the same annealing cycle of a quantum annealer, assuming enough physical qubits are available to embed more than one problem. Although the individual solution quality may be slightly decreased when solving several problems in parallel (as opposed to solving each problem separately), we demonstrate that our method may give dramatic speed-ups in terms of the Time-To-Solution (TTS) metric for solving instances of the Maximum Clique problem when compared to solving each problem sequentially on the quantum annealer. Additionally, we show that solving a single Maximum Clique problem using parallel quantum annealing reduces the TTS significantly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

Observation of second harmonic electron cyclotron resonance heating and current-drive transition during non-inductive plasma start-up experiment in QUEST

Here, noninductive plasma current start-up using 2nd harmonic electron cyclotron resonance heating (ECRH) with oblique radio frequency (RF) injection is demonstrated in a Q-shu University experiment with steady-state spherical tokamak. A strong transition was observed in the heating and plasma current ramp-up. The initial bulk electron heating regime exhibits T ebulk ~ 140 eV and no hard x-ray (HXR) emission with a low I p of ~15 kA; it abruptly transitions to a regime that exhibits a low T ebulk of ~10 eV and a strong HXR emission with a high I p of ~50 kA. This behavior is distinctly different from that observed in previous fundamental ECRH experiments. The mechanism of the heating and current drive transition are investigated considering wave power absorption and plasma power balance. The results indicate that the transition is caused by the favorable heating of tail electrons where the RF power absorption at the 2nd harmonic increases nearly linearly with T etail , while the power transfer from the tail electrons to the bulk electrons decreases with 1/T etail 0.5 . This causes a rapid transition to a state with high T etail while reducing T ebulk towards colder ion temperature. The understanding of the transition mechanism helps to consider plasma current start-up using 2nd harmonic ECRH for tokamak reactors such as JT-60 SA and ITER.

QUEST↗