SEARCH · Engineering Papers

Results for “kinetic simulation”

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 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

The Impact of Protonation on Early Translocation of Anthrax Lethal Factor: Kinetics from Molecular Dynamics Simulations and Milestoning Theory

We report atomically detailed molecular dynamics simulations of the permeation of the lethal factor (LF) N-terminal segment through the anthrax channel. The N-terminal chain is unstructured and leads the permeation process for the LF protein. The simulations were conducted in explicit solvent with milestoning theory, making it possible to extract kinetic information from nanosecond to millisecond time scales. Here, we illustrate that the initial event is strongly influenced by the protonation states of the permeating amino acids. While the N-terminal segment passes easily at high protonation state through the anthrax channel (and the Φ clamp), the initial permeation represents a critical step, which can be irreversible and establishes a hook in the channel mouth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Conversion kinetics during melting of simulated nuclear waste glass feeds measured by dissolution of silica

In this work, we experimentally investigated the batch-to-glass conversion kinetics of two low-activity and two high-level nuclear waste melter feeds that represent widely different behavior during melting; tests comprised thermogravimetric analysis, the feed expansion test, evolved gas analysis, and x-ray diffraction. Our study indicates that the conversion processes, which include gas-evolving reactions, dissolution of crystalline phases, and primary foam expansion and collapse, are interconnected and influenced by the feed’s thermal history. The fraction of dissolved silica was adopted as a measure of the extent of conversion inside the cold cap—the floating layer of feed material on the top surface of the melt in an electric melter. The interface between the cold cap and the free-flowing melt below was associated with a constant fraction of dissolved silica. The estimated viscosity of glass-forming melt at this interface lies within a narrow range for the feeds studied.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2014

Limitations of bootstrap current models

We assess the accuracy and limitations of two analytic models of the tokamak bootstrap current: (1) the well-known Sauter model and (2) a recent modification of the Sauter model by Koh et al. For this study, we use simulations from the first-principles kinetic code NEO as the baseline to which the models are compared. Tests are performed using both theoretical parameter scans as well as core- to-edge scans of real DIII-D and NSTX plasma profiles. The effects of extreme aspect ratio, large impurity fraction, energetic particles, and high collisionality are studied. In particular, the error in neglecting cross-species collisional coupling – an approximation inherent to both analytic models – is quantified. Moreover, the implications of the corrections from kinetic NEO simulations on MHD equilibrium reconstructions is studied via integrated modeling with kinetic EFIT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1996

Dissolution kinetics of a lunar glass simulant at 25 degrees C: the effect of pH and organic acids

The dissolution kinetics of a simulated lunar glass were examined at pH 3, 5, and 7. Additionally, the pH 7 experiments were conducted in the presence of citric and oxalic acid at concentrations of 2 and 20 mM. The organic acids were buffered at pH 7 to examine the effect of each molecule in their dissociated form. At pH 3, 5, and 7, the dissolution of the synthetic lunar glass was observed to proceed via a two-stage process. The first stage involved the parabolic release of Ca, Mg, Al, and Fe, and the linear release of Si. Dissolution was incongruent, creating a leached layer rich in Si and Ti which was verified by transmission electron microscopy (TEM). During the second stage the release of Ca, Mg, Al, and Fe was linear. A coupled diffusion/surface dissolution model was proposed for dissolution of the simulated lunar glass at pH 3, 5, and 7. During the first stage the initial release of mobile cations (i.e., Ca, Mg, Al, Fe) was limited by diffusion through the surface leached layer of the glass (parabolic release), while Si release was controlled by the hydrolysis of the Si-O-Al bonds at the glass surface (linear release). As dissolution continued, the mobile cations diffused from greater depths within the glass surface. A steady-state was then reached where the diffusion rate across the increased path lengths equalled the Si release rate from the surface. In the presence of the organic acids, the dissolution of the synthetic lunar glass proceeded by a one stage process. The release of Ca, Mg, Al, and Fe followed a parabolic relationship, while the release of Si was linear. The relative reactivity of the organic acids used in the experiments was citrate > oxalate. A thinner leached layer rich in Si/Ti, as compared to the pH experiments, was observed using TEM. Rate data suggest that the chemisorption of the organic anion to the surface silanol groups was responsible for enhanced dissolution in the presence of the organic acids. It is proposed that the increased rate of Si release is responsible for the one stage parabolic release of mobile cations and the relatively thin leached layer compared to experiments at pH 3 and 5.

NASA Discipline Life Support Systems↗

NASA NTRS1990

Simulation studies of plasma waves in the electron foreshock - The transition from reactive to kinetic instability

Particle simulation experiments were used to analyze the electron beam-plasma instability. It is shown that there is a transition from the reactive state of the electron beam-plasma instability to the kinetic instability of Langmuir waves. Quantitative tests, which include an evaluation of the dispersion relation for the evolving non-Maxwellian beam distribution, show that a quasi-linear theory describes the onset of this transition and applies again fully to the kinetic stage. This stage is practically identical to the late stage seen in simulations of plasma waves in the electron foreshock described by Dum (1990).

Dum, C. T.↗

DOE OSTI2012

Computer Simulations of Interstitial Loop Growth Kinetics in irradiated bcc Fe

The growth kinetics of (001) [001] interstitial loops in bcc Fe is studied by phase-field modeling. The effect of defect (vacancy/interstitial) concentration, generation, recombination, sink strength, and elastic interaction on the growth kinetics of interstitial loops is systematically simulated. Results show that the elastic interaction between the defects and interstitial loops speeds up the growth kinetics and affects the morphology of the interstitial loops. Linear growth rate, i.e., the loop average radius is linear to time, under both aging and irradiation are predicted, which is in agreement with experimental observation. The results also show that the interstitial loop growth rate, which is directly related to the sink strength of the interstitial loop for interstitials, increases linearly with the initial interstitial concentration during aging while changing logarithmically with the interstitial generation rate under irradiation.

growth kinetics↗

DOE OSTI2013

Pedestal fueling simulations with a coupled kinetic plasma–kinetic neutral transport code

The simplified neutral transport routine in the guiding center ion–electron-neutral neoclassical PIC code XGC0 has been replaced by one that calls the Monte Carlo neutral transport code DEGAS2. This provides a more realistic treatment of neutral atoms in the tokamak edge plasma and allows detailed atomic physics and plasma–material interaction processes to be incorporated into the simulations. The spatial profile of the neutral particle source used in DEGAS2 is determined from the fluxes of XGC0 ions to the material surfaces, and DEGAS2 is run in a time dependent manner synchronized with XGC0. The kinetic plasma–kinetic neutral transport capability is demonstrated with example pedestal fueling simulations.

Materials Science↗

NASA NTRS1995

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 2: Kinetic energy and resolution analysis

A kinetic energy (KE) analysis of the forcing of a mesoscale upper-tropospheric jet streak by organized diabatic processes within the simulated convective system (SCS) that was discussed in Part 1 is presented in this study. The relative contributions of the ageostrophic components of motion to the generation of KE of the convectively generated jet streak are compared, along with the KE generation by the rotational (nondivergent) and irrotational (divergent) mass transport. The sensitivity of the numerical simulations of SCS development to resolution is also briefly examined. Analysis within isentropic coordinates provides for an explicit determination of the influence of the diabatic processes on the generation of KE. The upper-level production of specific KE is due predominatly to the inertial advective ageostrophic component (IAD), and as such represents the primary process through which the KE of the convectively generated jet streak is realized. A secondary contribution by the inertial diabatic (IDI) term is observed. Partitioning the KE generation into its rotational and irrotational components reveals that the latter, which is directly linked to the diabatic heating within the SCS through isentropic continuity requirements, is the ultimate source of KE generation as the global area integral of generation by the rotational component vanishes. Comparison with an identical dry simulation reveals that the net generation of KE must be attributed to latent heating. Both the IAD and IDI ageostrophic components play important roles in this regard. Examination of results from simulations conducted at several resolutions supports the previous findings in that the effects of diabatic processes and ageostrophic motion on KE generation remain consistent. Resolution does impact the location and timing of SCS development, a result that has important implications in forecasting the onset of convection that develops from evolution of the large-scale flow and moisture transport. Marked differences are observed in the momentum field aloft subsequent to the life cycle of the SCS in the 1 deg, 30-level base case (MP130) simulation discussed in Part 1 versus its 2 deg counterparts in that the MP130 simulation with higher spatial resolution contains from 14% to 30% more total KE.

Wolf, Bart J.↗

DOE OSTI2013

Simulating the Heliosphere with Kinetic Hydrogen and Dynamic MHD Source Terms

The interaction between the ionized plasma of the solar wind (SW) emanating from the sun and the partially ionized plasma of the local interstellar medium (LISM) creates the heliosphere. The heliospheric interface is characterized by the tangential discontinuity known as the heliopause that separates the SW and LISM plasmas, and a termination shock on the SW side along with a possible bow shock on the LISM side. Neutral Hydrogen of interstellar origin plays a critical role in shaping the heliospheric interface, since it freely traverses the heliopause. Charge-exchange between H-atoms and plasma protons couples the ions and neutrals, but the mean free paths are large, resulting in non-equilibrated energetic ion and neutral components. In our model, source terms for the MHD equations are generated using a kinetic approach for hydrogen, and the key computational challenge is to resolve these sources with sufficient statistics. For steady-state simulations, statistics can accumulate over arbitrarily long time intervals. In this paper we discuss an approach for improving the statistics in time-dependent calculations, and present results from simulations of the heliosphere where the SW conditions at the inner boundary of the computation vary according to an idealized solar cycle.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

Parallel transport of long mean-free-path plasma along open magnetic field lines: Parallel heat flux

In a long mean-free-path plasma where temperature anisotropy can be sustained, the parallel heat flux has two components with one associated with the parallel thermal energy and the other the perpendicular thermal energy. Due to the large deviation of the distribution function from local Maxwellian in an open field line plasma with low collisionality, the conventional perturbative calculation of the parallel heat flux closure in its local or non-local form is no longer applicable. Here, a non-perturbative calculation is presented for a collisionless plasma in a two-dimensional flux expander bounded by absorbing walls. Specifically, closures of previously unfamiliar form are obtained for ions and electrons, which relate two distinct components of the species parallel heat flux to the lower order fluid moments such as density, parallel flow, parallel and perpendicular temperatures, and the field quantities such as the magnetic field strength and the electrostatic potential. The plasma source and boundary condition at the absorbing wall enter explicitly in the closure calculation. Although the closure calculation does not take into account wave-particle interactions, the results based on passing orbits from steady-state collisionless drift-kinetic equation show remarkable agreement with fully kinetic-Maxwell simulations. As an example of the physical implications of the theory, the parallel heat flux closures are found to predict a surprising observation in the kinetic-Maxwell simulation of the 2D magnetic flux expander problem, where the parallel heat flux of the parallel thermal energy flows from low to high parallel temperature region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Kinetics of Early Decomposition Stages in Diluted bcc Fe–Сu–Ni–Al Alloy: MC+MD Simulation

A combined approach including the Monte Carlo and molecular-dynamics simulation, decomposition kinetics and segregation formation in the multicomponent low-alloy Fe–1.5Cu–2Ni–1.5Al (at %) alloy is studied. It is shown that the formation of Cu nanoparticles surface-enriched with Ni and Al (coprecipitation mode) includes several stages: (i) the formation of clusters consisting of several Cu atoms, (ii) their enrichment with Ni and Al atoms, and (iii) redistribution of Ni and Al atoms with the formation of a surface layer providing stabilization of Cu nanoparticles. Observed structural features of segregations and their stability in Fe–Cu–Ni–Al alloys is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Spectral solver for multi-scale plasma physics simulations with dynamically adaptive number of moments

A spectral method for kinetic plasma simulations based on the expansion of the velocity distribution function in a variable number of Hermite polynomials is presented. The method is based on a set of non-linear equations that is solved to determine the coefficients of the Hermite expansion satisfying the Vlasov and Poisson equations. In this paper, we first show that this technique combines the fluid and kinetic approaches into one framework. Second, we present an adaptive strategy to increase and decrease the number of Hermite functions dynamically during the simulation. The technique is applied to the Landau damping and two-stream instability test problems. Performance results show 21% and 47% saving of total simulation time in the Landau and two-stream instability test cases, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Simulation of Nuclear Reactor Kinetics by the Monte Carlo Method

The KIR computer code intended for calculations of nuclear reactor kinetics using the Monte Carlo method is described. The algorithm implemented in the code is described in detail. Some results of test calculations are given.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2017

REVISITING A CLASSIC: THE PARKER–MOFFATT PROBLEM

The interaction of two colliding Alfvén wave packets is described here by means of magnetohydrodynamics (MHD) and hybrid kinetic numerical simulations. The MHD evolution revisits the theoretical insights described by Moffatt, Parker, Kraichnan, Chandrasekhar, and Elsässer in which the oppositely propagating large-amplitude wave packets interact for a finite time, initiating turbulence. However, the extension to include compressive and kinetic effects, while maintaining the gross characteristics of the simpler classic formulation, also reveals intriguing features that go beyond the pure MHD treatment.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Exceptional gettering response of epitaxially grown kerfless silicon

The bulk minority-carrier lifetime in p- and n-type kerfless epitaxial (epi) crystalline silicon wafers is shown to increase >500× during phosphorus gettering. We employ kinetic defect simulations and microstructural characterization techniques to elucidate the root cause of this exceptional gettering response. Simulations and deep-level transient spectroscopy (DLTS) indicate that a high concentration of point defects (likely Pt) is “locked in” during fast (60 °C/min) cooling during epi wafer growth. The fine dispersion of moderately fast-diffusing recombination-active point defects limits as-grown lifetime, but can also be removed during gettering, confirmed by DLTS measurements. Synchrotron-based X-ray fluorescence microscopy indicates metal agglomerates at structural defects, yet the structural defect density is sufficiently low to enable high lifetimes. Consequently, after phosphorus diffusion gettering, epi silicon exhibits a higher lifetime than materials with similar bulk impurity contents but higher densities of structural defects, including multicrystalline ingot and ribbon silicon materials. Device simulations suggest a solar-cell efficiency potential of this material >23%.

DLTS↗

DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗