Search NASA⌕ Search

SEARCH · Search NASA

Results for “MATHEMATICS - ITERATION”

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 325 records · Page 18

Feasibility of Formulating Ecosystem Biogeochemical Models From Established Physical Rules

Abstract To improve the predictive capability of ecosystem biogeochemical models (EBMs), we discuss the feasibility of formulating biogeochemical processes using physical rules that have underpinned the many successes in computational physics and chemistry. We argue that the currently popular empirically based approaches, such as multiplicative empirical response functions and the law of the minimum, will not lead to EBM formulations that can be continuously refined to incorporate improved mechanistic understanding and empirical observations of biogeochemical processes. Instead, we propose that EBM parameterizations, as a lossy data compression problem, can be better formulated using established physical rules widely used in computational physics and chemistry, and different biogeochemical processes can be more robustly integrated within a reactive‐transport framework. Through several examples, we demonstrate how mathematical representations derived from physical rules can improve understanding of relevant biogeochemical processes and enable more effective communication between modelers, observationalists, and experimentalists regarding essential questions, such as what measurements are needed to meaningfully inform models and how can models generate new process‐level hypotheses to test in empirical studies. Finally, while empirical models with more parameters are often less robust, physical rules‐based models can be more robust and show lower predictive equifinality, stemming from their enhanced consistency in representations of processes, interactions and spatial scaling.

54 ENVIRONMENTAL SCIENCES↗

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗

A Machine Learning Framework for Predicting Microphysical Properties of Ice Crystals From Cloud Particle Imagery

The microphysical properties of ice crystals are important because they significantly alter the radiative properties and spatiotemporal distributions of clouds, which in turn strongly affect Earth's climate. However, it is challenging to measure key properties of ice crystals, such as mass or morphological features. Here, we present a proof-of-concept framework for predicting three-dimensional (3D) microphysical properties of ice crystals from in situ two-dimensional (2D) imagery. First, we computationally generated synthetic ice crystals using 3D modeling software along with geometric parameters estimated from the 2021 Ice Cryo-Encapsulation Balloon (ICEBall) field campaign. Then, we used synthetic crystals to train machine learning (ML) models to predict effective density ($ρ_e$), effective surface area ($A_e$), and number of bullets ($N_b$) from synthetic rosette imagery. On unseen synthetic images, our ML models accurately predicted ice crystal properties. ResNet-18 performed best, achieving $R^2$ values of 0.99 and 0.98 for $ρ_e$ and $A_e$, respectively, and MAE of 0.10 for mathematical equation in single view tasks. Stereo view ResNet-18 further reduced RMSE by 40% for $ρ_e$ and $A_e$ and reduced MAE by 0.08 for $N_b$. This work provides a novel ML-driven framework for estimating ice microphysical properties from in situ imagery, which will allow for downstream constraints on microphysical parameterizations, such as the mass-size relationship.

Ko, J. [Columbia Univ., New York, NY (United State↗

Semi-empirical model for Henry’s law constant of noble gases in molten salts

Henry’s law constant, which describes the proportionality of dissolved gas to partial pressure of free gas in liquid–gas equilibrium systems, can also be applied to mass transport applications. In this work, we investigated an approach for determining the solubility of noble gases in a molten salt liquid utilizing the equilibrium concept of Henry’s gas constant. Henry’s gas constant is described as a mathematical function dependent on the van der Waals radius of the noble gas and the temperature of the molten salt. The alteration in Gibbs free energy encompasses contributions from both surface and volume energies. Enthalpy and entropy are deduced from these surface and volume energies in the Gibbs free energy formulation. A comparative analysis was conducted between the conventional method and our proposed model. Moreover, useful chemical properties can be determined from examination of surface and volume energies. Our findings provide an accurate and general theory of Gibbs free energy that can be validated experimentally based on the model proposed herein. This work unifies the prediction of Henry gas constant and subsequently the entropy and enthalpy calculation for noble gases in a molten salt solution to a single functional form using van der Waals radius of the gas and temperature of the system. This functional form is then used to perform a multiple regression method to find two parameters corresponding to the surface energy and volume energy. These two parameters are consistent between all combinations of noble gas and molten salt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical modeling of hepatitis C acute infection in liver-humanized mice support pre-clinical assessment of candidate viruses for controlled-human-infection studies

Designing and carrying out a controlled human infection (CHI) model for hepatitis C virus (HCV) is critical for vaccine development. However, key considerations for a CHI model protocol include understanding of the earliest viral-host kinetic events during the acute phase and susceptibility of the viral isolate under consideration for use in the CHI model to antiviral treatment before any infections in human volunteers can take place. Humanized mouse models lack adaptive immune responses but provide a unique opportunity to obtain quantitative understanding of early HCV kinetics and develop mathematical models to further understand viral and innate immune response dynamics during acute HCV infection. We show that the models reproduce the measured HCV kinetics in humanized mice, which are consistent with early acute HCV-host dynamics in immunocompetent chimpanzees. Our findings suggest that humanized mice are well-suited to support development of a CHI model. In-silico and in-vivo modeling estimates provide a starting point to characterize candidate viruses for testing in CHI model studies.

Agent-based modeling↗

Incorporating long-range dependence and fractal features in turbulence spectra

We introduce an advanced turbulence spectrum model developed from mathematical foundations from a covariance function class and empirically validated using extensive field data. This model captures the complex dynamics of long-range dependence, and fractal characteristics prevalent in riverine and atmospheric boundary layer (ABL) flows that are ignored by classical spectrum models, such as IEC (International Electrotechnical Commission) von Kármán and Kaimal model. The model delineates scaling behaviors across distinct frequency bands and offers substantial flexibility through five well-defined parameters each characterizing a distinct physical aspect of the velocity time series. A detailed procedure for obtaining each parameter from time series data is outlined. The comprehensive validations with field data from tidal currents and ABL flows substantiate the model’s fidelity in accurately replicating observed phenomena. This validation establishes the reliability of the proposed model and, when incorporated into stochastic full-field simulators such as TurbSim, demonstrates its potential to advance the predictive modeling and analysis of turbulent flows in environmental science and engineering contexts.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗

Spacetime pq theory for AC and DC electric power systems

The 50/60 Hz alternating current (AC) electric power has been the standard and most flexible energy source powering our modern societies for one and a half centuries since the war of the currents: AC versus direct current (DC). A reactive power concept that was introduced at the beginning of the AC power was very useful for circuit/system analysis, design, control, optimization, and ultimately for more efficient and stable generation, transmission, distribution, and consumption. The initial reactive power theory was based on single-phase sinusoidal AC power to capture inductive and capacitive power that yields to net-zero average power over one fundamental cycle. Soon it was expanded to non-sinusoidal AC power and finally to instantaneous three-phase AC power. However, these reactive power theories remain separate and limited to special cases and have never been consolidated and made valid to all cases. Today, more widespread adoption of power electronics and renewable energy is bringing back DC power into the electric grids. The reactive power concept has never been applied to DC power systems. There is no reactive power in DC power systems according to the existing reactive power theories. Do DC power systems really have no reactive power? Capacitors and inductors are widely used in DC just like in AC power systems. Are they not reactive power components? Why are they different from their AC counterparts? Furthermore, are batteries active or reactive power components? What about active devices like power converters (or inverters) with AC (or DC) on one side and DC (or AC) on the other? Do they generate or consume reactive power? Finally, what about AC and DC hybrid power systems? How to define reactive power in such a complex power system that has a multitude of loads, buses, and sources? Is there reactive power between any two loads, any two buses, or any two sources in a power system and what is the total reactive power in such a complex power system as a whole? As the motivation and goal of this paper to answer the above basic questions, to unify the existing AC reactive power theories and to ultimately provide theoretical and insightful guidance for system analysis, design, control, efficiency, optimization, and operation of complex power systems, a concept of spacetime (both spatial and temporal) active and reactive power (pq) theory—the spatiotemporal aspect of active and reactive power—is developed for both AC and DC power systems. The theoretical definitions and physical meanings of the spacetime reactive power will be developed, and real applications and thought experiments/cases/exercises will be explored and discussed. The developed mathematics to define the active (or real) and reactive (or imaginary) power— p and q respectively by dot (scalar) and cross (vector) products of multi-dimension spacetime vectors and time-space mapping principle/law can have some fundamental implications as well.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The design space of E(3)-equivariant atom-centred interatomic potentials

Abstract Molecular dynamics simulation is an important tool in computational materials science and chemistry, and in the past decade it has been revolutionized by machine learning. This rapid progress in machine learning interatomic potentials has produced a number of new architectures in just the past few years. Particularly notable among these are the atomic cluster expansion, which unified many of the earlier ideas around atom-density-based descriptors, and Neural Equivariant Interatomic Potentials (NequIP), a message-passing neural network with equivariant features that exhibited state-of-the-art accuracy at the time. Here we construct a mathematical framework that unifies these models: atomic cluster expansion is extended and recast as one layer of a multi-layer architecture, while the linearized version of NequIP is understood as a particular sparsification of a much larger polynomial model. Our framework also provides a practical tool for systematically probing different choices in this unified design space. An ablation study of NequIP, via a set of experiments looking at in- and out-of-domain accuracy and smooth extrapolation very far from the training data, sheds some light on which design choices are critical to achieving high accuracy. A much-simplified version of NequIP, which we call BOTnet (for body-ordered tensor network), has an interpretable architecture and maintains its accuracy on benchmark datasets.

Computer Science↗

Finding a natural rhythm

Paul J. Dauenhauer describes the mathematical basis for designing dynamic catalysts that are programmed to change with time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electrostatic approximation

The plasma weak turbulence theory is a perturbative nonlinear theory, which has been proven to be quite valid in a number of applications. However, the standard weak turbulence theory found in the literature is fully developed for highly idealized unmagnetized plasmas. As many plasmas found in nature and laboratory are immersed in a background static magnetic field, it is necessary to extend the existing discussions to include the effects of ambient magnetic field. Such a task is quite formidable, however, which has prevented fundamental and significant progresses in the subject matter. The central difficulty lies in the formulation of the complete nonlinear response functions for magnetized plasmas. The present paper derives the nonlinear susceptibilities for weakly turbulent magnetized plasmas up to the third order nonlinearity, but in doing so, a substantial reduction in mathematical complexity is achieved by the use of Bessel function addition theorem (or sum rule). The present paper also constructs the weak turbulence wave kinetic equation in a formal sense. For the sake of simplicity, however, the present paper assumes the electrostatic interaction among plasma particles. Fully electromagnetic generalization is a subject of a subsequent paper.

Physics↗

Nonlinear susceptibilities for weakly turbulent magnetized plasma: Electromagnetic formalism

This is a companion paper to the previous work [P. H. Yoon, Phys. Plasmas 31, 032309 (2024)] in which the nonlinear susceptibilities of weakly turbulent magnetized plasma are derived under a simplifying assumption of electrostatic interaction. The present paper extends the analysis to a general situation of electromagnetic interaction. The main novelty of the previous and present papers is that by employing the Bessel function addition theorem, the mathematical definitions for the susceptibilities are substantially simplified, a procedure that has not been discussed in the existing literature. In the present paper, a full set of Maxwell’s equations are considered in conjunction with the nonlinear Vlasov equation, which is solved by a perturbative method. The result is a fully general nonlinear susceptibility, given in tensorial form, which is applicable for weakly turbulent magnetized plasmas.

Materials Science↗

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↗