Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 235 records · Page 13

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of Pt-Metal Alloy Catalysts with Various Ionomers on Oxygen Reduction Reaction in Fuel Cell Application

Pt-M/C (M = Co, Ni, Mn, etc.) alloy catalysts exhibit superior oxygen reduction reaction (ORR) activity compared to pure Pt/C, leading to a high energy efficiency in hydrogen fuel cells. However, many Pt-M/C alloy catalysts were synthesized and evaluated at the lab scale in model test-bed systems like rotating disc electrodes, which don't always correlate to performance within a fuel cell system; there is a clear need to evaluate catalysts in electrodes that can be prepared at industrially relevant scales to evaluate how factors like ink formulation can greatly affect device-level of fuel cell performance. Herein, three commercial Pt-M/C alloy catalysts (two Pt-Co/C and one Pt-Ni/C) were comprehensively characterized by various techniques. The results show that the average particle sizes of the three catalysts are close to 5 nm; the atomic ratio of Pt/M is around 4; and the M was successfully embedded into Pt lattice, resulting in the positive shift of Pt 4f in XPS spectra and XRD patterns. These catalytic materials were incorporated into 9 different cathode catalyst layers (CCLs) with three kinds of ionomers (Nafion D2020, high oxygen permeability ionomer (HOPI), and Aquivion D79-25BS), and their performance in proton exchange membrane fuel cells (PEMFCs) were investigated. The results demonstrate that the Pt-Co/C catalysts possess a higher mass activity (MA) than Pt-Ni/C; the cathodes with Nafion ionomer provide the highest MA while electrodes with Aquivion ionomer showed the lowest activity, attributed to poor H+ conductivity resulting from suboptimal ionomer incorporation. Finally, these alloys were shown to exceed DOE targets for MA and H2/Air performance reported in the recent publications at beginning of life and after 90k cycle catalyst AST protocol. This study provides valuable performance benchmarks for these materials guiding future Pt-M/C catalyst design and material integration for heavy duty PEMFC applications.

08 HYDROGEN↗

GraMeR: Gra ph Me ta R einforcement learning for multi-objective influence maximization

Influence maximization (IM) is a combinatorial problem of identifying a subset of seed nodes in a network (graph), which when activated, provide a maximal spread of influence in the network for a given diffusion model and a budget for seed set size. IM has numerous applications such as viral marketing, epidemic control, sensor placement and other network-related tasks. However, its practical uses are limited due to the computational complexity of current algorithms. Recently, deep reinforcement learning has been leveraged to solve IM in order to ease the computational burden. However, there are serious limitations in current approaches, including narrow IM formulation that only consider influence via spread and ignore self-activation, low scalability to large graphs, and lack of generalizability across graph families leading to a large running time for every test network. In this work, we address these limitations through a unique approach that involves: (1) Formulating a generic IM problem as a Markov decision process that handles both intrinsic and influence activations; (2)incorporating generalizability via meta-learning across graph families. There are previous works that combine deep reinforcement learning with graph neural network, but this work solves a more realistic IM problem and incorporates generalizability across graphs via meta reinforcement learning. Extensive experiments are carried out in various standard networks to validate performance of the proposed Graph Meta Reinforcement learning (GraMeR) framework. Finally, the results indicate that GraMeR is multiple orders faster and generic than conventional approaches when applied on small to medium scale graphs.

97 MATHEMATICS AND COMPUTING↗

Jensen–Shannon divergence based novel loss functions for Bayesian neural networks

Bayesian neural networks (BNNs) are state-of-the-art machine learning methods that can naturally regularize and systematically quantify uncertainties using their stochastic parameters. Kullback–Leibler (KL) divergence-based variational inference used in BNNs suffer from unstable optimization and challenges in approximating light-tailed posteriors due to the unbounded nature of the KL divergence. To resolve these issues, we formulate a novel loss function for BNNs based on a new modification to the generalized Jensen–Shannon (JS) divergence, which is bounded. In addition, we propose a Geometric JS divergence-based loss, which is computationally efficient since it can be evaluated analytically. We found that the JS divergence-based variational inference is intractable, and hence employed a constrained optimization framework to formulate these losses. Our theoretical analysis and empirical experiments on multiple regression and classification data sets suggest that the proposed losses perform better than the KL divergence-based loss, especially when the data sets are noisy or biased. Specifically, there are approximately 5% and 8% improvements in accuracy for a noise-added CIFAR-10 dataset and a regression dataset, respectively. There is about 13% reduction in false negative predictions of a biased histopathology dataset. Additionally, we quantify and compare the uncertainty metrics for the regression and classification tasks.

97 MATHEMATICS AND COMPUTING↗

Comparative analysis of new crystal and plastic scintillators for fast and thermal neutron detection

This paper considers scintillation and physical properties of efficient organic scintillators tested as single crystals and as components of plastics with pulse shape discrimination (PSD). For the first time, single crystals of 9,9-dimethyl-2-phenyl-9H-fluorene (PhF) were grown for studies of the basic scintillation properties of this new compound. Comparison to classical organic crystals, like anthracene, trans-stilbene, and p-terphenyl, and to more recently introduced organic glass showed that the new crystal belongs to a group of the most efficient scintillators, with light output exceeding that of trans-stilbene and PSD comparable to that of p-terphenyl. Furthermore, additional studies were conducted to evaluate scintillation performance and physical properties, like hardness and dye leaching, of plastic scintillators prepared with PhF, organic glass, and liquid diisopropylnaphthalene (DIPN) that were considered as examples for potential replacement of PPO (2,5-diphenyloxazole) in current commercial PSD plastic scintillators. Comparison of the highest performing plastic scintillators prepared with these dyes shows that PhF formulations produced a record light output (LO) increase of 69 % relative to EJ-200. Similar improvements obtained with 6 Li-loaded formulations showed that future development of PSD plastics should not be limited to use of PPO but must involve the search and exploration of new efficient dyes that may lead to discovery of much brighter organic scintillators with improved physical properties required for fast and thermal neutron detection, fast neutron spectroscopy, and antineutrino detection applications.

9,9-dimethyl-2-phenyl-9H-fluorene↗

Variational neural network approach to QFT in the field basis

We present a variational neural network approach for solving quantum field theories in the field basis, focusing on the free Klein-Gordon model formulated in momentum space. While recent studies have explored neural-network-based variational methods for scalar field theory in position space, a systematic benchmark of the analytically solvable Klein-Gordon ground state—particularly in the momentum-space field basis—has been lacking. In this work, we represent the ground-state wavefunctional as a neural network defined on a discretized set of field configurations and train it by minimizing the Hamiltonian expectation value. This framework enables direct comparison to exact analytic results for a range of key observables, including the ground-state energy, two-point correlators, expectation value of the field, and the structure of the learned wavefunctional itself. Our results provide quantitative diagnostics of accuracy and establish a validated foundation for extending neural-network wavefunctional methods to interacting field theories and position-space formulations.

Klein-Gordon model↗

f-Element complexes with benzyl and cyclohexyl substituted trihydroborates

Actinide complexes containing the simplest borohydrides (BH 4 ) 1- and (MeBH 3 ) 1- can exhibit remarkably highly volatility, which creates unique hazards and handling challenges, especially when making measurements on solid samples under vacuum. Here we describe efforts to prepare new actinide borohydride complexes with attenuated volatility by adding bulkier benzyl (Bn) and cyclohexyl (Cy) substituents to boron. Reactions of ThCl 4 , UI 3 (thf) 4 , and NdI 3 with the mixed alkali metal salt Li/K(BnBH 3 )(thf) n yielded Th(BnBH 3 ) 4 (thf) 2 , U(BnBH 3 ) 4 (thf) 2 , and K[Nd(BnBH 3 ) 4 ], respectively. Notable amongst these, the reaction with UI 3 (thf) 4 proceeds via oxidation of U(III) to U(IV) despite the presence of reducing borohydride ligands. Similarly, reactions of the same metal halides with four equivalents of Li(CyBH 3 )(Et 2 O) n yielded Th(CyBH 3 ) 4 , U(CyBH 3 ) 4 (thf) 2 , and [Li(Et 2 O) 3 ][Nd(CyBH 3 ) 4 ]. Single crystal X-ray diffraction studies of the M(BnBH 3 ) 4 (thf) 2 complexes with M = Th and U confirmed their formulations. Furthermore, the complexes have approximate D 2d point group symmetry and adopt bicapped hexagonal antiprismatic coordination geometries with axial thf ligands and κ 3 -BnBH 3 ligands bound in the equatorial plane. K[Nd(BnBH 3 ) 4 ] and [Li(Et 2 O) 3 ][Nd(CyBH 3 ) 4 ], which were prepared for comparison to U(III) complexes that were unsuccessfully targeted, were also structurally characterized to reveal complex anions with tetrahedral arrangements of trihydroborate ligands bound to Nd(III). Crystals obtained for Th(CyBH 3 ) 4 and U(CyBH 3 ) 4 (thf) 2 were not suitable for XRD studies, but 1 H and 11 B NMR spectra were consistent with their formulations. Collectively, these complexes represent rare examples of structurally characterized f-element trihydroborate complexes with carbon substituents other than methyl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

JetGP: A derivative enhanced Gaussian process library

Derivative enhanced Gaussian Processes (DEGPs) can significantly improve surrogate model accuracy over standard Gaussian Process (GP) formulations by incorporating derivative information. However, standard implementations scale poorly with dimension, limiting their use in high dimensional engineering problems. JetGP is a Python framework that unifies existing derivative enhanced GP methodologies into a single library and extends them to support arbitrary order derivative information. The library implements four complementary formulations: standard derivative enhanced Gaussian Processes (DEGP), directional DEGP (DDEGP), generalized directional DEGP (GDDEGP), and weighted DEGP (WDEGP). By unifying these approaches in a consistent interface with robust numerical implementations, JetGP enables practitioners to balance predictive accuracy and computational efficiency for high dimensional optimization, uncertainty quantification, and sensitivity analysis in engineering design.

Derivative enhanced Gaussian process↗

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mean-field dynamo as a quantum-like modulational instability

Presented here is a novel formulation of the mean-field dynamo as a modulational instability of magnetohydrodynamic (MHD) turbulence. This formulation, termed mean-field wave kinetics (MFWK), is based on the Weyl symbol calculus and allows describing the interaction between the mean fields (magnetic field and fluid velocity) and turbulence without requiring scale separation that is commonly assumed in the literature. The turbulence is described by the Wigner–Moyal equation for the spectrum of the two-point correlation matrix (Wigner matrix) of magnetic-field and velocity fluctuations and depicts the turbulence as an effective plasma of quantum-like particles that interact via the mean fields. Eddy–eddy interactions, which serve as ‘collisions’ in this effective plasma, are modelled within the standard minimal tau approximation to aid comparison with existing theories. Using MFWK, the non-local electromotive force is calculated for generic turbulence from first principles, modulo the limitations of MFWK. This result is then used to study, both analytically and numerically, the modulational modes of MHD turbulence, which appear as linear instabilities of the said effective quantum-like plasma of fluctuations. The standard α 2 -dynamo and other known results are reproduced as special cases. A new dynamo effect is predicted that is driven by correlations between the turbulent flow velocity and the turbulent current.

astrophysical plasmas↗

Scale-Up Studies for the Dehydration of C 4+ Alcohols into Drop-In Diesel Fuel

Herein, we demonstrate the production of liter quantities of drop-in diesel-range ethers from biomass-derived alcohols. We report scale-up resultsfor the dehydration of a mixture of C 4+ alcohols using powder and pellet zeolite Ycatalyst in continuous flow and batch reactors. The activity of zeolite Y decreaseswith the introduction of an alumina binder. Large alcohols, as well as branchedand secondary alcohols, increase the coke content over the pellet Y catalyst. Thepellet formulation had lower carbon balances and selectivity to C 10+ ethers,suggesting that the pellet formulation increases alcohol absorption and cokeproduction. Crushing the pellet Y catalyst to smaller particle sizes does notrecover the activity of zeolite Y in its powder form. Cold flow experiments showthat particle agglomeration contributes to pressure buildup in the continuous flowreactor. C 10+ ether blends can be produced in batch reactors at high alcohol conversion regimes. One liter of a diesel #2 blend wasproduced by scaling up this reaction. The final blend consists of a substantial portion of C 10+ ethers (63.7 wt %), followed by heavyunknown products (27.4 wt %). Furthermore, the final alcohol and butyl ether concentration values were 7.3 and 1.5 wt %, respectively. Theblendstock reported in this paper can satisfy diesel #2 ASTM standards for density, cloud point, flashpoint, and cetane number.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrochemical CO 2 Conversion Commercialization Pathways: A Concise Review on Experimental Frontiers and Technoeconomic Analysis

Technoeconomic analysis (TEA) studies are vital for formulating guidelines that drive the commercialization of electrochemical CO 2 reduction (eCO 2 R) technologies. In this review, we first discuss the progress in the field of eCO 2 R processes by providing current state-of-the-art metrices (e.g., faradic efficiency, current density) based on the recent heterogeneous catalysts’ discovery, electrolytes, electrolyzers configuration, and electrolysis process designs. Next, we assessed the TEA studies for a wide range of eCO 2 R final products, different modes of eCO 2 R systems/processes, and discussed their relative competitiveness with relevant commercial products. Finally, we discuss challenges and future directions essential for eCO 2 R commercialization by linking suggestions from TEA studies. We believe that this review will catalyze innovation in formulating advanced eCO 2 R strategies to meet the TEA benchmarks for the conversion of CO 2 into valuable chemicals at the industrial scale.

54 ENVIRONMENTAL SCIENCES↗

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗

Multinuclear Solid-State NMR and NMR Crystallography of Solid Forms of Creatine and Creatinine

Creatine is a performance-enhancing supplement with two widely available commercial solid forms, namely, creatine monohydrate (creatine·H 2 O) and creatine HCl, the latter of which does not have a reported crystal structure. Moreover, commercial formulations of creatine may contain creatinine, an undesired impurity phase resulting from the self-cyclization of creatine during manufacturing. Therefore, reliable methods for characterizing the different solid forms of creatine and detecting the presence of creatinine are essential. Herein, we address these challenges using 13 C, 15 N, and 35 Cl solid-state NMR (SSNMR) spectroscopy to obtain distinct spectral fingerprints for creatine·H 2 O and creatine HCl, along with creatinine and creatinine HCl. The acquisition of these SSNMR spectra offers a robust approach for both the rapid characterization of each solid form and the detection of the impurity phases. Additionally, quadrupolar NMR crystallography-guided crystal structure prediction (QNMRX-CSP) was applied for the de novo crystal structure determination of creatine HCl, which was validated by the subsequently determined single-crystal X-ray diffraction (SCXRD) structure. Finally, to investigate the relationship between NMR parameters and structural features, 13 C and 15 N chemical shifts and 35 Cl electric field gradient (EFG) tensors were computed from geometry-optimized structures of the four solid forms by using dispersion-corrected DFT-D2* methods. Finally, this integrative approach offers a powerful framework for advancing the structural understanding and quality control of creatine-based supplements and next-generation formulations, as well as a wide range of other solid pharmaceuticals and nutraceuticals.

NMR↗

Dry Printing Pure Copper with High Conductivity and Adhesion for Flexible Electronics

Additive manufacturing of functional devices on various rigid and flexible substrates is rising rapidly due to their design flexibility, rapid manufacturing, and lower cost. Current printing technologies are ink-based and focused on printing silver (Ag) as conductive lines due to its matured ink formulation process, low sintering temperature, ease of printing, and low oxidation rate. However, Ag is the 68th most abundant element on Earth, while copper (Cu) is the 25th, making it much cheaper (>100×) while having a comparable conductivity to Ag. Therefore, printing Cu has become technologically and economically more attractive than Ag. Nevertheless, Cu printing is still a significant challenge in ink-based printing methods due to the higher sintering temperature relative to the glass-transition temperature of most flexible substrates, the higher oxidation rate, the challenging ink formulation process, and ink stability concerns. Here, we demonstrate printing highly conductive Cu on flexible polyimide substrates using a dry printing technique. Cu nanoparticles (~3–30 nm) are generated by on-demand laser ablation of a solid Cu target inside the printer head and under argon background gas. These Cu nanoparticles are then transported through a nozzle and onto the substrate, where they are laser-sintered in real time. The argon gas plays three critical roles in laser plume condensation for nanoparticle generation, transport, and sheath gas to avoid oxidation during sintering. The sintered nanoparticles thus show high electrical conductivity and mechanical stability under static and cyclic tests. Our dry printing technique can potentially revolutionize how electronic devices and sensors are additively manufactured for earth and space applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗