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 505 records · Page 28

Formation of active sites on transition metals through reaction-driven migration of surface atoms

Adopting low-index single-crystal surfaces as models for metal nanoparticle catalysts has been questioned by the experimental findings of adsorbate-induced formation of subnanometer clusters on several single-crystal surfaces. Here we used density functional theory calculations to elucidate the conditions that lead to cluster formation and show how adatom formation energies enable efficient screening of the conditions required for adsorbate-induced cluster formation. We studied a combination of eight face-centered cubic transition metals and 18 common surface intermediates and identified systems relevant to catalytic reactions, such as carbon monoxide (CO) oxidation and ammonia (NH 3 ) oxidation. We used kinetic Monte Carlo simulations to elucidate the CO-induced cluster formation process on a copper surface. Scanning tunneling microscopy of CO on a nickel (111) surface that contains steps and dislocations points to the structure sensitivity of this phenomenon. Metal-metal bond breaking that leads to the evolution of catalyst structures under realistic reaction conditions occurs much more broadly than previously thought.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultra-high energy density relativistic plasmas from nanostructures: scaling to ultra-high intensities (Final Report)

This project investigated the ultra-high energy density (UHED) regime of matter found in the center of stars, using a compact PW-class laser. These extreme conditions are typically only obtained in the laboratory in the central hot-spot of spherically imploded capsules in inertial confinement fusion experiments driven by the world’s largest lasers. We have shown that near-solid density arrays of aligned nanostructures can be volumetrically heated to multi-keV temperatures by irradiation at relativistic intensity with ultrafast laser pulses of modest energy, opening a path for the generation of UHED plasmas with compact lasers. This new UHED plasma generation approach promises to create an environment with extreme energy densities and degrees of ionization, record conversion of optical laser light into ultrafast x-ray pulses, gigantic magnetic fields and pressures, and directed beams of high energy particles . In this project we achieved record degree of ionization in volumetric heating solid density and near-solid density plasmas. Gold plasmas which spectra is characterized by L shell transition emission from ions with charge up to the Ne-like state, Au+ 72 were generated using ultrafast laser pulses of less than 10 J of energy from a compact laser focused to an intensity of ~ 3x10 21 Wcm -2 . We also conducted measurements to determine the heat penetration depth in Ni nanowire arrays as compared to Ni foil targets by monitoring the line emission of a Co buried tracer underneath a variable amount of Ni. The measurements revealed that the nanowire plasmas are roughly six times larger in depth than the solid density target plasmas. A result of this increased plasma in nanowire arrays is a greatly increased conversion of optical laser light into > 1 KeV x-rays, a record conversion efficiency of 20 %. Critical to the realization of the proposed experiments was the generation of ultrafast laser pulses with ultra-high contrast that can deposit the energy deep into the nanowire arrays before the nanowires explode to create a continuous plasma. Supporting the proposed experiments was the recent demonstration at Colorado State University of a Petawatt-class laser that emits 30 fs pulses at high repetition rates. The experiments combined this unique laser tool with a variety of tailored nanowire arrays fabricated in house and with an extensive suite of diagnostics. The experiments were accompanied by 3-dimensional particle-in-cell simulations and detailed atomic physics simulations with transient kinetics and radiation transport. The proposed research allowed us to continue training Ph.D students and post-docs with broad experience in HEDP Science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic effects of self-organization in low-temperature magnetized plasmas

This final report presents the progress made for the DOE-SC0020623 for the period of performance between January 1 2020 and September 30 2024. The research project is entitled "Kinetic effects of self-organization in low-temperature magnetized plasmas", in which kinetic theory and simulations are used to under stand self-organization patterns.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Substrate-explicit Modeling Tutorial

This tutorial Narrative demonstrates how to apply thermodynamic theory, also known as lambda theory, to convert molecular formulas of compounds (derived from FTICR-MS peaks using Formulatiry and R codes) into stoichiometric and kinetic forms of biogeochemical reactions and how to use the resulting kinetic equations to simulate dynamic conversion of compounds in batch and continuous stirred tank reactors. It will address under what conditions respiration rates are driven by thermodynamics, and how respiration rates respond to the variations in parameters and input variables and how to interpret the results.

54 ENVIRONMENTAL SCIENCES↗

Exploratory Studies of Mechanical Properties, Residual Stress, and Grain Evolution at the Local Regions near Pores in Additively Manufactured Metals

Directed energy deposition (DED) is gaining widespread acceptance in various industrial applications since its unique manufacturing features allow the DED to print metallic parts with very complex geometries. However, DED inevitably generates a lot of internal pores which can limit the widespread applications of the DED technique. The current studies on DED porosity are mostly focused on analyzing pores’ bulk-scale influences on mechanical properties and performances. Since DED pores have a micro-scale existence, with dimensions ranging from a few microns to several hundred microns, it is fundamental to explore the pores’ influences on the micro-scale, including local mechanical properties, residual stress, and grains near pores. However, this important research direction has been neglected. The objective of this work is to fill the above gap in DED porosity research and acquire a fundamental understanding of the role of porosity on a microscopic scale. The authors used nanoindentation approaches to investigate internal pores’ effects on mechanical properties and residual stress in local regions surrounding the pores. In addition, the grains near pores were observed through EBSD, and simulated with the Kinetic Monte Carlo model. The research findings can be provided for DED researchers and industrial practitioners as technical guidance. Most importantly, the research results can work as a good reference for tracing the source of bulk-scale mechanical performances and properties of DED parts with internal pores.

36 MATERIALS SCIENCE↗

DME-Propane Ignition Delay Time Measurements at Mixing Controlled Compression Ignition Engine-Relevant Conditions

The blend of dimethyl ether (DME, CH 3 OCH 3 ) and propane (C 3 H 8 ) is a potentially renewable fuel mixture that has the potential to replace diesel in compression ignition engines. The combination can potentially reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. However, detailed conceptual and simulation studies must be conducted before adopting a new fuel on a compression ignition engine. For these simulations, accurate chemical kinetic models are necessary. However, the validity of chemical kinetic mechanisms in the literature is unknown for mixing controlled compression ignition (MCCI) engine operating conditions. Hence, in this work, we studied the ignition of dimethyl ether (DME) and propane blends in a shock tube at MCCI engine conditions. Ignition delay time (IDT) data was collected behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Undiluted experiments spanning temperatures of 700 to 1100 K and pressures of 55 to 84 bar for various blends (100% CH 3 OCH 3 / 0% C 3 H 8, 100% C 3 H 8 / 0% CH 3 OCH 3 , 60% CH 3 OCH 3 / 40% C 3 H 8 ) of DME and propane were combusted in synthetic air (21% O 2 / 79% N 2 ). Some experiments were conducted at higher pressures (90-120 bar) to understand model performance at these conditions. Comparisons of IDT were made with the predictions of recent chemical kinetic mechanisms for DME-propane mixture, including the Aramco3.0, NUIG, and Dames et al. mechanisms. All mechanisms overpredicted IDT compared to experimental values. Sensitivity analysis was conducted with Dames et al. model, and critical reactions sensitive to IDT of DME-propane mixture near 100 bar are outlined.

Mohammed, Zuhayr Pasha↗

Toward closure between predicted and observed particle viscosity over a wide range of temperatures and relative humidity

Abstract. Atmospheric aerosols can exist in amorphous semi-solid or glassy phase states whose viscosity varies with atmospheric temperature and relative humidity. The temperature and humidity dependence of viscosity has been hypothesized to be predictable from the combination of a water–organic binary mixing rule of the glass transition temperature, a glass-transition-temperature-scaled viscosity fragility parameterization, and a water uptake parameterization. This work presents a closure study between predicted and observed viscosity for sucrose and citric acid. Viscosity and glass transition temperature as a function of water content are compiled from literature data and used to constrain the fragility parameterization. New measurements characterizing viscosity of sub-100 nm particles using the dimer relaxation method are presented. These measurements extend the available data of temperature- and humidity-dependent viscosity to −28 ∘C. Predicted relationships agree well with observations at room temperature and with measured isopleths of constant viscosity at ∼107 Pa s at temperatures warmer than −28 ∘C. Discrepancies at colder temperatures are observed for sucrose particles. Simulations with the kinetic multi-layer model of gas–particle interactions suggest that the observed deviations at colder temperature for sucrose can be attributed to kinetic limitations associated with water uptake at the timescales of the dimer relaxation experiments. Using the available information, updated equilibrium phase-state diagrams (-80∘C

54 ENVIRONMENTAL SCIENCES↗

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Combined TDLAS and chemiluminescence imaging in a flat flame burner operated with NH3/H2 blends

Ammonia is a promising hydrogen carrier due to its favorable storage and transport characteristics. However, its direct use in combustion systems is limited by low flammability and potential for high nitrogen oxide emissions. To better understand ammonia combustion, researchers conducted experiments using a flat flame burner and measured species profiles using tunable-diode-laser-absorption-spectroscopy and chemiluminescence imaging. They tested three flame conditions with different ammonia-hydrogen blends and oxygen levels. The results were compared to simulations using various kinetic mechanisms, including one that accounts for excited species chemistry. The goal is to provide direct information about species profiles in a simple system, isolating chemical kinetics from fluid dynamic effects, which can inform the development of more efficient and low-emission combustion systems using ammonia.

ammonia combustion↗

Bottom-Up Simulation, Reconstruction, and Quantification of Macromolecule Sequences from Experimental Polymerizations

Motivated by the canonical sequence–structure–function paradigm, tools to characterize chemical patterning in natural biomacromolecules, from proteins to nucleic acids, have grown exponentially in recent years. However, analogous strategies for synthetic macromolecules remain in nascent stages, complicated by sequence polydispersity and analytical limitations. To address this, we have developed a comprehensive and open-source Python package, PRISM (polymer rate insights and sequence modeling), an end-to-end workflow that provides a path from experimental kinetics measurements to quantitative and qualitative metrics for describing chemical patterning in stochastic polymers. First, a numerical integration strategy was constructed to simulate and fit experimental data from reversible addition–fragmentation chain transfer (RAFT) polymerization kinetics, enabling the facile estimation of relevant reactivity ratios. These ratios were then used in a mechanism-specific stochastic kinetic simulation strategy to simulate sequence ensembles corresponding to model systems spanning experimental copolymers, classes of statistical polymers (e.g., alternating, block, and gradient), and multiblock copolymers. Lastly, inspired by sequence homology metrics from bioinformatics, we introduce visualization strategies and quantitative metrics to facilitate comparisons of different sequence ensembles. As the sequence–structure–function paradigm becomes increasingly central in de novo design of synthetic macromolecules, this toolkit provides a first step toward accurate and representative sequence description and featurization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field

Magnetic reconnection is believed to be one of the most fundamental physical processes in laboratory and space plasmas. Fundamental physics of collisionless reconnection under general conditions of finite guide magnetic field, realistic electron-to-ion mass ratio, and in 3D geometries is still not adequately understood. The major goal of this project is to study the physics of collisionless magnetic reconnection in a current sheet under a finite guide field and with a realistic ion-to-electron mass ratio mi/me. Under the support of this DoE grant, we have made significant progress in (1) the investigation of linear and nonlinear current sheet instabilities under a wide range of guide field using the 3-D gyro kinetic electron and fully kinetic ion (GeFi) particle simulation model and (2) the investigation of fully nonlinear reconnection process with a strong guide field with the GeFi simulations. Furthermore, (3) we have extended the study to the realistic geometry of the magnetopause using 3D global hybrid simulations. This study, from local to global scales and for 3-D physics at plasma boundaries, has provides advanced theoretical understanding of the fundamental physics of magnetic reconnection in laboratory and space plasma regimes on multi-scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Combined TDLAS and chemiluminescence imaging in a flat flame burner operated with NH3/H2 blends

Ammonia is viewed as a viable hydrogen carrier due to favorable storage and transport characteristics. While it can be re-converted to hydrogen at point-of-use via thermal catalytic cracking, direct utilization in combustion systems can result in reduced costs and improved efficiency. A major barrier to this approach is the low flammability and potential for high nitrogen oxide emissions, driven by fuel-bound nitrogen and complex kinetic pathways. While a number of kinetic mechanisms currently exist for simulating ammonia combustion, a major need continues to be direct information about species profiles in easy-to-model systems capable of isolating chemical kinetics from multi-dimensional fluid dynamic effects. This paper reports on recent species measurements made in a flat flame burner using a combined tunable-diode-laser-absorption-spectroscopy (TDLAS) and chemiluminescence imaging approach. Three flame conditions were included representative of NH3/H2 blends, partially cracked NH3 (inc. N2), and 100% NH3 with enhanced air (30% oxygen). Two NIR distributed feedback (DFB) laser diodes were used to determine H2O concentration and temperature (via ratio thermometry) at various radial and axial positions, after which an inverse Abel transform was used to infer centerline values. Companion images were collected using a Princeton Instrument PI-MAX intensified camera equipped with a 105mm UV lens and multiple filter sets targeting OH*, NH*, and NH2* emission. Results were compared to companion Cantera burner-stabilized flame simulations using various kinetic mechanisms. A recent mechanism including excited species chemistry was also included, to investigate whether excited and ground state profiles exhibit significant differences.

ammonia combustion↗

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a synthetic phase contrast imaging diagnostic for turbulence studies at Wendelstein 7-X

Abstract We present a synthetic phase contrast imaging (PCI) diagnostic for studying turbulence at the Wendelstein 7-X (W7-X) stellarator. We first describe the implemented instrument response model, which captures diffraction effects, detector noise, and the long-wavelength cutoff due to the phase plate of the PCI system. To verify the instrument response model, we show that it is capable of reproducing the PCI signal generated by the sound wave speaker used for calibration at W7-X. Next, we discuss the calculation of synthetic PCI signals based on the global, nonlinear gyrokinetic codes GENE-3D and EUTERPE, including results from some of the first stellarator simulations of this type with kinetic electrons (KEs) in GENE-3D. While the simulations used in this work lack a neoclassical radial electric field, which is crucial for reproducing experimental PCI signals, they do indicate that the dominant rotation direction and velocities of the turbulent fluctuations can be inferred from the wave number-frequency spectra of the PCI signals, as expected. The synthetic PCI wave number spectra are further shown to be similar to those of the line-integrated fluctuating electron density, with distinct differences between adiabatic and KE simulations, explainable by previously published turbulence models. For example, the wave number spectra of all adiabatic electron simulations analyzed here follow a power law with an exponent close to −5 for sufficiently large wave numbers. This indicates that universal features of electron density turbulence at W7-X may be studied using the PCI system.

EUTERPE↗

A pyrolysis study of allylic hydrocarbon fuels

The pyrolysis of selected C 3 –C 5 allylic hydrocarbons has been studied using a single-pulse shock tube. A new single-pulse shock tube has been designed and constructed by recommissioning an existing conventional shock tube. This facility enables the investigation of high-temperature chemical kinetics with an emphasis on combustion chemistry. The modifications performed on the existing shock tube are described, and the details of the sampling system to analyze the species concentration using a gas chromatography-mass spectrometry-flame ionization detection (GC-MS with a flame ionization detector) system are also provided. This facility is characterized and validated by performing cyclohexene pyrolysis experiments. Furthermore, the performance of the shock tube is demonstrated by reproducing previous literature measurements on the pyrolysis of isobutene. Postvalidation, this setup is used to study the pyrolysis of trans-2-butene and 2-methyl-2-butene (2M2B). A newly developed mechanism, NUIGMech1.0, is used to simulate the experimental data of propene, isobutene, 2-butene, and 2M2B, allylic hydrocarbon fuels. A description using two different kinetic simulation approaches is provided using our isobutene experiments as a reference. We found no significant differences between the two methods. Additionally, the contribution of different reaction classes on fuel consumption is detailed and the influence of geometry on fuel consumption and first aromatic ring: benzene is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗