Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface interaction”

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 73 records · Page 4

Probing surface-adsorbate interactions through active particle dynamics

Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanoscale cooperative adsorption for materials control

Abstract Adsorption plays vital roles in many processes including catalysis, sensing, and nanomaterials design. However, quantifying molecular adsorption, especially at the nanoscale, is challenging, hindering the exploration of its utilization on nanomaterials that possess heterogeneity across different length scales. Here we map the adsorption of nonfluorescent small molecule/ion and polymer ligands on gold nanoparticles of various morphologies in situ under ambient solution conditions, in which these ligands are critical for the particles’ physiochemical properties. We differentiate at nanometer resolution their adsorption affinities among different sites on the same nanoparticle and uncover positive/negative adsorption cooperativity, both essential for understanding adsorbate-surface interactions. Considering the surface density of adsorbed ligands, we further discover crossover behaviors of ligand adsorption between different particle facets, leading to a strategy and its implementation in facet-controlled synthesis of colloidal metal nanoparticles by merely tuning the concentration of a single ligand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗

Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

alkaline solutions, IDREAM, gibbsite, molecular dy↗

Protein Adhesion on Semi-Fluorinated Polystyrene Surfaces in Static and Dynamic Measurements

Reducing protein adhesion is a critical strategy in fouling-resistant material innovation, with broad applications spanning biomedical and healthcare devices, biosensors, industrial and environmental systems, and other important technological domains. Here, in this study, we elucidated protein adhesion behavior on polystyrene-based thin films by neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), using both lysozyme and bovine serum albumin (BSA) as model proteins. To this end, semifluorinated polystyrene thin films with gradient wettability and surface energy were fabricated through dry processing using plasma oxidation and gas-phase deposition. Although it is believed that a fully fluorinated alkyl chain offers extremely low surface energy, thus rejecting foulants, and has been used in many fouling-resistant surface designs, enhanced protein–surface interactions were observed consistently in NR and QCM-D results, due to the combined effects of surface morphology and chemistry. On the contrary, depositing shorter fluorinated silane onto a hydrophilic PS surface contributed to a more homogeneous nanoscale fluorine coating, resulting in less initial protein adsorption and improved surface recovery. Comparative analysis of proteins with different sizes on the nanopatterned semifluorinated surface revealed the influence of molecular characteristics on surface interactions. Lysozyme, being smaller and more compact, showed faster adsorption kinetics and higher surface coverage but largely reversible binding, whereas BSA, with its larger and more flexible structure, formed broader and more stable interfacial layers. This study fills the gap in understanding protein adhesion within the range of hydrophobicity (water contact angle ∼90°), as current strategies often associate with extreme hydrophilic and superhydrophobic surfaces due to hydration or low-surface-energy rejection mechanisms, respectively. It also provides in-depth insights into current combinatorial fouling-resistant surface design.

Yuan, Yue [Oak Ridge National Laboratory (ORNL), O↗

Many‐body response of benzene at monolayer MoS 2 : Van der Waals interactions and spectral broadening

Abstract Models of surface enhancement of molecular electronic response properties are challenging for two reasons: (a) molecule‐surface interactions require a simultaneous solution of the molecular and the surface dynamic response (a daunting task), and (b) when solving for the electronic structure of the combined molecule + surface system, it is not trivial to single out the particular physical effects responsible for enhancement. To tackle this problem, in this work, we apply a formally exact decomposition of the system's response function into subsystem contributions by using subsystem density functional theory (DFT), which grants access to dynamic polarizabilities and optical spectra. In order to access information about the interactions between the subsystems, we extend a previously developed subsystem‐based adiabatic connection fluctuation‐dissipation theorem of DFT to separate the additive from the nonadditive correlation energy and identify the nonadditive correlation as the van der Waals interactions. As an example, we choose benzene adsorbed on monolayer MoS 2 . We isolate the contributions to benzene's dynamic response arising from the interaction with the surface, and for the first time, we evaluate the enhancements to the effectiveness of C 6 coefficients as a function of benzene‐MoS 2 distance and adsorption site. We also quantify the spectral broadening of the benzene's electronic excited states due to their interaction with the surface. We find that the broadening has a similar decay law with the molecule‐surface distance as the leading van der Waals interactions (ie, R −6 ) and that the surface enhancement of dispersion interactions between benzene molecules is less than 5% but is still large enough (0.5 kcal/mol) to likely play a role in the prediction of interface morphologies.

Umerbekova, Alina↗

Effect of Solvent on the Interaction of Lignin with a Zeolite Nanosheet in the Condensed Phase

Lignin is the essential building block of lignocellulosic biomass, an excellent renewable source of different aromatic monomers for the polymer and biofuel industry. The depolymerization of lignin into value-added chemicals and fuels through the catalytic process poses a significant challenge due to the complex structure of lignin. Understanding lignin’s conformational diversity and dynamics in the liquid phase is crucial to designing an effective depolymerization process. Here, we conducted all-atom molecular dynamics simulations to understand the conformation and dynamics of softwood lignin on the all-silica zeolite nanosheet based on the MFI topology in a binary mixture of water–methanol at three different molar compositions (0%, 50%, and 100% methanol). We observed that the methanol–surface interaction is stronger than the water–surface interaction, and methanol readily diffused into the MFI core. Lignin surface contacts decrease with increasing methanol composition due to higher solubility and dynamics. Lignin dynamics on the surface in neat water is an order of magnitude smaller than methanol. Here, we also found that lignin adopts a slightly extended conformation when it stays on the surface than in the bulk solution phase for the pure water case, whereas for pure methanol and the binary solution structures are statistically similar.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Affinity of small-molecule solutes to hydrophobic, hydrophilic, and chemically patterned interfaces in aqueous solution

Significance Increasing global population, urbanization, and climate change pose tremendous constraints on water availability. While membrane processes offer attractive, energy-efficient options, transformative advances in membrane technology are needed to enable reuse of poor-quality waters from diverse sources. This work brings fundamental understanding of solute–surface interactions to membrane phenomena, including selectivity and fouling. Molecular simulations reveal limitations of traditional predictive measures of solute–surface partitioning, suggesting alternative thermodynamic signatures for solute binding that highlight the contextual nature of hydrophobicity. An inverse design algorithm provides a tool to pattern membrane nanoscale chemistry to optimize solute affinity and selectivity. This work demonstrates how simulations can inform practical materials design for sustainable water treatment and, more broadly, many technologies involving water-mediated solute–surface interactions.

36 MATERIALS SCIENCE↗

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Interactions between atmospheric pressure plasmas and metallic catalyst particles in packed bed reactors

Atmospheric-pressure plasmas sustained in packed bed reactors (PBRs) are being investigated for chemical conversion of gases and pollution control. Metallic catalysts added to the surfaces of the dielectric beads of PBRs can increase the energy efficiency and selectivity of chemical processes by reducing operating temperature and providing additional reaction pathways. In this paper, results from a computational investigation of plasma surface interactions between micron-scale metallic catalysts and humid-air plasmas in PBRs are discussed. Here we found that high plasma density regions form in the proximity of the metallic catalysts. These higher-density plasma regions were confirmed experimentally using ICCD imaging. The intense plasmas result from geometrical electric field enhancement and redistribution of charges within the conductive particles, leading to further enhancement. The high electric field at the triple points of the catalysts can produce electric field emission of electrons, which provides a pre-ionization source or additional source of electrons. These regions of high electric field and sources of electrons guide discharges towards the catalysts and increases fluxes of excited species, ions, electrons and photons to their surfaces. These fluxes are focused primarily at the triple points between the metal, dielectric and gas. As a result, the catalyst is locally heated, which could lead to further increased rates of thermocatalytic reactions on the surface. Surface roughness of the metal inclusions can lead to additional electric field enhancement, which changes the character of the discharges in the vicinity of the catalysts while reducing breakdown voltage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The influence of surface heat fluxes on the growth of idealized monsoon depressions

This study explores the effect of surface sensible and latent heat fluxes on monsoon depressions using a series of idealized convection-permitting simulations. Each experiment is initialized with a small-amplitude wave that is allowed to grow within an environment representative of the South Asian monsoon. Comparing experiments with and without interactive surface heat fluxes, it is found that these fluxes enhance the growth of the simulated vortices. Without interactive surface fluxes, the strengthening period is short and the vortices fail to reach intensities characteristic of stronger monsoon depressions. Using a large set of experiments in which the vertical and meridional shear are systematically varied, it is found that surface heat fluxes enhance intensity the most when the upper-level shear is weak, the lower-level shear and associated moist static energy (MSE) gradient are sufficiently steep, and the lower-level meridional shear is strong. Here, these experiments reveal two different regimes of convection-coupled monsoon depression growth: one in which convection is driven by MSE advection and one in which it is driven by surface heat fluxes and quasigeostrophic forcing for ascent. Both regimes require sufficiently strong meridional shear to achieve initial growth by barotropic instability.

54 ENVIRONMENTAL SCIENCES↗

Establishing structure-sensitivity of ceria reducibility: real-time observations of surface-hydrogen interactions

The first layer of atoms on an oxide catalyst provides the first sites for adsorption of reactants and the last sites before products or oxygen are desorbed. We employ a unique combination of morphological, structural, and chemical analyses of a model ceria catalyst with different surface terminations under an H 2 environment to unequivocally establish the effect of the last layer of atoms on surface reduction. (111) and (100) terminated epitaxial islands of ceria are simultaneously studied in situ allowing for a direct investigation of the structure–reducibility relationship under identical conditions. Kinetic rate constants of Ce 4+ to Ce 3+ transformation and equilibrium concentrations are extracted for both surface terminations. Unlike the kinetic rate constants, which are practically the same for both types of islands, more pronounced oxygen release, and overall higher reducibility were observed for (100) islands compared to (111) ones. The findings are in agreement with coordination-limited oxygen vacancy formation energies calculated by density functional theory. The results point out the important aspect of surface terminations in redox processes, with particular impact on the catalytic reactions of a variety of catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Joint Bayesian Inference for Near-Surface Explosion Yield and Height-of-Burst

Forensic capabilities to understand chemical and nuclear explosions are greatly aided by an accurate estimate of explosive yield with uncertainty. The relationship between explosive size and geophysical observations of seismic, acoustic, and optical waves can be exploited to provide an estimate of yield. Any near-surface yield estimate is complicated by the surface interaction, so an estimate for the explosion height-of-burst is necessarily included in the relationship. Additionally, the relationship dictates a trade-off between estimates of yield and height-of-burst. Fortunately, the surface interaction for each type of observation is different, which breaks the trade-off, and the inclusion of height-of-burst with multiple data types improves yield estimation. We define simple parametric forward models to relate seismoäcoustoöptic observations from a data set of known explosive yields and height-of-bursts. The parameters of the models and a prediction for the yield and height-of-burst of a new event can then be estimated given new observations via Bayesian inference. We report posterior distribution estimates of the parametric models using a Markov chain Monte Carlo sampling technique. These models are then used to predict the yield and height-of-burst of SUGAR, a historical near-surface nuclear explosion, using its reported historical observations. The reported yield of 1.2 ktonne Trinitrotoluene (TNT)-equivalent (Department of Energy, 2015) is within the estimated posterior. Yield uncertainty can be estimated from the spread of the posterior, which is between 0.9 and 2.1 ktonne TNT-equivalent. The posterior for height-of-burst has a wider range between 10 m below and 8 m above ground that includes the true height-of-burst of 1 m.

58 GEOSCIENCES↗

Perspectives and challenges of ultra-high temperature ceramics for fusion plasma-facing applications

Ultra-high temperature ceramics (UHTCs) offer several potential advantages as plasma-facing components (PFCs) in fusion reactors due to their extreme melting points, tailorable thermal conductivity, and attractive unirradiated mechanical properties including fracture toughness comparable or superior to tungsten. Here, recent developments and material properties of UHTCs are briefly reviewed, along with an overview of limited studies on their responses to neutron irradiation and an evaluation of plasma-surface interactions. Five key research pathways, primarily focused on irradiation effects, for advancing UHTCs in PFC applications are discussed: (1) assessing irradiation effects on the coupled thermal–mechanical performance (2) addressing the lack of studies on irradiation, plasma-surface interactions, and their synergistic effects; (3) investigating high-temperature (>1000 °C) neutron irradiation effects critical for PFC performance; (4) optimizing multi-component UHTC compositions or composites to improve thermal or mechanical properties; (5) enhancing radiation resistance to mitigate microcracking and void swelling through strategies such as increasing sink strength by reducing grain size, introducing fine particles, and leveraging complex concentrated alloy concepts.

36 MATERIALS SCIENCE↗

Etching of Si 3 N 4 induced by electron beam plasma from hollow cathode plasma in a downstream reactive environment

An etching system based on the interaction of electrons extracted from a direct current (DC) hollow-cathode (HC) Ar plasma and injected towards a Si3N4 covered silicon substrate located in the downstream reactive environment created by an additional remote CF 4 /O 2 plasma source was developed and evaluated. By controlling the properties of the injected beam electrons, this approach allows to deliver energy to a surface functionalized by exposure to reactive species, and initiate surface etching. The energy of the primary beam electrons is controlled by the acceleration voltage relative to the HC discharge. Ar atoms flow from the high-pressure HC discharge into the low pressure downstream reactive environment in the process chamber. For an acceleration voltage greater than the ionization potential of Ar and/or process gas species, the energetic primary beam electrons produce a secondary plasma in the process chamber, and can also cause additional dissociation. We have characterized the properties of the secondary plasma and also surface etching of Si 3 N 4 as a function of process parameters, including acceleration voltage (0-80 V), discharge current of the HC discharge (1-2 A), pressure (3.5-20 mTorr), source to substrate distance (1.5 to 5 cm), and feed gas composition (20% and 80% O2 in CF 4 /O 2 ). The electron energy probability function (EEPF) measured with a Langmuir probe about 2.5 cm below the extraction ring suggests several major groups of electrons for this situation, including high energy primary beam electrons with an energy that varies as the acceleration voltage is changed and low energy electrons produced by beam electron-induced ionization of the Ar gas in the process chamber. When a remote CF 4 /O 2 plasma is additionally coupled to the process chamber, Si 3 N 4 surfaces can be functionalized, and by varying the energy of the beam electrons, Si 3 N 4 etching can be induced by electron-neutral synergy effect with plasma-surface interaction. For conditions without beam electron injection, the remote plasma etching rate (ER) of Si 3 N 4 depends strongly on the O 2 concentration in the CF 4 /O 2 processing gas mixture and can be suppressed for O 2 -rich process conditions by the formation of a SiONF passivation layer on the Si 3 N 4 surface. The combination of the HC electron beam source with the remote plasma source (HCEB-RP) makes it possible to induce Si 3 N 4 etching for O 2 -rich remote plasma conditions where remote plasma by itself produces negligible Si 3 N 4 etching. The electron enhanced etching of Si 3 N 4 depends strongly on O 2 /CF 4 mixing ratio reflecting changing arrival rates of O and F species at the surface. Optical emission spectroscopy was used to estimate the ratio of gas phase F and O densities and found to be controlled by the gas mixing ratio and independent of HC EB operating conditions. At this time the detailed sequence of events operative in the etching mechanism is unclear. While the increase of the electron energy is ultimately responsible for initiating surface etching, at present we cannot rule out a role of ions from the simultaneously produced secondary plasma in plasma surface interaction mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Role of Cooperative Interactions Among Surfaces, Solvents, and Reactive Intermediates on Catalysis at Liquid–Solid Interfaces (Final Report DE-SC0020224)

This project established quantitative links between inner‑sphere chemistry (active metal identity, coordination, and zeolite topology) and outer‑sphere organization (solvent identity, hydrogen‑bond networks, and pore condensation) that govern rates, activation barriers, and selectivities for alkene epoxidation and epoxide ring‑opening at solid–liquid and quasi‑liquid–solid interfaces. We deconvoluted contributions from covalent interactions at active sites and noncovalent, solvent‑mediated interactions within pores by pairing well‑defined metal substituted zeolites with controlled solvent environments. We then mapped those contributions onto measurable kinetics (ΔH‡, ΔS‡), adsorption thermodynamics (ITC), and in situ spectroscopy. The transferrable outcomes include a set of design rules that include the following understandings. First, tune silanol ((SiOH)x) density and pore topology to organize solvent networks that selectively stabilize transition states. Second, exploit activity‑coefficient‑normalized rates and adsorption– barrier correlations to diagnose when solvent reorganization rather than surface chemistry limits performance. Third, use partial pore condensation (e.g., acetonitrile, water but also generalizable to other solvents) to elicit liquid‑like stabilization effects even in nominally vapor‑phase reactors. Collectively, these results provide strategies to increase epoxidation rates, improve oxidant utilization (i.e., selectivities), and steer regioselectivity in zeolite‑based catalytic processes relevant to sustainable oxidation chemistry. These outcomes should be transferable to other classes of reactions that proceed in microporous materials and under confinement provided by organized solvents (e.g., electrochemical double layers).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

Measurements of He-Ar, He-N 2 , and He-air thermal creep slip and accommodation coefficients at high temperatures

Understanding gas flows in capillaries has many applications in modeling the transport of gases in nano-structured, porous, or fractured media. A network of capillaries can often approximate these media, and also information on gas-surface interactions obtained from capillary experiments can be used for modeling flows in these media. Experimental data on flows of different mixtures of He-Ar, He-N 2 , and He-air in the slip regime and in stainless steel capillaries at high temperatures were obtained by using a two-bulb apparatus. An accurate expression for the thermal creep slip coefficients with the Lennard-Jones potential parameters and diffuse-specular reflection gas-surface interaction conditions were then used to obtain the accommodation coefficients for the different gases. The experimental data are best described with values of accommodation coefficients in the range of 0.1–0.3 for He and 0.5–1.0 for Ar, N 2 , and air. Further, the use of values in this range is suggested for modeling gaseous flows in capillaries and nano-structured, porous, or fractured media if other direct measured values for a particular medium are unavailable.

42 ENGINEERING↗