Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Analytical Model”

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 811 records · Page 45

On a simplified solution of climate-carbon dynamics in idealized flat10MIP simulations

Abstract. Idealized experiments with coupled climate-carbon Earth system models (ESMs) provide a basis for understanding the response of the carbon cycle to external forcing and for quantifying climate-carbon feedbacks. Here, we analyze globally-averaged results from idealized esm-flat10 experiments and show that most models exhibit a quasi-linear relationship between cumulative carbon uptake on land and in the ocean during a period of constant fossil fuel emissions of 10 Pg C yr−1. We hypothesize that this relationship does not depend on emission pathways. Further, as a simplification, we quantify the relationship between cumulative ocean carbon uptake and changes in ocean heat content using a linear approximation. In this way, changes in oceanic heat content and atmospheric CO2 concentration become interdependent variables, reducing the coupled temperature-CO2 system to just one differential equation. The equation can be solved analytically or numerically for the atmospheric CO2 concentration as a function of fossil fuel emissions. This approach leads to a simplified description of global carbon and climate dynamics, which could be used for applications beyond existing analytical frameworks.

Brovkin, Victor

Concentration-Discharge Relationships in the Six Largest Arctic Rivers, 2003-2019

This dataset provides the results of the analysis of the relationship of dissolved analyte concentrations and river discharges in the six largest Arctic rivers across the global panarctic region (see Figure 1 in documentation file *.pdf). Long-term measurements of dissolved analyte concentrations and river discharge have been collected for each of the Kolyma, Lena, Mackenzie, Ob, Yenisey, and Yukon rivers by the Arctic Great Rivers Observatory (ArcticGRO) project from ~2003-present (Shiklomanov, 2021). The relationship of dissolved analyte concentrations and discharges in each river was characterized by statistical analysis of the slope of the log(concentration) vs log(discharge) (b), the coefficient of variation ratio (CVc/CVq), the 2.5% and 97.5% confidence intervals of b, and assigning a chemostatic, flushing, diluting, or non-systematic behavior category according to Koger (2018). The summary of these analyses for all six rivers is provided in one .csv file. The concentrations of 20 dissolved analytes and discharge measurement data for the individual Kolyma, Lena, Mackenzie, Ob, Yenisey, and Yukon rivers are also provided with this dataset. There are seven *.csv files; one for each river plus the statistical summary. These public ArcticGRO data at "https://www.arcticgreatrivers.org" (Shiklomanov, 2021) were downloaded on Feb 13, 2020, but each river has different measurement dates over the sampling and analysis period. The ArcticGRO metadata document (*.pdf) downloaded on Feb 13, 2020 is also included in this dataset. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES

Plutonium Oxidation State Distribution under WIPP Relevant Conditions (Rev. 2)

The oxidation state of plutonium in the Waste Isolation Pilot Plant (WIPP) environment has been a topic of interest since the initial compliance certification application. Plutonium (Pu) was initially expected to be present primarily as Pu(IV), but since the presence of Pu(III) could not be ruled out, it was also included in performance assessment (PA) calculations. The redox of the other variant actinides, uranium (U) and neptunium (Np), are also coupled to the plutonium redox. This led to a position in performance assessment calculations that is commonly referred to as “50/50”. In which half of the PA realizations assume solubility is dominated by the low oxidation state [U(IV), Np(IV), Pu(III)] and the other half assume solubility is dominated by the higher oxidation state [U(VI), Np(V), Pu(IV)]. Recently, a joint group of WIPP project participants from Los Alamos National Laboratory (LANL), Sandia National Laboratories (SNL), and the Department of Energy’s Carlsbad Field Office (DOE-CBFO) agreed on a new path forward for the oxidation state model in PA which will decouple the redox active actinides and will instead consider their redox as a function of E h . This will take the current binary approach and turn it into a system with four possible redox states.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Accelerating Combustion and Surface Chemistry Simulations

Design of modern combustion systems relies on computer models to predict how changes in design will affect performance. These models have largely displaced previous methods that rely on the designer’s intuition or costly and time-consuming physical testing. By using improved models, design cycles can be shortened, and cleaner and more efficient combustion devices can be created. This project aims to improve computer simulations of transportation fuels with the goal of making these simulations faster and more accurate for predicting combustion in vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Many-body interference in kagome crystals

When electrons in metals act collectively, they enable emergent phenomena and electronic functionalities that transcend the behaviour of individual particles. Coherent collective charge motion has so far been observed primarily in superconductors, in which it arises with the formation of Cooper pairs. Here we report experimental evidence for coherent charge transport in the normal state of the kagome metal CsV 3 Sb 5 , indicative of a distinct collective electronic state. The signature is a set of magnetoresistance oscillations in mesoscopic crystalline pillars under in-plane magnetic fields, with a periodicity determined by the number of magnetic flux quanta h/e threading between adjacent kagome layers—effectively forming an interlayer Aharonov–Bohm interferometer. The cooperative nature of this phenomenon is evidenced by a non-analytic angular dependence characterized by abrupt transitions between discrete oscillation frequencies and its persistence over length scales that exceed the single-particle mean free path. Notably, the oscillation amplitude matches other anomalous electronic responses reported in CsV 3 Sb 5 , pointing to an underlying mechanism that establishes intrinsic coherence. These findings shed new light on the debated nature of correlated order in kagome metals and establish CsV 3 Sb 5 as a platform for realizing long-range coherent charge transport in the absence of superconductivity—opening new directions for coherence in correlated electron systems beyond conventional models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Revealing systematic changes in the transcriptome during the transition from exponential growth to stationary phase

ABSTRACT The composition of bacterial transcriptomes is determined by the transcriptional regulatory network (TRN). The TRN regulates the transition from one physiological state to another. Here, we use independent component analysis to monitor the composition of the transcriptome during the transition from the exponential growth phase to the stationary phase. With Escherichia coli K-12 MG1655 as a model strain, we trigger the transition using carbon, nitrogen, and sulfur starvation. We find that (i) the transition to the stationary phase accompanies common transcriptome changes, including increased stringent responses and reduced production of cellular building blocks and energy regardless of the limiting element; (ii) condition-specific changes are strongly associated with transcriptional regulators ( e.g. , Crp, NtrC, CysB, Cbl) responsible for metabolizing the limiting element; and (iii) the shortage of each limiting element differentially affects the production of amino acids and extracellular polymers. This study demonstrates how the combination of genome-scale datasets and new data analytics reveals the fundamental characteristics of a key transition in the life cycle of bacteria. IMPORTANCE Nutrient limitations are critical environmental perturbations in bacterial physiology. Despite its importance, a detailed understanding of how bacterial transcriptomes are adjusted has been limited. By utilizing independent component analysis (ICA) to decompose transcriptome data, this study reveals key regulatory events that enable bacteria to adapt to nutrient limitations. The findings not only highlight common responses, such as the stringent response, but also condition-specific regulatory shifts associated with carbon, nitrogen, and sulfur starvation. The insights gained from this work advance our knowledge of bacterial physiology, gene regulation, and metabolic adaptation.

Lim, Hyun Gyu (ORCID:0000000204692388)

Verification of the REBUS Software

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. REBUS is central to this system. The driver for this effort requires the Triangular-Z and hexagonal-Z core geometry options of REBUS to be verified. Previous work identified the REBUS features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying REBUS’s ability to correctly intepret the user input model, verifying that the features identified yield the intended results, and verifying the correctness of the REBUS output tables. The REBUS software verification relies heavily upon the accuracy of the embedded DIF3D software, the verification of which was completed and documented elsewhere. Given that DIF3D produces an accurate solution, the primary focus of the verification in the REBUS software is to ensure that it properly uses the DIF3D solution and that the depletion system (Bateman equations) are correctly implemented. This manuscript reiterates the verification tasks and displays results with respect to the features needed for current design activities. Analytic solutions of the Batemen equations are displayed and the results calculated with REBUS are displayed demonstrating the accuracy. Since coupled Bateman and neutron diffusion/transport solutions are extremely difficult to obtain, much of the focus is placed on how REBUS uses a given DIF3D solution assuming the accuracy of the DIF3D solution. The verification effort identified no issues that are debilitating or otherwise impactful to the design usage of REBUS, and thus REBUS version 11.0, release 3012 is considered verified. It is important to note that several outputs of REBUS are identified to be inaccurate, such as burnup in MWD/MT. Most of the relevant ones for VTR are generally accurate with 10-20% errors which is not impactful as all regular REBUS users are aware of this issue and know how to hand calculate the results. The REBUS manual further makes it clear that these values are consistent with the methodology being used by REBUS and thus the “errors” are more of an inconsistent definition with respect to what a user would expect given a definition in literature. Other issues that were identified included unclear documentation and software bugs all of which were inconsequential to the final results.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Monte Carlo Simulations of Polarized Radiative Transfer in Neutron Star Atmospheres

Soft X-ray emission from neutron stars affords powerful diagnostic tools for uncovering their surface and interior properties, as well as their geometric configurations. In the atmospheres of neutron stars, the presence of magnetic fields alters the photon-electron scattering cross sections, resulting in nontrivial angular dependence of intensity and polarization of the emergent signals. This paper presents recent developments of our Monte Carlo simulation, MAGTHOMSCATT, which tracks the complex electric field vector for each photon during its transport. Our analysis encompasses the anisotropy and polarization characteristics of X-ray emission for field strengths ranging from nonmagnetic to extremely magnetized regimes that are germane to magnetars. In the very low field domain, we reproduced the numerical solution to the radiative transfer equation for nonmagnetic Thomson scattering, and provided analytical fits for the angular dependence of the intensity and the polarization degree. These fits can be useful for studies of millisecond pulsars and magnetic white dwarfs. By implementing a refined injection protocol, we show that, in the magnetar regime, the simulated intensity and polarization pulse profiles of emission from extended surface regions becomes invariant with respect to the ratio of photon (ω) and electron cyclotron (ω B ) frequencies once ω/ω B ≲ 0.01. This circumvents the need for simulations pertinent to really high magnetic field strengths, which are inherently slower. Our approach will be employed elsewhere to model observational data to constrain neutron star geometric parameters and properties of emitting hot spots on their surfaces.

79 ASTRONOMY AND ASTROPHYSICS

InterGraph-CPS: A Graph-Theoretic Approach to Characterize Cross-Domain Cyber-Physical Interdependencies and Uncertainties in Electric Grid Systems for Improved Decision-Making in Operation and Response

Critical infrastructure systems such as the electric grid are increasingly cyber-physical; yet, despite the cyber-physical characteristics of critical infrastructure systems, the physical process system and communication/control network system are traditionally analyzed in siloes. As these systems become more cyber-physical, it is crucial that models and methods are available to assess the cyber physical system (CPS) interdependencies, characteristics, and event propagation for improved planning, operation, and response. Thus, we proposed an integrated structural and temporal CPS interdependency analysis framework, InterGraph-CPS, that provides insight into the CPS function during normal operation as well as disturbances. This integrated structural and temporal interdependency framework is uniquely designed for assessing CPSs by account for the challenges of analyzing cyber and physical data streams together due to data availability, data type, and time scale differences. By leveraging both structural (e.g., graph analysis) and temporal (e.g., data analytics) techniques, different CPS behaviors and configurations can be accounted for.

24 POWER TRANSMISSION AND DISTRIBUTION

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Data driven investigation to understand the influence of total solids on biological biogas upgrading

In situ biogas upgrading achieves CO 2 conversion to CH 4 via hydrogenotrophic methanogenesis; however, gas-liquid mass transfer constraints limit the upgrading performance. Recognizing that optimization studies often underrepresent the effects of total solids (TS) and organic loading rate (OLR), this study undertook a holistic, statistics driven assessment of operating conditions for in situ H 2 assisted biogas upgrading, centering the analysis on TS and OLR. A dataset of 31 studies was compiled and comprised 99 observations. A rigorous analytical framework was employed, combining data standardization, fixed- and random-effects (REML) weighted regressions with cluster-robust errors, stratified analyses, and machine learning. Mixed-effects meta regression indicated that TS was the main factor explaining differences of methane fraction (CH 4 %) when considering the between studies heterogeneity. Focusing on a near-stoichiometric subset (H 2 /CO 2 ≈ 4:1), TS remained significant. Stratified results showed a stronger negative relationship between TS and CH 4 % in UASB reactors than in CSTRs, with a negative effect under mesophilic conditions and no significant effect under thermophilic conditions. A Random Forest model corroborated the statistical findings, consistently ranking H 2 /CO 2 ratio, OLR, TS, and hydrogen injection rate (HIR) as the most influential predictors. These findings delineate trends across increasing TS levels, particularly between 1% and 10%, and provide preliminary insights for TS above 15% in in situ biogas upgrading. They further provide insights for the influence of TS by reactor type and temperature, thereby advancing the evidence base for implementing biological CO 2 conversion to CH 4 in practice.

In situ biogas upgrading

A New Framework for the Attribution of Air‐Sea CO 2 Exchange

Abstract The air‐sea transfer of carbon dioxide can be viewed as a dynamic system through which atmospheric and oceanic processes push surface waters away from thermodynamic equilibrium, while diffusive gas transfer pulls them back toward local equilibrium. These push/pull processes drive significant sub‐seasonal, seasonal, and interannual variability in air‐sea carbon fluxes, the quantification of which is critical both for diagnosing the ocean response to fossil fuel emissions and for attempts to mitigate anthropogenic climate disruption through intentional modification of surface ocean biogeochemistry. In this study, we present a new approach for attributing air‐sea carbon fluxes to specific mechanisms. The new framework is first applied to a two‐box ocean nutrient and carbon cycle model as an illustrative example. Next, outputs from a regional eddy‐resolving model of the Southern Ocean are analyzed. The roles of multiple physical and biogeochemical processes are identified. The decomposition of the seasonal air‐sea carbon flux shows the dominant role of biological carbon pumps that are partially compensated by the transport convergence. Finally, the framework is used to diagnose the response to mesoscale iron and alkalinity release, explicitly quantifying transport feedback and eventual impacts on net air‐sea carbon flux. Ocean carbon transport has divergent influences between iron and alkalinity release, due to opposing near‐surface gradients of dissolved inorganic carbon. We suggest that our attribution framework may be a useful analytical technique for monitoring natural ocean carbon fluxes and quantifying the impacts of human intervention on the ocean carbon cycle.

Ito, Takamitsu [School of Earth and Atmospheric Sc

Positron stopping in multilayer materials

Positron annihilation spectroscopy provides a sensitive toolset for defect characterization. In beam based studies of single-layer targets, the form of implantation profiles is well established, depending on the kinetic energy and angle of incident positrons relative to the target surface and the density and average atomic number of the target. For multilayer systems, the difference in density and across the layers makes derivation of an analytical form difficult. To date, the determination of positron stopping profiles in multilayer targets has primarily involved Monte Carlo simulations. We present here an alternative approach that estimates the energy distribution dN/dE of those positrons transmitted past each layer boundary, by fitting the remaining tail of the stopping profile after each layer with a basis set comprised of calculated stopping profiles in the same material they are transmitted through. The stopping profile in the next layer is then found by summing a series of stopping profiles in the new medium in proportion to the determined distribution dN/dE. The results of our model are compared with simulation results in a system of alternating layers of Al and Au and find reasonable agreement in the predicted profile and excellent agreement in the predicted mean implantation depth. Finally, we derived a simple formula-based approach for the calculation of the mean implantation depth in two-layer systems that provides results in excellent agreement with the full model.

Jones, A. C. L. [Arizona State University, Tempe,

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation induced reactions is therefore key to innovating and optimizing next generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct dissolution–based reprocessing strategies. We will explore time resolved electron pulse radiolysis and alpha and gamma dose accumulation studies, integrated with multiscale computational modeling, to elucidate the molecular level roles of radiation driven, non equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity

Anion–Cation–Anion Ion Triplet Characterization by Computation and Photoelectron Spectroscopy

Ion triplets of the chloride salts of two commonly used weakly coordinating cations are reported (i.e., Cl – ·NMe 4 + Cl – (1 – ) and Cl – ·PPh 4 + Cl – (2 – )). Negative ion photoelectron spectra at 20 K afford vertical and adiabatic detachment energies of 5.18 and 5.00 eV (1 – ) and 5.03 and 4.70 eV (2 – ), respectively. These results are well reproduced by coupled cluster calculations with single, double, and perturbative triple excitations (CCSD(T)) whereas M06-2X is systematically too small by ~0.3 eV (i.e., 7 kcal mol –1 ). The structures of both 1 – and 2 – have five or six C–H···Cl – interactions that stabilize these cluster anions by 32 (1 – ) and 27 (2 – ) kcal mol –1 as given by their chloride dissociation enthalpies. In conclusion, these values drop to 7.4 and 3.8 kcal mol –1 in dichloromethane based up conductor-like polarizable continuum model calculations and suggest that X – ·M + X – ion triplets with a weakly coordinating cation maybe the reactive form of salts under some conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ Resonance Form Factors from 𝐵→ 𝜋⁢𝜋⁢ℓ⁢$\bar{v}$ in Lattice QCD

The decay 𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ is an attractive process for determining the magnitude of the smallest Cabibbo-Kobayashi-Maskawa matrix element, |𝑉 𝑢⁢𝑏 |, and can provide new insights into the origin of the long-standing exclusive-inclusive discrepancy in determinations of this standard-model parameter. This requires a nonperturbative QCD calculation of the 𝐵 → 𝜌 form factors 𝑉, 𝐴 0 , 𝐴 1 , and 𝐴 12 . The unstable nature of the 𝜌 resonance has prevented precise lattice QCD calculations of these form factors to date. Here, we present the first lattice QCD calculation of the 𝐵 → 𝜌 form factors in which the 𝜌 is treated properly as a resonance in 𝑃-wave 𝜋⁢𝜋 scattering. To this end, we use the Lellouch-Lüscher finite-volume formalism to compute the 𝐵 → 𝜋⁢𝜋 form factors as a function of both momentum transfer and 𝜋⁢𝜋 invariant mass, and then analytically continue to the 𝜌 resonance pole. This calculation is performed with 2 + 1 dynamical quark flavors at a pion mass of approximately 320 MeV, and demonstrates a clear path toward results at the physical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS