Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of Scalable Reactive Transport Framework for PFAS

Research on PFAS chemicals is extensive and covers all aspects, including analytical quantification, determination of properties, toxicology, ex situ treatment, and in situ remediation. PFAS chemicals are very stable and thus persistent/recalcitrant in the environment. Although there are many unknowns about PFAS chemicals, degradation pathways, reaction rates, etc., many different sorption, oxidation, reduction, biological, and innovative treatment approaches are being developed. Sorption with activated carbon is currently the only fully available in situ treatment technology for PFAS-impacted groundwater. Given the wide array of PFAS chemicals and transformation products, remediation may need multi-step treatment trains to fully address the PFAS contamination. The work here provides kinetic reaction modules that represent an initial cut at functionality representing PFAS migration and reaction in groundwater aquifer flow and transport models. One reaction kinetics module provides a method to model kinetically limited adsorption using a mass transfer model. The second reaction module represents biological transformation of 8:2 FTOH and daughter species, illustrating how a complex reaction pathway network can be represented. Both reaction modules allow for spatially variable parameter values so that a variety of remediation approaches (e.g., a PRB or volumetric treatment or variations in geochemical conditions) can be simulated. The intent with these PFAS reaction modules is to provide tools for practitioners to aid in the selection, design, and assessment of potential in situ PFAS remediation strategies. It is anticipated that, as PFAS remediation technologies and scientific understanding advances, these modules would be refined or replaced to match new knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Simplified half-life methods for the analysis of kinetic data

The analysis of reaction rate data has as its goal the determination of the order rate constant which characterize the data. Chemical reactions with one reactant and present simplified methods for accomplishing this goal are considered. The approaches presented involve the use of half lives or other fractional lives. These methods are particularly useful for the more elementary discussions of kinetics found in general and physical chemistry courses.

Eberhart, J. G.↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2024

Comparison of second-order spectral line widths formulae

Frequently used second-order spectral line width formulae from the projector operator and kinetic theory methods have not been formally compared previously. Here it is shown that a systematic second-order expansion of the projection operator expression including initial correlations agrees with the second-order kinetic theory result. The agreement assumes a common approximation in the projector operator method that introduces a screened radiator–perturber interaction to account for neglected electron–electron correlations. On the other hand, it is shown that the usual width expression from the projection operator approach neglecting initial correlations differs from kinetic theory. The differences, however, are at least third order in the radiator–perturber interaction. The comparisons suggest using the more compact width expression from kinetic theory, which includes initial correlations and a systematic screening of the radiator–perturber interactions, as the starting point for second-order width calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Identification of local modes using modal mass distributions

An approach for local modal identification using a modal mass spatial distribution matrix is proposed. This modal mass matrix is obtained by overlay multiplying the modal matrix to the product of the mass and the modal matrices of the structure. Each column in the spatial modal mass matrix can be interpreted as a measure of the kinetic energy distribution over the physical degrees of freedom for the corresponding mode. Local modes are identified due to the highly concentrated kinetic energy distribution in the localized regions. This spatial modal mass approach is more sensitive than the modal amplitude ratio method, since kinetic energy is in the quadratic form of modal amplitudes. Numerical examples involving identification of local modes using the proposed modal mass spatial distribution matrix approach are presented and discussed.

Chou, Chaur-Ming↗
NASA NTRS1995

Grain formation around carbon stars. 1: Stationary outflow models

Asymptotic giant branch (AGB) stars are known to be sites of dust formation and undergo significant mass loss. The outflow is believed to be driven by radiation pressure on grains and momentum coupling between the grains and gas. While the physics of shell dynamics and grain formation are closely coupled, most previous models of circumstellar shells have treated the problem separately. Studies of shell dynamics typically assume the existence of grains needed to drive the outflow, while most grain formation models assume a constant veolcity wind in which grains form. Furthermore, models of grain formation have relied primarily on classical nucleation theory instead of using a more realistic approach based on chemical kinetics. To model grain formation in carbon-rich AGB stars, we have coupled the kinetic equations governing small cluster growth to moment equations which determine the growth of large particles. Phenomenological models assuming stationary outflow are presented to demonstrate the differences between the classical nucleation approach and the kinetic equation method. It is found that classical nucleation theory predicts nucleation at a lower supersaturation ratio than is predicted by the kinetic equations, resulting in significant differences in grain properties. Coagulation of clusters larger than monomers is unimportant for grain formation in high mass-loss models but becomes more important to grain growth in low mass-loss situations. The properties of the dust grains are altered considerably if differential drift velocities are ignored in modeling grain formation. The effect of stellar temperature, stellar luminosity, and different outflow velocities are investigated. The models indicate that changing the stellar temperature while keeping the stellar luminosity constant has little effect on the physical parameters of the dust shell formed. Increasing the stellar luminosity while keeping the stellar temperature constant results in large differences in grain properties. For small outflow velocities, grains form at lower supersaturation ratios and close to the stellar photosphere, resulting in larger but fewer grains. The reverse is true when grains form under high outflow velocities, i.e., they form at higher supersaturation ratios, farther from the star, and are much smaller but at larger quantities.

Egan, Michael P.↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Study of Kinetics of Solid Phase Transition in Tetracosane С{sub 24}Н{sub 50} by High-Resolution Synchrotron X-Ray Powder Diffraction

Using the small-angle synchrotron X-ray diffraction method at the research facility “BELOK” of the Kurchatov synchrotron radiation source at the NRC Kurchatov Institute, it was shown that the first-order solid phase transition in tetracosane C{sub 24}H{sub 50} passes according to a heterogeneous mechanism in a narrow temperature range ΔT ~ 1 K in accordance with the theory of diffuse phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Time-dependent numerical modeling of dust halo formation at comets

The evolution of gas and dust distributions following a spatially and temporally localized comet outburst was calculated using a hybrid kinetic-hydrodynamic method. In the inner coma the time-dependent continuity, momentum, and energy equations of the dusty gas flow were solved simultaneously using 12 dust sizes. Beyond 300 km a three-dimensional kinetic model was used to calculate the trajectory of each individual dust grain. It was found that following the onset of the comet outburst a gas-dust blast wave propagates outward in the inner coma. About 15 minutes after the increased gas and dust production was initiated at the nucleus, a new equilibrium was reached in the inner coma. The most important feature of this new steady state was the significant increase of the dust terminal velocities. These higher terminal velocity values resulted in larger apex distances for dust particles emitted during the outburst. The dust particles spend a relatively long time near their apex points; therefore, the outburst generates long-lasting distinct dust envelopes in front of the regular dust coma.

Gombosi, T. I.↗
DOE OSTI2020

Trace Chemical Detection Using Intercalated MXenes as a Signal Enhancing Substrate in Optical Probes

MXenes are 2D materials composed of layered transition metal nitrides or carbides. These materials are synthesized by HF exfoliation from MAX phases (Ti{sub 3}AlC{sub 2}). The 2D nanomaterial was synthesized by the removal of the 'A' element, resulting in a Mxene product (Ti{sub 3}C{sub 2}). MXenes have the general formula M{sub n+1}X{sub n}T{sub x}, where M is an early transition metal, X is Carbon and/or Nitrogen, and T accounts for surface terminated functional groups such as Fluoride, hydroxyl, and oxygen. These materials have very unique properties, similar to graphene, that allows them to be applied in a variety of trace detection techniques including surface-enhanced Raman spectroscopy (SERS). MXenes have also been demonstrated to selectively uptake uranyl ion, UO{sub 2}{sup 2+}. If this property can be combined with SERS or fluorescence detection, it may be possible to use MXenes as the basis for an alternative method to kinetic phosphorescence analysis (KPA) for trace uranyl measurements. Objectives: To confirm that MXene Nano materials are suitable substrates for SERS and sensor development by enhancing Raman signaling. To determine if certain MXene preparation methods yield materials that are more suitable for trace sensing methods. To determine uranyl uptake properties of these MXene materials and test them for analytical signals. Sample Preparation: Preparation of Ti{sub 3}C{sub 2}MXene (at FSU). MXenes were prepared by etching Al from Ti{sub 3}AlC{sub 2} (MAX phase)material. Two etching techniques yield different MXene products: LiF/HCl: Milder reaction, larger MXene flakes. HF: Harsher reaction, smaller flakes, larger layer separation. Products washed to remove etchant, vacuum filtered, and dried. Dried MXene flakes are air-stable. Film preparation for sensor testing (at SRNL): Suspend powder in diH{sub 2}O, purge with Ar, sonicate for 30 min. Centrifuge and collect supernate with suspended particles. Observed LiF-etched Mxene yielded a higher density of particles and darker collected solution. Drop-cast (4 ml) supernate onto slides and dried with Ar. For Rhodamine B (RhB) testing, drop-cast 4 ml drops onto Mxene spots and dried with Ar. Scanning Electron Microscopy conditions: 10 kV Beam energy, high vacuum; Working distance of 8 mm; beam penetration depth appx. 4 microns, beam spot size appx. 2 nanometers. Results: Detection of aluminum correlates with bright spots on image. Presence of aluminum shows that LiF/HCl etching was less thorough than HF etching. Trace Cl detection in LiF images suggests incomplete rinsing. HF has smaller feature size, more layer structure, and increased homogeneity, consistent with expectations. Macroscopic Raman spectroscopy measurements: 532 nm excitation, ∼50 mW with a ∼100 micron spot size (InPhotonics RPB probe). Kaiser Optical Holospec f/1.8 spectrometer with cooled (-60 deg.C) Andor iDus OE420 CCD. LiF 1x supernate showed good signal for trace measurements of Rhodamine B. HF and 1/4x LiF supernates showed little Mxene or Rhodamine B signal. Low deposition densities led to excess background signal from glass slides. For LiF film, response is linear with Rhodamine B concentration over range tested. Will retest with Raman microscope (∼1 micron spot size) to characterize SERS of more dilute LiF and HF etched Mxenes. Conclusions: The LiF etched material was more suitable for macroscopic SERS measurements because it was more concentrated, resulting in a thicker film than the HF etched Mxene and diluted LiF sample. However, the other materials may give greater SERS enhancements, which we hope to determine from measurements with the Raman microscope. From characterization with SEM we concluded that the HF etched Mxene is more uniform/homogenous and has smaller particle size than the LiF etched Mxene. There is still aluminum present in both samples indicating that etching wasn't complete, but the removal of the aluminum was more efficient in the HF method. Path Forward: Observe SERS with Raman microscopy, to obtain better signals for the more diluted samples and be able to compare enhancement effects for the different MXenes. Characterize uranyl sorption into MXene films and test Raman and fluorescence signals. Revisit the etching conditions to improve removal of aluminum. FSU and SRNL will continue to collaborate to create and characterize different Mxene materials and test their usefulness for sensor applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Predicting Asteroid Material Properties from a DART-like Kinetic Impact

Abstract NASA’s Double Asteroid Redirection Test (DART) mission is the first full-scale test of the kinetic impactor method for asteroid deflection, in which a spacecraft intentionally impacts an asteroid to change its trajectory. DART represents an important first step for planetary defense technology demonstration, providing a realistic assessment of the effectiveness of the kinetic impact approach on a near-Earth asteroid. The momentum imparted to the asteroid is transferred from the impacting spacecraft and enhanced by the momentum of material ejected from the impact site. However, the magnitude of the ejecta contribution is dependent on the material properties of the target. These properties, such as strength and shear modulus, are unknown for the DART target asteroid, Dimorphos, as well as most asteroids since such properties are difficult to characterize remotely. This study examines how hydrocode simulations can be used to estimate material properties from information available post-impact, specifically the asteroid size and shape, the velocity and properties of the impacting spacecraft, and the final velocity change imparted to the asteroid. Across >300 three-dimensional simulations varying seven material parameters describing the asteroid, we found many combinations of properties could reproduce a particular asteroid velocity. Additional observations, such as asteroid mass or crater size, are required to further constrain properties like asteroid strength or outcomes like the momentum enhancement provided by impact ejecta. Our results demonstrate the vital importance of having as much knowledge as possible prior to an impact mission, with key material parameters being the asteroid’s mass, porosity, strength, and elastic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

KiT-RT: An Extendable Framework for Radiative Transfer and Therapy

Here, in this article, we present Kinetic Transport Solver for Radiation Therapy (KiT-RT), an open-source C++-based framework for solving kinetic equations in therapy applications available at https://github.com/CSMMLab/KiT-RT . This software framework aims to provide a collection of classical deterministic solvers for unstructured meshes that allow for easy extendability. Therefore, KiT-RT is a convenient base to test new numerical methods in various applications and compare them against conventional solvers. The implementation includes spherical harmonics, minimal entropy, neural minimal entropy, and discrete ordinates methods. Solution characteristics and efficiency are presented through several test cases ranging from radiation transport to electron radiation therapy. Due to the variety of included numerical methods and easy extendability, the presented open-source code is attractive for both developers, who want a basis to build their numerical solvers, and users or application engineers, who want to gain experimental insights without directly interfering with the codebase.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Linking stem rehydration kinetics to hydraulic traits using a novel method and mechanistic model

Background: Despite the recognized importance of hydraulic capacitance as a mechanism used by plants to maintain hydraulic functioning during high transpiration, characterizing the dynamics of capacitance remains a challenge. Methods: For this work, we used a novel ‘two-balance method’ to investigate relationships between stem rehydration kinetics and other hydraulic traits in multiple tree species, and we developed a model to explore stem rehydration kinetics further. Key Results: We found that: (1) rehydration time constants and the amount of water uptake occurring during rehydration differed significantly across species; (2) time constants did not change with declining water potential (Ψ), while water uptake increased at lower Ψ in some species; (3) longer time constants were associated with lower wood density, higher capacitance and less negative stem pressures causing 50 % loss of hydraulic conductivity (P 50 ); (4) greater water uptake occurred in stems with lower wood density and less negative P 50 values; and (5) the model could estimate the total hydraulic resistance of the rehydration path, which cannot be measured directly. Conclusions: Overall, the two-balance method can be used to examine rehydration dynamics quickly and thoroughly in detached woody stems. This method has the potential to improve our understanding of how capacitance functions across tree species, which is an often-overlooked component of whole-plant hydraulics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Kinetically controlled metal-elastomer nanophases for environmentally resilient stretchable electronics

Nanophase mixtures, leveraging the complementary strengths of each component, are vital for composites to overcome limitations posed by single elemental materials. Among these, metal-elastomer nanophases are particularly important, holding various practical applications for stretchable electronics. However, the methodology and understanding of nanophase mixing metals and elastomers are limited due to difficulties in blending caused by thermodynamic incompatibility. Here, we present a controlled method using kinetics to mix metal atoms with elastomeric chains on the nanoscale. We find that the chain migration flux and metal deposition rate are key factors, allowing the formation of reticular nanophases when kinetically in-phase. Moreover, we observe spontaneous structural evolution, resulting in gyrified structures akin to the human brain. The hybridized gyrified reticular nanophases exhibit strain-invariant metallic electrical conductivity up to 156% areal strain, unparalleled durability in organic solvents and aqueous environments with pH 2–13, and high mechanical robustness, a prerequisite for environmentally resilient devices.

36 MATERIALS SCIENCE↗