Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic separation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dual phase patterning during a congruent grain boundary phase transition in elemental copper

The phase behavior of grain boundaries can have a strong influence on interfacial properties. Little is known about the emergence of grain boundary phases in elemental metal systems and how they transform. Here, we observe the nanoscale patterning of a grain boundary by two alternating grain boundary phases with distinct atomic structures in elemental copper by atomic resolution imaging. The same grain boundary phases are found by computational grain boundary structure search indicating a first-order transformation. Finite temperature atomistic simulations reveal a congruent, diffusionless transition between these phases under ambient pressure. The patterning of the grain boundary at room temperature is dominated by the grain boundary phase junctions separating the phase segments. Our analysis suggests that the reduced mobility of the phase junctions at low temperatures kinetically limits the transformation, but repulsive elastic interactions between them and disconnections could additionally stabilize the pattern formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
DOE OSTI2023

A Primer on Nuclear "Recoil" Data

For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗
DOE OSTI2020

Achieving Balanced Crystallization Kinetics of Donor and Acceptor by Sequential-Blade Coated Double Bulk Heterojunction Organic Solar Cells

Sequential deposition has great potential to achieve high performance in organic solar cells due to the resulting well-controlled vertical phase separation. In this work, double bulk heterojunction organic solar cells are fabricated by sequential-blade cast in ambient conditions. Probed by the in situ grazing incidence X-ray diffraction and in situ UV–vis absorption measurements, the seq-blade system exhibits a different tendency from each of the binary films during the film formation process. Due to the extensive aggregation of FOIC, the binary PBDB-T:FOIC film displays a strong and large phase separation, resulting in low current density (J sc) and unsatisfactory power conversion efficiency. In the seq-blade cast system, the bottom layer PBDB-T:IT-M produces many crystal nuclei for the top layer PBDB-T:FOIC, so the PBDB-T molecules are able to crystallize easily and quickly. Balanced crystallization kinetics between polymer and small molecule and an ideal percolation network in the film are observed. In addition, the balanced crystallization kinetics are favorable toward realizing lower recombination loss through charge transport processes.

14 SOLAR ENERGY↗
DOE OSTI2024

Turning Normal to Abnormal: Reversing CO 2 /C2–Hydrocarbon Selectivity in HKUST–1

Metal–organic frameworks (MOFs) can efficiently purify hydrocarbons from CO 2 , but their rapid saturation, driven by preferential hydrocarbon adsorption, requires energy-intensive adsorption–desorption processes. To address these challenges, an innovative approach is developed, enabling control over MOF flexibility through densification and defect engineering, resulting in an intriguing inverse CO 2 /C2 hydrocarbon selectivity. Here, in this study, the densification process induces the shearing of the crystal lattice and contraction of pores in a defective CuBTC MOF. These changes have led to a remarkable transformation in selectivity, where the originally hydrocarbon-selective CuBTC MOF becomes CO 2 -selective. The selectivity values for densified CuBTC are significantly reversed when compared to its powder form, with notable improvements observed in CO 2 /C 2 H 6 (4416 vs 0.61), CO 2 /C 2 H 4 (15 vs 0.28), and CO 2 /C 2 H 2 (4 vs 0.2). The densified material shows impressive separation, regeneration, and recyclability during dynamic breakthrough experiments with complex quinary gas mixtures. Simulation studies indicate faster CO 2 passage through the tetragonal structure of densified CuBTC compared to C 2 H 2 . Experimental kinetic diffusion studies confirm accelerated CO 2 diffusion over hydrocarbons in the densified MOF, attributed to its small pore window and minimal interparticle voids. This research introduces a promising strategy for refining existing and future MOF materials, enhancing their separation performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kinetic versus Thermodynamic Stability of LLZO in Contact with Lithium Metal

Li 7 La 3 Zr 2 O 12 (LLZO) garnet-based oxides are a promising class of solid electrolytes used as the separator in all-solid-state batteries (ASSBs). While LLZO is considered to have a wide electrochemical stability window, its intrinsic stability in contact with lithium metal is not sufficiently well understood, and there is still a debate on the key question of whether LLZO does or does not form passivation layers before and during cycling. Utilizing both in situ and operando X-ray photoelectron spectroscopy techniques, in this work we reveal the presence of a kinetic barrier to the reduction of LLZO by Li metal, with the extent of oxygen-deficient interphase (ODI) formation depending sensitively on the energetics of Li metal arriving at the Li vertical bar LLZO interface. Despite the clear presence of a kinetic barrier to reduction, the electrochemical response of the Li vertical bar LLZO interface is unchanged by the presence of the ODI, indicating that ODI formation during electrochemical cycling does not hinder charge transfer across the Li vertical bar LLZO interface. Overall, these results reveal that the reactivity of LLZO with Li metal depends not only on the material properties of the adjoining phases (i.e., surface purity and active contact) and their resulting thermodynamic stability but also on the energy input at the interface and the resulting reaction kinetics. Furthermore, the presence of a kinetic barrier to reduction highlights the additional complexities governing the reactivity of solid-state interfaces in ASSBs and underscores the importance of operando characterization of interfacial stability to design more robust, high-performance protection strategies for solid electrolytes in contact with reactive electrodes.

25 ENERGY STORAGE↗
DOE OSTI2022

bolt: A GPU-Accelerated Solver for Marginally Collisional Kinetic Physics Using a High- Resolution Constrained Transport Scheme

Collisionless and marginally collisional plasmas are very expensive to simulate due to high dimensionality and large scale separations, especially in global problems with additional length scales. bolt addresses this challenge by achieving high performance on GPUs using the arrayfire performance portability library. We have demonstrated accuracy on a range of test problems, and the bolt code paper is currently in preparation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evolution of the role of molybdenum in duplex stainless steels during thermal aging: From enhancing spinodal decomposition to forming heterogeneous precipitates

Thermal aging of ferrite in duplex stainless steels leads to a complex microstructural evolution that can give way to large deleterious changes in mechanical properties. The microstructural evolution mechanisms that happen concurrently are strongly linked to minor changes in composition. By characterizing duplex stainless steels with various compositions over a range of aging conditions, it was revealed that adding a modest amount of Mo to cast austenitic stainless steel (such as changing from grade CF3 to grade CF3M) not only increases the kinetics of spinodal decomposition and Ni-Si-Mn-rich G-phase precipitation, but it also forms a separate Mo-rich phase. At shorter equivalent aging times, Mo segregates from Fe and towards Cr accelerating spinodal decomposition. At longer aging times, Mo enriches in the Ni-Si-Mn clusters before phase separating into an adjoined Mo-rich cluster. We found that these phase changes and precipitation events result in a rather steep drop in fracture and impact toughness. These results will inform the design of similar alloys for use at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fast charging of Li-ion cells: Effect of separator membranes and mapping of “safe lines” to avoid Li plating

An experimental method is introduced that uses microprobe Li/Cu reference electrodes to map "safe" charging regimes to avoid Li metal plating on graphite anodes at rates to 8C. With this method, a Li-ion cell closely approaches the Li plating condition but avoids descending into this regime. This "safe" approach is used to demonstrate effects of the microporous separator on Li plating in cells consisting of a layered nickel-manganese -cobalt oxide cathode and graphite anode. We argue that most of the difference in behaviors for different separators arises from time delayed Li + ion percolation through the membrane. In conclusion, our study suggests that explicitly treating kinetics of phase transitions in the lithiated graphite anode is essential for electrochemical models predicting "safe" lines for fast charging cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Physics of beam-driven ion cyclotron emission in the large plasma device

Abstract Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v ⊥ / v A , is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas ‘might as well’ be occurring in a cylinder.

Samant, O. (ORCID:0000000226055363)↗
DOE OSTI2025

Towards modelling AR Sco: calibration – reproducing high-energy pulsar emission and testing convergence to Aristotelian electrodynamics

In recent years, kinetic simulations have been crucial to further our understanding of pulsar electrodynamics. Yet, due to the large-scale separation between the gyro-period and the stellar rotation period, resolving the particle gyration has been computationally unfeasible for realistic pulsar parameters. The main aim of this work is comparing our gyro-phase-resolved model with a gyro-centric pulsar model, where our model solves the general equations of motion with included radiation reaction using a higher order numerical solver with adaptive time-steps. Specifically, we aim to (i) reproduce a pulsar’s high-energy emission maps, namely one with 10 per cent of the surface B-field strength of Vela, and the spectra produced by an independent gyro-centric pulsar emission model; and (ii) test convergence of these results to the radiation-reaction limit of Aristotelian electrodynamics. (iii) Additionally, we identify the effect that a large $E_{\parallel }$-field has on the trajectories and radiation calculations. We find that we can reproduce the curvature radiation emission maps and spectra well, using 10 per cent field strengths of the Vela pulsar and injecting our particles at a higher altitude in the magnetosphere. Using sufficiently large $E_{\parallel }$-fields, our numeric results converge to the analytic radiation-reaction limit trajectories. Additionally, we illustrate the importance of accounting for the $\mathbf {E}\times \mathbf {B}$-drift in the particle trajectories and radiation calculations, validating the Harding and collaborators’ model approach. Lastly, we found that our model deals very well with the high-radiation-reaction and high-field regimes present in pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Preliminary Thermoluminescent Dosimeter Glow Curve Analysis with Automated Glow Peak Identification for LiF Mg,Ti

When appropriately analyzed, thermoluminescent dosimeter glow curve analysis allows for improved quantification of thermoluminescent material behavior while flagging abnormalities. The mathematical separation of a glow curve into contributions from energetically unique trap states, or glow curve analysis, may be used to remove undesired effects of signal fading for complex materials. A generalized glow curve analysis software for the separation of glow curves is presented in this paper. Written in C ++ , the software uses the first-order kinetics model with automatic peak identification. The automatic identification of peaks is achieved through a unique peak-finding algorithm. Here, the program was performance tested using experimental glow curve data from LiF:Mg,Ti, and comparative results are presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Multifunctional electrochemical memory stabilized by phase coexistence

Our growing computing needs, especially in applications that heavily rely on artificial intelligence (AI), motivate a search for new components that could substantially augment the performance of general-purpose digital computers. Beyond ON/OFF switching, new components with linear multistate analog resistive tuning, nonlinear volatile switching, spiking, oscillatory, stochastic and other complex functionalities could enable highly efficient neuromorphic computing schemes for AI information processing. Compared to the extreme multifunctionality of biological neurons, realizing all the above characteristics in a single, scalable analog component remains a grand challenge. Here we investigate electrochemical gating combined with localized thermal activation to program and switch a single, vertically integrated and dimensionally scaled electrothermal chemical random access memory (ETCRAM) with a channel and reservoir composed of phase-separated vanadium oxide. Closely related to electrochemical RAM (ECRAM), ETCRAM uses an integrated gate-heater electrode to overcome kinetic barriers that help retain states at ambient temperatures. In addition to synapse-like stable and programmable analog resistance states arising from redox-tunable phase coexistence, a single component exhibits neuron-like nonlinear conductance switching with a tunable threshold and self-driven dynamics owing to the thermally driven metal-insulator phase transition in vanadium dioxide. More broadly, we demonstrate that electrochemically stabilized phase coexistence could unlock analog electronics with novel functionality, stability, reconfigurability, and scalability.

Oh, Sangheon [Sandia National Lab. (SNL-CA), Liver↗
DOE OSTI2024

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Primer terminal ribonucleotide alters the active site dynamics of DNA polymerase η and reduces DNA synthesis fidelity

DNA polymerases catalyze DNA synthesis with high efficiency, which is essential for all life. Extensive kinetic and structural efforts have been executed in exploring mechanisms of DNA polymerases, surrounding their kinetic pathway, catalytic mechanisms, and factors that dictate polymerase fidelity. Recent time-resolved crystallography studies on DNA polymerase η (Pol η) and β have revealed essential transient events during the DNA synthesis reaction, such as mechanisms of primer deprotonation, separated roles of the three metal ions, and conformational changes that disfavor incorporation of the incorrect substrate. DNA-embedded ribonucleotides (rNs) are the most common lesion on DNA and a major threat to genome integrity. While kinetics of rN incorporation has been explored and structural studies have revealed that DNA polymerases have a steric gate that destabilizes ribonucleotide triphosphate binding, the mechanism of extension upon rN addition remains poorly characterized. Using steady-state kinetics, static and time-resolved X-ray crystallography with Pol η as a model system, we showed that the extra hydroxyl group on the primer terminus does alter the dynamics of the polymerase active site as well as the catalysis and fidelity of DNA synthesis. During rN extension, Pol η error incorporation efficiency increases significantly across different sequence contexts. Finally, our systematic structural studies suggest that the rN at the primer end improves primer alignment and reduces barriers in C2'-endo to C3'-endo sugar conformational change. Overall, our work provides further mechanistic insights into the effects of rN incorporation on DNA synthesis.

59 BASIC BIOLOGICAL SCIENCES↗