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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A general metriplectic framework with application to dissipative extended magnetohydrodynamics

General equations for conservative yet dissipative (entropy producing) extended magnetohydrodynamics are derived from two-fluid theory. Keeping all terms generates unusual cross-effects, such as thermophoresis and a current viscosity that mixes with the usual velocity viscosity. While the Poisson bracket of the ideal version of this model has already been discovered, we determine its metriplectic counterpart that describes the dissipation. This is done using a new and general thermodynamic point of view to derive dissipative brackets, a means of derivation that is natural for understanding and creating dissipative dynamics without appealing to underlying kinetic theory orderings. Finally, the formalism is used to study dissipation in the Lagrangian variable picture where, in the context of extended magnetohydrodynamics, non-local dissipative brackets naturally emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1976

Effect of anisotropy and rotation on turbulence production

The dynamical equations for the Reynolds stresses in a form which is suitable for incompressible homogeneous flows are used as a starting point in the investigation. It is assumed that the mean flow is two-dimensional. The equation for the rate of change of the kinetic energy of turbulence is derived. If the turbulence is isotropic, there is no turbulence production. However, the isotropic condition is not stable. The production mechanism is found to be autocatalytic. Rotation always inhibits production because it reduces anisotropy.

Ferziger, J. H.↗

NASA NTRS1999

The Impact of Buoyancy and Flame Structure on Soot, Radiation and NOx Emissions from a Turbulent Diffusion Flame

It is hypothesized that the spatial structure of a turbulent diffusion flame plays an important role in determining the emissions of radiative energy, soot and NO, from a combustor. This structure, manifested in the two point statistics, is influenced by buoyancy. Radiation, soot and NOx emissions are the cumulative result of processes that occur throughout a flame. For example, radiation fluxes along a line of sight can be found from summing up the contributions from sources in individual pockets of hot soot that emit, and from sinks in cold soot that absorb. Soot and NOx are both the results of slow chemistry and are not equilibrium products. The time that is available for production and burnout is crucial in determining the eventual emissions of these pollutants. Turbulence models generally rely on a single point closure of the appropriate time averaged equations. Hence, spatial information is lost and needs to be modeled using solution variables such as turbulence kinetic energy and dissipation rate, often with the assumption of isotropy. However, buoyancy can affect the physical structure of turbulent flames and can change the spatial extent of soot bearing regions. Theoretical comparisons with models are best done in the limit of infinite Froude number because the inclusion of buoyancy in flow models introduces significant uncertainties. Hence, LII measurements of soot, measurements of radiation fluxes from soot, Particle Imaging Velocimetry (PIV) of the flow field and measurements of post flame NOX will be carried out on the NASA Lewis 2.2 sec drop tower and eventually on the parabolic flight aircraft. The drop rig will be a modified version of a unit that has been successfully used at Lewis in the past.

Kennedy, I. M.↗

DOE OSTI2025

Particle acceleration during classical phase transitions on a spherical lattice

Abstract When compressed, certain lattices undergo phase transitions that may allow nuclei to gain significant kinetic energy. To explore the dynamics of this phenomenon, we develop a methodology to study Coulomb coupled N -body systems constrained to a sphere, as in the Thomson problem. We initialize N total Boron nuclei as point particles on the surface of the sphere, allowing them to equilibrate via Coulomb scattering with a viscous damping term. To simulate a phase transition, we remove N rm particles, forcing the system to rearrange into a new equilibrium. With this model, we consider the Thomson problem as a dynamical system, providing a framework to explore how non-zero temperature affects structural imperfections in Thomson minima. We develop a scaling relation for the average peak kinetic energy attained by a single particle as a function of N and N rm . For certain values of N , we find an order of magnitude energy gain when increasing N rm from 1 to 6. The model may help to design a lattice that maximizes the energy output.

Bachmann, Aidan M. (ORCID:0000000294706404)↗

DOE OSTI2020

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗

DOE OSTI2024

Heterogeneous Nucleating Agent for High‐Boiling‐Point Nonhalogenated Solvent‐Processed Organic Solar Cells and Modules

Abstract High‐boiling‐point nonhalogenated solvents are superior solvents to produce large‐area organic solar cells (OSCs) in industry because of their wide processing window and low toxicity; while, these solvents with slow evaporation kinetics will lead excessive aggregation of state‐of‐the‐art small molecule acceptors (e.g. L8‐BO), delivering serious efficiency losses. Here, a heterogeneous nucleating agent strategy is developed by grafting oligo (ethylene glycol) side‐chains on L8‐BO (BTO‐BO). The formation energy of the obtained BTO‐BO; while, changing from liquid in a solvent to a crystalline phase, is lower than that of L8‐BO irrespective of the solvent type. When BTO‐BO is added as the third component into the active layer (e.g. PM6:L8‐BO), it easily assembles to form numerous seed crystals, which serve as nucleation sites to trigger heterogeneous nucleation and increase nucleation density of L8‐BO through strong hydrogen bonding interactions even in high‐boiling‐point nonhalogenated solvents. Therefore, it can effectively suppress excessive aggregation during growth, achieving ideal phase‐separation active layer with small domain sizes and high crystallinity. The resultant toluene‐processed OSCs exhibit a record power conversion efficiency (PCE) of 19.42% (certificated 19.12%) with excellent operational stability. The strategy also has superior advantages in large‐scale devices, showing a 15.03‐cm 2 module with a record PCE of 16.35% (certificated 15.97%).

Chemistry↗

DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗

DOE OSTI2025

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1985

The magnetohydrodynamics of current sheets

Examples of current sheets are summarized and their formation is described. A universal phenomenon in cosmic plasmas is the creation of sheets off intense current near X-type neutral points (where the magnetic field vanishes). These sheets are important as sites where the magnetic-field energy is converted efficiently into heat and bulk kinetic energy and where particles can be accelerated to high energies. Examples include disruptions in laboratory tokamaks, substorms in the earth's magnetosphere, and flares on the sun. The basic behavior of a one-dimensional sheet is presented, together with an account of the linear tearing-mode instability that can cause the field lines in such a sheet to reconnect. Such reconnection may develop in different ways: it may arise from a spontaneous instability or it may be driven, either from outside by motions or locally by a resistivity enhancement. Various processes are described that may occur during the nonlinear development of tearing, along with the many numerical and laboratory experiments that are aiding our understanding of this intriguing cosmical process.

Priest, E. R.↗

NASA NTRSDate not supplied

Decoding the energy inputs that drive the upper atmosphere with NASA’s Geospace Dynamics Constellation (GDC)

NASA’s Geospace Dynamics Constellation (GDC) mission will study the effect of space weather phenomena on Earth’s upper atmosphere ion and neutral gases and, in turn, their impact on human assets. GDC will provide multi-point observations to determine how the Earth’s upper atmosphere is driven to dynamic changes and extremes by the solar wind. Magnetospheric energy input, both kinetic and electromagnetic, will be captured by the Comprehensive Auroral Precipitation Experiment (CAPE), the Thermal Plasma Sensor (TPS), and the Near Earth Magnetometer Instrument in a Small Integrated System (NEMISIS). CAPE captures the complete kinetic energy inputs with observations of the precipitating and upgoing electrons and precipitating ions. TPS captures part of the electromagnetic energy inputs with observations of the thermal plasma velocity that map E-field inputs, and NEMISIS captures the other part of the electromagnetic energy input with measurements of magnetic field variations due to field aligned currents. We use similar observations from the Defense Meteorological Space Platform (DMSP) satellites to demonstrate the balance between the different energy inputs (electron and ion Joule heating and Poynting flux) during the different phases of storms. We explore the local time and spatial structure and the impact they have on atmospheric density and winds with simulations of the Global Ionosphere and Thermosphere Model (GITM) and demonstrate how GDC will enable us to reconstruct those energy inputs locally and regionally

Eftyhia Zesta↗

NASA NTRS1994

(abstract) Odd Hydrogen in the Atmospheres of Earth and Mars

The Martian atmosphere has many features in common with the terrestrial mesosphere. Both share similar pressure and temperature ranges, and much of the same chemistry operates in each. For example, the radical species H, OH, and H(sub 2)O, which comprise the odd hydrogen family, are of central importance in the catalytic destruction of CO and O(sub 3) in both atmospheres. The inclusion of recent chemical kinetics data, specifically temperature dependent CO(sub 2) absorption cross-sections, into our one dimensional photochemical model of the Martial atmosphere shows that oxidation of CO by odd hydrogen is too efficient. The incorporation of smaller cross sections for CO(sub 2) leads to an enhanced photolysis rate of water vapor, increasing odd hydrogen to the point where the predicted mixing ratio of CO in our model is substantially less than the observed value of 6.5 x 10(sup -4). Interestingly, most photochemical models of the terrestrial mesosphere underestimate the CO and O(sub 3) densities using currently accepted photodissociation and kinetic rate coefficients. This has also been attributed to an overabundance of odd hydrogen in the models. We shall show that agreement between models and observations of CO in the Martian atmosphere as well as of CO and O(sub 3) in the terrestrial mesosphere can be achieved by revising the rate constants for the reactions OH + HO(sub 2) and CO + OH within their published uncertainties. The fact that similar revisions alleviate discrepancies in both the terrestrial and Martian atmospheres warrants a re-evaluation of these key rate constants at the appropriate temperatures and pressures.

dioxide ozone odd hydrogen family pressure tempera↗

DOE OSTI2021

Evaluating the benefits of kinetic Monte Carlo and microkinetic modeling for catalyst design studies in the presence of lateral interactions

Popular computational catalyst design strategies rely on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends. The main benefit of this approach is related to the inherent dimensionality reduction of the large material space to just a few catalyst descriptors. Conversely, it is well documented that a small set of descriptors is insufficient to capture the intricacies and complexities of a real catalytic system. The inclusion of coverage effects through lateral adsorbate-adsorbate interactions can narrow the gap between simplified descriptor predictions and real systems, but mean-field MKMs cannot properly account for local coverage effects. This shortcoming of the mean-field approximation can be rectified by switching to a lattice-based kinetic Monte Carlo (kMC) method using cluster expansion representation of adsorbate–adsorbate lateral interactions. Using the prototypical CO oxidation reaction as an example, we critically evaluate the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. Here, the mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, the MKM approach does provide a better entry point for computational catalyst design.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Crystallization Mechanisms of Poly(vinylidene Fluoride- co -chlorotrifluoroethylene): Nucleation Transitions, Growth Kinetics, and Microstructure Evolution

The crystallization kinetics of FK-800, a commercial semicrystalline copolymer of chlorotrifluoroethylene and vinylidene fluoride (poly(CTFE-co-VDF)), dictate its microstructure and functional performance. Crystallization occurs between the glass transition (T g ≈ 31 °C) and melting point (T m ≈ 110 °C), with grain boundaries playing a key role in applications such as memristors, where they regulate filament growth. This study integrates in situ atomic force microscopy (AFM), grazing-incidence wide-angle X-ray scattering (GIWAXS), and differential scanning calorimetry (DSC) to characterize nucleation, growth, and morphology across the full crystallization window. Hot-stage AFM reveals a transition from homogeneous to heterogeneous nucleation near 45 °C and a shift from reaction-limited to diffusion-limited growth above 60 °C. Kinetics are analyzed using Lauritzen–Hoffman, Turnbull–Fisher, and Avrami models to extract surface free energies and activation barriers. A multiscale modeling framework is developed in which nucleation densities and growth rates measured by AFM are used to reconstruct Avrami kinetics, which are then validated against coverage data, and extrapolated to predict bulk crystallization behavior observed by DSC. GIWAXS confirms a primarily flat-on chain orientation at all temperatures, linking orientation to domain morphology. This integrative approach quantitatively connects nanoscale crystallization dynamics with macroscopic phase evolution, establishing a quantitative framework for modeling crystallization across time, temperature, and length scales.

Crystallization↗

NASA NTRS1998

Joining of correct and incorrect DNA ends at double-strand breaks produced by high-linear energy transfer radiation in human fibroblasts

DNA double-strand breaks (DSBs) were measured within a 3.2-Mbp NotI fragment on chromosome 21 of cells of a normal human fibroblast cell line. Correct rejoining of DSBs was followed by measuring reconstitution of the original-size NotI fragment, and this was compared to total rejoining as measured by a conventional pulsed-field gel electrophoresis technique (FAR assay). After 80 Gy of particle irradiations with LETs in the range of 7-150 keV/microm, it was found that the repair kinetics was generally slower after irradiation with high-LET particles compared to X irradiation and that a larger proportion of the breaks remained unrepaired after 24 h. On the other hand, the misrejoining frequency as measured by the difference between correct and total rejoining after 24 h did not change with LET, but was approximately the same for all radiations at this dose, equal to 25-30% of the initial breaks. This result is discussed in relation to formation of chromosomal aberrations, deletion mutations and other biological end points.

NASA Discipline Radiation Health↗

NASA NTRS1992

Comments on the present state and future directions of PDF methods

The one point probability density function (PDF) method is examined in light of its use in actual engineering problems. The PDF method, although relatively complicated, appears to be the only format available to handle the nonlinear stochastic difficulties caused by typical reaction kinetics. Turbulence modeling, if it is to play a central role in combustion modeling, has to be integrated with the chemistry in a way which produces accurate numerical solutions to combustion problems. It is questionable whether the development of turbulent models in isolation from the peculiar statistics of reactant concentrations is a fruitful line of development as far as propulsion is concerned. There are three issues for which additional viewgraphs are prepared: the one point pdf method; the amplitude mapping closure; and a hybrid strategy for replacing a full two point pdf treatment of reacting flows by a single point pdf and correlation functions. An appeal is made for the establishment of an adequate data base for compressible flow with reactions for Mach numbers of unity or higher.

Obrien, E. E.↗

DOE OSTI2023

Energetics of high temperature degradation of fentanyl into primary and secondary products

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50–100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. In conclusion, we also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Poudel, Bharat↗

DOE OSTI2022

Chloride Salt Purification by Reaction With Thionyl Chloride Vapors to Remove Oxygen, Oxygenated Compounds, and Hydroxides

Molten chloride salts (including MgCl 2 , KCl, NaCl, and ZnCl 2 ) are being considered for heat transfer media for renewable (solar) and nuclear power generators, as fuel carrier for nuclear reactors, and as thermal energy storage media. Impurities such as oxygen, hydroxides, moisture, and sulfur are known to negatively influence the corrosion of materials in contact with the salt (e.g., structural metals). Commercially available chloride salts come with a range of impurities. Before using the chloride salts at high temperature, it is desirable to remove the impurities to increase the performance of the salt and reduce corrosion. In this study, we tested the use of thionyl chloride vaporized into a stream of argon to react with oxygenated impurities in a mixture of MgCl 2 -KCl-NaCl, removing them as HCl and SO 2 . The reagent was bubbled through the salt when both above and below the melting point. The reaction was followed using thermocouple data from the salt and by Fourier transform infrared (FTIR) spectroscopy on the exhaust of the reactor. The reaction kinetics were followed by comparing the peaks from SO 2 product to SOCl 2 reagent in the FTIR spectra. The purity of the salt was assessed at the end of the purification process by x-ray diffraction and inductively coupled plasma analysis. Although the process was effective in removing the oxygen content of the mixture, ternary compounds were formed in the process, including KNiCl 3 and KMgCl 3 . The nickel in KNiCl 3 came from the reaction between the salt and the nickel vessel. Thus, these experiments suggest that improvements to the process must be made before using SOCl 2 vapors for the purification of chloride salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗