Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric model dependence”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Fantômas unconfined: global QCD fits with Bézier parameterizations

Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.

Bézier curves↗
DOE OSTI2025

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
DOE OSTI2022

Impact of Irradiation-induced Lattice Defects and Dopants on Thermal Transport in Oxide Nuclear Fuels

Advanced nuclear reactor concepts aim to use fuels that must withstand unprecedented temperature and radiation extremes. In these fuels, thermal energy transport under irradiation is directly related to fuel longevity, reactor safety, and is arguably one of the most important performance metrics. Here we provide a comprehensive, first-principles-informed treatment of phonon mediated thermal transport in a defect-bearing actinide oxide with direct comparison to experimental measurements. Pristine and proton irradiated thorium dioxide was chosen as a model system to treat the complexity of thermal transport in the presence of lattice defects. A thermal transport model is implemented using the linearized Boltzmann transport equation (LBTE) with input from first principles calculations and defect evolution models. The output of the LBTE is compared directly to mesoscopic measurements of thermal conductivity on length scales commensurate with defect accumulation. Parametric measurements of conductivity with irradiation dose and temperature reveal a monotonic decrease in conductivity with irradiation dose, in the range of 0.001 dpa to 0.1 dpa. Besides radiation-induced defects, doping the thorium dioxide crystal with small quantities of uranium atoms also results in a drastic reduction in thermal conductivity. Temperature-dependent measurements of thermal conductivity in uranium-doped thorium dioxide single crystals suggests a resonant scattering mechanism is responsible for the observed reduction. This comprehensive, atomistic- to meso-scale treatment provides the necessary basis to investigate thermal transport under irradiation in more complex systems that exhibit strong electron correlation. This is a poster presentation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Investigation of methane and ethane pyrolysis with highly active and durable iron-alumina catalyst to produce hydrogen and valuable nano carbons: Continuous fluidized bed tests and reaction rate analysis

The goal of this research work was to evaluate a novel iron-alumina catalyst for methane and ethane pyrolysis with high conversion rates to H 2 with a long life time to address two major issues reported with previous catalysts. Fluidized bed tests demonstrated continuous hydrogen production while collecting carbon containing fines. During the 230 h fluid bed test the methane conversion rate to hydrogen at 700 °C remained above 80% for 160 h while creating valuable nano carbons. Similarly high catalyst lifetimes with high methane conversions have not been reported before. Ethane was more active than methane for pyrolysis with 100% conversion to H 2 . Transmission electron microscopy and Raman analysis of the carbon showed single and multiwall nanotubes and conglomerate nano-fibers. Ethane did not impact the allotrope of carbon formed. TGA was used to create a data array that was used for kinetic analysis and reaction model parametric regressions. Ethane had a higher rate of pyrolysis that was 2.3–2.8 times that of methane. Ethane pyrolysis was found to be more sensitive to temperature with a higher activation energy of 62 kJ/mol as compared to methane at 43 kJ/mol. The frequency factor, A i , was significantly higher for ethane at 0.73 [[kmol/m 3 -sec]/[N/m2] n ] compared to 0.03 for methane on the same basis. Both gases had a functional partial pressure order dependence of ~0.6 based on the regression analysis.

08 HYDROGEN↗
DOE OSTI2021

Influence of gas expansion on the propagation of a premixed flame in a spatially periodic shear flow

It has been previously demonstrated that thermal gas expansion might have a role in boundary layer flashback of premixed turbulent flames, inducing local flow-reversal in the boundary layer's low-velocity streaks on the reactants’ side of the flame and facilitating its upstream propagation. Here, we perform a two-dimensional numerical investigation of the interaction between a periodic shear flow and a laminar premixed flame. The periodic shear is a simplified model for the oncoming prolonged streamwise velocity streaks with alternating regions of high and low velocities found in turbulent boundary layers in the vicinity of the walls. The parametric study focuses on the amplitude and wavelength of the periodic shear flow and on the gas expansion ratio (unburnt-to-burnt density ratio). With the increase of the amplitudes of the periodic shear flow and of the gas expansion, the curved flame velocity increases monotonically. The flame velocity dependence on the periodic shear wavelength is non-monotonic, which is consistent with previous theoretical studies of curved premixed flame velocity. The flame shape that is initially formed by the oncoming periodic shear appears to be metastable. At a later stage of the flame propagation, the flame shape transforms into the stationary one dominated by the Darrieus-Landau instability.

42 ENGINEERING↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Evidence of Quadrupole and Octupole Deformations in Zr 96 + Zr 96 and Ru 96 + Ru 96 Collisions at Ultrarelativistic Energies

In the hydrodynamic model description of heavy-ion collisions, the elliptic flow v 2 and triangular flow v 3 are sensitive to the quadrupole deformation β 2 and octupole deformation β 3 of the colliding nuclei. The relations between v n and β n have recently been clarified and were found to follow a simple parametric form. The STAR Collaboration has just published precision v n data from isobaric Ru 96 + Ru 96 and Zr 96 + Zr 96 collisions, where they observe large differences in central collisions v 2 , Ru > v 2 , Zr and v 3 , Ru < v 3 , Zr . Using a transport model simulation, we show that these orderings are a natural consequence of $\beta_{2,Ru} \gg \beta_{2,Zr}$ and $\beta_{3,Ru} \ll \beta_{3,Zr}$ We reproduce the centrality dependence of the v 2 ratio qualitatively and v 3 ratio quantitatively and extract values of β 2 and β 3 that are consistent with those measured at low-energy nuclear structure experiments. STAR data provide the first direct evidence of strong octupole correlations in the ground state of Zr 96 in heavy-ion collisions. Our analysis demonstrates that flow measurements in high-energy, heavy-ion collisions, especially using isobaric systems, are a new precision tool to study nuclear structure physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

CMB-S4: Iterative Internal Delensing and r Constraints

Abstract The tightest constraints on the tensor-to-scalar ratio r can only be obtained after removing a substantial fraction of the lensing B -mode sample variance. The planned Cosmic Microwave Background (CMB)-S4 experiment ( cmb-s4.org ) will remove the lensing B -mode signal internally by reconstructing the gravitational lenses from high-resolution observations. We document here a first lensing reconstruction pipeline able to achieve this optimally for arbitrary sky coverage. We make it part of a map-based framework to test CMB-S4 delensing performance and its constraining power on r , including inhomogeneous noise and two non-Gaussian Galactic polarized foreground models. The framework performs component separation of the high-resolution maps, followed by the construction of lensing B -mode templates, which are then included in a parametric small-aperture map cross-spectra-based likelihood for r . We find that the lensing reconstruction and framework achieve the expected performance, compatible with the target σ ( r ) ≃ 5 · 10 −4 in the absence of a tensor signal, after an effective removal of 92%–93% of the lensing B -mode variance, depending on the simulation set. The code for the lensing reconstruction can also be used for cross-correlation studies with large-scale structures, lensing spectrum reconstruction, cluster lensing, or other CMB lensing-related purposes. As part of our tests, we also demonstrate the joint optimal reconstruction of the lensing potential with the lensing curl potential mode at second order in the density fluctuations.

Astronomy & Astrophysics↗
DOE OSTI2020

Planck 2018 results

We report on the implications for cosmic inflation of the 2018 release of the Planck cosmic microwave background (CMB) anisotropy measurements. The results are fully consistent with those reported using the data from the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be n s = 0.9649 ± 0.0042 at 68% CL. We find no evidence for a scale dependence of n s , either as a running or as a running of the running. The Universe is found to be consistent with spatial flatness with a precision of 0.4% at 95% CL by combining Planck with a compilation of baryon acoustic oscillation data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, r0.002 < 0.10, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain r 0.002 < 0.056. In the framework of standard single-field inflationary models with Einstein gravity, these results imply that: (a) the predictions of slow-roll models with a concave potential, V"(Φ) < 0, are increasingly favoured by the data; and (b) based on two different methods for reconstructing the inflaton potential, we find no evidence for dynamics beyond slow roll. Three different methods for the non-parametric reconstruction of the primordial power spectrum consistently confirm a pure power law in the range of comoving scales 0.005 Mpc -1 ≲ k ≲ 0.2 Mpc -1 . A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectra. For the case of oscillatory features that are logarithmic or linear in k, this result is further strengthened by a new combined analysis including the Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions for the cosmological fluctuations. In correlated, mixed adiabatic and isocurvature models, the non-adiabatic contribution to the observed CMB temperature variance is constrained to 1.3%, 1.7%, and 1.7% at 95% CL for cold dark matter, neutrino density, and neutrino velocity, respectively. Planck power spectra plus lensing set constraints on the amplitude of compensated cold dark matter-baryon isocurvature perturbations that are consistent with current complementary measurements. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadupolar modulation of the primordial fluctuations. However, the polarization data do not support physical models for a scale-dependent dipolar modulation. All these findings support the key predictions of the standard single-field inflationary models, which will be further tested by future cosmological observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Linear two-pool models are insufficient to infer soil organic matter decomposition temperature sensitivity from incubations

Abstract Terrestrial carbon (C)-climate feedbacks depend strongly on how soil organic matter (SOM) decomposition responds to temperature. This dependency is often represented in land models by the parameter Q 10 , which quantifies the relative increase of microbial soil respiration per 10 °C temperature increase. Many studies have conducted paired laboratory soil incubations and inferred “active” and “slow” pool Q 10 values by fitting linear two-pool models to measured respiration time series. Using a recently published incubation study (Qin et al. in Sci Adv 5(7):eaau1218, 2019) as an example, here we first show that the very high parametric equifinality of the linear two-pool models may render such incubation-based Q 10 estimates unreliable. In particular, we show that, accompanied by the uncertain initial active pool size, the slow pool Q 10 can span a very wide range, including values as high as 100, although all parameter combinations are producing almost equally good model fit with respect to the observations. This result is robust whether or not interactions between the active and slow pools are considered (typically these interactions are not considered when interpreting incubation data, but are part of the predictive soil carbon models). This very large parametric equifinality in the context of interpreting incubation data is consistent with the poor temporal extrapolation capability of linear multi-pool models identified in recent studies. Next, using a microbe-explicit SOM model (RESOM), we show that the inferred two pools and their associated parameters (e.g., Q 10 ) could be artificial constructs and are therefore unreliable concepts for integration into predictive models. We finally discuss uncertainties in applying linear two-pool (or more generally multiple-pool) models to estimate SOM decomposition parameters such as temperature sensitivities from laboratory incubations. We also propose new observations and model structures that could enable better process understanding and more robust predictive capabilities of soil carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Higher-order correlations between different moments of two flow amplitudes in Pb-Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

The correlations between different moments of two flow amplitudes, extracted with the recently developed asymmetric cumulants, are measured in Pb-Pb collisions at $\sqrt {s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the CERN Large Hadron Collider. The magnitudes of the measured observables show a dependence on the different moments as well as on the collision centrality, indicating the presence of nonlinear response in all even moments up to the eighth. Furthermore, the higher-order asymmetric cumulants show different signatures than the symmetric and lower-order asymmetric cumulants. Comparisons with state-of-the-art event generators using two different parametrizations obtained from Bayesian optimization show differences between data and simulations in many of the studied observables, indicating a need for further tuning of the models behind those event generators. These results provide new and independent constraints on the initial conditions and transport properties of the system created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗
DOE OSTI2025

An analysis of parameter compression and Full-Modeling techniques with Velocileptors for DESI 2024 and beyond

In anticipation of forthcoming data releases of current and future spectroscopic surveys, we present the validation tests and analysis of systematic effects within velocileptors modeling pipeline when fitting mock data from the AbacusSummit N-body simulations. We compare the constraints obtained from parameter compression methods to the direct fitting (Full-Modeling) approaches of modeling the galaxy power spectra, and show that the ShapeFit extension to the traditional template method is consistent with the Full-Modeling method within the standard ΛCDM parameter space. We show the dependence on scale cuts when fitting the different redshift bins using the ShapeFit and Full-Modeling methods. We test the ability to jointly fit data from multiple redshift bins as well as joint analysis of the pre-reconstruction power spectrum with the post-reconstruction BAO correlation function signal. We further demonstrate the behavior of the model when opening up the parameter space beyond ΛCDM and also when combining likelihoods with external datasets, namely the Planck CMB priors. Finally, we describe different parametrization options for the galaxy bias, counterterm, and stochastic parameters, and employ the halo model in order to physically motivate suitable priors that are necessary to ensure the stability of the perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

36 MATERIALS SCIENCE↗