Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis sets”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Machine-Learning of Nonlocal Kernels for Anomalous Subsurface Transport from Breakthrough Curves

Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evaluation of global teleconnections in CMIP6 climate projections using complex networks

In climatological research, the evaluation of climate models is one of the central research subjects. As an expression of large-scale dynamical processes, global teleconnections play a major role in interannual to decadal climate variability. Their realistic representation is an indispensable requirement for the simulation of climate change, both natural and anthropogenic. Therefore, the evaluation of global teleconnections is of utmost importance when assessing the physical plausibility of climate projections. We present an application of the graph-theoretical analysis tool δ-MAPS, which constructs complex networks on the basis of spatio-temporal gridded data sets, here sea surface temperature and geopotential height at 500 hPa. Complex networks complement more traditional methods in the analysis of climate variability, like the classification of circulation regimes or empirical orthogonal functions, assuming a new non-linear perspective. While doing so, a number of technical tools and metrics, borrowed from different fields of data science, are implemented into the δ-MAPS framework in order to overcome specific challenges posed by our target problem. Those are trend empirical orthogonal functions (EOFs), distance correlation and distance multicorrelation, and the structural similarity index. δ-MAPS is a two-stage algorithm. In the first place, it assembles grid cells with highly coherent temporal evolution into so-called domains. In a second step, the teleconnections between the domains are inferred by means of the non-linear distance correlation. We construct 2 unipartite and 1 bipartite network for 22 historical CMIP6 climate projections and 2 century-long coupled reanalyses (CERA-20C and 20CRv3). Potential non-stationarity is taken into account by the use of moving time windows. The networks derived from projection data are compared to those from reanalyses. Our results indicate that no single climate projection outperforms all others in every aspect of the evaluation. But there are indeed models which tend to perform better/worse in many aspects. Differences in model performance are generally low within the geopotential height unipartite networks but higher in sea surface temperature and most pronounced in the bipartite network representing the interaction between ocean and atmosphere.

58 GEOSCIENCES↗
DOE OSTI2020

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗
DOE OSTI2026

Operationally induced preferred basis in unitary quantum mechanics

The preferred-basis problem and the definite-outcome aspect of the measurement problem persist even if the detector is modeled unitarily, because experimental data are necessarily represented in a Boolean event algebra of mutually exclusive records whereas the theoretical description is naturally formulated in a noncommutative operator algebra with continuous unitary symmetry. This change of mathematical type constitutes the core of the 'cut': a structurally necessary interface from group-based kinematics to set-based counting. In the presented view the basis relevant for recorded outcomes is not determined by the system Hamiltonian alone; it is induced by the measurement mapping, i.e., by the detector channel together with the coarse-grained readout that defines an instrument. The probabilistic mapping is anchored in symmetry and measure theory: by Gleason-type uniqueness (Gleason for projections in $d>2$ and Busch's extension for Positive Operator-Valued Measures (POVMs) including $d=2$), the trace rule is the unique probability measure consistent with additivity over exclusive events and basis-independence of the unitary sector. A compact qubit--pointer model yields an induced unsharp POVM $E_\pm=\tfrac12(\id\pm η\,σ_z)$ with $η$ fixed by pointer resolution, displaying explicitly how the detector induces the relevant basis. Finally, nested-observer paradoxes are tightened into a non-composability lemma: joint assignment of outcome propositions is obstructed unless a joint instrument exists. This relocates the origin of randomness to the stochasticity of the transition rules.

Pronskikh, Vitaly [Fermilab] (ORCID:00000002518174↗
DOE OSTI2024

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2020

Characterizing The Mechanical Properties of Polyurethane Foams

The Department of Energy has many radioactive facilities that are on the path of deactivation and decommissioning (D and D). These facilities can sit cold and dark for many years before final disposition, and must be maintained to ensure no radiological release occurs in the interim. Improvements and additions to D and D tool sets can greatly: Save time and money, Reduce worker risk. Conventional fixatives widely used often take the form of paints or films that are not readily applicable to 3-dimensional void spaces. Foams are one promising platform that may offer solutions to a number of contaminated problem sets such as: Gloveboxes, Pipes, Tanks. SRNL is working to characterize various commercial foams that would encapsulate the interior volume of a given space and are capable of immobilizing any remaining contamination. One key performance metric of these foams is how well they will perform in an accident scenario. To this end, SRNL is researching the mechanical properties of these foams to ensure that the material can withstand the environment of application while maintaining structural integrity. ASTM standard E3191-18 served as a guiding document for this project, outlining the requirements that foaming fixatives used for the mitigation of radioactive contamination need to meet before being implemented. Objective: Quantification of the mechanical properties of 6 commercial polyurethane foams was performed using multiple ASTM standards to record measurements for: Compression testing: Flexible Foams, Force required to produce 50% compression, Rigid Foams, Compressive and apparent modulus, Point of 10% core deformation, The 'Zero Deformation' point, Compressive strength, Yield point. Tensile testing: Tensile strength, Tensile stress, Percent elongation. Experiment 1: An electromechanical compressive tester (MTS Criterion Series 43) was used to evaluate 6 foams (4 flexible, 2 rigid). Per ASTM D1621, each rigid foam was compressed at 10% of the measured thickness per minute until the sample was 13% of it's original thickness. Per ASTM D35/4, each flexible foam was pre-flexed twice to 80% original thickness at a rate of 250 mm/min, then compressed to 80% original thickness at a rate of 50 mm/min. Experiment 2: An electromechanical tensile tester (MTS Criterion Series 43) was used to evaluate 6 foams (4 flexible, 2 rigid). Per ASTM D1623, rigid foams were pulled apart at a rate of 1.27 mm/min until the sample broke. Per ASTM D35/4, flexible foams were pulled apart at a rate of 500 mm/min until the sample broke. The strongest material in both compression and tensile testing scenarios was found to be the rigid intumescent polyurethane Hilti foam. The experiments revealed that the Hilti foam in a tensile scenario had a peak stress value that was larger than the closest competitor by a factor of 2.3 and a compressive yield point that was larger than the closest competitor by a factor of 1.4, indicating that the Hilti foam is the best choice for implementation in mechanically harsh environments. The performance metrics measured can serve as a basis for future mechanical tests that would help set relevant ASTM standards (E3191) for intumescent polyurethane foams in fixating applications. Tests like surface adhesion, impact, and flame tests would serve as a better indicator as to how this material would perform in environmentally harsh scenarios often found in decommissioned nuclear facilities. Further tests of the foams' intumescent properties would also be important should these foams be implemented in environmentally harsh scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
DOE OSTI2020

Data-driven selection of stiff chemistry ODE solver in operator-splitting schemes

Most computational fluid dynamics simulations of practical combustion applications employ operator-splitting schemes, where chemistry and transport are separated and integrated with distinct numerical methods. The changes in composition due to chemistry are evaluated by solving ordinary differential equations (ODE) in each cell of the computational domain, which typically dominates the computational cost when detailed chemistry is considered. In this work, a data-driven approach for the selection of chemistry ODE solvers in operator-splitting schemes is presented. Neural networks are used to predict the ODE solvers CPU times and errors for a given thermochemical state. This allows the selection of an optimal ODE solver on a cell-by-cell, timestep-by-timestep basis. The models are trained using a wide set of thermochemical states generated through partially-stirred reactors and flames simulations. The methodology is validated by quantifying the prediction errors, the classification accuracy, and the computational speedup. The model predicts the optimal ODE solver for 70 to 95% of the validation cases and decreases the computional cost by a factor of 3 or more. The generalizability of the methodology to different chemical mechanisms and different fuels is assessed and it is shown that the model’s performance is only slightly degraded and its applicability is significantly enhanced if the inputs to the neural networks are restricted to a small set of thermochemical state variables present in most chemical mechanisms. In conclusion, the models are used in an homogeneous reactor case and a multi-dimensional CFD simulation of a diesel spray at high pressure where a speedup of more than 3 is achieved.

42 ENGINEERING↗
DOE OSTI2020

Cluster dynamics modeling of irradiation growth in single crystal Zr

Irradiation growth of single-crystal Zr is modeled with a reduced set of cluster dynamics equations. Nucleation and growth of basal and prismatic dislocation loops are both accounted for in the model. Equations are developed for the time evolution of single point defects, a limited number of point defect clusters, interstitial and vacancy loop nucleation rates, as well as for the growth of vacancy loops on basal planes and interstitial loops on prismatic planes. Reduction of the usual infinite hierarchy of cluster dynamics equations to the simple set studied here is justified on the physical basis of the stability of small loops once nucleated. This simplified cluster dynamics model with a small number of adjustable parameters avoids the complexity of explicit representation of higher order point defect clusters. The model shows consistency with experimental observations of the following aspects: (1) the growth rates of Zr crystals along the a- and c-axes; (2) the onset dose for breakaway irradiation growth; (3) the saturation dislocation loop densities of vacancy and interstitial loops; (4) the effects of cold work and temperature on irradiation growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Programmable Heisenberg interactions between Floquet qubits

Abstract The trade-off between robustness and tunability is a central challenge in the pursuit of quantum simulation and fault-tolerant quantum computation. In particular, quantum architectures are often designed to achieve high coherence at the expense of tunability. Many current qubit designs have fixed energy levels and consequently limited types of controllable interactions. Here by adiabatically transforming fixed-frequency superconducting circuits into modifiable Floquet qubits, we demonstrate an XXZ Heisenberg interaction with fully adjustable anisotropy. This interaction model can act as the primitive for an expressive set of quantum operations, but is also the basis for quantum simulations of spin systems. To illustrate the robustness and versatility of our Floquet protocol, we tailor the Heisenberg Hamiltonian and implement two-qubit iSWAP, CZ and SWAP gates with good estimated fidelities. In addition, we implement a Heisenberg interaction between higher energy levels and employ it to construct a three-qubit CCZ gate, also with a competitive fidelity. Our protocol applies to multiple fixed-frequency high-coherence platforms, providing a collection of interactions for high-performance quantum information processing. It also establishes the potential of the Floquet framework as a tool for exploring quantum electrodynamics and optimal control.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

LANSCE Science Overview [Slides]

LANSCE’s combination of beam power, flexibility, and authorization basis uniquely positions it to address a broad set of NNSA (and beyond) science. We have a programmatically motivated science vision for the facility that extends through the next several decades.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Microbial partner (MiPner) analysis

Introduction Although a few bacteria have been studied in great depth, relatively little is known about the characteristics of microbe-microbe interactions that occur within ecosystems on a daily basis. A simple, robust technique was developed to set up the foundation for investigating pairwise bacterial-bacterial interactions, using cell-cell binding as a self-selective mechanism to identify interesting bacterial species pairs. Methods Using a Serratia marcescens strain (SMC43) isolated from Georgia soil as a “bait”, specific bacteria were purified by their specificity in binding SMC43 bacteria that were themselves attached to a wooden applicator stick. Results The isolated Microbial Partners (MiPners) were greatly enriched for members of the genera Sphingobium and Caulobacter . Two streaked MiPners were unable to grow on the plates employed after separation from SMC43to be separated from, and grow on the plate type tested without, SMC43. Discussion This suggests that the MiPner technology will be one strategy for purifying bacteria that were previously recalcitrant to culturing.

Bennetzen, Jeffrey L.↗
DOE OSTI2022

Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry

Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d 1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d 1 MOPs to be suitable for nearly eliminating all energetic delocalization and static correlation errors. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization and static correlation errors within a non-empirical framework.

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

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs 2 BB'Cl 6 ) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (t) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs 2 BB'Cl 6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs 2 [Alk] + [TM] 3+ Cl 6 , where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs 2 BB'Cl 6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

14 SOLAR ENERGY↗
DOE OSTI2022

U-splines: Splines over unstructured meshes

U-splines are a novel approach to the construction of a spline basis for representing smooth objects in Computer-Aided Design (CAD) and Computer-Aided Engineering (CAE). A spline is a piecewise-defined function that satisfies continuity constraints between adjacent cells in a mesh. U-splines differ from existing spline constructions, such as Non-Uniform Rational B-splines (NURBS), subdivision surfaces, T-splines, and hierarchical B-splines, in that they can accommodate local variation in cell size, polynomial degree, and smoothness simultaneously over more varied mesh configurations. Mixed cell types (e.g., triangle and quadrilateral cells in the same mesh) and T-junctions are also supported, although the continuity of interfaces with triangle and tetrahedral cells is limited in the present work. The U-spline algorithm introduces a new technique for using local null space solutions to construct basis functions for the global spline null space problem. The U-spline construction is presented for curves, surfaces, and volumes with higher dimensional generalizations possible. Lastly, a set of requirements are given to ensure that the U-spline basis is positive, forms a partition of unity, is complete, and is locally linearly independent.

42 ENGINEERING↗
DOE OSTI2023

Implementation of a Frenkel exciton-based controlled phase shifter

Frenkel excitons are of interest for ultrafast switching applications due to the femtosecond timescale at which coherent exciton transfer between chromophores takes place. The Frenkel Hamiltonian, which governs the dynamics of Frenkel excitons, belongs to a class of Hamiltonians that enable universal quantum computation. It is thus of interest to determine how a complete set of exciton-based gates for quantum computing could be constructed as aggregates of chromophores. Here, we demonstrate that a controlled phase shifter can be constructed as a pair of exciton transmission lines that employ two types of chromophores and that are cross coupled by a two-exciton interaction. This element facilitates the construction of controlled basis-change gates, thereby enabling the implementation of a complete set of exciton-based gates for universal quantum computation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Bill Savings vs. Backup Power: Evaluating operational tradeoffs for home solar+storage systems [Slides]

This study explores tradeoffs between the use of home solar+storage systems for backup power versus day-to-day utility bill savings. The study focuses specifically on the “reserve setting” available with most home battery storage systems, which allow the customer to maintain some minimum level of storage in reserve in case of an unforeseen power interruption. The more capacity that is held in reserve, the greater the customer’s ability to ride-through possible power interruptions, but less capacity is then available to manage utility bills on a day-to-day basis. This study evaluates this operational tradeoff across a diverse set of locations and residential electricity tariff structures, relying on Berkeley Lab’s PRESTO model to stochastically simulate power interruption events, and exploring a range of sensitivities, including variations in customer value of lost load (VoLL), interruption frequency, and other key drivers. The results show that, in most circumstances, the opportunity cost of holding storage capacity in reserve, in terms of foregone bill saving, tends to outweigh any gains in reliability value associated with mitigated power interruptions. This finding is robust across tariff structures and across most of the sensitivities considered, including those related to rate level, customer load level, and storage sizing. There are a limited set of circumstances where raising the reserve setting improves the overall customer value (comprised of bill savings plus reliability value). Specifically, that exception occurs when all of the following conditions apply: (a) the customer resides in a location with exceptionally poor reliability, (b) the customer has exceptionally high VoLL; (c) the customer is on a net billing rate or on a TOU rate that allows grid discharging but not grid charging; and (d), depending on the location, the price arbitrage differential on that rate is relatively small. In all other circumstances analyzed, total customer value declines with reserve level.

14 SOLAR ENERGY↗