SEARCH · Engineering Papers

Results for “Weight function”

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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Reducing Mass of Steel Auto Bodies using Thin Advanced High Strength Steel with Carbon-Fiber Reinforced Epoxy

Diversitak, a company based in Detroit, MI, has developed a proprietary, low specific gravity, carbon fiber-reinforced epoxy (CFRE) under U.S. patent number 9,963,58832. Preliminary testing on this new material conducted in collaboration with ArcelorMittal Steel Company proved out the CFRE concept. A thin layer of this CFRE was applied to a stamped sheet of steel with residual stamping oils from a mill, in a time corresponding to automotive processing (e.g., ~15 seconds), and processed following automotive e-coat procedures (phosphating + 175–200°C heating), to complete the curing. No problems with adherence or performance were noted. While the CFRE does add weight to a thin gauge steel panel, it weighs much less than what is displaced by using thicker conventional mild steel gauges. The application of the coating showed a significant increased dent resistance, oil canning resistance, and part stiffness.This current two-year project was designed to mature this new technology to near manufacturing readiness to reduce the weight of a vehicle and lower the cost of weight reduction. The process involves the use of thinner gauge steels than are currently used. The collaborative development team included two industrial manufacturers: Diversitak and ArcelorMittal Steel Company; LightMAT; and two National Laboratories: Oak Ridge National Laboratory (ORNL) and Idaho National Laboratory (INL). The team developed a new manufacturing process to reduce the weight of a vehicle and lower the cost of weight reduction, as well as a better understanding of how to apply the coating so that it will perform to a high standard in-service. The team also performed an in-depth study to determine the long-term durability of the materials manufactured using this technology and well-known automotive industry standard tests.The overall process involved stiffening the thinner gauge steel by applying the CFRE on only one side. To accomplish this goal, the optimal reinforcement fiber length and fiber concentration was first determined. This was followed by measuring the coefficient of thermal expansion (CTE) in all three directions, so it could be fed into manufacturing models and methods for rapidly and inexpensively applying the coating. This was followed by panel level evaluations of the coating and steel combination, and then by full part demonstration of the technology on door panels. The final step was corrosion testing of the parts.ArcelorMittal characterized the advanced high strength steel (AHSS) (e.g., metallurgy-heat treatment for required AHSS properties as a function of the sheet thickness, state of internal stress) and quantified CFRE adhesion to the steel as a function of sheet preparation (e.g., rolling and stamping).ORNL optimized the fiber length, fiber concentration, and coating thickness for best vehicle function and performance at the least cost. Along with the suppliers, ORNL developed a durable CFRE application process (e.g., gun material, design, robotic dispensing process) and identified the adhesion stability of the CFRE during process holding. An approach to ensure that application/curing timing conforms to conventional assembly line speed and plant cycle times was determined. ORNL also determined the CTE of the material in all three directions and performed material scanning electron microscopy (SEM) analyses.INL characterized the corrosion properties of the steel panels coated with CFRE. The panels were investigated for corrosion resistance and stability as replacement materials used in automotive body panels to reduce mass. The coupons tested at INL were supplied by Diversitak after an optimized CFRE formulation was achieved in the already coated form for corrosion testing.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Analyzing and Exploring Training Recipes for Large-Scale Transformer-Based Weather Prediction

Abstract The rapid rise of deep learning (DL) in numerical weather prediction (NWP) has led to a proliferation of models which forecast atmospheric variables with comparable or superior skill than traditional physics-based NWP. However, among these leading DL models, there is a wide variance in both the training settings and architecture used. Further, the lack of thorough ablation studies makes it hard to discern which components are most critical to success. In this work, we show that it is possible to attain high forecast skill even with relatively off-the-shelf architectures, simple training procedures, and moderate compute budgets. Specifically, we train a minimally modified Swin Transformer V2 (SwinV2) on ERA5 data and find that it attains superior skill in terms of mean-square errors of deterministic forecasts when compared against the European Centre for Medium-Range Weather Forecasts’ Integrated Forecasting System (IFS). Almost all DL–NWP systems share a core set of hyperparameters and design decisions. To aid and expedite future DL–NWP research, we present an in-depth, systematic exploration of different loss functions, model sizes and depths, patch sizes, and multistep training objectives. We also examine the model performance with metrics beyond the typical accuracy (ACC) and RMSE and investigate how the performance scales with model size. Through our open-source code, scoring pipelines, and models, we share our findings on key aspects of the training pipeline. These ablations reduce the necessity for expensive hyperparameter tuning and lower the barrier to entry for future DL–NWP research. Significance Statement This study investigates the potential of using large-scale transformer-based models for weather prediction, showing that it is possible to achieve high forecast accuracy with simpler, off-the-shelf architectures. By training a minimally modified SwinV2 transformer on ERA5 data, we show that the model achieves competitive forecast skill in terms of mean-square error for key variables, outperforming the European Centre for Medium-Range Weather Forecasts’ Integrated Forecasting System (IFS) at all lead times. Our findings suggest that effective training strategies, such as multistep fine-tuning and channel-weighted losses, significantly enhance the model’s performance. However, we also highlight that these improvements come with trade-offs in other areas, such as ensemble spread and high-frequency spatial detail. This work highlights the promise of deep learning in improving weather forecasts, which could lead to better preparedness and response to weather events, ultimately benefiting society by providing more reliable weather predictions.

Willard, Jared D. [Lawrence Berkeley National Labo↗

DOE OSTI2026

Measurement of the Neutron Elastic Electromagnetic Form Factor Ratio at Large Momentum Transfer

Exploring nucleon structure is vital both for understanding its origin and existence as well as for the advancement of the sciences. It helps us answer key questions such as how quark and gluon dynamics create 99% of the nucleon mass. Electron- nucleon scattering has been widely used for precision studies of the nucleon and nuclear structure since the Nobel Prize winning investigations by Robert Hofstadter and collaborators in the 1950s. These studies provide information about the spatial charge and current densities of the nucleon in terms of the electromagnetic form factors. The form factors are functions of four momentum transfer squared (Q2). Extending the electromagnetic form factor measurements to higher Q2 plays a critical role in furthering the understanding of nucleon structure. This motivated the Super BigBite Spectrometer (SBS) program at Jefferson Lab. The open nature of the spectrometers and the direct line of sight from the target to the tracking detector locations in experimental setups such as SBS creates high levels of background at the detectors. This necessitates the use of tracking detectors with high rate capability and good position resolution. Gas Electron Multiplier (GEM) detectors are an excellent choice for tracking detectors in such experiments. Understanding the performance of the GEM detectors is important not just for SBS experiments but also for future high-luminosity experiments. This thesis reports the exploratory results from the measurement of the neutron elastic electromagnetic form factor ratio (Gn E/Gn M) at high momentum transfer. A longitudinally polarized electron beam was scattered off a polarized 3He target, used as an effective polarized neutron target. In this experiment, the polarized 3He target achieved a world record polarization weighted luminosity at a beam current of 45 µA. Double spin asymmetry of the scattered neutron events is used to extract the neutron form factor ratio. Measurements were taken at Q2 = 3.0, 6.8, 9.8 (GeV/c)2. The lowest Q2 measurement is in good agreement with the existing world data, and the higher-Q2 measurements extend the Q2 reach well beyond the existing world data and are expected to remain unmatched for a long time.

Gamage, Vimukthi Haththotuwa [Univ. of Virginia, C↗

DOE OSTI2022

Additive manufacturing of LiNi 1/3 Mn 1/3 Co 1/3 O 2 battery electrode material via vat photopolymerization precursor approach

Additive manufacturing, also called 3D printing, has the potential to enable the development of flexible, wearable and customizable batteries of any shape, maximizing energy storage while also reducing dead-weight and volume. In this work, for the first time, three-dimensional complex electrode structures of high-energy density LiNi 1/3 Mn 1/3 Co 1/3 O 2 (NMC 111) material are developed by means of a vat photopolymerization (VPP) process combined with an innovative precursor approach. This innovative approach involves the solubilization of metal precursor salts into a UV-photopolymerizable resin, so that detrimental light scattering and increased viscosity are minimized, followed by the in-situ synthesis of NMC 111 during thermal post-processing of the printed item. The absence of solid particles within the initial resin allows the production of smaller printed features that are crucial for 3D battery design. The formulation of the UV-photopolymerizable composite resin and 3D printing of complex structures, followed by an optimization of the thermal post-processing yielding NMC 111 is thoroughly described in this study. Based on these results, this work addresses one of the key aspects for 3D printed batteries via a precursor approach: the need for a compromise between electrochemical and mechanical performance in order to obtain fully functional 3D printed electrodes. In addition, it discusses the gaps that limit the multi-material 3D printing of batteries via the VPP process.

25 ENERGY STORAGE↗

DOE OSTI2022

Star formation and AGN feedback in the local Universe: Combining LOFAR and MaNGA

The effect of active galactic nuclei (AGN) on their host galaxies – in particular their levels of star formation – remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. Here, in this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey (LoTSS) with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), to probe the star-forming properties of 307 local (z < 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad star-formation rate (SFR) distributions, typically lie below the SFMS, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies’ star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Identification and characterization of a small-molecule metallophore involved in lanthanide metabolism

Many bacteria secrete metallophores, low-molecular-weight organic compounds that bind ions with high selectivity and affinity, in order to access essential metals from the environment. Previous work has elucidated the structures and biosynthetic machinery of metallophores specific for iron, zinc, nickel, molybdenum, and copper. No physiologically relevant lanthanide-binding metallophore has been discovered despite the knowledge that lanthanide metals (Ln) have been revealed to be essential cofactors for certain alcohol dehydrogenases across a diverse range of phyla. Here, we report the biosynthetic machinery, the structure, and the physiological relevance of a lanthanophore, methylolanthanin. The structure of methylolanthanin exhibits a unique 4-hydroxybenzoate moiety which has not previously been described in other metallophores. We find that production of methylolanthanin is required for normal levels of Ln accumulation in the methylotrophic bacterium Methylobacterium extorquens AM1, while overexpression of the molecule greatly increases bioaccumulation and adsorption. Our results provide a clearer understanding of how Ln-utilizing bacteria sense, scavenge, and store Ln; essential processes in the environment where Ln are poorly bioavailable. More broadly, the identification of this lanthanophore opens doors for study of how biosynthetic gene clusters are repurposed for additional functions and the complex relationship between metal homeostasis and fitness.

Science & Technology - Other Topics↗

DOE OSTI2024

Role of the likelihood for elastic scattering uncertainty quantification

In the last decade, uncertainty quantification (UQ) for optical model potentials (OMPs) has become a focal point for nuclear reaction theory, and several competing approaches for OMP UQ have recently been developed. Here, we clarify recent efforts to compare frequentist and Bayesian approaches in the context of OMP UQ [G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019)]. We replicate a portion of that OMP UQ study but use independent statistical tools. Specifically, we compare two methods for OMP parameter inference from elastic scattering data: the Levenberg-Marquardt algorithm for χ 2 minimization on one hand and Markov chain Monte Carlo (MCMC) sampling on the other. Separately, we assess the common practice of using a renormalized likelihood (χ 2 /N), N being the number of data points, instead of the canonical weighted-least-squares likelihood (χ 2 ), as a way of accounting for unknown data correlations. Here, we show that for a generic linear model and for a five-parameter OMP analysis, frequentist and uniform-prior Bayesian approaches recover the same optimum and uncertainty estimates—not systematically larger uncertainties for the Bayesian approach, as was concluded in G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019). Further, we show that if an additional, near-degenerate parameter is introduced into the same OMP analysis such that the parameter posterior becomes non-Gaussian, then covariance-based estimates of uncertainty become unreliable. Finally, we show that regardless of optimization approach, if χ 2 /N is used for the likelihood, the resulting parametric uncertainties increase by $\sqrt{N}$, and that this is responsible for the conclusions drawn in the revisited study. Based on our replication results, we find that a fortuitous cancellation of unreported errors and the renormalization factor can lead to improvement in empirical coverages, as was the case in the original comparative study. We emphasize that developing and applying a realistic likelihood function is an essential task in a UQ analysis, and that several recent UQ studies that employed a renormalized likelihood (i.e., including a 1/N factor) may have yielded unrealistically large uncertainties for elastic-scattering observables. If the parameter posterior deviates from multivariate-normal, a sampling-based approach like MCMC has a clear advantage over methods that assume the Laplace approximation holds. We note that empirical coverage can serve as an important internal check for the analyst whose model or data may have additional, unaccounted-for uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Electronic coherence evolution at the nearly commensurate-incommensurate CDW boundary of 1⁢T−TaS 2

Transition-metal dichalcogenides host a variety of charge-density-wave phases that couple lattice, charge, and correlation effects. In 1⁢𝑇−TaS 2 , the commensurate and nearly commensurate states are well characterized, yet the transition near 350 K into the incommensurate phase has lacked direct momentum-resolved insight. In this work, we use temperature-dependent angle-resolved photoemission spectroscopy to track the electronic structure across this transition. We observe a suppression of quasiparticle spectral weight at the Brillouin-zone center, coincident with the transport anomaly, but without clear evidence of a full band-gap opening. The transition appears to involve momentum-dependent redistribution of spectral weight, consistent with a loss of coherence that reshapes the Fermi surface while leaving conduction dispersions largely intact. These results suggest that the nearly commensurate–incommensurate transition may not align with a conventional metal-insulator transition picture, but rather as an electronic reconstruction driven by loss of coherence. Our work provides new microscopic insight into the resistivity anomaly near room temperature and may guide design principles for collective electronic switching in transition-metal dichalcogenides.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗

DOE OSTI2023

High-Strength, High-Ductility, High Entropy Alloys with High-Efficiency Native Oxide Solar Absorbers for Concentrating Solar Power Systems

This EPSCoR Project has been investigating the synergy between the excellent high-temperature mechanical behavior of FeMnNiAlCr high entropy alloys (HEA) and the high solar absorptance of their native oxides for high efficiency concentrated solar thermal power (CSP) systems working at >700°C. While HEAs have attracted substantial interest in recent years, most investigations have focused on their applications as structural materials rather than functional materials. This EPSCoR project discovered that FeMnNiAlCr HEAs can potentially be applied synergistically as both a structural and functional material for high-efficiency concentrating solar thermal power (CSP) systems working at >700°C. The HEA itself would be used in high-temperature tubing to carry molten salts or supercritical CO 2 , while its surface oxide would act as a high-efficiency solar thermal absorber. With Fe and Mn being the major components in these HEAs (adding up to ~70 at.% of the alloy), these materials are much more cost-effective than the Ni-based superalloys currently being investigated for high-temperature CSP systems. Through this research, these Fe-Mn based HEAs have demonstrated yield strengths 2-3x greater than that of stainless steel at 700°C and a creep lifetime >800 h at 700ºC under a typical CSP tubing mechanical load of 35 MPa. Their Mn-rich surface oxides maintain a high optical-to-thermal conversion efficiency of ~87% under 1000x solar concentration ratio for 20 simulated day-night thermal cycles between 750ºC and room temperature. In preliminary corrosion studies, these HEAs have sustained immersion in unpurified bromide molten salts for 14 days at 750°C with <2% weight loss, in contrast to 70% weight loss from a 316 stainless steel reference. The simultaneous achievement of promising mechanical, optical, and thermochemical properties in this FeMnNiAlCr system opens the door to new applications of HEAs in solar energy harvesting. Partnerships with Ames Laboratory and Oak Ridge National Laboratory (ORNL) also advanced our understanding of the fundamental structure-property relationships through atomic scale material characterization and first-principles computational modeling. The key research results in this project can potentially be extended to other HEAs and their native oxides. In terms of applications, the proposed FeMnNiAlCr HEA/native oxide system could potentially exceed the mechanical and the optical performance of existing tubing and solar coating materials under EERE’s CSP program at lower cost, which also aligns well with the EPSCoR Science and Technology strategies of New Hampshire in boosting the deployment of renewable energy.

14 SOLAR ENERGY↗

DOE OSTI2024

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗

DOE OSTI2025

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Screening a knowledge‐based library of low molecular weight compounds against the proline biosynthetic enzyme 1‐pyrroline‐5‐carboxylate 1 ( PYCR1)

Abstract Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 1 (PYCR1) is the last enzyme of proline biosynthesis and catalyzes the NAD(P)H‐dependent reduction of Δ 1 ‐pyrroline‐5‐carboxylate toL‐proline. High PYCR1 gene expression is observed in many cancers and linked to poor patient outcomes and tumor aggressiveness. The knockdown of thePYCR1gene or the inhibition of PYCR1 enzyme has been shown to inhibit tumorigenesis in cancer cells and animal models of cancer, motivating inhibitor discovery. We screened a library of 71 low molecular weight compounds (average MW of 131 Da) against PYCR1 using an enzyme activity assay. Hit compounds were validated with X‐ray crystallography and kinetic assays to determine affinity parameters. The library was counter‐screened against human Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 3 and proline dehydrogenase (PRODH) to assess specificity/promiscuity. Twelve PYCR1 and one PRODH inhibitor crystal structures were determined. Three compounds inhibit PYCR1 with competitive inhibition parameter of 100 μM or lower. Among these, (S)‐tetrahydro‐2H‐pyran‐2‐carboxylic acid (70 μM) has higher affinity than the current best tool compoundN‐formyl‐l‐proline, is 30 times more specific for PYCR1 over human Δ 1 ‐pyrroline‐5‐carboxylate reductase isoform 3, and negligibly inhibits PRODH. Structure‐affinity relationships suggest that hydrogen bonding of the heteroatom of this compound is important for binding to PYCR1. The structures of PYCR1 and PRODH complexed with 1‐hydroxyethane‐1‐sulfonate demonstrate that the sulfonate group is a suitable replacement for the carboxylate anchor. This result suggests that the exploration of carboxylic acid isosteres may be a promising strategy for discovering new classes of PYCR1 and PRODH inhibitors. The structure of PYCR1 complexed withl‐pipecolate and NADH supports the hypothesis that PYCR1 has an alternative function in lysine metabolism.

Biochemistry & Molecular Biology↗

DOE OSTI2022

Multiobjective optimization of nuclear microreactor reactivity control system operation with swarm and evolutionary algorithms

To improve the marketability of novel microreactor designs, there is a need for automated and optimal control of these reactors. This paper presents a methodology for performing multiobjective optimization of control drum operation for a microreactor under normal and off-nominal conditions. Here, two different case studies are used where the control drum configuration is optimized for the reactor to be critical with some desired power distribution that would satisfy peaking limits. A surrogate model for power distribution is developed based on a feedforward neural network. The process for determining weights for scalarization of the multiobjective optimization problem is also detailed. Six optimization algorithms: evolutionary strategies, differential evolution, grey wolf optimization, Harris hawks optimization, moth flame optimization and particle swarm optimization, are all applied to these cases and the results analyzed. Although all these algorithms will demonstrate optima-seeking behavior, for real-time control it is necessary to identify the best algorithm to efficiently provide reasonable optima without operator interference. The moth flame optimization algorithm was found to perform particularly well on both cases. Overall, it was found that the algorithms capable of supplying the best optima were also the most consistent. Finally, the found optima were verified with the original model used to train surrogates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Photochemical spin-state control of binding configuration for tailoring organic color center emission in carbon nanotubes

Abstract Incorporating fluorescent quantum defects in the sidewalls of semiconducting single-wall carbon nanotubes (SWCNTs) through chemical reaction is an emerging route to predictably modify nanotube electronic structures and develop advanced photonic functionality. Applications such as room-temperature single-photon emission and high-contrast bio-imaging have been advanced through aryl-functionalized SWCNTs, in which the binding configurations of the aryl group define the energies of the emitting states. However, the chemistry of binding with atomic precision at the single-bond level and tunable control over the binding configurations are yet to be achieved. Here, we explore recently reported photosynthetic protocol and find that it can control chemical binding configurations of quantum defects, which are often referred to as organic color centers, through the spin multiplicity of photoexcited intermediates. Specifically, photoexcited aromatics react with SWCNT sidewalls to undergo a singlet-state pathway in the presence of dissolved oxygen, leading to ortho binding configurations of the aryl group on the nanotube. In contrast, the oxygen-free photoreaction activates previously inaccessible para configurations through a triplet-state mechanism. These experimental results are corroborated by first principles simulations. Such spin-selective photochemistry diversifies SWCNT emission tunability by controlling the morphology of the emitting sites.

36 MATERIALS SCIENCE↗