SEARCH · Engineering Papers

Results for “Point Kinetics”

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 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS2003

The Mass of a Solar Quiescent Prominence

This paper follows up on our recent paper on the role of prominence mass in the storage of magnetic energy for driving a coronal mass ejection (CME). The previous paper erroneously rejected a set of sheet- prominence solutions, the recovery of which allows for a simple theoretical estimate of the mass of a quiescent prominence. For coronal fields of 5-10 G, these hydromagnetic solutions suggest that a prominence mass of (1-26) x 10(exp 6) g is needed to hold detached magnetic fields of intensity comparable to the coronal fields in an unbounded atmosphere such that the global magnetic field is energetically able to spontaneously open up and still have enough energy to account for the kinetic and gravitational potential energies carried away in a CME. This simple result is discussed in relation to observed prominence magnetic field intensities, densities, and masses, pointing to the relevance of such observations to the question of magnetic energy storage in the solar corona.

Low, B. C.↗

DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Enhanced steam oxidation resistance of uranium nitride nuclear fuel pellets

Here, the steam oxidation resistance of UN and UN-(20 vol%)ZrN fuel pellets is evaluated to enhance understanding of steam corrosion mechanisms in advanced nuclear fuel materials. In situ neutron diffraction shows the modified UN fuel pellets form a (U 0.77 ,Zr 0.23 )N solid-solution and the sole crystalline oxidation product detected in bulk is (U 0.77 ,Zr 0.23 )O 2 . U 2 N 3 is not detected in significant quantities during the steam oxidation of UN or (U 0.77 ,Zr 0.23 )N and stable lattice parameters show that hydriding does not take place. Steam oxidation rates, obtained via sequential Rietveld refinement show how (U 0.77 ,Zr 0.23 )N has a higher activation energy (79 ± 1 kJmol -1 vs. 50 ± 5 kJmol -1 ), higher onset temperature (430 °C vs. 400 °C) and slower reaction rates for steam oxidation up to 616 °C, than pure UN. Throughout, both UN and (U 0.77 ,Zr 0.23 )N exhibit linear (non-protective) oxidation kinetics, signifying that degradation of the fuel pellets is caused by the evolution of gaseous products at the interface followed by oxide scale spallation. This quantitative and mechanistic understanding of material degradation enables better defined operating regimes and points towards (U,Zr)N solid solutions as a promising strategy for the design of advanced nuclear fuel materials with enhanced steam corrosion resistance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2005

Thermal and Kinetic Properties of Motions in a Prominence Activation and Nearby Loop

We perform a quantitative analysis of the thermal properties of a prominence activation and motions in a nearby loop. In order to make measurements of the quickly moving features seen in prominences in the UV we use the SOHO/SUMER spectrograph to take a time series of exposures from a single pointing position, providing a measurement of spectral line properties as a function of time and position along the slit. The lines observed cover a broad range of temperatures from 80,000 - 1.6 million K. These measurements are combined with TRACE movies in transition region and coronal temperature bands to obtain more complete information concerning prominence structure and motions. The resulting observations allow us to analyze the thermal and kinetic energy of the moving sources as functions of time. The loop and prominence are most apparent in lines formed at temperatures below 250,000 K. We find that in most cases the temperature distribution of plasma in a moving feature changes relatively little over time periods of about 20 minutes.

Kucera, Therese↗

DOE OSTI2022

Crystallization kinetics and nanoparticle ordering in semicrystalline polymer nanocomposites

There has been considerable interest in the nucleation and crystallization of polymers in the presence of nanoparticles (NPs, or nanofillers in general, NFs). Most of the extensive work in this area has focused on anisotropic, non-Brownian NFs (e.g., clay sheets, carbon nanotubes) whose spatial dispersion state in these nanocomposites is controlled by the process by which they are formed. Hence, NF spatial dispersion is generally limited and often remains poorly characterized. Thermodynamic handles that can be used to control NF dispersion state in the polymer melt include (a) favorable interactions between the polymer chains and the bare NP surfaces, or (b) the density and length of the chains, with the same chemistry as the matrix, grafted to the NP surface. These relatively large NFs merely act as stationary objects that affect the kinetics of nucleation by providing heterogeneous sites, and the crystallization rate by confining the polymer in the melt state. The dispersion state of the NFs can dramatically affect the nucleation and crystallization of the matrix, but in most cases reported, the NFs increase nucleation efficiency relative to the neat polymer. At higher NF loadings, the effect of polymer confinement by the NFs dominates, leading to a decrease in crystal growth rates. This review describes the most important lessons learned from these commonly studied systems and then extends to polymer composite systems containing small, mobile spherical NPs (typically smaller than 100 nm in size). The role of NP mobility, which provides for dynamic confinement of the polymer melt, on the kinetics of polymer crystallization (nucleation, growth, and overall crystallization) and how this behavior is mostly consistent with the case of immobile NF is a second important focus of this review. In addition to the role of NFs on crystallization kinetics, recently reported nanoparticle ordering phenomena such as the effect of matrix crystallization on the organization of small spherical NPs within the amorphous regions of the semicrystalline morphology are discussed. In conclusion, such phenomena are clearly not observed for large NFs and hence provide a point of departure from past works in this area.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1996

Single residue substitutions that change the gating properties of a mechanosensitive channel in Escherichia coli

MscL is a channel that opens a large pore in the Escherichia coli cytoplasmic membrane in response to mechanical stress. Previously, we highly enriched the MscL protein by using patch clamp as a functional assay and cloned the corresponding gene. The predicted protein contains a largely hydrophobic core spanning two-thirds of the molecule and a more hydrophilic carboxyl terminal tail. Because MscL had no homology to characterized proteins, it was impossible to predict functional regions of the protein by simple inspection. Here, by mutagenesis, we have searched for functionally important regions of this molecule. We show that a short deletion from the amino terminus (3 amino acids), and a larger deletion of 27 amino acids from the carboxyl terminus of this protein, had little if any effect in channel properties. We have thus narrowed the search of the core mechanosensitive mechanism to 106 residues of this 136-amino acid protein. In contrast, single residue substitutions of a lysine in the putative first transmembrane domain or a glutamine in the periplasmic loop caused pronounced shifts in the mechano-sensitivity curves and/or large changes in the kinetics of channel gating, suggesting that the conformational structure in these regions is critical for normal mechanosensitive channel gating.

Non-NASA Center↗

NASA NTRS1993

A Mechanism for Bulk Energization in the Impulsive Phase of Solar Flares: MHD Turbulent Cascade

We propose that the large production rate (approximately 10(exp 36)/s) of energetic electrons (greater than or approximately equal to 25 keV) required to account for the impulsive-phase hard X-ray burst in large flares is achieved through MHD turbulent cascade of the bulk kinetic energy of the outflows from many separate reconnection events. Focusing on large two- ribbon eruptive flares as representative of most large flares, we envision the reconnection events to be the driven reconnection of oppositely directed elementary flux tubes pressing into the flare-length current-sheet interface that forms in the wake of the eruption of the sheared core of the preflare bipolar field configuration. We point out that, because the outflows from these driven reconnection events have speeds of order the Alfven speed and because the magnetic field reduces the shear viscosity of the plasma, it is reasonable that the outflows are unstable and turbulent, so that the kinetic energy of an outflow is rapidly dissipated through turbulent cascade. If the largest eddies in the turbulence have diameters of order the expected widths of the outflows (10(exp 7)-10(exp 8)cm), then the cascade dissipation of each of these eddies could produce approximately 10(exp 26) erg burst of energized electrons (approximately 3 x (10(exp 33) 25 keV electrons) in approximately 0.3 s, which agrees well with hard X-ray and radio sub-bursts commonly observed during the impulsive phase. Of order 10(exp 2) simultaneous reconnection events with turbulent outflow would produce the observed rate of impulsive-phase plasma energization in the most powerful flares (approximately 10(exp 36) 25 keV electrons/ s); this number of reconnection sites can easily fit within the estimated 3 x 10(exp 9) cm span of the overall current-sheet dissipation region formed in these large flares. We therefore conclude that MHD turbulent cascade is a promising mechanism for the plasma energization observed in the impulsive phase of solar flares.

LaRosa, T. N.↗

DOE OSTI2026

PyJMAK: An Open-Source Python Toolkit for Modeling Solid-State Metallurgical Phase Transformations

Accurate prediction of metallurgical phase transformations is an essential basis for autonomous optimization and rapid part qualification. Several methods can be used to estimate the evolution of phase fractions such as JMAK kinetics-based models, phase-field models, thermodynamic models, and data-driven machine learning models. Thermodynamic and phase-field-based methodologies solve multiphysics equations requiring numerous calibration parameters and significant computational resources. As a result, the computation domain is limited to a point or on order of micron-meters. The data-driven models rely on large datasets from experiments and simulations. While the JMAK model only provides information about phase fraction evolution, it can predict this evolution in near real-time using thermal history and thermodynamic data without restriction on the domain. JMAK models have been popularly used by researchers to model phase transformations occuring during additive manufacturing or over arbitrary temperature profiles. Commercial proprietary software such as Abaqus and Ansys or closed-source in-house implementations offer the ability to model JMAK based kinetics to predict phase transformation. However, these software packages are not open-source or freely available for use and development in conjunction with manufacturing machines, sensors, and machine learning algorithms. In addition, the use of the model is restricted by a license token. In contrast, given temperature profiles at multiple points in the domain, this Python-based PyJMAK model can compute phase evolution in parallel due to its stand-alone modular, voxel-based structure, and it can be executed on high-performance computing resources without any license restrictions.

Prabhune, Bhagya [Oak Ridge National Laboratory (O↗

DOE OSTI2025

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The reaction NH 2 + CH 2 O: Kinetic measurements and computational studies

The reaction of amidogen with formaldehyde is relevant to astrochemistry and couples nitrogen and carbon chemistry in flames. The first measurements of the temperature dependence of the rate constant k were made, over 409–643 K, and are summarized as k = 7.1 × 10 −12 exp(−18.2 kJ mol −1 /RT) cm 3 molecule −1 s −1 with a confidence interval of ±20%. These results address a controversy over the presence of a significant barrier in the addition/elimination pathway leading to the production of formamide between theoretical models that include the zero-point energy of all modes at the transition state and a model that excludes transitional modes. The latter would lead to a negligible barrier, which is inconsistent with the experiments. Contrary to earlier claims, CCSDT(Q)-corrected energies combined with transition state theory yield quantitative accord with the measurements.

Astrochemistry↗

DOE OSTI2021

A Kinetic Monte Carlo Approach for Simulating Cascading Transmission Line Failure

In this work, cascading transmission line failures are studied through a dynamical model of the power system operating under fixed conditions. The power grid is modeled as a stochastic dynamical system where first-principles electromechanical dynamics are excited by small Gaussian disturbances in demand and generation around a specified operating point. In this context, a single line failure is interpreted in a large deviation context as a first escape event across a surface in phase space defined by line security constraints. The resulting system of stochastic differential equations admits a transverse decomposition of the drift, which leads to considerable simplification in evaluating the quasipotential (rate function) and, consequently, computation of exit rates. Tractable expressions for the rate of transmission line failure in a restricted network are derived from large deviation theory arguments and validated against numerical simulations. Extensions to realistic settings are considered, and individual line failure models are aggregated into a Markov model of cascading failure inspired by chemical kinetics. Cascades are generated by traversing a graph composed of weighted edges representing transitions to degraded network topologies. Numerical results indicate that the Markov model can produce cascades with qualitative power-law properties similar to those observed in empirical cascades.

Kinetic Monte Carlo, rare events, Cascading Failur↗

DOE OSTI2024

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗

NASA NTRS1987

Effect of shaft rotation on the incompressible flow in a labyrinth seal

The incompressible flow in a labyrinth seal at low leakage rates was computationally and experimentally investigated over a wide range of seal rotation rates. QUICK differencing was employed in the finite difference code to reduce the effects of false diffusion. The use of measured inlet boundary conditions for the axial and swirl velocity components and for the turbulent kinetic energy resulted in good agreement between velocity predictions and hot-film measurements. It was found that when the rotation rate is increased beyond a certain point, a second recirculation zone forms inside the seal cavity, altering the flow field in the cavity and resulting in a substantial increase in the pressure drop across it.

Demko, J. A.↗

NASA NTRS1990

Numerical methods for TVD transport and coupled relaxing processes in gases and plasmas

The construction of second-order upwind schemes for nonequilibrium plasmas, for both one- and two-fluid formulations is demonstrated. Coupled relaxation processes, including ionization kinetics and radiative processes and their algorithms for nonequilibrium, multiple temperature conditions are described as well. The paper applies the numerical techniques on some simple test cases, points out critical problems and their solutions, and makes qualitative comparisons with known results, whenever possible.

Cambier, Jean-Luc↗

NASA NTRS1993

Linear theory of fast reconnection at an X-type neutral point

A complete analytic theory of dynamic linear reconnection at an X-type neutral point is developed. An eigenmode analysis, using cylindrical coordinates centered on the neutral point, extends the work of Craig & McClymont (1991) to include nonazimuthally symmetric perturbations on the flux function. It is demonstrated that all physically significant disturbances, both reconnective and nonreconnective, decay resistively on a 'fast' time scale. The significance of the linear theory is discussed within the general context of steady state and dynamic reconnection studies. Disturbances can be expected to focus explosively in the vicinity of the neutral point. This suggests the formation of a 'flux pile-up' current layer in which the bulk of the magnetic energy is released as heat rather than kinetic energy of mass motion.

Craig, I. J. D.↗

DOE OSTI2025

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗

DOE OSTI2022

A Gauss’s law analysis of redox active adsorbates on semiconductor electrodes: The charging and faradaic currents are not independent

A detailed framework for modeling and interpreting the data in totality from a cyclic voltammetric measurement of adsorbed redox monolayers on semiconductor electrodes has been developed. A three-layer model consisting of the semiconductor space-charge layer, a surface layer, and an electrolyte layer is presented that articulates the interplay between electrostatic, thermodynamic, and kinetic factors in the electrochemistry of a redox adsorbate on a semiconductor. Expressions are derived that describe the charging and faradaic current densities individually, and an algorithm is demonstrated that allows for the calculation of the total current density in a cyclic voltammetry measurement as a function of changes in the physical properties of the system (e.g., surface recombination, dielectric property of the surface layer, and electrolyte concentration). The most profound point from this analysis is that the faradaic and charging current densities can be coupled. That is, the common assumption that these contributions to the total current are always independent is not accurate. Their interrelation can influence the interpretation of the charge-transfer kinetics under certain experimental conditions. More generally, this work not only fills a long-standing knowledge gap in electrochemistry but also aids practitioners advancing energy conversion/storage strategies based on redox adsorbates on semiconductor electrodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗