Atomistic Multi-Lattice Kinetic Monte Carlo (Kmc) Modeling of Hyperthermal Oxidation of Multi-Layer Graphene
Explore the source record for details and available documents.
Sharon Edward↗
SEARCH · Search NASA
Results for “Kinetic Modelling”
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.
Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)
Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.
Convex Optimization↗
Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)
Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.
GN&C↗
Development of Kinetic PIC-DSMC Model for Breakdown in the Presence of a Dielectric.
Abstract not provided.
Moore, Christopher Hudson↗
Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields
We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann-Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.
79 ASTRONOMY AND ASTROPHYSICS↗
A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications
Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.
ab initio↗
Non-Fermi Liquids from Kinetic Constraints in Tilted Optical Lattices
Here we study Fermi-Hubbard models with kinetically constrained dynamics that conserves both total particle number and total center of mass, a situation that arises when interacting fermions are placed in strongly tilted optical lattices. Through a combination of analytics and numerics, we show how the kinetic constraints stabilize an exotic non-Fermi liquid phase described by fermions coupled to a gapless bosonic field, which in several respects mimics a dynamical gauge field. This offers a novel route towards the study of non-Fermi liquid phases in the precision environments afforded by ultracold atom platforms.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Predictive Modeling of Fast-Curing Thermosets in Nozzle-Based Extrusion
This work presents an approach to modeling the dynamic spreading and curing behavior of thermosets in nozzle-based extrusions. Thermosets cover a wide range of materials, some of which permit low-temperature processing with subsequent high-temperature and high-strength working properties. Extruding thermosets may overcome the limited working temperatures and strengths of conventional thermoplastic materials used in additive manufacturing. This project aims to produce technology for the fabrication of thermoset-based structures leveraging advances made in nozzle-based extrusion, such as fused deposition modeling (FDM), material jetting, and direct writing. Understanding the synergistic interactions between spreading and fast curing of extruded thermosetting materials will provide essential insights for applications that require accurate dimensional controls, such as additive manufacturing [1], [2] and centrifugal coating/forming [3]. Two types of thermally curing thermosets -- one being a soft silicone (Ecoflex 0050) and the other being a toughened epoxy (G/Flex) -- served as the test materials in this work to obtain models for cure kinetics and viscosity. The developed models align with extensive measurements made with differential scanning calorimetry (DSC) and rheology. DSC monitors the change in the heat of reaction, which reflects the rate and degree of cure at different crosslinking stages. Rheology measures the change in complex viscosity, shear moduli, yield stress, and other properties dictated by chemical composition. By combining DSC and rheological measurements, it is possible to establish a set of models profiling the cure kinetics and chemorheology without prior knowledge of chemical composition, which is usually necessary for sophisticated mechanistic modeling. In this work, we conducted both isothermal and dynamic measurements with both DSC and rheology. With the developed models, numerical simulations yielded predictions of diameter and height of droplets, along with width and height of extruded lines cured at varied temperatures. Experimental results carried out on a goniometric platform and a nozzle-based 3D printer showed agreement with the numerical simulations. Finally, this presentation will show how the models are adaptable to the planning of tool paths and designs in additive manufacturing.
Xie, Jingjin↗
The theory of kinetic effects on resistive wall mode stability in tokamaks
Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Novel Approach to Modeling Biomass Pyrolysis in a Fluidized Bed Reactor
A novel approach to simulating biomass pyrolysis in a fluidized bed of mostly inert material is presented. The bed is assumed to be externally heated, although simulation of autothermal operation with partial oxidation of the products is a future objective. Combining pyrolysis reaction kinetics developed by the CRECK Modeling Group and a previously published component properties model, material and energy balances are closed by tracking the residence time of biomass particles without regard for their exact spatial distribution. The model is used to simulate a pilot-scale fluidized bed pyrolysis process being developed at Iowa State University, and model predictions are compared with experimental results for red oak and corn stover feedstocks. The results are in general agreement, with the model typically predicting more noncondensable gas and less low-boiling liquid products. The differences are ascribed to the limitations of the kinetics model. This model provides a rigorous energy balance for fluidized bed pyrolysis processes that can be incorporated into commercial flowsheet models with the capacity for future addition of partial oxidation reactions and implementation of improved kinetics models.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes
We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.
97 MATHEMATICS AND COMPUTING↗
Experimental and modeling study of chemical-based strategies for mitigating dust formation in fusion reactors [Plasma Flow Reactor Chemical Kinetics: A Method Toward Optimizing Nano-Manufacturing and Laser Inertial Fusion Energy Strategies]
We studied carbon/hydrogen/oxygen chemical kinetics at time scales and thermal conditions relevant to fusion energy applications using a custom-built plasma flow reactor to investigate chemical-based strategies for eliminating carbon dust formation in fusion reactors. Acetylene and oxygen gases under varying conditions of initial concentrations are injected into an inductively coupled argon plasma where complete molecular dissociation occurs. The evolution of chemical species is investigated along the plasma flow reactor as a function of temperature and residence time. Atomized species of C, H, and O cool from 5000 K to 1000 K within 30 ms at atmospheric pressures. We employed optical emission and infrared absorption spectroscopy to measure the reaction intermediates (e.g. C 2 ) and products (e.g. C 2 H 2 ). Chemical equilibrium models are inadequate to describe the evolution of carbon molecular products, and thus a chemical kinetics model is developed. In both experiments and kinetic modelling, we find that the addition of oxygen in 1:1 proportion to carbon strongly favors the formation of CO, preventing the formation of acetylene (an important soot precursor) in less than 10 milliseconds. The kinetics model is also used to perform a reaction sensitivity and a rate of production analyses to identify the rate determining steps and the major chemical pathways that control the acetylene production/consumption. The results demonstrate the feasibility of chemical-based strategies for eliminating the formation of carbonaceous particles in fusion energy reactors.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Kinetic inductance photodetectors based on nonequilibrium response in superconducting thin-film structures
We modeled the kinetic inductance response of detectors made from disordered superconducting Nb, NbC, and MoRe films.
kinetics↗
Three-dimensional stochastic model of actin–myosin binding in the sarcomere lattice
The effect of molecule tethering in three-dimensional (3-D) space on bimolecular binding kinetics is rarely addressed and only occasionally incorporated into models of cell motility. The simplest system that can quantitatively determine this effect is the 3-D sarcomere lattice of the striated muscle, where tethered myosin in thick filaments can only bind to a relatively small number of available sites on the actin filament, positioned within a limited range of thermal movement of the myosin head. Here we implement spatially explicit actomyosin interactions into the multiscale Monte Carlo platform MUSICO, specifically defining how geometrical constraints on tethered myosins can modulate state transition rates in the actomyosin cycle. The simulations provide the distribution of myosin bound to sites on actin, ensure conservation of the number of interacting myosins and actin monomers, and most importantly, the departure in behavior of tethered myosin molecules from unconstrained myosin interactions with actin. In addition, MUSICO determines the number of cross-bridges in each actomyosin cycle state, the force and number of attached cross-bridges per myosin filament, the range of cross-bridge forces and accounts for energy consumption. At the macroscopic scale, MUSICO simulations show large differences in predicted force-velocity curves and in the response during early force recovery phase after a step change in length comparing to the two simplest mass action kinetic models. The origin of these differences is rooted in the different fluxes of myosin binding and corresponding instantaneous cross-bridge distributions and quantitatively reflects a major flaw of the mathematical description in all mass action kinetic models. Consequently, this new approach shows that accurate recapitulation of experimental data requires significantly different binding rates, number of actomyosin states, and cross-bridge elasticity than typically used in mass action kinetic models to correctly describe the biochemical reactions of tethered molecules and their interaction energetics.
60 APPLIED LIFE SCIENCES↗
Modeling nitrogen chemistry in combustion
Understanding of the chemical processes that govern formation and destruction of nitrogen oxides (NO x ) in combustion processes continues to be a challenge. Even though this area has been the subject of extensive research over the last four decades, there are still unresolved issues that may limit the accuracy of engineering calculations and thereby the potential of primary measures for NO x control. In this review our current understanding of the mechanisms that are responsible for combustion-generated nitrogen-containing air pollutants is discussed. The thermochemistry of the relevant nitrogen compounds is updated, using the Active Thermochemical Tables (ATcT) approach. Rate parameters for the key gas-phase reactions of the nitrogen species are surveyed, based on available information from experiments and high-level theory. The mechanisms for thermal and prompt-NO, for fuel-NO, and NO formation via NNH or N 2 O are discussed, along with the chemistry of NO removal processes such as reburning and Selective Non-Catalytic Reduction of NO. Each subset of the mechanism is evaluated against experimental data and the accuracy of modeling predictions is discussed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Resolving Discrepancies Between Theory and Experiment for the NCN + H Reaction
Flame modeling studies have highlighted the role of branching in the kinetics for the prompt NO switch reaction, NCN + H. In this reaction, there is strong competition between the CH + N2 (R1) and HCN + N (R2) product channels. Increased branching towards reaction R2 promotes the subsequent formation of Fenimore or Prompt NO. Recent direct shock tube studies on this complex reaction conclude that branching to reaction R2 is much more favorable than the predicted branching from theoretical studies. The experimentally predicted prompt NO switch temperature (TS), i.e., the temperature at which k2/(k1+k2) = 0.5, is at 1670 K, in contrast to theory that predicts Ts > 3200 K. The present theory/modeling work was initiated to highlight potential causes for this significant discrepancy between experiment and theory. Our analysis of the simulations reported in the experimental studies indicate that including a fall-off representation for C2H5 radical dissociation (with C2H5 produced from C2H5I dissociation and acting as an H-atom source) has a noticeable influence on the simulations for NCN decay and HCN formation in the shock tube studies. With the inclusion of a theory-based pressure fall-off representation for C2H5 dissociation, we show that this radical persists for much longer timescales even at the high-temperatures in the shock tube studies. This persistence then facilitates the rapid conversion of reactive H-atoms to CH3 radicals at the shock tube conditions via the H + C2H5 → CH3 + CH3 reaction. An updated theoretical analysis for this classic addition-elimination reaction is also provided in this work. The rapid formation of CH3 radicals necessitates the inclusion of additional reactions to adequately simulate the NCN and HCN data. Our simulations conclude that, with the inclusion of these additional reactions, the NCN and HCN profiles can be reasonably well simulated by the most recent theoretical predictions.
Ethyl Radicals↗