Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

The balance of kinetic and total energy simulated by the OSU two-level atmospheric general circulation model for January and July

A three-year simulation which prescribes seasonally varying solar radiation and sea surface temperature is the basis of the present study of the horizontal structure of the balances of kinetic and total energy simulated by Oregon State University's two-level atmospheric general circulation model. Mechanisms responsible for the local energy changes are identified, and the energy balance requirement's fulfilment is examined. In January, the vertical integral of the total energy shows large amounts of external heating over the North Pacific and Atlantic, together with cooling over most of the land area of the Northern Hemisphere. In July, an overall seasonal reversal is found. Both seasons are also characterized by strong energy flux divergence in the tropics, in association with the poleward transport of heat and momentum.

Wang, J.-T.↗
DOE OSTI2016

Hybrid framework for the simulation of stochastic chemical kinetics

Stochasticity plays a fundamental role in various biochemical processes, such as cell regulatory networks and enzyme cascades. Isothermal, well-mixed systems can be modelled as Markov processes, typically simulated using the Gillespie Stochastic Simulation Algorithm (SSA) [25]. While easy to implement and exact, the computational cost of using the Gillespie SSA to simulate such systems can become prohibitive as the frequency of reaction events increases. This has motivated numerous coarse-grained schemes, where the “fast” reactions are approximated either using Langevin dynamics or deterministically. While such approaches provide a good approximation when all reactants are abundant, the approximation breaks down when one or more species exist only in small concentrations and the fluctuations arising from the discrete nature of the reactions become significant. This is particularly problematic when using such methods to compute statistics of extinction times for chemical species, as well as simulating non-equilibrium systems such as cell-cycle models in which a single species can cycle between abundance and scarcity. In this paper, a hybrid jump-diffusion model for simulating well-mixed stochastic kinetics is derived. It acts as a bridge between the Gillespie SSA and the chemical Langevin equation. For low reactant reactions the underlying behaviour is purely discrete, while purely diffusive when the concentrations of all species are large, with the two different behaviours coexisting in the intermediate region. A bound on the weak error in the classical large volume scaling limit is obtained, and three different numerical discretisations of the jump-diffusion model are described. The benefits of such a formalism are illustrated using computational examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Electron Distributions in Kinetic Scale Field Line Resonances: A Comparison of Simulations and Observations

Observations in kinetic scale field line resonances, or eigenmodes of the geomagnetic field, reveal highly field-aligned plateaued electron distributions. By combining observations from the Van Allen Probes and Cluster spacecraft with a hybrid kinetic gyrofluid simulation we show how these distributions arise from the nonlocal self-consistent interaction of electrons with the wavefield. This interaction is manifested as electron trapping in the standing wave potential. The process operates along most of the field line and qualitatively accounts for electron observations near the equatorial plane and at higher latitudes. In conjunction with the highly field-aligned plateaus, loss cone features are also evident, which result from the action of the upward-directed wave parallel electric field on the untrapped electron populations.

42 ENGINEERING↗
NASA NTRS1984

Threshold kinetics of a solar-simulator-pumped iodine laser

A model of the chemical kinetics of the n-C3F7I solar-simulator-pumped iodine laser is utilized to study the major kinetic processes associated with the threshold behavior of this experimental system. Excited-state diffusion to the cell wall is the dominant limiting factor below 5 torr. Excited-state diffusion to the cell wall is the dominant limiting factor below 5 torr. Excited-state recombination with the alkyl radical and quenching by the parent gas control threshold at higher pressures. Treatment of the hyperfine splitting and uncertainty in the pressure broadening are important factors in fixing the threshold level. In spite of scatter in the experimental data caused by instabilities in the simulator high-pressure high-pressure arc, reasonable agreement is achieved between the model and experiment. Model parameters arrived at are within the uncertainty range of values found in the literature.

Wilson, J. W.↗
DOE OSTI2020

Kinetic entropy-based measures of distribution function non-Maxwellianity: theory and simulations

We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson (J. Geophys. Res., vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Gyrokinetic and kinetic particle-in-cell simulations of guide-field reconnection. I. Macroscopic effects of the electron flows

In this work, we compare gyrokinetic (GK) with fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields (b{sub g}). Here, we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma (β{sub i} = 0.01), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field (b{sub g} ≳ 30) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive (b{sub g} ≳ 5). Kinetic PIC simulations using guide fields b{sub g} ≲ 30 reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma (β{sub i} = 1.0) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields (b{sub g} ≲ 3)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Solar-simulator-pumped atomic iodine laser kinetics

The literature contains broad ranges of disagreement in kinetic data for the atomic iodine laser. A kinetic model of a solar-simulator-pumped iodine laser is used to select those kinetic data consistent with recent laser experiments at the Langley Research Center. Analysis of the solar-simulator-pumped laser experiments resulted in the following estimates of rate coefficients: for alkyl radical (n-C3F7) and atomic iodine (I) recombination, 4.3 x 10 to the 11th power (1.9) + or - cu cm/s; for n-C3F7I stabilized atomic iodine recombination (I + I) 3.7 x 10 to the -32nd power (2.3) + or -1 cm to the 6th power/s; and for molecular iodine (I2) quenching, 3.1 x 10 to the -11th power (1.6) + or - 1 cu cm/s. These rates are consistent with the recent measurements.

Wilson, H. W.↗
DOE OSTI2016

TURBULENCE AND PROTON–ELECTRON HEATING IN KINETIC PLASMA

Analysis of particle-in-cell simulations of kinetic plasma turbulence reveals a connection between the strength of cascade, the total heating rate, and the partitioning of dissipated energy into proton heating and electron heating. A von Karman scaling of the cascade rate explains the total heating across several families of simulations. The proton to electron heating ratio increases in proportion to total heating. We argue that the ratio of gyroperiod to nonlinear turnover time at the ion kinetic scales controls the ratio of proton and electron heating. The proposed scaling is consistent with simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Verification of electromagnetic fluid-kinetic hybrid electron model in global gyrokinetic particle simulation

The fluid-kinetic hybrid electron model is verified in global gyrokinetic particle simulation of linear electromagnetic drift-Alfvenic instabilities in tokamak. In particular, we have recovered the {beta}-stabilization of the ion temperature gradient mode, transition to collisionless trapped electron mode, and the onset of kinetic ballooning mode as {beta}{sub e} (ratio of electron kinetic pressure to magnetic pressure) increases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Simulation and assessment of ion kinetic effects in a direct-drive capsule implosion experiment

The first simulations employing a kinetic treatment of both fuel and shell ions to model inertial confinement fusion experiments are presented, including results showing the importance of kinetic physics processes in altering fusion burn. A pair of direct drive capsule implosions performed at the OMEGA facility with two different gas fills of deuterium, tritium, and helium-3 are analyzed. During implosion shock convergence, highly non-Maxwellian ion velocity distributions and separations in the density and temperature amongst the ion species are observed. Finally, diffusion of fuel into the capsule shell is identified as a principal process that degrades fusion burn performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

High-order continuum kinetic method for modeling plasma dynamics in phase space

Continuum methods offer a high-fidelity means of simulating plasma kinetics. While computationally intensive, these methods are advantageous because they can be cast in conservation-law form, are not susceptible to noise, and can be implemented using high-order numerical methods. Advances in continuum method capabilities for modeling kinetic phenomena in plasmas require the development of validation tools in higher dimensional phase space and an ability to handle non-cartesian geometries. To that end, a new benchmark for validating Vlasov-Poisson simulations in 3D (x,v x ,v y ) is presented. The benchmark is based on the Dory-Guest-Harris instability and is successfully used to validate a continuum finite volume algorithm. To address challenges associated with non-cartesian geometries, unique features of cylindrical phase space coordinates are described. Preliminary results of continuum kinetic simulations in 4D (r,z,v r ,v z ) phase space are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Studying protein assembly with reversible Brownian dynamics of patchy particles

Assembly of protein complexes like virus shells, the centriole, the nuclear pore complex, or the actin cytoskeleton is strongly determined by their spatial structure. Moreover, it is becoming increasingly clear that the reversible nature of protein assembly is also an essential element for their biological function. Here we introduce a computational approach for the Brownian dynamics of patchy particles with anisotropic assemblies and fully reversible reactions. Different particles stochastically associate and dissociate with microscopic reaction rates depending on their relative spatial positions. The translational and rotational diffusive properties of all protein complexes are evaluated on-the-fly. Because we focus on reversible assembly, we introduce a scheme which ensures detailed balance for patchy particles. We then show how the macroscopic rates follow from the microscopic ones. As an instructive example, we study the assembly of a pentameric ring structure, for which we find excellent agreement between simulation results and a macroscopic kinetic description without any adjustable parameters. This demonstrates that our approach correctly accounts for both the diffusive and reactive processes involved in protein assembly.

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

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
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↗