Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of Low-Carbon Energy Adoption on Airborne Particulate Matter Concentrations With Feedbacks to Future Climate Over California

California plans to reduce emissions of long-lived greenhouse gases (GHGs) through adoption of new energy systems that will also lower concentrations of short-lived absorbing soot contained in airborne particulate matter (PM). Here we examine the direct and indirect effects of reduced PM concentrations under a low-carbon energy (GHG-Step) scenario on radiative forcing in California. Simulations were carried out using the source-oriented WRF/Chem (SOWC) model with 12 km spatial resolution for the year 2054. Furthermore, the avoided aerosol emissions due to technology advances in the GHG-Step scenario reduce ground level PM concentrations by ~8.85% over land compared to the Business-as-Usual (BAU) scenario, but changes to meteorological parameters are more modest. Top-of-atmospheric forcing predicted by the SOWC model increased by 0.15 W m -2 , surface temperature warmed by 0.001 K, and planetary boundary layer height (PBLH) increased by 2.20 cm in the GHG-Step scenario compared to the BAU scenario. PM climate feedbacks are small because the significant changes in ground level PM concentrations associated with the GHG-Step scenario are limited to the first few hundred meters of the atmosphere, with little change for the majority of the vertical column above that level. As an order-of-magnitude comparison, the long-term effects of global reductions in GHG emissions (RCP8.5–RCP4.5) lowered average surface temperature over the California study domain by approximately 0.76 K. The effects of long-lived climate pollutants such as CO 2 are much stronger than the effects of short-lived climate pollutants such as PM soot over California in the year 2054.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Hybrid Pathway to Biojet Fuel via 2,3-Butanediol

Production of biomass-derived sustainable alternative jet fuels (SAJF) has been considered as an important approach to decarbonize the aviation industry but still possesses various challenges in technology advancement, particularly in achieving high carbon efficiency. Here we report a hybrid pathway to SAJF from 2,3-butanediol (2,3-BDO), integrating biologically converting biomass to 2,3-BDO with catalytically upgrading of 2,3-BDO to jet-range hydrocarbons. This pathway is demonstrated to have a high carbon recovery to liquid hydrocarbons from corn stover (25–28%) (74–82% of the theoretical maximum efficiency). The catalytic conversion steps involve 2,3-BDO to C3+ olefins, oligomerization, and hydrogenation where the first two steps are the focus of this study. Under optimum reaction conditions (523 K, 115 kPa, 1.0 h-1 weight hourly space velocity), 2,3-BDO conversion and C3+ olefin selectivity are >97% and 94–98% during 40 h time on stream, respectively. To demonstrate the adaptability of this technology with bio-derived 2,3-BDO, we also investigated the impact of water and other organic coproducts (acetoin and acetic acid) inherited from the fermentation broth on the catalyst performance and product selectivities. We have shown that the catalyst can handle a significant amount of water in the liquid feed (40 wt% water/60 wt% 2,3-BDO) and maintain catalyst stability for ~40 h. Acetoin can be converted to similar C3+ olefins as 2,3-BDO with complete conversion of acetoin. Co-feeding 10 wt% acetoin with 2,3-BDO is found to have no impact on 2,3-BDO conversion, C3+ olefin selectivity, and catalyst stability. The utilization of organic coproduct like acetoin can help to improve overall carbon conversion efficiency when using real biomass-derived 2,3-BDO. On the other hand, the presence of 10 wt% acetic acid is shown to drastically inhibit methyl ethyl ketone (MEK) hydrogenation and butene oligomerization as revealed by the increased MEK and butene selectivities, implying the importance of separating organic acids when feeding bio-derived 2,3-BDO. The formed C3–C6 olefins from 2,3-BDO are further oligomerized over Amberlyst-36 catalyst to longer-chain hydrocarbons with >70 wt% jet-range hydrocarbons including predominantly iso-olefins/iso-paraffins. The overall carbon efficiency for the jet-range hydrocarbons is 19–22%, exceeding most of the reported biojet pathways, which makes it a promising approach for SAJF production.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Controlled hydrogenation of a biomass-derived platform chemical formed by aldol-condensation of 5-hydroxymethyl furfural (HMF) and acetone over Ru, Pd, and Cu catalysts

Her we report that we studied the hydrogenation at temperatures from 313–393 K of a biomass-derived platform molecule, 5-hydroxymethyl furfural (HMF)-acetone-HMF (HAH) over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Lumped reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. The 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, the 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Batch conditions are determined for each hydrogenated product that give the highest yields in both batch and plug flow reactors.

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

Understanding the origin of structure sensitivity in hydrodechlorination of trichloroethylene on a palladium catalyst

Periodic density functional theory (DFT) calculations are employed to understand the origin of structure sensitivity in hydrodechlorination (HDC) of trichloroethylene (TCE) over different facets of a palladium catalyst. The HDC reaction is simulated on the terrace (Pd (111) and Pd (100)) and undercoordinated (Pd (211) and Pd (110)) sites of the Pd catalyst. The most stable binding configuration of TCE on the Pd surfaces is observed to be through the di-σ mode of binding, wherein each carbon atom of the TCE molecule is adsorbed atop of the Pd atom. On comparing TCE adsorption over different facets of Pd, a maximum binding energy of –178 kJ mol –1 is calculated over the Pd (110) surface. TCE, upon adsorption on Pd catalyst, undergoes dechlorination followed by hydrogenation of the hydrocarbon intermediates. The activation energies for C–Cl bond dissociation steps are significantly low when compared to the hydrogenation steps. The chlorine released from dechlorination tends to block the active sites, thereby poisoning the surface with high binding energies (B.E > –160 kJ mol –1 ) on all the surfaces. The trend in chlorine binding energies on Pd facets follows: Pd (110) > Pd (211) > Pd (100) > Pd (111). The removal of surface chlorine is facilitated by its reaction with surface hydrogen to form hydrogen chloride. The activation energy for hydrogen chloride formation is calculated to be 90 kJ mol –1 and 88 kJ mol –1 on Pd (111) and Pd (100) terrace sites, respectively and 109 kJ mol –1 and 118 kJ mol –1 on the step Pd (211) and corrugated Pd (110) facets, respectively. This suggests the ease of removal of Cl as HCl from the terrace sites as compared to the step and corrugated sites. The structure sensitivity in the TCE HDC reaction could possibly arise due to the differences in the energetics of Cl removal on different Pd facets. As a result, this mechanistic understanding could provide a rationale for designing suitable catalysts for the HDC of TCE.

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

Multivariate zeolitic imidazolate frameworks with an inverting trend in flexibility

Here, through systematic linker substitution in a flexible zeolitic imidazolate framework (ZIF) with step-shaped adsorption–desorption, structural intermediates between the known open and closed phases were isolated. Reflecting this, modulative sorption behaviour with an inverting adsorption pressure trend—in which the step pressure decreases and then increases again with increasing mixed linker concentration—is observed, highlighting how linker substitution modifies the energetic landscape of framework flexibility

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

Structural dynamics of a thermally silent triiron( II ) spin crossover defect grid complex

The structural evolution of spin crossover (SCO) complexes during their spin transition at equilibrium and out-of-equilibrium conditions needs to be understood to enable their successful utilisation in displays, actuators and memory components. In this study, diffraction techniques were employed to study the structural changes accompanying the temperature increase and the light irradiation of a defect [2 × 2] triiron(II) metallogrid of the form [Fe II 3 L H 2 (HL H ) 2 ](BF 4 ) 4 ·4MeCN (FE3), L H = 3,5-bis{6-(2,2'-bipyridyl)}pyrazole. Although a multi-temperature crystallographic investigation on single crystals evidenced that the compound does not exhibit a thermal spin transition, the structural analysis of the defect grid suggests that the flexibility of the grid, provided by a metal-devoid vertex, leads to interesting characteristics that can be used for intermolecular cooperativity in related thermally responsive systems. Time-resolved photocrystallography results reveal that upon excitation with a ps laser pulse, the defect grid shows the first two steps of the out-of-equilibrium process, namely the photoinduced and elastic steps, occurring at the ps and ns time scales, respectively. Similar to a previously reported [2 × 2] tetrairon(II) metallogrid, FE3 exhibits a local distortion of the entire grid during the photoinduced step and a long-range distortion of the lattice during the elastic step. Although the lifetime of the pure photoinduced high spin (HS) state is longer in the tetranuclear grid than in the defect grid, suggesting that the global nuclearity plays a crucial role for the lifetime of the photoinduced species, the influence of the co-crystalising solvent on the lifetime of the photoinduced HS state remains unknown. This study sheds light on the out-of-equilibrium dynamics of a thermally silent defect triiron SCO metallogrid.

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

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Emerging edge states on the surface of the epitaxial semimetal CuMnAs thin film

Epitaxial thin films of CuMnAs have recently attracted attention due to their potential to host relativistic antiferromagnetic spintronics and exotic topological physics. Here, we report on the structural and electronic properties of a tetragonal CuMnAs thin film studied using scanning tunneling microscopy (STM) and density functional theory (DFT). STM reveals a surface terminated by As atoms, with the expected semi-metallic behavior. An unexpected zigzag step edge surface reconstruction is observed with emerging electronic states below the Fermi energy. DFT calculations indicate that the step edge reconstruction can be attributed to an As deficiency that results in changes in the density of states of the remaining As atoms at the step edge. Furthermore, this understanding of the surface structure and step edges on the CuMnAs thin film will enable in-depth studies of its topological properties and magnetism.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

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

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

LaFeSi–LaFe 13-x Si x composites: Modulating magnetic and magnetocaloric properties through inherent stress manipulation

We examine structural and magnetic properties of a series of La–Fe–Si alloys in the region of concentrations where they naturally form two-phase LaFeSi–LaFe 13-x Si x composites with variable content and connectivity of LaFe 13-x Si x grains distributed within the LaFeSi matrix. Theoretical calculations confirm that the LaFeSi constituent is magnetically and structurally inert below room temperature and at pressures between -10 and 10 GPa. The LaFe 13-x Si x constituent, on the other hand, is magnetically and structurally active: it exhibits first-order magnetostructural transformations that, in addition to x Si , can be controlled with temperature, magnetic field, and pressure. In composites where the concentration of the inactive constituent is ~70 wt. % or greater, the standard, single-step, LaFe 13-x Si x first-order phase transformation proceeds in two steps separated by over 30 K in a zero magnetic field. Increasing the magnetic field recouples the two steps and restores the single-step phase transformation pathway. We analyze the roles of stresses caused by both thermal expansion mismatch and the first-order magnetic phase transition in LaFe 13-x Si x to rationalize the observed physical behaviors that emerge as the temperature or/and magnetic field vary.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗
DOE OSTI2022

The pathway for coenzyme M biosynthesis in bacteria

Mercaptoethane sulfonate or coenzyme M (CoM) is the smallest known organic cofactor and is most commonly associated with the methane-forming step in all methanogenic archaea but is also associated with the anaerobic oxidation of methane to CO 2 in anaerobic methanotrophic archaea and the oxidation of short-chain alkanes in Syntrophoarchaeum species. It has also been found in a small number of bacteria capable of the metabolism of small organics. Although many of the steps for CoM biosynthesis in methanogenic archaea have been elucidated, a complete pathway for the biosynthesis of CoM in archaea or bacteria has not been reported. Here, we present the complete CoM biosynthesis pathway in bacteria, revealing distinct chemical steps relative to CoM biosynthesis in methanogenic archaea. The existence of different pathways represents a profound instance of convergent evolution. The five-step pathway involves the addition of sulfite, the elimination of phosphate, decarboxylation, thiolation, and the reduction to affect the sequential conversion of phosphoenolpyruvate to CoM. The salient features of the pathway demonstrate reactivities for members of large aspartase/fumarase and pyridoxal 5'-phosphate–dependent enzyme families.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Cross-Code Verification of Neutronics Analysis Tools at INL Applied for 238 Pu Production in the Advanced Test Reactor

Here, analyses are completed for experiments prior to experiment irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Various codes are used to qualify all experiments planned for insertion in the reactor, thereby ensuring that all safety and programmatic requirements are satisfied preirradiation. Among the common experiment analysis tools at INL are MCNP5 coupled to ORIGEN2 (MOPY) and MC21. MOPY uses MCNP5 for transport calculations along with calculations for fluxes and select reaction rates, and then ORIGEN2 handles the step-by-step and postirradiation depletion. MC21 handles all in-reactor transport and step-by-step, during-irradiation, depletion calculations, and then ORIGEN (SCALE 6.2.3) is used for decay and dose calculations postirradiation. The MOPY results, along with those obtained via two variations of the MC21 model, were compared in terms of 238 Pu production in the ATR’s H10 position. For the MOPY model, the MC21 model utilizing the HELIOS-based fission product (FP) library, and the MC21 model utilizing the expanded 1300 FP library, the during-cycle irradiation in-core heating results were sufficiently equivalent; however, the MOPY model and the MC21 model with the HELIOS library showed some differences relating to the respective FP libraries. Ultimately, the MC21 model with a 1300 FP library produced the most consistent results throughout the cycle, whereas the MC21 model that utilized the (smaller) HELIOS library was able to handle during-irradiation analysis but lacked certain short-lived FPs that significantly contributed to the total decay heat at shutdown. MOPY, on the other hand, was found to overpredict fission gas production, as a result of limitations in the ORIGEN2 code.

ATR↗
DOE OSTI2021

Methanol carbonylation to acetaldehyde on Au particles supported by single-layer MoS 2 grown on silica

Homogenous single-layer MoS 2 films coated with sub-single layer amounts of gold are found to isolate the reaction of methanol with carbon monoxide, the fundamental step toward higher alcohols, from an array of possible surface reactions. Active surfaces were prepared from homogenous single-layer MoS 2 films coated with sub-single layer amounts of gold. These gold atoms formed clusters on the MoS2 surface. A gas mixture of carbon monoxide (CO) and methanol (CH 3 OH) was partially converted to acetaldehyde (CH 3 CHO) under mild process conditions (308 kPa and 393 K). This carbonylation of methanol to a C 2 species is a critical step toward the formation of higher alcohols. Density functional theory modeling of critical steps of the catalytic process identify a viable reaction pathway. Imaging and spectroscopic methods revealed that the single layer of MoS 2 facilitated formation of nanoscale gold islands, which appear to sinter through Ostwald ripening. Here, the formation of acetaldehyde by the catalytic carbonylation of methanol over supported gold clusters is an important step toward realizing controlled production of useful molecules from low carbon-count precursors.

2D catalyst↗
DOE OSTI2023

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗