SEARCH · Search NASA

Results for “Closed form solution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Boson star normal modes

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schrödinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schrödinger approximation for the high-energy excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Large-Area Periodic Arrays of Atomically Flat Single-Crystal Gold Nanotriangles Formed Directly on Substrate Surfaces

The advancement of nanoenabled wafer-based devices requires the establishment of core competencies related to the deterministic positioning of nanometric building blocks over large areas. Within this realm, plasmonic single-crystal gold nanotriangles represent one of the most attractive nanoscale components but where the formation of addressable arrays at scale has heretofore proven impracticable. In this work, a benchtop process is presented for the formation of large-area periodic arrays of gold nanotriangles. The devised growth pathway sees the formation of an array of defect-laden seeds using lithographic and vapor-phase assembly processes followed by their placement in a growth solution promoting planar growth and threefold symmetric side-faceting. The nanotriangles formed in this high-yield synthesis distinguish themselves in that they are epitaxially aligned with the underlying substrate, grown to thicknesses that are not readily obtainable in colloidal syntheses, and present atomically flat pristine surfaces exhibiting gold atoms with a close-packed structure. As such, they express crisp and unambiguous plasmonic modes and form photoactive surfaces with highly tunable and readily modeled plasmon resonances. The devised methods, hence, advance the integration of single-crystal gold nanotriangles into device platforms and provide an overall fabrication strategy that is adaptable to other nanomaterials.

42 ENGINEERING↗

DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Iodate interactions with calcite: implications for natural attenuation

Solid-phase interactions and speciation are important to radioiodine transport in groundwater. At the Hanford Site in Southeastern Washington State, iodate (IO 3 - ) is the main aqueous species in dilute radioiodine groundwater plumes. Like other oxyanions, IO 3 - may be incorporated into and/or adsorbed onto calcite, a common mineral at Hanford, decreasing its mobility in the environment. A series of macroscale batch experiments combined with solid phase characterization were conducted to identify variables impacting time-dependent aqueous IO 3 - removal via calcite precipitation and determine the location of IO 3 - within the calcite crystal structure. Results demonstrated 11.5-97% aqueous IO 3 - removal during initial rapid calcite precipitation. Incorporation was apparently the main removal mechanism, although later slower precipitation and/or adsorption may have also contributed to IO 3 - removal. Using a higher concentration of the calcite forming solutions (i.e., using 1M vs. 0.1M concentrations) resulted in an increase in the amount of precipitated calcite and a greater percentage of IO 3 - removed; however, calcite formed with lower molarity solutions resulted in higher IO 3 - mass (µg/g) removal. Solubility testing of laboratory produced calcites showed only small differences in solubility for calcite with and without IO 3 - incorporated into the structure. Evidence collected from SEM/FIB and TEM/SAED suggested that the IO 3 - incorporated into calcite was present in regions close to surface (implying potential easy release upon calcite dissolution).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Incremental Interval Assignment by Integer Linear Algebra with Improvements

Interval Assignment (IA) is the problem of selecting the number of mesh edges (intervals) for each curve for conforming quad and hex meshing. The intervals x is fundamentally integer-valued. Many other approaches perform numerical optimization then convert a floating-point solution into an integer solution, which is slow and error prone. We avoid such steps: we start integer, and stay integer. Incremental Interval Assignment (IIA) uses integer linear algebra (Hermite normal form) to find an initial solution to the meshing constraints, satisfying the integer matrix equation Solving for reduced row echelon form provides integer vectors spanning the nullspace of A. Here we add vectors from the nullspace to improve the initial solution, maintaining Ax = b Heuristics find good integer linear combinations of nullspace vectors that provide strict improvement towards variable bounds or goals. IIA always produces an integer solution if one exists. In practice we usually achieve solutions close to the user goals, but there is no guarantee that the solution is optimal, nor even satisfies variable bounds, e.g. has positive intervals. We describe several algorithmic changes since first publication that tend to improve the final solution. The software is freely available.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗

DOE OSTI2023

Non-Fermi Liquid Behavior in the Three-Dimensional Hubbard Model

We present a numerical study on the non-Fermi liquid behavior of a three-dimensional strongly correlated system. The Hubbard model in a simple cubic lattice is simulated by the dynamical cluster approximation; in particular, the quasi-particle weight is calculated at finite dopings for a range of temperatures. By fitting the quasi-particle weight to the marginal Fermi liquid form at finite doping near the putative quantum critical point, we find evidence of a separatrix between Fermi liquid and non-Fermi liquid regions. Our results suggest that a marginal Fermi liquid and possibly a quantum critical point exist in the non-symmetry broken solution of the three-dimensional interacting electron systems. We also calculate the spectral function, close to the half-filling, and we obtain evidence of pseudogap.

Kellar, Samuel↗

DOE OSTI2024

Synthetically Reversible, Proton-Mediated Nitrite N–O Bond Cleavage at a Dicopper Site

A monocationic dicopper(I,I) nitrite complex [Cu 2 (μ-κ 1 :κ 1 -O 2 N)DPFN][NTf 2 ] (2) (DPFN = 2,7-bis(fluoro-di(2-pyridyl)methyl)-1,8-naphthyridine, NTf 2 – = N(SO 2 CF 3 ) 2 – ), was synthesized by treatment of a dicopper acetonitrile complex, [Cu 2 (μ-MeCN)DPFN][NTf 2 ] 2 (1), with tetrabutylammonium nitrite ([nBu 4 N][NO 2 ]). DFT calculations indicate that 2 is one of three linkage isomers that are close in energy and presumably accessible in solution. Reaction of the μ-κ 1 :κ 1 -O 2 N complex with p-TolSH produces nitrous acid (HONO) and the corresponding dicopper thiolate species via an acid–base exchange reaction. Notably, treatment of 2 with HNTf 2 results in N–O bond cleavage in the putative, HONO-ligated complex to form the more thermodynamically favorable nitrosyl-bridged dicopper complex [Cu 2 (μ-NO)(μ-OH)DPFN][NTf 2 ] 2 (4). This scission can be reversed via deprotonation of the hydroxy ligand with KO t Bu. X-ray diffraction studies confirmed the solid-state molecular structures of 2 and 4. DFT calculations were used to construct a reaction coordinate diagram detailing formation of the μ-NO complex and to describe its electronic structure. The nitrosyl ligand in 4 is chemically labile, as demonstrated by its ready displacement in reactions with CO or NO 2 – .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

An Analytical Solution for the Initiation and Early Progression of Fretting Wear in Spherical Contacts

Abstract This article derives analytical solutions to calculate the wear volume at the initiation of fretting motion and its early progression over the first few oscillation cycles. The Archard-based model considers a deformable hemisphere that is contact with a deformable flat block. The material pairs investigated are special alloys, the Inconel 617/Incoloy 800H, and Inconel 617/Inconel 617. The analytical study begins with a unidirectional frictional sliding contact, where the local interfacial sliding distance and the nominal sliding distance at the initiation of gross slip are derived. The obtained analytical expressions for unidirectional sliding are then used to derive the corresponding wear volume for the initiation and early progression of gross slip and the wear volume for a general fretting cycle under elastic conditions. These analytical derivations are all verified by the finite element analysis (FEA). The FEA method and the analytical solutions render virtually identical results for both similar and dissimilar material pairs. The effects of plasticity on the wear volume under elastic–plastic conditions are also investigated. It is found that the fretting wear volumes obtained from the FEA simulations, which include plasticity, are close to those obtained from the analytical expressions for purely elastic regimes. All the results are presented in normalized forms, which can be easily generalized and applied to three-dimensional fretting wear of other material pairs.

Engineering↗

DOE OSTI2021

Solidification dynamics in metal additive manufacturing: analysis of model assumptions

Solidification characteristics in the meltpool drive the process-microstructure relationship which helps dictate the material properties of as-built parts in additive manufacturing; therefore, being able to accurately and quickly simulate the size, shape, and solidification characteristics in the melt pool is of great interest to the field. In this study, we investigate various important physical phenomena (dynamic material properties, fluid-flow, radiation and vaporization) which can either be included or neglected in a continuum finite volume model (FVM) and their effect on the solidification conditions. Additionally, since the simplest form of such a model (conduction only) has an analytic solution which is much faster, its viability is also considered. Since the inclusion of some of these physical phenomena will inherently change the net energy input as well as the amount of energy needed to achieve melting of a control volume, each set of included phenomena had an effective absorption efficiency which was calibrated to closely match the dimensions of the melt pool to that of the ground truth data. The ground truth data for this study was defined to be the output of the FVM which included all the physical phenomena (OF). This study then goes on to compare the effects on solidification conditions each of these calibrated models has. It was found that most of the change in solidification conditions comes from the inclusion of latent heat. A posterior correlation factor (PCF) is then introduced to enable an analytic model to predict similar solidification conditions to OF model.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗

DOE OSTI2020

Assembly of Disordered Cocontinuous Morphologies by Multiblock Copolymers with Random Block Sequence and Length Dispersity

Among the wide variety of nanostructures accessible by self-assembly of block copolymers, cocontinuous morphologies are of great interest for separations, catalysis, sensing, and energy storage due to their ability to integrate distinct properties of the constituent polymers. Furthermore, we investigate the tendency to form disordered cocontinuous phases of multiblock copolymers (MBCs) with random block sequence and length dispersity prepared by step-growth polymerization of pre-made telechelic polymers, specifically using thiol-Michael addition between diacrylate functionalized polystyrene (PS) and poly (D,L - lactide) (PLA) chains and a small molecule dithiol. In the phase diagram so established, disordered cocontinuous nanostructures occupy a wide window (≈ 25 % vol) at a modest segregation strength of χN = 15, close to the critical value for microphase separation. At a higher segregation strength of χN = 30, only a moderately wide (≈ 10 vol %) cocontinuous window is found, with MBCs instead tending to form ordered nanostructures. Although the breadth of the cocontinuous widow is less than that found in prior work on analogous randomly end-linked copolymer networks (RECNs), the MBC architecture offers offsetting advantages in terms of solution and melt processability.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Machine learning of superconducting critical temperature from Eliashberg theory

Abstract The Eliashberg theory of superconductivity accounts for the fundamental physics of conventional superconductors, including the retardation of the interaction and the Coulomb pseudopotential, to predict the critical temperature T c . McMillan, Allen, and Dynes derived approximate closed-form expressions for the critical temperature within this theory, which depends on the electron–phonon spectral function α 2 F ( ω ). Here we show that modern machine-learning techniques can substantially improve these formulae, accounting for more general shapes of the α 2 F function. Using symbolic regression and the SISSO framework, together with a database of artificially generated α 2 F functions and numerical solutions of the Eliashberg equations, we derive a formula for T c that performs as well as Allen–Dynes for low- T c superconductors and substantially better for higher- T c ones. This corrects the systematic underestimation of T c while reproducing the physical constraints originally outlined by Allen and Dynes. This equation should replace the Allen–Dynes formula for the prediction of higher-temperature superconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗

DOE OSTI2025

In situ monitoring of lanthanide reactions with oxide species via combined absorption spectroscopy and electrochemical methods

Molten salts for engineering scale applications of spent nuclear fuel pyrochemical processing will inevitably have some level of oxygen impurities which can form insoluble oxide and oxychloride species with fission products. This work demonstrates real-time concentration monitoring of two lanthanide (Ln3+) fission products, Nd3+ and Pr3+, and their reactions with oxygen (O2-) impurities to form insoluble products in LiCl-NaCl-KCl eutectic salt. Combined absorption spectroscopy and electrochemical testing were used to track lanthanide concentrations. O2- impurity levels were controlled by adding Li2O to lanthanide-salt solutions. After the introduction of O2- impurities, Ln3+ concentrations were monitored via time-resolved absorption spectroscopy. Concentrations of both Ln3+ species in solution decreased with time as insoluble products formed. The initial impurity concentration controlled whether insoluble products were predominantly oxychlorides or mixtures of oxychlorides and oxides. However, absorption spectroscopy is limited for weakly absorbing species, such as Pr3+, and under conditions of high impurity concentrations where solutions can be turbid. To circumvent this limitation, the concentrations of Pr3+ were monitored with cyclic voltammetry (CV). Reaction rates and extent of Pr3+ removal from solution as monitored by CV agreed closely with results found from spectroscopic monitoring. This demonstrates that simultaneous electrochemical testing complements the capabilities of absorption spectroscopy to monitor reactions of fission products in molten salts.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Silicon Anodes with Improved Calendar Life Enabled By Multivalent Additives

Silicon is widely recognized as the most promising upgrade for graphite anodes due to its much higher capacity, natural abundance, and ability to be directly applied in the slurry-based, roll-to-roll production lines. However, in addition to the fast capacity decay, silicon anodes also suffer from inferior calendar life in practical applications due to the unstable solid-electrode interface (SEI). Until now, strategies to effectively improve the calendar life by tailored SEIs remain largely unclear, especially in high-Si content, zero-graphite anodes. Here, silicon anodes with superior calendar life are developed by adding small concentrations of multivalent salts into the baseline electrolyte. The Ca additive reacts with the F ions in the electrolyte, forming a layer of nanocrystalline CaF 2 that is closely coated around the silicon particles. The CaF 2 -enabled new SEI is strong and dense, which effectively protects the silicon core from side reactions, leading to lower capacity decay after calendar aging at high voltage. More importantly, the Ca additive is effective universally for all available commercial silicon or SiO sources. This study provides a feasible and low-cost solution for developing silicon anodes with long calendar life, paving the way towards commercially viable silicon anodes.

25 ENERGY STORAGE↗

DOE OSTI2025

In situ Monitoring of Lanthanide Reactions with Oxide Species via Combined Absorption Spectroscopy and Electrochemical Methods

Molten salts for engineering scale applications of spent nuclear fuel pyrochemical processing will inevitably have some level of oxygen impurities which can form insoluble oxide and oxychloride species with fission products. This work demonstrates real-time concentration monitoring of two lanthanide (Ln3+) fission products, Nd3+ and Pr3+, and their reactions with oxygen (O2-) impurities to form insoluble products in LiCl-NaCl-KCl eutectic salt. Combined absorption spectroscopy and electrochemical testing were used to track lanthanide concentrations. O2- impurity levels were controlled by adding Li2O to lanthanide-salt solutions. After the introduction of O2- impurities, Ln3+ concentrations were monitored via time-resolved absorption spectroscopy. Concentrations of both Ln3+ species in solution decreased with time as insoluble products formed. The initial impurity concentration controlled whether insoluble products were predominantly oxychlorides or mixtures of oxychlorides and oxides. However, absorption spectroscopy is limited for weakly absorbing species, such as Pr3+, and under conditions of high impurity concentrations where solutions can be turbid. To circumvent this limitation, the concentrations of Pr3+ were monitored with square wave voltammetry (SWV). Reaction rates and extent of Pr3+ removal from solution as monitored by SWV agreed closely with results found from spectroscopic monitoring. This demonstrates that simultaneous electrochemical testing complements the capabilities of absorption spectroscopy to monitor reactions of fission products in molten salts.

absorption spectroscopy↗

DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗