Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixture models”

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 307 records · Page 17

Automatic Volume Balancing for Online Refueling in Molten Salt Reactor Simulations with SCALE

This work introduces recently implemented capabilities in the SCALE code system’s TRITON reactor physics sequence that improve molten salt reactor (MSR) modeling: (1) a volume balancing option, which automatically balances the volumes of the fed material with a corresponding material removal to maintain fixed mixture volumes in the neutron transport model and accurate densities, and (2) continuous feed from mixtures, which enables users to define material feed streams directly from salt mixture definitions rather than individual nuclides. The new capabilities were demonstrated via SCALE/TRITON simulations of two representative MSR concepts. Simulations of the 180 MWth molten chloride fast reactor—a system with a fixed fuel salt volume in the reactor core—applied continuous refueling with fresh fuel salt. The results confirm approximately constant reactivity when the new volume balancing capability is used. Additionally, excellent agreement with a manual feed-and-drain method for volume balancing further verified the new implementation. Simulations of the 400 MWth EIRENE reactor, an integral MSR in which the salt volume may grow over time within the reactor core, demonstrated the capability to represent growing salt volume within the TRITON depletion calculation. Compared with reference solutions, the results show consistent trends in reactivity and isotopic evolution, confirming that these new user-friendly capabilities provide accurate, physically consistent approaches for modeling MSRs in SCALE.

Elzohery, Rabab [ORNL] (ORCID:0000000160043633)↗

High-throughput screening tools facilitate calculation of a combined exposure-bioactivity index for chemicals with endocrine activity

High-throughput and computational tools provide a new opportunity to calculate combined bioactivity of exposure to diverse chemicals acting through a common mechanism. We used high throughput in vitro bioactivity data and exposure predictions from the U.S. EPA’s Toxicity and Exposure Forecaster (ToxCast and ExpoCast) to estimate combined estrogen receptor (ER) agonist activity of non-pharmaceutical chemical exposures for the general U.S. population. High-throughput toxicokinetic (HTTK) data provide conversion factors that relate bioactive concentrations measured in vitro (µM), to predicted population geometric mean exposure rates (mg/kg/day). These data were available for 22 chemicals with ER agonist activity and were estimated for other ER bioactive chemicals based on the geometric mean of HTTK values across chemicals. For each chemical, ER bioactivity across ToxCast assays was compared to predicted population geometric mean exposure at different levels of in vitro potency and model certainty. Dose additivity was assumed in calculating a Combined Exposure-Bioactivity Index (CEBI), the sum of exposure/bioactivity ratios. Combined estrogen bioactivity was also calculated in terms of the percent maximum bioactivity of chemical mixtures in human plasma using a concentration-addition model. Estimated CEBIs vary greatly depending on assumptions used for exposure and bioactivity. In general, CEBI values were <1 when using median of the estimated general population chemical intake rates, while CEBI were ≥1 when using the upper 95th confidence bound for those same intake rates for all chemicals. Concentration-addition model predictions of mixture bioactivity yield comparable results. Based on current in vitro bioactivity data, HTTK methods, and exposure models, combined exposure scenarios sufficient to influence estrogen bioactivity in the general population cannot be ruled out. Future improvements in screening methods and computational models could reduce uncertainty and better inform the potential combined effects of estrogenic chemicals.

54 ENVIRONMENTAL SCIENCES↗

On multidomain multiscale modeling and simulation of a novel partial pressure and temperature swing adsorptive reactor ($\mathrm{PPTSAR}$) with application to the water gas shift reaction

Here in this work, a dynamic, multidomain, multiscale model of a novel Partial Pressure and Temperature Swing Adsorptive Reactor (PPTSAR) Process is developed and simulated for a Water Gas Shift Reaction (WGSR) application. The proposed model couples the reactor wall, catalyst/adsorbent pellet, and fluid mixture domains at their respective scales to comprehensively model mass, momentum, and energy transport phenomena occurring within each domain. The resulting model elucidates the PPTSAR’s start-up, transient, and long-term, periodic, dynamic behavior, during both adsorption/reaction and desorption operating modes. In the WGSR application considered, the use of two PPTSARs is shown to allow continuous operation that delivers CO conversion greater than 95% with simultaneous complete CO 2 capture.

42 ENGINEERING↗

Optimization of Refrigerant Compositions for Low-GWP Refrigerant Mixtures Using Segment-by-segment Heat Exchanger and Detailed System Models

The recently introduced hydrofluoroolefin (HFO) refrigerants, including R1234yf and R1234ze(E), have significantly lower global warming potentials (GWPs) than traditional hydrofluorocarbon (HFC) refrigerants like R410A. However, prior tests show that direct drop-in of pure R1234yf or R1234ze(E) into equipment designed for R410A results in a decrease in heat exchanger capacity and the system coefficient of performance. The primary reason is the lower in-tube heat transfer performance of R1234yf and R1234ze(E) compared with that of R410A. To address this issue, previous studies have mixed the mildly flammable HFC R32 with HFOs to improve system performance, with HFC R125 also added to suppress flammability. Previous studies selected compositions based on simple cycle analyses and did not consider modifications of the heat exchanger circuitry configuration to adapt to the new refrigerants. This study presents a novel multi-objective optimization approach to design a refrigerant composition that maximizes energy efficiency within flammability and GWP limits. The approach in this work simultaneously optimizes mixture composition and heat exchanger circuitry configuration. A case study on a rooftop unit indicates that, compared with mixture-only optimization, simultaneous optimization of mixture and heat exchanger circuitry yields a 5.9% improvement in cycle efficiency and a 48.6% reduction in refrigerant flammability with a GWP of 268. Circuitry optimization using refrigerants with different temperature glides shows that the larger the temperature glide is, the larger EER improvement is obtained. The results show that zeotropic blends with a large temperature glide are more sensitive to the refrigerant circuitry than pure refrigerants and may suffer significant performance degradation with subpar heat exchanger circuitry design. The proposed optimization approach is generally applicable to mixtures with any number of components. Using this approach to design a HVAC system can yield higher system efficiency within flammability and GWP constraints.

Li, Zhenning↗

Graph neural networks for CO 2 solubility predictions in Deep Eutectic Solvents

Deep Eutectic Solvents (DESs) are a promising class of solvents for CO 2 capture. DESs are complex mixtures that can be designed to optimize CO solubility and overall capture process efficiency. However, the vast design landscape of DES mixtures makes experimental investigation prohibitive; as such, there is a need for computational models that can quickly and efficiently navigate the design space and inform data collection efforts. In this work, we propose Graph Neural Network (GNN) models for predicting CO 2 solubility for DESs; the GNN leverages a mixture graph representation that captures the molecular structure of the DES components as well as their intermolecular interactions. Here, we compare the GNN framework against alternative architectures (neural networks, graph convolution networks, and random forests) and data representations (molecular fingerprints, sigma profiles, and graphs). We show that the proposed approach offers superior predictive performance; specifically, we show that solubility can be predicted reliably directly from molecular structure (without the need of using sigma profiles as proposed in previous studies). This result is important, as obtaining sigma profiles requires expensive density functional theory computations. We also explored the ability of GNNs to predict solubility for new DES mixtures and operating conditions. We found that the model extrapolates across temperature reliably. However, we also found deficiencies in the ability of the model to predict solubility for DES mixtures, pressures, and molar ratio not included in the training sets; we show that this is due to an inherent lack of chemical diversity in datasets available in the literature. The proposed computational capabilities can thus help navigate the design space of DES and inform data collection efforts. Our models, data, and benchmarks are shared as Python code implemented in Jupyter notebooks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A simulation-based model studying monoethanolamine and aprotic heterocyclic anion ionic liquid (AHA-IL) mixtures for carbon capture

Capturing of waste CO 2 , particularly from large industrial point sources, is necessary for either permanent storage or CO 2 utilization. Amine solvents used for carbon capture have certain drawbacks, particularly the high energy required by the reboiler for solvent regeneration and degradation of the amine with absorption/regeneration cycling. These challenges may be overcome using ionic liquids (ILs) that can be tailored for operating under typical flue gas stream conditions (8–10% CO 2 , 18–20% H 2 O, 2–3% O 2 , and 67–72% N 2 ). ILs, however, oftentimes suffer from high viscosities, leading to diffusion limitations when compared to amine solvents. This work models the combinations of aprotic heterocyclic anion ILs (AHA-IL) with monoethanolamine (MEA) and studies the carbon capture effectiveness of these mixtures. Furthermore, this study found that all AHA-IL/MEA mixtures exhibited 32 % lower enthalpy of reaction (compared to MEA) and 64 % lower viscosity (compared to AHA-IL), which are desirable characteristics for CO 2 capture solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving Mixtures of Soot Characterized by SP-AMS Spectra Using a Latent Dirichlet Allocation Model

Soot produced by detonation or combustion events exhibits different chemical properties depending on the fuel, device construction, and environmental conditions in which the event occurs. These properties can be useful for defining relevant signatures for probabilistically identifying the different types of events that occurred, based on the soot that is produced from these events. However, it is rare to observe samples of soot from a detonation or combustion that are not contaminated by outside particles. In this paper, we present a method for resolving mixtures of soot to determine the contributions of sources that may be present in samples of recovered soot. We use Latent Dirichlet Allocation to describe the generative process for a sample of recovered soot, and use Variational Bayesian Inference to learn about the parameters associated with the generative model. We demonstrate the utility of this method by considering real samples of mixtures of soot under various frameworks to show that the model is able to identify the different components present in a sample of soot as well as their mixing proportions.

54 ENVIRONMENTAL SCIENCES↗

Experimental and Kinetic Modeling Study of Laminar Burning Velocities of Cyclopentanone and Its Binary Mixtures with Ethanol and n-Propanol

Cyclopentanone is a promising biofuel that can enable more efficient engine operation and increase the fuel economy of the light duty fleet over current and planned technology developments. While the ignition of cyclopentanone has been investigated in detail, more studies on the laminar burning velocities of cyclopentanone are called for. In this work, the laminar burning velocities of cyclopentanone (C 5 H 8 O) have been measured using the heat flux and spherical flame methods at 1 atm, equivalence ratios from 0.7 to 1.6, and initial temperatures of 328, 353, and 428 K. To further investigate the relationship between the molecular structure and laminar burning velocity, identical experiments were also performed for binary mixtures of cyclopentanone with ethanol and n-propanol at 1:1 (mol). The consistency between the experimental data sets obtained in this work and literature data sets has been evaluated. Here, a recently published mechanism of cyclopentanone was used for simulation after adopting the submechanism of n-propanol. Good agreement has been seen between experimental and simulated results for all flames. To qualitatively explain the characteristics of the laminar burning velocity of cyclopentanone and the differences with those of ethanol and n-propanol, sensitivity analysis and reaction pathway analysis have been performed to compare the chemistry of the fuels under flame conditions, which revealed how the molecular structure of cyclopentanone could affect its laminar burning velocity. Compared to ethanol and n-propanol, cyclopentanone does not have primary carbon atoms in its molecule, leading to lower production of methyl radicals. Meanwhile, the carbonyl group in the cyclopentanone molecule is mostly released as CO in the decomposition of multiple intermediates accompanied by the production of unsaturated C 2 and C 4 species, especially C 2 H 4 and C 2 H 3 . Both features contribute to the high laminar burning velocity of cyclopentanone.

09 BIOMASS FUELS↗

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗

Estimation of mean electron energy in helium surface ionization waves on dielectric substrates

Abstract The determination of basic plasma parameters in atmospheric pressure discharges is critical to advancing their use in applications. Atmospheric pressure plasma jets have found use in the fields of medicine, agriculture, material modification and others. Atmospheric pressure plasma jets often generate plasma surface ionization waves (SIW) which interact with and propagate over surfaces. Electrical diagnostics are challenging in SIW due to high collision frequencies and small scale of the plasma discharge. This work employs a passive optical emission line ratio technique to estimate the mean electron energy in SIW over planar dielectric substrates. The method uses an intensity ratio of two helium triplet lines: He( 3 3 S ) at 706.5 nm and He( 3 3 D ) at 587.56 nm. A collisional-radiative model is used to correlate line ratio to mean electron energy and determine dependencies on electron density and He/air gas mixture. Mean electron energies ranging from 3–8 eV are determined in He/air mixtures and are found to remain constant as the surface wave propagates radially. This work provides a 2D, time-resolved, mean electron energy diagnostic for surface ionization wave propagation and validation of numerical modeling in atmospheric pressure systems with spatially varying He/air gas mixtures. The model in question is designed for use with any He line ratio in the n = 3 excitation level.

Morsell, Joshua (ORCID:0000000331073002)↗

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗

Competitive Metabolism of Polycyclic Aromatic Hydrocarbons (PAHs): An Assessment Using In Vitro Metabolism and Physiologically Based Pharmacokinetic (PBPK) Modeling

Humans are routinely exposed to complex mixtures such as polycyclic aromatic hydrocarbons (PAHs) rather than to single compounds, as are often assessed for hazards. Cytochrome P450 enzymes (CYPs) metabolize PAHs, and multiple PAHs found in mixtures can compete as substrates for individual CYPs (e.g., CYP1A1, CYP1B1, etc.). The objective of this study was to assess competitive inhibition of metabolism of PAH mixtures in humans and evaluate a key assumption of the Relative Potency Factor approach that common human exposures will not cause interactions among mixture components. To test this objective, we co-incubated binary mixtures of benzo[a]pyrene (BaP) and dibenzo[def,p]chrysene (DBC) in human hepatic microsomes and measured rates of enzymatic BaP and DBC disappearance. We observed competitive inhibition of BaP and DBC metabolism and measured inhibition coefficients (K i ), observing that BaP inhibited DBC metabolism more potently than DBC inhibited BaP metabolism (0.061 vs. 0.44 µM K i , respectively). We developed a physiologically based pharmacokinetic (PBPK) interaction model by integrating PBPK models of DBC and BaP and incorporating measured metabolism inhibition coefficients. The PBPK model predicts significant increases in BaP and DBC concentrations in blood AUCs following high oral doses of PAHs (≥100 mg), five orders of magnitude higher than typical human exposures. We also measured inhibition coefficients of Supermix-10, a mixture of the most abundant PAHs measured at the Portland Harbor Superfund Site, on BaP and DBC metabolism. We observed similar potencies of inhibition coefficients of Supermix-10 compared to BaP and DBC. Overall, results of this study demonstrate that these PAHs compete for the same enzymes and, at high doses, inhibit metabolism and alter internal dosimetry of exposed PAHs. This approach predicts that BaP and DBC exposures required to observe metabolic interaction are much higher than typical human exposures, consistent with assumptions used when applying the Relative Potency Factor approach for PAH mixture risk assessment.

54 ENVIRONMENTAL SCIENCES↗

TOUGH3/EOS7MG (EOS7MG) v1.1

EOS7MG is a new TOUGH3 module for modeling water, brine, and multiple gas mixtures in the subsurface multiphase flow system. The module allows mixture of as more as three gas components selecting from a list of predefined gases (CH4, C2H6, C3H8, H2S, CO2, N2, NO2, H2, and AIR). EOS7MG was developed based on EOS7C and inherited the general modeling capabilities of the TOUGH3 platform. It uses cubic equation of states for gas density, viscosity, and enthalpy calculation. The module provides the selection options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state. The subsurface flow and transport of gas and aqueous phase mixtures over a wide range of pressures and temperatures is modeled with a multiphase Darcy's Law and an accurate chemical equilibrium approach as originally implemented in EOS7C is used for calculation of the solubilities of different mass components. EOS7MG was designed to allow flexible number of components based on the number of gases included in the model. The potential application of this module may include subsurface gas storage, energy storage, CCS, enhanced gas recovery, and more.

Zhang, Keni↗

Neutron Absorber Plate Characterization Plan for Criticality Experiments Design

After being used in nuclear installations, depleted fuel can still be highly reactive and must be handled securely to prevent any radiological or criticality concerns. In particular, spent fuel from use in nuclear power reactors must be stored and transported in specifically designed containers using neutron absorber materials to prevent criticality. Various neutron absorber material types exist and are manufactured by various entities, as thoroughly described in the Handbook of Neutron Absorber Materials for Spent Nuclear Fuel Storage and Transportation Applications written by EPRI. Presently, one of the most modern and most widely used types of neutron absorber material contains particles of boron carbide, or B 4 C, embedded in aluminum matrix: Boralcan, manufactured by Rio Tinto. It is very important for the community to know as much as possible about such neutron absorber materials. Therefore, in the recent years, a US Department of Energy National Nuclear Security Administration–Nuclear Criticality Safety Program funded project initiated design of an experiment that places Boralcan neutron-absorbing plates in an established critical assembly using low-enriched uranium fuel at the Sandia Pulsed Reactor Facility/Critical Experiments (SPRF/CX) apparatus at Sandia National Laboratories. The goal of the experiment is to produce high-quality benchmark data to submit to the International Criticality Safety Benchmark Evaluation Project (ICSBEP), for use in validating calculational tools and nuclear data by criticality safety analysts. The project, named IER-554, is currently in its final design stage, following a successful preliminary design. In the work documented in the design study, ten critical configurations using Boralcan neutron absorber plates were designed, and the experiment was proven to be feasible, with a predicted low k eff uncertainty around 100 pcm. An overview of the modeled cutout of the critical assembly with a Boralcan plate is shown in Figure 1, representing one of the configurations planned for the critical experiments. Before the plates are inserted in the critical assembly, it is necessary to know more about their composition and uniformity. This summary focuses on the plate characterization plans. Each plate will undergo (1) neutron transmission measurements at different locations to determine the 10 B areal density and (2) an in-depth x-ray computed tomography (XCT) examination to obtain the exact Sizes and distribution of the B4C powder particles inside the plates. In parallel, plate modeling studies are performed with a goal to determine the validity of the currently used approximation of modeling the neutron absorber plates as a homogeneous mixture of Aluminum 1100 alloy and B4C— instead of explicitly modeling the B4C particles. By using the experimental 10 B areal density measurements, and the exact size and location of the B4C particles obtained by XCT, a plate model can theoretically be built that reproduces the plate with extremely high fidelity. The results of this modeling study could increase the confidence of the criticality safety community in its modeling methods when using this type of neutron absorber material, and the industry could use these validations to change the boron loading credit limits from the U.S. Nuclear Regulatory Commission standard review plan for dry cask storage of spent nuclear fuel. The modeling calculations are performed with SCALE 6.3.0 using the KENO V.a sequence for criticality calculations with the ENDF/B-VIII.0 continuous-energy cross section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Planetesimal Impact Vapor Plumes and Nebular Shocks Form Chondritic Mixtures

The origin of chondrules and the chondritic sedimentary rocks that dominate the meteoritic record is a long-standing problem in planetary science. Here, we develop a physical model for the formation of chondritic mixtures as an outcome of vaporizing collisions between planetesimals that were dynamically excited by the growth and migration of planets. We present calculations of nebular shock waves generated by impact vapor plumes and focus on aspects of the plume interaction with the nebular gas and dust that have been neglected in previous studies of impact ejecta. We find that when water dominates the vapor, the plumes are relatively cool. However, the plume expansion is supersonic and can drive strong shock waves in the dusty nebular gas. Portions of these nebular shock fronts initially melt nebular dust, forming chondrules that are coupled to the moving front. As the shock front expands and cools, the chondrules solidify, while the shock front entrains additional dust. Eventually, the plume expansion stalls and then hydrodynamically collapses, turbulently mixing variably processed dust and size-sorted chondrules. For probable impact parameters and nebular conditions during giant planet growth and migration, the impact-generated mixtures have characteristics that span the range observed in chondritic meteorites, providing an environment for rapid formation of chondritic assemblages after chondrule formation. Our impact vapor and nebular shocks model, named the IVANS model, links chondrule formation to the overall context of planet formation and provides a framework for interpreting the detailed chronological and geochemical record contained in chondritic meteorites.

79 ASTRONOMY AND ASTROPHYSICS↗

Data-driven predictions of complex organic mixture permeation in polymer membranes

Membrane-based organic solvent separations are rapidly emerging as a promising class of technologies for enhancing the energy efficiency of existing separation and purification systems. Polymeric membranes have shown promise in the fractionation or splitting of complex mixtures of organic molecules such as crude oil. Determining the separation performance of a polymer membrane when challenged with a complex mixture has thus far occurred in an ad hoc manner, and methods to predict the performance based on mixture composition and polymer chemistry are unavailable. Here, we combine physics-informed machine learning algorithms (ML) and mass transport simulations to create an integrated predictive model for the separation of complex mixtures containing up to 400 components via any arbitrary linear polymer membrane. We experimentally demonstrate the effectiveness of the model by predicting the separation of two crude oils within 6-7% of the measurements. Integration of ML predictors of diffusion and sorption properties of molecules with transport simulators enables for the rapid screening of polymer membranes prior to physical experimentation for the separation of complex liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shatter Thresholds and Fragment Size Distributions of Deuterium–Neon Mixture Cryogenic Pellets for Tokamak Thermal Mitigation

Reliable mitigation is necessary to eliminate the detrimental effects of a disruption event in large high-current tokamaks such as ITER. To avoid serious damage to plasma-facing components during the thermal quench phase of a disruption, material is injected to radiate the plasma energy over the inner surface of the machine. The most promising method of material injection is a process known as shattered pellet injection (SPI). SPI utilizes cryogenic cooling to desublimate gas into the barrel of a pipe gun to form a solid pellet. High-pressure gas or a mechanical punch is used to dislodge the pellet and accelerate it into a bent tube to intentionally fracture it. Pellets made of a mixture of deuterium and neon are likely candidates for thermal mitigation. The survivability of these pellets throughout their flight path, before striking the shatter tube, is essential for reliable SPI operation. Experiments were conducted to determine intact speed limits for various mixtures. In this paper, we outline the details of brittle fracture theory and compare a theory-based model to experimental results from various mixtures of deuterium and neon pellets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

1D modeling of plasma streamers at ammonia-air flame conditions

Abstract Self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures has been performed in this work. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of O ( 1 D ) in the production of key radicals, such as O, OH, and NH 2 . Furthermore, the contributions of the dissociative electronic excitation of NH 3 towards the production of H and NH 2 radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗