Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICS - MATRICES”

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 307 records · Page 17

Acid–base concentration swing for direct air capture of carbon dioxide

This work demonstrates the first experimental evidence of the acid–base concentration swing (ABCS) for direct air capture of CO 2 . This process is based on the effect that concentrating particular acid–base chemical reactants will strongly acidify solution, through Le Chatelier's principle, and result in outgassing absorbed CO 2 . After collecting the outgassed CO 2 , diluting the solution will result in a reversal of the acid–base reaction, basifying the solution and allowing for atmospheric CO 2 absorption. The experimental study examines a system that includes sodium cation as the alkalinity carrier, boric acid, and a polyol complexing agent that reversibly reacts with boric acid to strongly acidify solution upon concentration. Though the tested experimental system faces absorption rate and water capacity limitations, the ABCS process described here provides a basis for further process optimization. A generalized theoretical ABCS reaction framework is developed and different reaction orders and conditions are studied mathematically. Higher order reactions yield favorable cycle output results, reaching volumetric cycle capacity above 50 mM for third-order and 80 mM for fourth-order reactions. Optimal equilibrium constants are determined in order to guide alternative chemical searches and synthetic chemistry design targets. There is a substantial energetic benefit for reaction orders above the first, with second- and third-order ABCS cycles exhibiting a thermodynamic minimum work for the concentrating and outgassing steps around 150 kJ per mole of CO 2 . A significant advantage of the ABCS is that it can be driven through well-developed and widely-deployed desalination technologies, such as reverse osmosis, with opportunities for energy recovery when recombining the concentrated and diluted streams, and extraction can occur directly from the liquid phase upon vacuum application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Enceladus dust plume based on in situ measurements performed with the Cassini Cosmic Dust Analyzer

We analyzed data recorded by the Cosmic Dust Analyzer on board the Cassini spacecraft during Enceladus dust plume traversals. Our focus was on profiles of relative abundances of grains of different compositional types derived from mass spectra recorded with the Dust Analyzer subsystem during the Cassini flybys E5 and E17. The E5 profile, corresponding to a steep and fast traversal of the plume, has already been analyzed. In this paper, we included a second profile from the E17 flyby involving a nearly horizontal traversal of the south polar terrain at a significantly lower velocity. Additionally, we incorporated dust detection rates from the High Rate Detector subsystem during flybys E7 and E21. We derived grain size ranges in the different observational data sets and used these data to constrain parameters for a new dust plume model. This model was constructed using a mathematical description of dust ejection implemented in the software package DUDI. Further constraints included published velocities of gas ejection, positions of gas and dust jets, and the mass production rate of the plume. Our model employs two different types of sources: diffuse sources of dust ejected with a lower velocity and jets with a faster and more colimated emission. From our model, we derived dust mass production rates for different compositional grain types, amounting to at least 28 kg s –1 . Previously, salt-rich dust was believed to dominate the plume mass based on E5 data alone. The E17 profile shows a dominance of organic-enriched grains over the south polar terrain, a region not well constrained by E5 data. By including both E5 and E17 profiles, we find the salt-rich dust contribution to be at most 1% by mass. This revision also results from an improved understanding of grain masses of various compositional types that implies smaller sizes for salt-rich grains. Our new model can predict grain numbers and masses for future mission detectors during plume traversals.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear causality and strong hyperbolicity of baryon-rich Israel-Stewart hydrodynamics

We present the first set of fully-nonlinear, necessary and sufficient conditions guaranteeing causal evolution of the initial data for the Israel-Stewart hydrodynamic equations with shear and bulk viscosity coupled to a nonzero baryon current. These constraints not only provide nonlinear causality: they also (a) guarantee the existence of a locally well-posed evolution of the initial data (they enforce strong hyperbolicity) when excluding the endpoints of the bounds, (b) arise from purely algebraic constraints that make no underlying symmetry assumptions on the degrees of freedom and (c) propagate the relevant symmetries of the degrees of freedom over the entire evolution of the problem. Our work enforces a mathematically rigorous foundation for future studies of viscous relativistic hydrodynamics with baryon-rich matter including neutron star mergers and heavy-ion collisions.

Cordeiro, Ian [University of Illinois Urbana-Champ↗

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗

Investigating the Impact of Temporal and Directional Traffic Distribution on Crash Frequencies

Safety Performance Functions (SPFs) are mathematical models that establish relationships between the frequency of various crash types and site-specific characteristics, serving as essential tools for traffic safety analysis and roadway design. Traditional SPFs, however, often overlook the temporal fluctuations in traffic flow (such as peak-hour surges) and directional imbalances between opposing traffic streams. These traffic patterns can exacerbate congestion, disrupt driver behavior, and create unexpected conflict points, potentially leading to increased crash frequencies and more severe accidents. In light of this gap, this study aims to explore the potential of incorporating K-factors (representing peak-hour traffic proportions) and D-factors (reflecting the imbalance of directional traffic) into the development of SPFs to assess whether these factors can effectively represent the impact of temporal and spatial traffic distribution on roadway safety. Using crash data from Pennsylvania urban-suburban collector roadways, it is found that the D-factor plays a significant role in predicting the frequency of total crashes, fatal + injury crashes, and angle crashes, with positive coefficient signs indicating that higher directional imbalances correspond to increased crash risks. Similarly, the K-factor emerges as a critical predictor for fatal + injury crashes and rear-end crashes, with negative coefficients suggesting that a more pronounced traffic peak is associated with a reduction in expected crash frequencies. These results highlight the importance of accounting for uneven traffic distribution in both time and direction when developing SPFs, offering deeper insights into crash patterns and supporting more effective safety interventions and roadway designs.

Xu, Guanhao [ORNL] (ORCID:0000000214326357)↗

Deconvolution of dynamic heterogeneity in protein structure

Heterogeneity is intrinsic to the dynamic process of a chemical reaction. As reactants are converted to products via intermediates, the nature and extent of heterogeneity vary temporally throughout the duration of the reaction and spatially across the molecular ensemble. The goal of many biophysical techniques, including crystallography and spectroscopy, is to establish a reaction trajectory that follows an experimentally provoked dynamic process. It is essential to properly analyze and resolve heterogeneity inevitably embedded in experimental datasets. We have developed a deconvolution technique based on singular value decomposition (SVD), which we have rigorously practiced in diverse research projects. In this review, we recapitulate the motivation and challenges in addressing the heterogeneity problem and lay out the mathematical foundation of our methodology that enables isolation of chemically sensible structural signals. We also present a few case studies to demonstrate the concept and outcome of the SVD-based deconvolution. Finally, we highlight a few recent studies with mechanistic insights made possible by heterogeneity deconvolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms

Non-thermal advanced fuel fusion trades the requirement of a large amount of recirculating tritium in the system for that of large recirculating power. Phase space engineering technologies utilizing externally injected electromagnetic fields can be applied to meet the challenge of maintaining non-thermal particle distributions at a reasonable cost. The physical processes of the phase space engineering are studied from a theoretical and algorithmic perspective. It is emphasized that the operational space of phase space engineering is limited by the underpinning symplectic dynamics of charged particles. The phase space incompressibility according to the Liouville theorem is just one of many constraints, and Gromov's non-squeezing theorem determines the minimum footprint of the charged particles on every conjugate phase space plane. In this sense and level of sophistication, the mathematical abstraction of phase space engineering is symplectic topology. To simulate the processes of phase space engineering, such as the Maxwell demon and electromagnetic energy extraction, and to accurately calculate the minimum footprints of charged particles, recently developed structure-preserving geometric algorithms can be used. The family of algorithms conserves exactly, on discretized spacetime, symplecticity and thus incompressibility, non-squeezability, and symplectic capacities. The algorithms apply to the dynamics of charged particles under the influence of external electromagnetic fields as well as the charged particle–electromagnetic field system governed by the Vlasov–Maxwell equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the statistical theory of self-gravitating collisionless dark matter flow: High order kinematic and dynamic relations

Dark matter, if it exists, accounts for five times as much as ordinary baryonic matter. To better understand the self-gravitating collisionless dark matter flow on different scales, a statistical theory involving kinematic and dynamic relations must be developed for different types of flow, e.g., incompressible, constant divergence, and irrotational flow. This is mathematically challenging because of the intrinsic complexity of dark matter flow and the lack of a self-closed description of flow velocity. Here, this paper extends our previous work on second-order statistics Xu to kinematic relations of any order for any type of flow. Dynamic relations were also developed to relate statistical measures of different orders. The results were validated by N-body simulations. On large scales, we found that (i) third-order velocity correlations can be related to density correlation or pairwise velocity; (ii) the pth-order velocity correlations follow ∝ a (p+2)/2 for odd p and ∝ a p/2 for even p, where a is the scale factor; (iii) the overdensity δ is proportional to density correlation on the same scale, $\langle$δ$\rangle$∝$\langle$δδ'$\rangle$; (iv) velocity dispersion on a given scale r is proportional to the overdensity on the same scale. On small scales, (i) a self-closed velocity evolution is developed by decomposing the velocity into motion in haloes and motion of haloes; (ii) the evolution of vorticity and enstrophy are derived from the evolution of velocity; (iii) dynamic relations are derived to relate second- and third-order correlations; (iv) while the first moment of pairwise velocity follows $\langle$Δu L $\rangle$=-Har (H is the Hubble parameter), the third moment follows $\langle$(Δu L ) 3 $\rangle$ ∝ ε u ar that can be directly compared with simulations and observations, where ε u ≈ 10 -7 m 2 /s 3 is the constant rate for energy cascade; (v) the pth order velocity correlations follow ∝ a (3p-5)/4 for odd p and ∝ a 3p/4 for even p. Finally, the combined kinematic and dynamic relations lead to exponential and one-fourth power-law velocity correlations on large and small scales, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Current potential patches

A novel form of the current potential, a mathematical tool for the design of stellarators and stellarator coils, is developed. Specifically, these are current potentials with a finite-element basis, called current potential patches. Current potential patches leverage the relationship between distributions of magnetic dipoles and current potentials to explore limits of the access properties of stellarator coil sets. An example calculation is shown using the Helically Symmetric Experiment (HSX) equilibrium, demonstrating the method's use in coil design and understanding the limits of the access properties of coil sets. Current potential patches have additional desirable properties such as of promoting sparse current sheet solutions and identifying crucial locations of shaping current placement. A result is found for the HSX equilibrium that shaping currents covering only 22% of the winding surface are sufficient to produce the equilibrium to a good accuracy, provided a toroidal field is generated by an exterior coil set.

Elder, Todd↗

Measurement of Thomson-scattering spectra with continuous angular resolution (invited)

Here, a novel Thomson-scattering diagnostic with continuous angular resolution over a span of 120° was developed for the characterization of plasmas produced at the Omega Laser Facility. Spectrally resolving light scattered from electron plasma wave features as a function of emission angle provides a means to efficiently probe a large range of plasma frequencies and k vectors. Together, these spectra contain critical constraints on the plasma-physics models used to interpret the data and enable experimental measurements of the electron-velocity distribution function over several orders of magnitude without assumptions about its mathematical form. Major components of the instrument include (1) a reflective collection objective that gathers light over a range of 120° × 12°; (2) a spatial-filter image relay for measurement localization; (3) cylindrical optics for producing a line image of the collection aperture; (4) a transmission grating spectrometer; and (5) a time-gated, image-intensified camera. Thomson-scattered light collected from an ~(50 – μM) 3 volume of plasma is recorded with 0.8-nm spectral and 1° angular resolution. Initial experiments examined the properties of the electron-velocity distribution in gas-jet-produced plasmas in the presence of heating via inverse bremsstrahlung absorption.

47 OTHER INSTRUMENTATION↗

A review of fluids under nanoconfinement: Reactivity, geomechanics, phase transitions, and flow

Due to its many industrial applications, the physics associated with the nanoconfinement of fluids has garnered great attention from the scientific community in roughly the last fifteen years. One pioneer of that research has been the Multi-Scale Fluid-Solid Interactions in Architected and Natural Materials (MUSE) Energy Frontier Research Center, which produced more than 50 journal articles from 2018 to 2023. We structure the findings of MUSE into four categories and contextualize the work with the broader literature. Here, the four categories are (1) the dependence of chemical reactions on the degree of confinement, which given the lower freedom of motion of molecules yet improved availability of reaction sites generally lead to enhancement from bulk reactions; (2) the mechanical properties of the porous matrix, which are dependent on the heterogeneity of the system in terms of both composition and structure; (3) the phase transitions under nanoconfinement, which are influenced by the pore size and morphology, surface chemistry, and substrate composition; and (4) the physical controls of slip flow, which demonstrate the importance of several factors, such as wettability. These aspects of nanoconfinement are becoming clearer, but there are other elements such as mathematical theory that are still being developed.

58 GEOSCIENCES↗

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Near-ideal relaxed MHD in slab geometry

We investigate the solutions of the relaxed magnetohydrodynamic (MHD) model (RxMHD) of R. Dewar and Z. Qu. This model generalizes Taylor relaxation by including the ideal Ohm's law constraint using an augmented Lagrangian method, providing a pathway to extend the multi-region relaxed MHD (MRxMHD) model. We present the first numerical solution of the RxMHD model by Dewar and Qu, demonstrating that it is mathematically well-defined and computationally feasible for constructing MHD equilibria in slab geometry. We also show that a cross-field flow can exist without enforcing an arbitrary constraint on the angular momentum, as is done in the case of MRxMHD with flow. Our results also demonstrate the self-organization of fully relaxed regions during the optimization, which was an important motivation behind developing this model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scattering theory in noncanonical phase space: A Drift-Kinetic collision operator for weakly collisional plasmas

After developing a scattering theory for grazing collisions in general noncanonical phase spaces, we introduce a guiding center collision operator in five-dimensional phase space designed for plasma regimes characterized by long wavelengths (relative to the Larmor radius), low frequencies (relative to the cyclotron frequency), and weak collisionality (where repeated Coulomb collisions induce cumulatively small changes in particle magnetic moment). The collision operator is fully determined by the noncanonical Hamiltonian structure of guiding center dynamics and exhibits a metriplectic structure, ensuring the conservation of particle number, momentum, energy, and interior Casimir invariants. It also satisfies an H-theorem, allowing for deviations from an equilibrium Maxwellian distribution due to the nontrivial kernel of the noncanonical guiding center Poisson tensor, spanned by the magnetic moment. We propose that this collision operator and its underlying mathematical structure may offer valuable insight into the study of turbulence, transport, and self-organizing phenomena in both laboratory and astrophysical plasmas.

Hamiltonian mechanics↗

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗

Engineering tertiary chirality in helical biopolymers

Tertiary chirality describes the handedness of supramolecular assemblies and relies not only on the primary and secondary structures of the building blocks but also on topological driving forces that have been sparsely characterized. Helical biopolymers, especially DNA, have been extensively investigated as they possess intrinsic chirality that determines the optical, mechanical, and physical properties of the ensuing material. Here, we employ the DNA tensegrity triangle as a model system to locate the tipping points in chirality inversion at the tertiary level by X-ray diffraction. We engineer tensegrity triangle crystals with incremental rotational steps between immobile junctions from 3 to 28 base pairs (bp). We construct a mathematical model that accurately predicts and explains the molecular configurations in both this work and previous studies. Our design framework is extendable to other supramolecular assemblies of helical biopolymers and can be used in the design of chiral nanomaterials, optically active molecules, and mesoporous frameworks, all of which are of interest to physical, biological, and chemical nanoscience.

Science & Technology - Other Topics↗

Glyceraldehyde-3-phosphate dehydrogenase homologs as bifunctional gatekeepers of metabolic segregation in Pseudomonas putida

Metabolically versatile Pseudomonas species can assimilate various glycolytic and gluconeogenic substrates. Simultaneous assimilation is known to segregate carbons from each substrate type into different metabolic pathways. However, the mechanisms of this metabolic segregation remain unresolved. Here, we investigate Pseudomonas putida KT2440 during processing of the sugar glucose through glycolysis versus the phenolic acid ferulate through gluconeogenesis. Metabolome profiling reveals up to twofold less tricarboxylic acid cycle metabolites but up to 10-fold higher metabolites of upper glycolysis, pentose-phosphate, and Entner–Doudoroff pathways in glucose-grown cells compared to ferulate-grown cells. After 13C-substrate switching, kinetic isotopic profiling captures rapid assimilation of new substrate carbons into initial catabolic pathways, but incorporation into downstream pathways is absent or incomplete. Proteomics identifies a 22-fold higher abundance of one homolog of glyceraldehyde-3-phosphate dehydrogenase (GAPDH, GapA) in cells fed on glucose relative to ferulate, while abundance of another homolog (GapB) remains unchanged. Growth phenotypes and quantitative metabolomics for single and double knockout mutants of these GAPDH homologs indicate only GapA involvement in glycolytic flux, which can be compensated by the Entner–Doudoroff pathway, and distinct preference of GapB with minimal role of GapA for gluconeogenic flux. Accordingly, growth of triple knockout mutant with deletion of gapA, gapB, and edd is possible only when glycolytic and gluconeogenic substrates are provided together to meet metabolic demands in a segregated fashion, but metabolic tradeoffs lead to slow growth. A mathematical, experimentally constrained, model of the GAPDH node shows that tuning of GapA and GapB concentrations enables transition between flux regimes for nutritional adaptability.

Zhou, Nanqing [Northwestern University, Evanston]↗

Closed-form expressions for unitaries of spin-adapted fermionic operators

One of the open challenges in quantum computing simulations of problems of chemical interest is the proper enforcement of spin symmetry. Efficient quantum circuits implementing unitaries generated by spin-adapted operators remain elusive, while naïve Trotterization schemes break spin symmetry. Here, in this work, we analyze the mathematical structure of spin-adapted operators and derive closed-form expressions for unitaries generated by singlet spin-adapted generalised single and double excitations. These results represent significant progress toward the economical enforcement of spin symmetry in quantum simulations.

fermionic algebra↗