Search NASA⌕ Search

SEARCH · Search NASA

Results for “RECOMBINATION COEFFICIENT”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electron ITB formation in EAST high poloidal beta plasmas under dominant electron heating

Abstract Plasma confinement and transport in tokamaks play a crucial role in the development of high poloidal beta steady-state operation scenarios. Therefore, it is very important to study the relevant mechanism of internal transport barriers (ITBs), which can help plasmas to obtain better confinement in order to achieve higher fusion gain. This paper mainly introduces the analysis of the characteristics of electron heat transport of discharges with ITB in high β P operation regime in Experimental Advanced Superconducting Tokamak (EAST). Based on the statistical analysis of stable discharges with β P > 1.5, it is found that there is an obvious bifurcation of the normalised electron temperature gradient (ETG) ( R / L T e ) in the range of β P = 2–2.2. Then the discharges of lower β P ( β P < 2, where the value of β P is below the bifurcation threshold) and of higher β P ( β P > 2.2, where the value of β P is above the threshold) were selected for analysis. The diagnostic data provided by Thomson scattering, x-ray crystal spectrometry and charge exchange recombination spectroscopy are used to provide reliable parameter profiles and then combined with the data of external magnetic probe measurements and the polarisation interferometer diagnosis system to fully reconstruct the balance. A relevant plasma current calculation model is used to calculate and analyse the current density profiles and power deposition, and then the transport analysis is carried out. Interestingly, in the higher β P discharges, it is found that the turbulence intensity provided by the CO 2 laser collective scattering system gradually decreases and the normalised ETG R / L T e gradually increases. It is also found that the electron heat transport coefficient decreases in the discharge with higher β P and the growth rate of the electron-scale turbulence calculated by transport gyro-Landau fluid (TGLF) is significantly reduced. Meanwhile, similar conclusions are also obtained in the discharges when the β P is further increased.

Wang, Z. H. (ORCID:0009000203431036)↗
NASA NTRS1995

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗
DOE OSTI2022

Modeling of deuterium and carbon radiation transport in MAST-U tokamak advanced divertors

Modest effects of deuterium and carbon radiation opacity in the Super-X and snowflake divertor plasmas are predicted for MAST Upgrade tokamak with core plasma input power 2.5–5 MW and plasma current 1 MA. The radiation transport modeling is based on the SOLPS-EIRENE and UEDGE code divertor plasma predictions. Two radiation transport models are used: one is based on a full radiation transport equation implemented in the radiation transport and collisional-radiative code CRETIN (without feedback on the background plasma), and another is an internal self-consistent UEDGE model with ionization, recombination, and heating rates corrected for Ly α line trapping based on the escape probability model implemented in CRETIN. In MAST-U, the Super-X and snowflake divertor plasmas are predicted to reach detached regimes at lower upstream densities than the standard divertor, and the conclusion still holds with radiation transport effects included. At neutral densities m -3 , modest Ly α deuterium line trapping with optical depths 10–15 is predicted in the Super-X divertor. Divertor plasmas are optically thin to other Lyman and Balmer lines, as well as to strong C III and C IV lines that are responsible for most of divertor radiated power. Insignificant changes (within a few percent) to divertor deuterium ionization and recombination rates are found. Radiation fluxes on outer divertor target are modified within a factor of 2–3 when the radiation transport is accounted for, and a similar variation is found due to the line shape models that define the absorption and emission line profiles in the radiation transport modeling. The predicted Lyman and Balmer spectral intensities are significantly modified due to radiation trapping. A measurement of divertor radiation transport effects is discussed using the Ly β /Ba α line ratio. In the snowflake divertor configuration, divertor plasmas are found to be optically thin to Lyman series lines within a large range of parameter variations that include upstream density, divertor transport coefficients, and magnetic configurations. Finally, modest radiation transport effects are only found in a few cases with strongest divertor transport and magnetic configurations closest to the ideal snowflake configuration, however, the plasma background models that were used are yet to be validated with an experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2021

Comprehensive Laboratory Measurements Resolving the LMM Dielectronic Recombination Satellite Lines in Ne-like Fe XVII Ions

We investigated experimentally and theoretically dielectronic recombination (DR) populating doubly excited configurations 3l3l¢ (LMM) in Fe XVII, the strongest channel for soft X-ray line formation in this ubiquitous species. We used two different electron beam ion traps and two complementary measurement schemes for preparing the Fe XVII samples and evaluating their purity, observing negligible contamination effects. This allowed us to diagnose the electron density in both EBITs. We compared our experimental resonant energies and strengths with those of previous independent work at a storage ring as well as those of configuration interaction, multiconfiguration Dirac–Fock calculations, and many-body perturbation theory. This last approach showed outstanding predictive power in the comparison with the combined independent experimental results. From these we also inferred DR rate coefficients, unveiling discrepancies from those compiled in the OPEN-ADAS and AtomDB databases.

Filipe Grilo↗
DOE OSTI2022

PolyDODT: a macrocyclic elastomer with unusual properties

The effect of reaction conditions on the structure of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by Reversible Radical Recombination Polymerization (R3P) using triethylamine (TEA), H 2 O 2 and air was investigated. 800 MHz (1 and 2D) NMR was used to investigate the polymer structures. Sensitivity analysis provided direct evidence for high purity cyclic polyDODT up to M n ~ 100 000 g mol –1 . Here, comparative analysis by High Resolution Multidetector Size Exclusion Chromatography (SEC) using integrated data showed that the cycles had lower viscosity and were more compact (both for R g and R h ) than linear samples of similar molecular weight. However, differential data revealed unusual behavior. While lower molecular weight cyclic polymers eluted later and had lower intrinsic viscosity than their linear counterparts at the same molecular weight, at higher molecular weights the polymers showed strange behavior: both the diffusion coefficient measured by Quasielastic Light Scattering (QELS) and Mark–Houwink–Sakurada plots of intrinsic viscosity for linear and cyclic polyDODT were found to converge. R3P, an aqueous based “green” method is capable of producing polymers at the 10–100 g scale in the lab, which will allow more detailed studies of this new class of biodegradable elastomers so further experimentation can be performed to elucidate the reasons for the unusual findings.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Quantum Mechanical Evaluation of CZTS Properties for Photovoltaic Applications

As a quaternary semiconductor, CZTS (Cu 2 ZnSnS 4 ) shows promise as a thin-film solar cell material. The constituent elements in CZTS are abundant and non-toxic, and CZTS has an ideal band gap (~1.5 eV) for solar cell applications, which produces a large absorption coefficient. Unfortunately, the low efficiencies and difficulties in synthesizing the desired crystallographic phase (kesterite) and the associated formation of antisite defects limit its deployment. Using theoretical and computational tools, we investigated these current limitations, explored strategies to resolve these issues, and developed design guidelines to improve the efficiency of CZTS-based sunlight absorbers. We first sought to understand the challenges associated with synthesizing the kesterite phase and proposed ways to favor kesterite growth, thus adding to the understanding of CZTS photovoltaics. We also analyzed the major limitations at the atomic scale lowering the efficiency of CZTS: the low open-circuit voltage (V oc ) and the short carrier lifetime, thus addressing the technical effectiveness of CZTS photovoltaics. We believe the latter is due to Shockley-Read-Hall (SRH) recombination centers associated with point defects. We also studied the electronic structure of interfaces and defects introducing such recombination centers, and explored doping strategies to reduce SRH recombination. Ultimately, these efforts will benefit the public by contributing to the development of cheap and efficient solar energy conversion technologies made from abundant, non-toxic materials.

14 SOLAR ENERGY↗
NASA NTRS1987

Calculations of the displacement damage and short-circuit current degradation in proton irradiated (AlGa)As-GaAs solar cells

A theoretical model for computing the displacement damage defect density and the short-circuit current (I sub sc) degradation in proton-irradiated (AlGa)As-GaAs p-n junction solar cells is presented. Assumptions were made with justification that the radiation induced displacement defects form an effective recombination center which controls the electron and hole lifetimes in the junction space charge region and in the n-GaAs active layer of the irradiated GaAs p-n junction cells. The degradation of I sub sc in the (AlGa)As layer was found to be negligible compared to the total degradation. In order to determine the I sub sc degradation, the displacement defect density, path length, range, reduced energy after penetrating a distance x, and the average number of displacements formed by one proton scattering event were first calculated. The I sub sc degradation was calculated by using the electron capture cross section in the p-diffused layer and the hole capture cross section in the n-base layer as well as the wavelength dependent absorption coefficients. Excellent agreement was found between the researchers calculated values and the measured I sub sc in the proton irradiated GaAs solar cells for proton energies of 100 KeV to 10 MeV and fluences from 10 to the 10th power p/square cm to 10 to the 12th power p/square cm.

Yeh, C. S.↗
NASA NTRS1985

Spectral identification/elimination of molecular species in spacecraft glow

Computer models of molecular electronic and vibrational emission intensities were developed. Known radiative emission rates (Einstein coefficients) permit the determination of relative excited state densities from spectral intensities. These codes were applied to the published spectra of glow above shuttle surface and to the Spacelab 1 results of Torr and Torr. The theoretical high-resolution spectra were convolved with the appropriate instrumental slit functions to allow accurate comparison with data. The published spacelab spectrum is complex but N2+ Meinel emission can be clearly identified in the ram spectrum. M2 First Positive emission does not correlate well with observed features, nor does the CN Red System. Spectral overlay comparisons are presented. The spectrum of glow above shuttle surfaces, in contrast to the ISO data, is not highly structured. Diatomic molecular emission was matched to the observed spectral shape. Source excitation mechanisms such as (oxygen atom)-(surface species) reaction product chemiluminescence, surface recombination, or resonance fluorescent re-emission will be discussed for each tentative assignment. These assignments are the necessary first analytical step toward mechanism identification. Different glow mechanisms will occur above surfaces under different orbital conditions.

Green, B. D.↗
NASA NTRS1998

Theoretical Studies of Important Processes in Planetary and Comet Atmospheres

Using theoretical quantum chemical calculations, I have successfully described the dissociative recombination (DR) of O2(+) leading to the excited S-1 state of atomic oxygen, the upper state of the well known green line emission. The process is described by O2(+) + e(-) yields O(S-1)) + O(D-1) (1) where e(-) is an electron and the product oxygen atoms are both excited. This process is important in the atmospheres of Venus, Mars and Earth. I have shown in prior calculations that only one repulsive potential curve of O2, f(sup l)Sigma(sub u, sup +) can generate O(S-1) from DR of the lowest vibrational levels of O2(+). However, in the prior results, the calculated quantum yield (i.e. the number of O(S-1) atoms produced for every two product atoms) from the O2(+) v = 0 level was smaller than the laboratory and atmospheric measurements by more than an order of magnitude. Including only direct recombination, the calculated quantum yield for O(S-1) is only 0.0016. In a calculation that accounts for both direct and indirect recombination, the quantum yield is 0.0012. The range of experimentally determined quantum yields is between 0.01 and 0.23. Because of this large difference between the theoretical and experimental quantum yields, it was thought for some time that the ionospheric and laboratory O2(+) must be vibrationally excited since for excited levels, theory gave quantum yields that are similar to experimental yields. It was also suggested that some other process was generating O(S-1) but none could be identified. Under current NASA support, I have found that reaction (1) proceeds via an unusual mechanism. The f(sup 1)Sigma(sub u, sup,+) state does not cross the ion between the turning points of the v = 0 level of the O2(+) ground state. The lack of a favorable crossing leads to a very small calculated DR rate coefficient. However, this mechanism assumes that initial electron capture must occur into the repulsive state that leads to )(S-1). Instead, I have found that initial electron capture occurs mostly into the B(sup 3)Sigma(sup u, sub -) state which crosses the ion between the turning points of the v = 0 ion level and has a large DR rate coefficient. The B state dissociates to O(D-1) and )(P-3). After capture, some of the flux is transferred to the f(sup 1)Sigma(sup u, sub +) state via symmetry mixed intermediate Rydberg states. The neutral Rydberg states are a mixture of Sigma(sub u, sup +)-1 and Sigma(sub u, sup -)-3 symmetry.

Guberman, Steven L.↗
DOE OSTI2024

Reaction Pathways and Energy Consumption in NH 3 Decomposition for H 2 Production by Low Temperature, Atmospheric Pressure Plasma

Pathways for NH 3 decomposition to N 2 and N 2 H 4 by atmospheric pressure nonthermal plasma are analyzed using a combination of molecular beam mass spectrometry measurements and zero-dimensional kinetic modeling. Experimental measurements show that NH 3 conversion and selectivity towards N 2 formation scale monotonically with the specific energy input into the plasma with ~ 100% selectivity to N 2 formation achieved at specific energy inputs above 0.12 J cm −3 (3.1 eV (molecule NH 3 ) −1 ). The kinetic model recovers these trends, although it underpredicts N 2 selectivity at low specific energy input. These discrepancies can be explained by the underestimation of reaction rate coefficients for reactions that consume N 2 H x species in collisions with H radicals and/or radial nonuniformities in power deposition, gas temperature, and species concentrations that are not represented by the plug flow approximation used in the model. The kinetic model shows that N 2 formation proceeds through N 2 H x decomposition pathways rather than NH x decomposition pathways in low temperature, atmospheric pressure plasma. Higher selectivity toward N 2 production can be achieved by operating at higher NH 3 conversion and with a higher gas temperature. Furthermore, the high energy cost of NH 3 decomposition by atmospheric pressure nonthermal plasma found in this work (25–50 eV (molecule NH 3 converted) −1 ; 17–33 eV (molecule H 2 formed) −1 ) is a result of the energy requirement for electron-impact dissociation of NH 3 and the significant re-formation of NH 3 by three-body recombination reactions between NH 2 and H.

Nonthermal plasma↗
DOE OSTI2025

Ultrafast and Highly Mobile Photocarriers in Monolayer WSe 2 Doped by α-RuCl 3

We report an experimental investigation of photocarrier dynamics in a heterostructure formed by a WSe 2 monolayer stacked on a multilayer α-RuCl 3 flake. The samples were fabricated using mechanical exfoliation and dry transfer techniques. Photoluminescence measurements showed that the excitonic photoluminescence of WSe 2 is quenched by more than 3 orders of magnitude upon contact with α-RuCl 3 . Transient absorption measurements revealed an ultrashort photocarrier lifetime of 0.5 ps in the heterostructure. Spatially resolved transient absorption microscopy measurements were employed to probe the transport properties of these photocarriers, resulting in a room-temperature diffusion coefficient of 370 cm 2 s –1 , which is higher than the exciton diffusion coefficients of most previously studied 2D semiconductors. Furthermore, these results suggest that α-RuCl 3 -doped WSe 2 could be utilized in high-speed optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ultra-wide bandgap β-Ga 2 O 3 films: Optical, phonon, and temperature response properties

Optical and phonon interactions of Ga 2 O 3 thin films with nanocrystalline morphology were studied at extreme temperatures. The films were grown using a sputtering technique and analyzed via temperature response transmission, Raman scattering, and high-resolution deep-UV photoluminescence (PL). Raman modes indicated that the structure corresponds to the β-phase. The optical-gap at the range of 77–620 K exhibited a redshift of ~200 meV, with a temperature coefficient of ~0.4 meV/K. The optical-gap at room-temperature is 4.85 eV. The electron–phonon interaction model at that temperature range pointed to a low energy phonon, ~31 meV, that is involved in the thermal properties of the optical-gap. Detailed Urbach energy analysis indicated that defects are the dominant mechanism controlling the band-edge characteristics even at an elevated temperature regime where phonon dominance is usually expected. Defects are attributed to the disordered forms of graphite that were detected via Raman scattering and to the granular morphology of the film. A deep-UV laser with an above bandgap exaction line of 5.1 eV was employed to map the PL of the films. The highly resolved spectra, even at room-temperature, show a strong emission of ~3.56 eV attributed to self-trapped holes (STHs). The STH is discussed and modeled in terms of the self-trapped exciton. Moreover, a very distinct but low-intensity emission was found at 4.85 eV that agrees with the value of the optical-gap and is attributed to bandgap recombination. The intensity ratio between the STH and that of the bandgap was found to be 6:1.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CdTe-based thin film photovoltaics: Recent advances, current challenges and future prospects

Cadmium telluride (CdTe)-based cells have emerged as the leading commercialized thin film photovoltaic technology and has intrinsically better temperature coefficients, energy yield, and degradation rates than Si technologies. More than 30 GW peak (GWp) of CdTe-based modules are installed worldwide, multiple companies are in production, modules are shipping at up to 18.6% efficiency, and lab cell efficiency is above 22%. We review developments in the science and technology that have occurred over approximately the past decade. These achievements were enabled by manufacturing innovations and scaling module production, as well as maximizing photocurrent through window layer optimization and alloyed CdSe x Te 1-x (CST) absorbers. Improved chlorine passivation processes, film microstructure, and serendipitous Se defect passivation significantly increased minority carrier lifetime. Efficiencies >22% have been realized for both Cu and As doped CST-based cells. The path to further efficiency gains hinges primarily on increasing open circuit voltage (V oc ) and fill factor (FF) through innovations in materials, fabrication methods, and device stacks. Replacing the longstanding Cu doping with As doping is resulting in better module stability and is being translated to large-scale production. To realize 25% efficiency and >1 V V oc , research and development is needed to increase the minority carrier lifetime beyond 100 ns, reduce grain boundary and interface recombination, and tailor band diagrams at the front and back interfaces. Many of these goals have been realized separately however combining them together using scalable manufacturing approaches has been elusive to date. We review these achievements and outstanding opportunities for this remarkable photovoltaic technology.

14 SOLAR ENERGY↗
NASA NTRS2002

Investigation of N2O Production from 266 and 532 nm Laser Flash Photolysis of O3/N2/O2 Mixtures

Tunable diode laser absorption spectroscopy has been employed to measure the amount of N2O produced from laser flash photolysis of O3/N2/O2 mixtures at 266 and 532 nm. In the 532 nm photolysis experiments very little N2O is observed, thus allowing an upper limit yield of 7 x 10(exp -8) to be established for the process O3 + N2 yield N2O + O2, where O3 is nascent O3 that is newly formed via O(3P(sub J)) + O2 recombination (with vibrational excitation near the dissociation energy of O3). The measured upper limit yield is a factor of approx. 600 smaller than a previous literature value and is approximately a factor of 10 below the threshold for atmospheric importance. In the 266 nm photolysis experiments, significant N2O production is observed and the N2O quantum yield is found to increase linearly with pressure over the range 100 - 900 Torr in air bath gas. The source of N2O in the 266 nm photolysis experiments is believed to be the addition reaction O(1D(sub 2)) + N2 + M yields (k(sub sigma)) N2O + M, although reaction of (very short-lived) electronically excited O3 with N2 cannot be ruled out by the available data. Assuming that all observed N2O comes from the O(1D(sub 2)) + N2 + M reaction, the following expression describes the temperature dependence of k(sub sigma) (in its third-order low-pressure limit) that is consistent with the N2O yield data: k(sub sigma) = (2.8 +/- 0.1) x 10(exp -36)(T/300)(sup -(0-88+0.36)) cm(sup 6) molecule(sup -2)/s, where the uncertainties are 2(sigma) and represent precision only. The accuracy of the reported rate coefficients at the 95% confidence level is estimated to be 30 - 40% depending on the temperature. Model calculations suggest that gas phase processes initiated by ozone absorption of a UV photon represent about 1.4% of the currently estimated global source strength of atmospheric N2O. However, these processes could account for a significant fraction of the oxygen mass-independent enrichment observed in atmospheric N2O, and they appear to be the first suggested photochemical mechanism that is capable of explaining the altitude dependence of the observed mass -independent isotopic signature.

Estupinan, E. G.↗
DOE OSTI2020

Transverse-momentum and event-shape dependence of D-meson flow harmonics in Pb–Pb collisions at s NN = 5.02 TeV

The elliptic and triangular flow coefficients v 2 and v 3 of prompt D 0 , D + , and D * + mesons were measured at midrapidity ( | y | < 0.8 ) in Pb–Pb collisions at the centre-of-mass energy per nucleon pair of s NN = 5.02 TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decays in the transverse momentum interval 1 < p T < 36 GeV / c in central (0–10%) and semi-central (30–50%) collisions. Compared to pions, protons, and J / ψ mesons, the average D-meson v n harmonics are compatible within uncertainties with a mass hierarchy for p T ≲ 3 GeV / c , and are similar to those of charged pions for higher p T . The coupling of the charm quark to the light quarks in the underlying medium is further investigated with the application of the event-shape engineering (ESE) technique to the D-meson v 2 and p T -differential yields. The D-meson v 2 is correlated with average bulk elliptic flow in both central and semi-central collisions. Within the current precision, the ratios of per-event D-meson yields in the ESE-selected and unbiased samples are found to be compatible with unity. All the measurements are found to be reasonably well described by theoretical calculations including the effects of charm-quark transport and the recombination of charm quarks with light quarks in a hydrodynamically expanding medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2017

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗