Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion processes”

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

Escape and finite-size scaling in diffusion-controlled annihilation

In this paper, we study diffusion-controlled single-species annihilation with a finite number of particles. In this reaction-diffusion process, each particle undergoes ordinary diffusion, and when two particles meet, they annihilate. We focus on spatial dimensions d>2 where a finite number of particles typically survive the annihilation process. Using scaling techniques we investigate the average number of surviving particles, M, as a function of the initial number of particles, N. In three dimensions, for instance, we find the scaling law M ~ N 1/3 in the asymptotic regime N»1. We show that two time scales govern the reaction kinetics: the diffusion time scale, T ~ N 2/3 , and the escape time scale, τ ~ N 4/3 . The vast majority of annihilation events occur on the diffusion time scale, while no annihilation events occur beyond the escape time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Graphite Lithiation under Fast Charging Conditions: Atomistic Modeling Insights

The charging of lithium ion batteries in a fast and safe manner is critical for promoting the mass adoption of electric vehicles. Li intercalation in graphite electrodes is known to be one of the bottlenecks during the fast charging process. The mechanism of Li diffusion in highly polarized graphite anode at high current rates remains, however, not well understood. In this work, Density Functional Theory (DFT) calculations are used to gain insights into the Li diffusion process in graphite when it is far from equilibrium under fast charging conditions. The effect of uncompensated charges on Li mobility is determined in the highly polarized regions of the anode close to the interfaces. The extra charge was found to increase the interlayer spacing in the diffusion layer and adjacent channels, increasing the diffusivity and promoting the formation of Li clusters. A concerted diffusion mechanism at the edge of high-concentration Li domains is proposed to enhance the diffusion of Li.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast Advective Water Flow through Nanochannels in Clay Interlayers: Implications for Moisture Transport in Soils and Unconventional Oil/Gas Production

Water flow in nanometer or sub-nanometer hydrophilic channels bears special importance in diverse fields of science and engineering. However, the nature of such water flow remains elusive. In this work, we report our molecular-modeling results on water flow in a sub-nanometer clay interlayer between two montmorillonite layers. We show that a fast advective flow can be induced by evaporation at one end of the interlayer channel, that is, a large suction pressure created by evaporation (~818 MPa) is able to drive the fast water flow through the channel (~0.88 m/s for a 46 Å-long channel). Scaled up for the pressure gradient to a 2 μm particle, the velocity of water is estimated to be about 95 μm/s, indicating that water can quickly flow through a μm-sized clay particle within seconds. The prediction seems to be confirmed by our thermogravimetric analysis of bentonite hydration and dehydration processes, which indicates that water transport at the early stage of the dehydration is a fast advective process, followed by a slow diffusion process. The possible occurrence of a fast advective water flow in clay interlayers prompts us to reassess water transport in a broad set of natural and engineered systems such as clay swelling/shrinking, moisture transport in soils, water uptake by plants, water imbibition/release in unconventional hydrocarbon reservoirs, and cap rock integrity of supercritical CO 2 storage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

(Photo)Electrocatalytic CO 2 Reduction at the Defective Anatase TiO 2 (101) Surface

Excessive carbon dioxide (CO 2 ) emissions by combustion of fossil fuels are linked to global warming and rapid climate change. One promising route to lowering the concentration of CO 2 in the atmosphere is to reduce it to useful small molecules via photoelectrocatalytic hydrogenation, which would enable solar energy storage with a zero-carbon emission cycle and perform a more efficient separation of the photogenerated electron and hole pair than pure photocatalysis. Indeed, photoelectrocatalytic CO 2 reduction has been an intense focus of research. Using the density functional theory (DFT), we studied the CO 2 reduction reaction on the defective anatase TiO 2 (101) surface, at both the solvent/catalyst and the electrolyte/catalyst interfaces. The analysis of the electronic structure of the surface shows a contrast between the solvent/catalyst and the electrolyte/catalyst interfaces, which results in the two corresponding catalytic cycles being distinct. Furthermore, our study explains at the electronic and mechanistic levels why methanol is the main product in the presence of the electrolyte and why the overpotential is not only controlled by the reaction process but also by the diffusion process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗

Mass dependence of overshooting beneath convective envelopes

State of the art stellar evolution codes [Paxton et al., 2010, Demarque et al., 2004, Weiss and Schlattl, 2008, Siess et al., 2013, Christensen-Dalsgaard, 2008] evoke a diffusive process to model convective overshooting. The diffusion coefficient [Freytag et al., 1996, Pratt et al., 2017] used in the model can be set to change with the classification of the convective zone as non-burning, H-burning, He-burning, or metal-burning. It can also be set to change at defined evolutionary points, such as the bottom of the asymptotic giant branch, or during the third dredge up [Herwig, 2000, Lugaro et al., 2003]. Aside from these abrupt changes, the diffusion coefficient is typically locked to a percentage of the pressure scale height measured at the convective boundary. However, there is no theoretical reason for convective overshooting, or indeed other convective properties, to change in the same way that the pressure scale height changes as a star evolves.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding the Mechanism of Electronic Defect Suppression Enabled by Nonidealities in Atomic Layer Deposition

Silicon germanium (SiGe) is a multifunctional material considered for quantum computing, neuromorphic devices, and CMOS transistors. However, implementation of SiGe in nanoscale electronic devices necessitates suppression of surface states dominating the electronic properties. The absence of a stable and passive surface oxide for SiGe results in the formation of charge traps at the SiGe-oxide interface induced by GeOx. In an ideal ALD process in which oxide is grown layer by layer, the GeO x formation should be prevented with selective surface oxidation (i.e., formation of an SiO x interface) by controlling the oxidant dose in the first few ALD cycles of the oxide deposition on SiGe. However, in a real ALD process, the interface evolves during the entire ALD oxide deposition due to diffusion of reactant species through the gate oxide. In this work, this diffusion process in nonideal ALD is investigated and exploited: the diffusion through the oxide during ALD is utilized to passivate the interfacial defects by employing ozone as a secondary oxidant. Periodic ozone exposure during gate oxide ALD on SiGe is shown to reduce the integrated trap density (D it ) across the band gap by nearly 1 order of magnitude in Al 2 O 3 (<6 × 10 10 cm -2 ) and in HfO 2 (<3.9 × 10 11 cm -2 ) by forming a SiO x -rich interface on SiGe. Depletion of Ge from the interfacial layer (IL) by enhancement of volatile GeO x formation and consequent desorption from the SiGe with ozone insertion during the ALD growth process is confirmed by electron energy loss spectroscopy (STEM-EELS) and hypothesized to be the mechanism for reduction of the interfacial defects. In this work, the nanoscale mechanism for defect suppression at the SiGe-oxide interface is demonstrated, which is engineering of diffusion species in the ALD process due to facile diffusion of reactant species in nonideal ALD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian Static Parameter Estimation for Partially Observed Diffusions via Multilevel Monte Carlo

In this article we consider static Bayesian parameter estimation for partially observed diffusions that are discretely observed. We work under the assumption that one must resort to discretizing the underlying diffusion process, for instance, using the Euler--Maruyama method. Given this assumption, we show how one can use Markov chain Monte Carlo (MCMC) and particularly particle MCMC [C. Andrieu, A. Doucet, and R. Holenstein, J. R. Stat. Soc. Ser. B Stat. Methodol., 72 (2010), 269--342] to implement a new approximation of the multilevel (ML) Monte Carlo (MC) collapsing sum identity. Our approach comprises constructing an approximate coupling of the posterior density of the joint distribution over parameter and hidden variables at two different discretization levels and then correcting by an importance sampling method. The variance of the weights are independent of the length of the observed data set. The utility of such a method is that, for a prescribed level of mean square error, the cost of this MLMC method is provably less than i.i.d. sampling from the posterior associated to the most precise discretization. However the method here comprises using only known and efficient simulation methodologies. Finally, the theoretical results are illustrated by inference of the parameters of two prototypical processes given noisy partial observations of the process: the first is an Ornstein--Uhlenbeck process and the second is a more general Langevin equation.

97 MATHEMATICS AND COMPUTING↗

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The leaching kinetics of cadmium from hazardous Cu-Cd zinc plant residues

Highlights: • A cleaner cadmium, zinc and copper recovery method from zinc plant waste residues was investigated. • The mechanisms and kinetics of cadmium recovery were discussed. • The effects of leaching variables were discussed. •The XRD and morphology analysis of raw residue and leaching residue were conducted. - Abstract: A large amount of Cu-Cd zinc plant residues (CZPR) are produced from the hydrometallurgical zinc plant operations. Since these residues contain substantial amount of heavy metals including Cd, Zn and Cu, therefore, they are considered as hazardous wastes. In order to realize decontamination treatment and efficient extraction of the valuable metals from the CZPR, a comprehensive recovery process using sulfuric acid as the leaching reagent and air as the oxidizing reagent has been proposed. The effect of temperature, sulfuric acid concentration, particle size, solid/liquid ratio and stirring speed on the cadmium extraction efficiency was investigated. The leaching kinetics of cadmium was also studied. It was concluded that the cadmium leaching process was controlled by the solid film diffusion process. Moreover, the order of the reaction rate constant versus H{sub 2}SO{sub 4} concentration, particle size, solid/liquid ratio and stirring speed was calculated. The XRD and SEM-EDS analysis results showed that the main phases of the secondary sulfuric acid leaching residues were lead sulfate and calcium sulfate.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Direct Visualization of Defect‐Controlled Diffusion in van der Waals Gaps

Abstract Diffusion processes govern fundamental phenomena such as phase transformations, doping, and intercalation in van der Waals (vdW) bonded materials. Here, the diffusion dynamics of W atoms by visualizing the motion of individual atoms at three different vdW interfaces: hexagonal boron nitride (BN)/vacuum, BN/BN, and BN/WSe 2 , by recording scanning transmission electron microscopy movies is quantified. Supported by density functional theory (DFT) calculations, it is inferred that in all cases diffusion is governed by intermittent trapping at electron beam‐generated defect sites. This leads to diffusion properties that depend strongly on the number of defects. These results suggest that diffusion and intercalation processes in vdW materials are highly tunable and sensitive to crystal quality. The demonstration of imaging, with high spatial and temporal resolution, of layers and individual atoms inside vdW heterostructures offers possibilities for direct visualization of diffusion and atomic interactions, as well as for experiments exploring atomic structures, their in situ modification, and electrical property measurements of active devices combined with atomic resolution imaging.

Chemistry↗

Tritium species diffusion on and desorption from γ-LiAlO 2 (100) surface: A first-principles investigation

γ-LiAlO 2 enriched in the 6 Li isotope is a good candidate for tritium (T, 3 H) production in nuclear reactor. In this work, to better understand its 3 H production performance, first-principles calculations are performed to study the 3 H-related species diffusion on and desorption from the γ-LiAlO 2 (100) surface. For the diffusion process, we investigate the pathways and energy barriers for the T and OT diffusion on the pristine and defective γ-LiAlO 2 (100) surfaces. Our results reveal that the O and Li/Al atoms of the surface layer play important roles for T and OT diffusion, respectively. By comparing with its diffusion barrier in bulk γ-LiAlO 2 , the order of energy barrier for T diffusion is: bulk < surface < defective surface, which is reversed for OT diffusion. Space is the prime factor for OT diffusion, while diffusion-mediating atom instead of space becomes the key factor for T diffusion. For the desorption process, we investigate the 3 H-related species (T, OT, T 2 , and T 2 O) desorption from the γ-LiAlO 2 (100) surface. The desorption energies we obtained are so high that the desorption behavior of all the 3 H-related species is nearly prohibited. We expect the possible T 2 or T 2 O release from the large accumulation of 3 H atoms on the surface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sub-cloud moist entropy curvature as a predictor for changes in the seasonal cycle of tropical precipitation

Convective Quasi-Equilibrium (CQE) may be a useful framework for understanding the precipitation minus evaporation (P-E) response to CO2-induced warming. To explore this proposition, a suite of aquaplanet simulations with a slab ocean from TRACMIP (the Tropical Rain belts with an Annual cycle and a Continent Model Intercomparison Project) is analyzed. A linear relationship between P-E and the curvature of sub-cloud moist entropy, a marker for the spatial distribution of sub-cloud moist energy and onset of a tropical direct overturning circulation under CQE conditions, is shown to exist across many of the TRACMIP simulations. Furthermore, this linear relationship is a skillful predictor of changes in P-E in response to CO2-induced warming. The curvature metric also shows improvement in predicting P-E changes compared to the simpler method of relating P-E directly to the sub-cloud moist entropy field or a simple ‘wet-get-wetter’ type null hypothesis, especially on seasonal and shorter timescales. Using fixed relative humidity in the curvature metric and sub-cloud moist entropy degrades their ability to predict P-E changes, implying that both temperature and relative humidity changes in the boundary layer are important for characterizing future precipitation changes. To understand why the curvature metric is a skillful predictor of hydrological changes, a moist static energy (MSE) budget analysis is performed, which shows that MSE divergence on sub-daily timescales is well parameterized as a downgradient diffusive process. Additionally, this transient MSE divergence has a similar structure to the curvature term. Assuming that the transient MSE divergence is related to convective activity, these findings suggest that the theoretical underpinning of the linear relationship between the curvature term and P-E resides in cloud processes that both remove spatial gradients of sub-cloud moist entropy and generate precipitation at sub-seasonal and shorter timescales.

aquaplanet simulations↗

Sorption Modeling and Verification for Off-Gas Treatment

This project was successfully executed to provide valuable adsorption data and improve a comprehensive model developed in previous work by the authors. Data obtained were used in an integrated computer program to predict the behavior of adsorption columns. The model is supported by experimental data and has been shown to predict capture of off gas similar to that evolving during the reprocessing of nuclear waste. The computer program structure contains (a) equilibrium models of off-gases with the adsorbate; (b) mass-transfer models to describe off-gas mass transfer to a particle, diffusion through the pores of the particle, and adsorption on the active sites of the particle; and (c) incorporation of these models into fixed bed adsorption modeling, which includes advection through the bed. These models are being connected with the MOOSE (Multiphysics Object-Oriented Simulation Environment) software developed at the Idaho National Laboratory through DGOSPREY (Discontinuous Galerkin Off-gas SeParation and REcoverY) computer codes developed in this project. Experiments for iodine and water adsorption have been conducted on reduced silver mordenite (Ag0Z) for single layered particles. Adsorption apparatuses have been constructed to execute these experiments over a useful range of conditions for temperatures ranging from ambient to 250°C and water dew points ranging from -69 to 19°C. Experimental results were analyzed to determine mass transfer and diffusion of these gases into the particles and to determine which models best describe the single and binary component mass transfer and diffusion processes. The experimental results were also used to demonstrate the capabilities of the comprehensive models developed to predict single-particle adsorption and transients of the adsorption-desorption processes in fixed beds. Models for adsorption and mass transfer have been developed to mathematically describe adsorption kinetics and transport via diffusion and advection processes. These models were built on a numerical framework for solving conservation law problems in one-dimensional geometries such as spheres, cylinders, and lines. Coupled with the framework are specific models for adsorption in commercial adsorbents, such as zeolites and mordenites. Utilizing this modeling approach, the authors were able to accurately describe and predict adsorption kinetic data obtained from experiments at a variety of different temperatures and gas phase concentrations. A demonstration of how these models, and framework, can be used to simulate adsorption in fixed- bed columns is provided. The CO 2 absorption work involved modeling with supportive experimental information. A dynamic model was developed to simulate CO 2 absorption using high alkaline content water solutions. The model is based upon transient mass and energy balances for chemical species commonly present in CO 2 absorption. A computer code was developed to implement CO 2 absorption with a chemical reaction model. Experiments were conducted in a laboratory scale column to determine the model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. Continuing work could employ the model to control column operation and predict the absorption behavior under various input conditions and other prescribed experimental perturbations. The value of the validated models and numerical frameworks developed in this project is that they can be used to predict the sorption behavior of off-gas evolved during the reprocessing of nuclear waste and thus reduce the cost of the experiments. They can also be used to design sorption processes based on concentration limits and flow-rates determined at the plant level.

42 ENGINEERING↗

Measurement of moisture-dependent ion diffusion constants in wood cell wall layers using time-lapse micro X-ray fluorescence microscopy

Our future bioeconomy depends on increased utilization of renewable lignocellulosic biomass. Controlling the diffusion of chemicals, such as inorganic ions, within secondary plant cell walls is central to many biomass applications. However, insufficient understanding of intra-cell-wall diffusion within secondary plant cell walls is hindering the advancement of many lignocellulosic biomass applications. In this work, X-ray fluorescence microscopy was used to measure diffusion constants of K + , Cu 2+ , and Cl - diffusing through loblolly pine (Pinus taeda) cell wall layers under 70%, 75%, or 80% relative humidity (RH). Results revealed that diffusion constants increased with RH, the larger Cu 2+ diffused more slowly than the K + , and the Cl - diffusion constant was the same as that for the counter cation, indicating cations and anions diffused together to maintain charge neutrality. Comparison with electrical conductivity measurements showed that conductivity is being controlled by ion mobility over these RH. The results further support that intra-cell-wall diffusion of inorganic ions is a Fickian diffusion process occurring through rubbery amorphous polysaccharides, which contradicts previous assertions that intra-cell-wall diffusion is an aqueous process occurring through water pathways. Researchers can now utilize polymer science approaches to engineer the molecular architecture of lignocellulosic biomass to optimize properties for specific end uses.

59 BASIC BIOLOGICAL SCIENCES↗

Growth and Properties of Vapor Diffused Nb3Sn Coating for Superconducting Radiofrequency Accelerator Cavity Applications

The desire to reduce construction and operating costs of future SRF accelerators motivates the search for higher-performing alternative materials. Nb3Sn (Tc ~ 18.3 K and Hsh ~ 425 mT) is the front runner. The tin vapor diffusion process is currently the technique of choice to produce promising Nb3Sn-coated cavities. Understanding Nb3Sn nucleation and growth in this process is essential to progress. Samples representing different stages of Nb3Sn formation have been produced and studied to elucidate the effects of nucleation, growth, process conditions, and impurities. Nb3Sn films with thickness from a few hundred nm up to ~15 µm were grown and characterized. The microscopic examination of samples suggests the mechanisms of thin film nucleation and growth. Broadly, nucleation deposits tin as a thin surface phase and, under some conditions, a few hundred nanometer sized particles as well. Conditions that impair nucleation promote the formation of defects, such as ?patches?, in subsequent coating growth. Analysis of coated samples is consistent with the model of Nb3Sn grown in which tin diffuses via grain boundaries to Nb3Sn-Nb interface, where the growth Nb3Sn into the niobium bulk takes place. Similar scaling laws are found for grain growth and layer thickness. Non-parabolic layer growth is consistent with significant grain growth, which reduces the number of Sn-transport channels. Examination of patchy region in Nb3Sn coating revealed it to be large single crystalline grains, pointing to impeded Nb3Sn layer growth due to low grain boundary density, resulting in a significantly thin coating in those areas. Examination of RF loss regions from a coated cavity, identified with a thermometry mapping system showed patchy regions and carbonus defects were associated with strong local field-dependent surface resistance. RF measurements of coated cavities were combined with material characterization of witness samples and coated-cavity cutouts to improve the coating process. Understanding obtained and applied to cavity coatings, resulted in single-cell Nb3Sn cavities with a quality factor of ~2 ×1010 up to 15 MV/m accelerating gradient at 4 K, without "Wuppertal" Q-slope. We have also produced Nb3Sn-coated CEBAF 5-cell cavities with accelerating gradients useful for accelerator cryomodules. This dissertation will discuss the genesis of the Nb3Sn coating in a typical tin vapor diffusion process, effects of different process parameters, and its consequences to the coating of single-cell and multi-cell SRF cavities.

Pudasaini, Uttar↗

First-passage time statistics on surfaces of general shape: Surface PDE solvers using Generalized Moving Least Squares (GMLS)

Here, we develop numerical methods for computing statistics of stochastic processes on surfaces of general shape with drift-diffusion dynamics d X t = a (X t ) dt + b(X t ) d W t . We formulate descriptions of Brownian motion and general drift-diffusion processes on surfaces. We consider statistics of the form u (x) = E x [$∫^{τ}_{0}$ g (X t ) dt ] + E x [ f (X τ )] for a domain Ω and the exit stopping time τ = inf t { t >0 | X i Ω}, where f , g are general smooth functions. For computing these statistics, we develop high-order Generalized Moving Least Squares (GMLS) solvers for associated surface PDE boundary-value problems based on Backward- Kolmogorov equations. We focus particularly on the mean First Passage Times (FPTs) given by the case f = 0, g = 1 where u (x) = E x [τ]. We perform studies for a variety of shapes showing our methods converge with high-order accuracy both in capturing the geometry and the surface PDE solutions. We then perform studies showing how statistics are influenced by the surface geometry, drift dynamics, and spatially dependent diffusivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗