Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational formulation”

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

Next Generation Solvent Vapor Pressure Testing

Previous work by Savannah River National Laboratory (SRNL) indicated that the actual Next Generation Solvent vapor pressure would be higher than the Original Caustic Side Solvent Extraction (CSSX) solvent vapor pressure, but below a bounding vapor pressure using Raoult’s Law. Since solvent vapor pressure has a significant impact on Composite Lower Flammability Limits (CLFL) and attendant accident analyses, obtaining an additional margin from the bounding NGS vapor pressure would be valuable. In order to quantify how much margin might be gained from the NGS bounding solvent, SRNL researchers were requested to perform vapor pressure testing with the Next Generation Solvent (NGS) by Savannah River Mission Completion (SRMC). The vapor pressure curve for the NGS formulation set to be deployed at the Salt Waste Processing Facility (SWPF) has been determined by SRNL up to 55°C (131°F) using headspace Gas Chromatography (GC). It was expected that NGS would have a higher vapor pressure than the Original CSSX solvent; however, experimental results indicate a lower vapor pressure. At this time, it is uncertain if this difference is due to the 7x increase in concentration of the large calixarene in the solvent (0.007M BOBCalix in Original CSSX solvent vs. 0.05M MaxCalix in NGS) or to minor batch-to-batch variations in Isopar-L constituents. The NGS and the Original CSSX solvent vapor pressure data were fitted to the Antoine Equation. The Antoine equation gives a more accurate representation of the expected vapor pressure of the solvents outside of the temperature range tested. The Antoine equation fitting for NGS is given below: $P=10^{6.364⁻\frac{1788.4}{T+219.3}}$. Where, p is the vapor pressure (partial pressure) of NGS in mmHg and T is the temperature in °C. It is recommended that SRMC either continue using the more conservative Isopar-L vapor pressure calculations at SWPF with the equation developed for the Original CSSX solvent, or the equation presented above for the NGS solvent. Additionally, it is recommended to study the variability in Isopar-L vapor pressure between lots.

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

A review of crosstalk polymorphic circuits and their scalability

Using a control variable, the functionality of Polymorphic circuits can be modified, making them adaptable and useful for reconfiguring circuit behavior — all the way from gate level to system level. State-of-the art polymorphic circuits are based on custom non-linear circuit design or emerging devices such as ambipolar FET, configurable magnetic devices etc. While some of these approaches are inefficient in performance, others involve exotic devices. The Crosstalk computing based polymorphic circuits offer a fresh perspective. In Crosstalk, the interconnect interference between nanoscale metal lines is intentionally engineered to exhibit the programmable Boolean logic behavior. This approach relies on the coupling between metal lines and not on the transistors for computing, resulting in better scalability, security by obscurity, and fault tolerance by reconfiguration. Our novel approach is backed by the mathematical formulation that conveys the rationale to generalize and achieve a wide variety of polymorphic circuits. Our experiments, including design, simulation, and Power Performance Area (PPA) characterization results indicate that crosstalk circuits provide significant improvement in transistor count (about 3x), switching energy (2x), and speed (1.5x) for polymorphic logic circuits. In the best-case scenario, the transistor count reduction is 5x. This paper presents Crosstalk computing’s fundamentals, polymorphism and the scalability aspects to compete/co-exist with CMOS for digital logic implementations below 10 nm. Our scalability study uses Open Source 7 nm PDK, considers all process variation aspects and accommodates worst-case scenarios. The study results for various benchmark circuits show that the Crosstalk technology is a viable alternative to CMOS for digital logic implementations below 10 nm, having 48% density, 57% power, and 10% performance gains over equivalent CMOS counterparts. Finally, we compare Crosstalk Polymorphic Circuit design technique with similar approaches described in related works and discuss its features and constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enhancing Cyber Resilience of Networked Microgrids using Vertical Federated Reinforcement Learning

This paper presents a novel federated reinforcement learning (Fed-RL) methodology to inject sufficient resiliency into the operations of the network of microgrids. We consider adversarial actions to the voltage and power control loop reference signals at the grid forming (GFM) inverters in the microgrids which are essential to integrate renewable resources. Therefore, we formulate a resilient reinforcement learning training setup that uses these adversarial injections to generate episodic trajectories and train the RL agents to alleviate their impact on performance. To circumvent the concerns about data-sharing and privacy for different owners of the microgrids in the networked setting, we bring in the aspects of the federated operation to propose novel Fed-RL algorithms. As the dynamics of each microgrid are coupled due to electrical interlinks, the conventional federated RL approaches using decoupled independent environments are not applicable, which leads us to propose a multi-agent vertically federated variation of actor-critic algorithms, namely federated soft actor-critic (FedSAC). We have performed numerical simulations on an IEEE 123-bus benchmark test feeder with three microgrids by creating a customized simulation setup by encapsulating the microgrid dynamic simulations in GridLAB-D/HELICS co-simulation platform with the OpenAI Gym environment and validated the proposed resilient and secured learning methodology.

Artificial Intelligence (AI), reinforcement learni↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ionic Tunability of Conjugated Polyelectrolyte Solutions

Conjugated polyelectrolytes (CPEs), which combine a π-conjugated polymer backbone with pendant ionic functionalities, offer an opportunity for electrostatic control of materials properties. In this work, the mesoscale morphology and physical properties of a high-mobility conjugated polyelectrolyte are tuned by the addition of salt, variation of the charge-compensating counterion, and complexation with an oppositely charged polyelectrolyte containing the same π-conjugated backbone. In systems with a single polyelectrolyte species, added ions screen the electrostatic repulsions stabilizing the gel-phase, resulting in the dissolution of ionic cross-links and hydrophobic collapse. Further, exchanging the charge compensating counterion is found to enable finer structural control of the solution behavior and modulation of the self-doping behavior. Finally, novel CPE–CPE complexes resulting in dense solutions and gels of semiconductive material are produced by combination of oppositely charged polyelectrolytes. Such concentrated CPE formulations should be useful materials for mixed electronic–ionic conduction and pseudocapacitative energy storage.

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

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗
DOE OSTI2021

Variational approach to relative entropies with an application to QFT

Abstract We define a new divergence of von Neumann algebras using a variational expression similar in nature to Kosaki’s formula for Umegaki’s relative entropy. Our divergence satisfies several of the usual desirable properties, upper bounds the sandwiched Renyi entropy and reduces to the fidelity in a limit. As an illustration, we use the formula in quantum field theory to compute our divergence between the vacuum in a bipartite system and an “orbifolded”—in the sense of a conditional expectation—system in terms of the Jones index. We take the opportunity to point out an entropic certainty relation associated with an inclusion of von Neumann factors related to the relative entropy. This certainty relation has an equivalent formulation in terms of error correcting codes.

Hollands, Stefan↗
DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗
DOE OSTI2021

Application of molecular convergent close-coupling cross sections in a collisional radiative model for the triplet system of molecular hydrogen

Collisional radiative (CR) models for molecular hydrogen are of high relevance for performing qualitative and quantitative analysis of excited-state population densities measured in plasmas or predicting the dependence of plasma emission on parameter variations. Although the development of such models for H 2 started decades ago, major uncertainties still exist regarding the most important set of input parameters, namely the cross sections for electron-impact excitation. The deviations between cross sections from different datasets are particularly pronounced in the energy region close to the threshold energy, strongly increasing the uncertainty of CR models applied to low-temperature plasmas. This paper presents experimental validation of a set of newly calculated non ro-vibrationally resolved electron-impact cross sections calculated for the triplet system of H 2 using the molecular convergent close-coupling method in the adiabatic-nuclei formulation. These cross sections are implemented into a CR model based on the flexible solver Yacora. A first comparison of CR calculations with the different datasets to experimentally-determined population densities is performed at a planar ICP discharge for varying pressure (between 1 and 10 Pa) and RF power (between 700 and 1100 W). For the experimentally-accessible electron temperature and density range (2.5–10 eV and 1.8–3.3 × 10 16 m –3 , respectively), very good agreement between the model and experiment is obtained using the new data set, in contrast to previously used cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling and simulation of earthquake soil structure interaction excited by inclined seismic waves

Presented is an application of wave potential formulation (WPF) together with domain reduction method (DRM) to modeling earthquake soil structure interaction (ESSI) behavior in horizontally layered ground under inclined incident seismic waves. Wave potential formulation is used to develop a spatially varying, inclined seismic wave field from incident Primary (P) and Secondary (S) waves that propagate through layered ground. Developed seismic wave field is then used to develop effective forces for Domain Reduction Method that are then used for analyzing ESSI response of a soil structure system. Developed methodology, called WPF-DRM, is verified using analytic solution for a free field response of layered ground subjected to inclined incident waves. Developed WPF-DRM methodology is illustrated through analysis of an ESSI response of a deeply embedded structure, a small modular reactor (SMR) subjected to incident S wave polarized in vertical plane (SV) with variation in inclinations and frequencies. Presented example highlights the influences of incident wave inclination and frequency on ESSI response of analyzed SMR.

58 GEOSCIENCES↗
DOE OSTI2024

Identification of Suitable Vacuum Gas Oils as Plasticizers Using HT-GC × GC-HRMS

Vacuum gas oils (VGOs) have long been heralded as effective plasticizers due to their high boiling-points and lubricating properties. One of these VGOs, HyVac Oil 93050, is a paraffinic plasticizer that is currently utilized in polymer–plasticizer explosive formulations, while this oil effectively increases the elasticity and decreases the sensitivity for current formulations. We seek to identify alternative VGOs with similar density, viscosity, molecular composition, and impurities for future formulations. In this study, 18 VGOs, including HyVac Oil 93050, were initially evaluated for their density and viscosity. The oils were then ranked based on physical characteristics and analyzed for molecular composition using high-temperature comprehensive two-dimensional gas chromatography with high-resolution time-of-flight mass spectrometry (HT-GC × GC-HRMS). Further, the HT-GC × GC-HRMS chromatograms of the top 10 most similar VGOs to HyVac Oil 93050, according to density, were compared utilizing the Pearson correlation coefficient. Three of the oils with densities similar to those of HyVac Oil 93050 had a Pearson correlation within the lot-to-lot variation of HyVac Oil 93050. For the top 3 candidates, Pearson correlation was then utilized as a feature selection technique to discover significant chemical differences. Positive chemical ionization (PCI) and negative chemical ionization (NCI) HT-GC × GC-HRMS chromatograms were also evaluated and aided in the discovery and identification of several key compounds including additives that acted as stabilizers for the VGOs. After this further chemical analysis, two of the original 17 potential VGOs were determined to be physically and chemically similar to HyVac Oil 93050.

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

Variation of Aliphatic Diisocyanates in Bio-Based TPUs

The utilization of bio-based materials for polymer production poses a challenge to both the industrial and academic sectors due to the availability and production costs of the necessary raw materials. Diisocyanates necessary for polyurethane synthesis have posed a particular challenge, given the paucity of natural diamine precursors for traditional phosgenation. We recently developed a phosgene-free flow chemistry methodology that allows for the safe, efficient, and scalable preparation of diisocyanates from naturally produced diacids. This chemistry broadens the potential for the development of renewable diisocyanates for a broad variety of applications. Here, we expand upon this work by demonstrating the scaled production of a panel of bio-based linear, aliphatic diisocyanates, with chain lengths ranging from four to eight carbons in good yields of 57–81% and high purity (~97%). Furthermore, these are applied to the synthesis of thermoplastic polyurethanes (TPUs) containing up to a 100% renewable carbon content. TPUs formulated using shorter carbon-chained diisocyanates displayed a higher tensile strength compared to those formulated using longer chains. Ready access to diisocyanates of varying chain length affords the ability to tailor TPU properties through careful selection of the diisocyanate, polyol, and chain extender, offering hard and soft TPUs for multiple end-use applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Noise-induced barren plateaus in variational quantum algorithms

Abstract Variational Quantum Algorithms (VQAs) may be a path to quantum advantage on Noisy Intermediate-Scale Quantum (NISQ) computers. A natural question is whether noise on NISQ devices places fundamental limitations on VQA performance. We rigorously prove a serious limitation for noisy VQAs, in that the noise causes the training landscape to have a barren plateau (i.e., vanishing gradient). Specifically, for the local Pauli noise considered, we prove that the gradient vanishes exponentially in the number of qubits n if the depth of the ansatz grows linearly with n . These noise-induced barren plateaus (NIBPs) are conceptually different from noise-free barren plateaus, which are linked to random parameter initialization. Our result is formulated for a generic ansatz that includes as special cases the Quantum Alternating Operator Ansatz and the Unitary Coupled Cluster Ansatz, among others. For the former, our numerical heuristics demonstrate the NIBP phenomenon for a realistic hardware noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Probabilistic simulation of quantum circuits using a deep-learning architecture

The fundamental question of how to best simulate quantum systems using conventional computational resources lies at the forefront of condensed matter and quantum computation. It impacts both our understanding of quantum materials and our ability to emulate quantum circuits. Here we present an exact formulation of quantum dynamics via factorized generalized measurements which maps quantum states to probability distributions with the advantage that local unitary dynamics and quantum channels map to local quasistochastic matrices. This representation provides a general framework for using state-of-the-art probabilistic models in machine learning for the simulation of quantum many-body dynamics. Using this framework, we have developed a practical algorithm to simulate quantum circuits using an attention network based on a powerful neural network ansatz responsible for the most recent breakthroughs in natural language processing. We demonstrate our approach by simulating circuits that build Greenberger-Horne-Zeilinger and linear graph states of up to 60 qubits, as well as a variational quantum eigensolver circuit for preparing the ground state of the transverse field Ising model on several system sizes. Our methodology constitutes a modern machine learning approach to the simulation of quantum physics with applicability both to quantum circuits as well as other quantum many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mathematical modeling and parameter estimation of MFI membranes for para/ortho‐xylene separation

Abstract A mathematical model of a membrane with a thin, oriented, and selective MFI layer, which also includes contributions from defects, pore blockages, support layer, and external mass transfer, was formulated based on the Maxwell‐Stefan equations. It was validated using reported (Kim et al., Angew Chemie Int Ed, 2018, 57:480–485; Jeon et al., Nature, 2017, 543:690–694) para/ortho‐xylene separation data from five MFI membranes. The diffusivities of the xylenes were considered to be the same for all membranes, while the thickness and the density of defects and pore blockage were treated as unique properties of each membrane, and their contributions were estimated by fitting the model to the corresponding separation data. The effects of these properties and the role of permeate pressure on the separation performance were subsequently analyzed. The proposed modeling, parameter estimation, and analysis framework allow one to quantitatively interpret the variation of separation performance, to understand separation bottlenecks, and to provide guidance for designing membranes with desired performance.

Liu, Jingjun↗
DOE OSTI2026

Extrusion‐Spheronization of Energetic Materials

The prevailing method to produce plastic‐bonded explosive (PBX) molding powder, or “prills”, is a complex, multiphase, and bespoke process that was developed by the high explosives (HEs) manufacturing industry several decades ago. This work demonstrates the utility of a simpler, widely‐used mechanical process—extrusion‐spheronization—to produce PBX prills. We begin by detailing precautions taken to enable safe remote operation of extrusion‐spheronization equipment with HE. We then perform a study investigating the effect of lacquer solvent composition on the particle packing, pressed density, and compressive strength properties of a 95 wt.% TATB/5 wt.% polymer binder formulation akin to PBX 9502. It was found that increased composition of low vapor pressure solvents caused prolonged retention of the solvent, resulting in tackier materials that would agglomerate and form larger prills. The larger prills also led to lower poured density and tapped density of HE prills and compressive strength of pressed PBX articles. The samples prepared with a 75% propyl acetate/25% butyl acetate lacquer solvent composition exhibited the highest compressive strength. However, it is believed that the prill packing and compressive strength properties are primarily driven by the prill size rather than the chemical composition of the lacquer itself. Extrusion‐spheronization remains a promising method to reliably and repeatably produce HE prills that is less sensitive to feedstock or process variation than traditional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗