SEARCH · Engineering Papers
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.
A culture is a way of life of a group of people - the behaviors, beliefs, values, and symbols that they accept, generally without thinking about them, and that are passed along by communication and imitation from one generation to the next. We have different cultures between Technical Areas, facilities, divisions, work groups, etc.
PDQ is a tool for the management of the input and execution of batch jobs for simulation codes that use a text based input system. It accomplishes this goal by operating at two levels. First, it takes input file templates (commonly known at LANL as input deck templates) and creates multiple instantiations by performing substitutions of data from table files into symbols (variables) found in the template. Second, it provides commands to submit the created files to the SLURM batch system for execution. These two activities taken together produce a whole that is greater than the sum of its parts and provides an elegant way of executing studies across multiple similar simulations while minimizing the risk of typographical errors in the input files. PDQ was originally developed as a job management system called XVS by Jeff McAninch while he was at LANL. Besides the capabilities described here, XVS had many other features specific for interactions with particular simulation codes. After Jeff’s departure, maintenance of XVS was taken over by Rendell Carver; he added some new features as well as kept it functioning as the batch system at LANL was changed from LSF to MOAB to SLURM. In 2017, Rob Pelak decided to develop a different version that removed the additional features (many of which were rendered obsolete with the retirement of the simulation code or batch system that they supported) and produced a cleaner “bare bones” version of XVS. A few other behaviors of XVS that Rob found irksome were altered. Rob gave the resulting code a new name: PDQ. In 2022 Danielle McDermott developed a version that runs under Python 3.X. As suggested by Rob, she used the python2to3 utility to identify most changes. Given that PDQ continues to operate with Python version 2.7 we have advanced the version number to 1.20.
This summer project started out to quantify the error associated with an approximate model of a superconducting circuit, where two qubits are dispersively coupled through a bus resonator. Understanding this error helps us decide if the approximate model is sufficient for the purpose of optimal control. We find out that the error of the approximate model covers up a nontrivial problem of the underlying circuit, which limits the fidelity of two-qubit gates under 99% for practical parameters. This suggests that high fidelity gates of this circuit cannot be achieved through optimal control alone, and changes of the parameters/design of the circuits appear necessary. Additionally, to carry out the algebra involved, a symbolic commutator algebra package HamMaker is developed in Julia. The package is capable of automating some difficult algebra in optimal control that previously had to be performed by hand. Within this work, HamMaker has been used to derived some simple but nontrivial corrections to the existing effective Hamiltonian.
Collections management staff from the National Security Research Center (NSRC) recently digitized more than 10,000 McKibbin Cards to make them accessible on the Lab’s unclassified network, said NSRC collections management team leader Patricia Cote (WRS-NSRCMS). The cards, named after Dorothy McKibbin, who was known as the gatekeeper of Los Alamos because she was often the first point of contact for new hires, have become symbolic of the Lab when the world’s greatest minds secretly gathered to create the first atomic bomb and end World War II.
The damage mechanisms that lead to failure in engineering alloys have been studied extensively, but converting this knowledge into constitutive models that are suitable for engineering-scale analysis remains a challenge. Evolution laws for continuum damage have been developed in the past and have proven effective but suffer from many non-physical assumptions that inhibit the overall accuracy of the model. Further, the assumptions inherent in these existing models prevent them from being applicable to a broad class of materials. At the same time, computational models of fine-scale damage mechanisms continue to advance making it tractable to generate large training data sets through computer simulation. Data-driven machine learning approaches can leverage these data sets to avoid making limiting assumptions, and instead produce models directly from the results of microstructural simulations and/or experiments. Many of these machine learning approaches are rapid and accurate, but they offer little to no insight into the underlying relationships among state variables being discovered. Conversely, genetic programming symbolic regression (GPSR) is a machine learning method that produces analytic expressions relating the state variables, allowing maximal insight and interpretability. To that end, we propose using GPSR as a data-driven method of obtaining microstructurally informed continuum damage models. Data is generated using microstructural simulations of damage evolution, parameterized over microstructural statistics (i.e., pore shape) and nominally applied deformations. Analytic expressions for damage evolution are obtained from the data using GPSR, and these expressions are then utilized within a continuum constitutive model. Overall, this approach is a promising method of automatically obtaining analytic relations describing constitutive phenomena in a material.
This technical report summarized the contribution of the DADAIST project funded by the Data Model Convergence Initiative via the Laboratory Directed Research and Development (LDRD) investments at Pacific Northwest National Laboratory (PNNL). Specifically, we report the development of the NeuroMANCER (Neural Modules with Adaptive Nonlinear Constraints and Efficient Regularizations), a new open-source Scientific Machine Learning library for formulating and solving parametric constrained optimization problems, physics-informed system identification, and parametric optimal control problems. NeuroMANCER is using differentiable programming to combine modern data-driven models and optimization modeling language into a coherent algorithmic and software framework. NeuroMANCER is a Pytorch-based framework and adopts much of its philosophy focused on research and development, rapid prototyping, and streamlined deployment. Strong emphasis is given to extensibility, interoperability with the PyTorch ecosystem, and quick adaptability to custom domain problems. Neuromancer repository contains a comprehensive library of differentiable modules, including custom activation functions, matrix factorizations, deep learning architectures, neural differential equations, differential equation solvers, implicit layers such as iterative solvers, high-level API for symbolic expressions, API for modeling and control of dynamical systems, and extensive set of tutorial code examples in the form of python scripts and jupyter notebooks.
Computer programs are often viewed as collections of functions – each function has parameters (inputs) and computes a return value, and each has potential side effects that modify program state (outputs). In this research, a Sandia symbolic execution tool designed to support “human-in-the-loop” analysis was modified to automatically create “function summaries,” and a new tool, “diaboolical,” was created to support enhancing readability of the summary using a novel approach to bit-vector simplification that leverages logic synthesis and tautologies. For this effort, students at Auburn University created several finite state machines (FSMs) to serve as exemplars for high-consequence systems. Function summaries for each of the machines were obtained, and then portions of the summaries were simplified using both diaboolical and the simplification procedure of a popular SMT solver. A comparison of the results shows that diaboolical can often produce smaller function summaries, with expression length improvements over the unsimplified function summaries ranging from 0% to 90% for diaboolical and 0% to 65% for the SMT solver, though diaboolical had a significantly greater cost in time. Diaboolical was evaluated against a collection of “arbitrary” C-code as well as FSM exemplars, and for both datasets it achieved an approximately 10% improvement in expression length compared to simplifications that could be obtained using existing techniques. Function summaries can assist assurance efforts that evaluate existing systems and their executable code. A smaller function summary is likely easier for humans to understand and could thus increase the ability and efficacy of assurance practices centered around the analysis of executable artifacts.
The objective is to better predict the bulk-scale mechanical behavior of porous metals that have over pressurized cavities (e.g., irradiated metals with helium bubbles) by quantifying the complex coupling among cavity aspects (e.g., size distribution, inhomogeneous overpressure values, spatial arrangement) and metal properties (e.g., rate-dependency, crystallographic lattice). This requires up-scaling local mechanical fields from the single crystal scale and will be accomplished using a homogenization approach that combines full-field numerical simulations, analytical formalisms, and physics-informed machine learning to produce symbolically-defined constitutive equations (e.g., gauge functions). These equations will satisfy the objective because they enable computationally efficient predictions that approach the accuracy of computationally expensive full-field numerical simulations, abide by theoretical requirements (e.g., conservation of energy, work conjugacy), and retain the transparency of analytical models.
A national symbol of the United States, the Bald Eagle (Haliaeetus leucocephalus) plays an important role in the riverine ecosystem at the U. S. Department of Energy’s (DOE) Hanford Site. Historically, Bald Eagles have occupied the Hanford Site during the winter and early spring and have more recently remained through the nesting season. Monitoring is essential to maintain current biological information about Bald Eagle abundance and distribution on the Hanford Site, to ensure compliance with protection regulations, and to inform future protection and management efforts. This monitoring report provides an overview of Bald Eagle activity on the Hanford Site, along with Bald Eagle management guidelines and monitoring objectives between November 2021 and June 2022.
A national symbol of the United States, the Bald Eagle (Haliaeetus leucocephalus) plays an important role in the riverine ecosystem at the U. S. Department of Energy (DOE) Hanford Site. Historically, Bald Eagles have occupied the Hanford Site during the winter and early spring and have more recently remained through the nesting season. Monitoring is essential to maintain current biological information about Bald Eagle abundance and distribution on the Hanford Site, to ensure compliance with protection regulations, and to inform future protection and management efforts. This monitoring report provides an overview of Bald Eagle activity on the Hanford Site, along with Bald Eagle management guidelines and monitoring objectives between November 2020 and June 2021.
A national symbol of the United States, the Bald Eagle ( Haliaeetus leucocephalus ) plays an important role in the riverine ecosystem at the U. S. Department of Energy (DOE) Hanford Site. Historically, Bald Eagles have occupied the Hanford Site during the winter and early spring and have more recently remained through the nesting season. Monitoring is essential to maintain current biological information about Bald Eagle abundance and distribution on the Hanford Site, to ensure compliance with protection regulations, and to inform future protection and management efforts. This monitoring report provides an overview of Bald Eagle activity on the Hanford Site, along with Bald Eagle management guidelines and monitoring objectives between November 2019 and June 2020.
A national symbol of the United States, the Bald Eagle ( Haliaeetus leucocephalus ) plays an important role in the riverine ecosystem on the U. S. Department of Energy’s (DOE) Hanford Site. Historically, Bald Eagles have occupied the Hanford Site during the winter and early spring, and have more recently remained through the nesting season. Monitoring is essential to maintain current biological information about Bald Eagle abundance and distribution on the Hanford Site, to ensure compliance with protection regulations, and to inform future protection and management efforts. This monitoring report provides an overview of Bald Eagle activity on the Hanford Site, along with Bald Eagle management guidelines and monitoring objectives between November 2018 and July 2019.
The deployment of artificial intelligence systems in critical applications requires higher levels of assurance for safety, security, and interpretability. While neurosymbolic (NESY) approaches combining neural networks with symbolic reasoning offer potential advantages for assured AI, existing differentiable neurosymbolic frameworks face significant limitations including computational overhead and performance constraints. This report investigates the ISED (InferSampleEstimateDescend) framework as an alternative approach that enables neurosymbolic learning without requiring endtoend differentiability. We evaluate ISED’s utility for geointelligence applications by comparing neurosymbolic models against standard neural networks on aircraft classification tasks using the RarePlanes and MTARSI imagery datasets. Our results demonstrate that while ISEDbased models achieve slightly lower accuracy (89.7% vs 92.1% on RarePlanes; 91.1% vs 92.5% on MTARSI), they provide critical explainability capabilities that enable tracing incorrect predictions back to specific attribute misclassifications. We also present an automated pipeline that generates both attributeclass mappings and neurosymbolic model architectures from natural language descriptions, significantly reducing the manual effort required for NESY model deployment. These findings suggest that ISED offers a promising direction for developing assured AI systems where interpretability and reasoning transparency are prioritized alongside performance.
The primary objective of this project is to strengthen the trustworthiness of AI systems by designing algorithms that make their internal decision-making processes more understandable to human users. This involves creating clear, interpretable explanations for AI decisions and developing metrics to assess these explanations' validity and reliability. Significant progress has been achieved through (i) developing symbolic explanations, (ii) generating meaningful interpretive insights, (iii) establishing accuracy and confidence metrics, and (iv) devising methods to evaluate the knowledge boundaries of AI models. To date, the research findings have been shared in peer-reviewed publications, with accompanying scientific and technical information (STI) detailed below.
The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.
Motivation: Predicting microbial gene fitness across environmental conditions remains a central challenge for predictive phenomics and autonomous experimentation. Fitness assays generate large volumes of genotype–phenotype measurements difficult to integrate with experimental metadata and biological function in a form that supports mechanistic reasoning. Knowledge graphs offer a semantic framework for unifying modalities and enabling context-aware inference. Results: We build GIMME (Graph Inference for Microbial Metabolism Exploration), a semantically grounded knowledge graph that unifies gene fitness measurements spanning 10 Pseudomonas species with experimental metadata and biological context. Media are decomposed into chemical components and experiments carry structured links to natural-language descriptions. The resulting graph supports two inference modes: (1) symbolic graph traversal to surface candidate gene–environment and gene–chemical associations, and (2) learned inference using heterogeneous graph neural networks that propagate information across neighborhoods. We formulate link regression over (gene, media, experiment) triplets, combining learned gene embeddings with pretrained LLM sourced text embeddings of node descriptions to predict gene fitness. We then augment a baseline MLP with an auxiliary message-passing encoder (GraphSAGE/GAT) that propagates information over gene–protein–function and media–chemical subgraphs, and fuse the two pathways with a gated residual connection. This approach produces strong agreement with held-out fitness measurements (GraphSAGE Pearson r 0.74) while also highlighting inference challenges in extreme-fitness regimes. We aggregate GAT edge-attention weights by relation type and layer to estimate which biological and environmental relations most influence fitness predictions. Conclusion: This work explores using knowledge graphs as “context graphs” for microbial phenotype prediction. They provide a rich substrate which enables explainable retrieval of supporting evidence, and provides a natural bridge to autonomous workflows that prioritize the next experiment.
The new quaternary phases Eu5Zn2As5O and Eu5Cd2As5O have been synthesized by metal flux reactions and their structures have been established through single-crystal X-ray diffraction. Both compounds crystallize in the centrosymmetric space group Cmcm (No. 63, Z = 4; Pearson symbol oC52), with unit cell parameters a = 4.3457(11) Å, b = 20.897(5) Å, c = 13.571(3) Å; and a = 4.4597(9) Å, b = 21.112(4) Å, c = 13.848(3) Å, for Eu5Zn2As5O and Eu5Cd2As5O, respectively. The crystal structures include one-dimensional double-strands of corner-shared MAs4 tetrahedra (M = Zn, Cd) and As–As bonds that connect the tetrahedra to form pentagonal channels. Four of the five Eu atoms fill the space between the pentagonal channels and one Eu atom is contained within the channels. An isolated oxide anion O2– is located in a tetrahedral hole formed by four Eu cations. Applying the valence rules and the Zintl concept to rationalize the chemical bonding in Eu5M2As5O (M = Zn, Cd) reveals that the valence electrons can be counted as follows: 5 × [Eu2+] + 2 × [M2+] + 3 × [As3–] + 2 × [As2–] + O2–, which suggests an electron-deficient configuration. The presumed h+ hole is confirmed by electronic band structure calculations, where a fully optimized bonding will be attained if an additional valence electron is added to move the Fermi level up to a narrow band gap (Eu5Zn2As5O) or pseudo-gap (Eu5Cd2As5O). In order to achieve such a formal charge balance, and hence, narrow-gap semiconducting behavior in Eu5M2As5O (M = Zn, Cd), europium is theorized to be in a mixed-valent Eu2+/ Eu3+ state.
Crystals of a new ternary compound in the Ca-In-As family, Ca3InAs3, have been successfully synthesized via flux growth techniques. This is only the third known compound between the respective elements. As elucidated by single-crystal X-ray diffraction measurements, Ca3InAs3 crystallizes in the orthorhombic space group Pnma (No. 62, Pearson symbol oP28) with unit cell parameters a = 12.296(2) Å, b = 4.2553(7) Å, and c = 13.735(2) Å. The smallest building motifs of the structure are InAs4 tetrahedra, which are connected to one another by shared As corners, forming infinite [InAs2As2/2] chains. The latter propagate along the crystallographic b-axis. The As-In-As bond angles within the InAs4 tetrahedra deviate from the ideal 109.5° value and range from 98.12(2)° to 116.53(2)°, attesting to a small distortion from the regular tetrahedral geometry. Electronic structure calculations indicate the opening of a bandgap, consistent with the expected (Ca2+)3(In3+)(As3–)3 formula breakdown based on conventional oxidation numbers. The calculations also show that the Ca–As interactions are an intermediate between covalent and ionic, while providing evidence of strong covalent features of the In–As interactions. Weak s-p hybridization of In states was observed, supporting the experimentally found deviation of the InAs4 moiety from the ideal tetrahedral symmetry.