Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROUGHNESS”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Potentials of mean force fail to describe chemical bond-breaking in solution

Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.

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

Saturation and decay of an induced scattering instability driven by static field asymmetries in a non-neutral plasma

Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A point-particle-based hydride shell-shedding model for ejecta particle transport in reactive environments

A shock wave passing over a rough or perturbed metal surface will induce a limiting case of Richtmyer–Meshkov instability and will cause small particles to eject from the surface and transport into the surrounding medium. These particles are known as ejecta and can be either solid or liquid in nature. Recent experiments have shown that liquid cerium ejecta clouds exhibit unexpected non-monotonic acceleration behaviors as well as temperature plateaus after a brief temperature rise if they are transporting in a chemically reactive, hydrogen-based medium while they act as expected in an inert medium. This work details a point-particle model developed for reactive cerium ejecta transport, which attempts to account for these new physics through the behavior of a developing solid hydride shell, which is believed to form as a product of the reaction. The overall model incorporates the effects of the reaction on the particle properties as well as the effects of potential shedding of the shell into sub-micrometer scale flakes and potential phase change of the hydride if the ejecta particles reach the melt point of the hydride layer. The model is tested by performing simulations of the original motivating experiments and comparing quantities, such as ejected mass, velocimetry, and temperature profiles, against the experimental data. While the model is able to capture many general features of the observed anomalies, some inaccuracies still exist. These point to both missing physics in the model (such as a deuterium adsorption mechanism on the hydride layer) as well as a lack of knowledge of certain material properties (such as the strength of cerium hydride to determine dynamic fracture thicknesses) needed to fully reduce the uncertainties in the model by up to an order of magnitude and perform a true attempt at model validation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Bremsstrahlung constraints on proton-boron 11 inertial fusion

Proton-boron 11 (pB11) fusion is relatively safe and clean, but difficult to use for net power production, since bremsstrahlung radiation tends to radiate away power more quickly than it can be generated by fusion power, particularly once poisoning by alpha particles is taken into account. While in magnetic confinement fusion, this problem can be addressed by deconfining the alphas, in inertial confinement fusion (ICF), the alphas that heat the plasma linger for the duration of the reaction. Thus, it becomes essential to trap the bremsstrahlung radiation in the hotspot. Through burn simulations incorporating bremsstrahlung emission and reabsorption, we infer the necessary conditions to capture enough radiation to produce scientific breakeven in a pB11 ICF plasma. We find that breakeven requires a stagnation areal density roughly two orders of magnitude higher than the current state-of-the-art, at pressures three orders of magnitude higher.

Alpha particles↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2025

Surface and sub-surface porosity effects on the initial free surface expansion in shocked single crystal aluminum

Defects and roughness at a material’s surface can allow for the formation of a material jet when a shockwave reaches a surface. The formation of a jet can lead to significant material damage of the shocked material, as well as the downstream material that the jet impacts. Here, we utilize planar shocks in single crystal aluminum to assess initial jetting formation for surface notches and sub-surface porosity. In the case of surface notches, we find that changing the width of the notch can cause the shock breakout response to transition from a standard jetting mechanism to an atomic atomization at high velocities. This transition occurs due to the mechanisms associated with lateral relaxations within the notch and how this interacts with the shock focusing event that induces jetting. For sub-surface pores, the pore collapse induces a re-shock, localized near the pore, that is stronger than the initial shock, but unsupported. Changing the pore dimensions and location relative to the surface influences the strength of the shock that emanates from its collapse when it reaches the surface, leading to a variety of different breakout conditions including jetting and localized material failure.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Low-energy ion beam sputtering of polycrystalline diamond

Ion beam milling is a highly versatile and precise technique for the microfabrication of diamond films, enabling controlled material removal at both micrometer and sub-micrometer scales while accommodating complex geometries and surface features. Comprehensive sputtering yield data are essential for process optimization. In this work, we systematically measure the sputter yield of polycrystalline diamond bombarded with $\textrm{N}_2^+$, Ne + , Ar + , Kr + , and Xe + ions over an energy range of 300–900 eV and beam incidence angles of 0°−80°, providing data critical for precision machining applications. Results show that the sputter yield increases with both ion energy and the incidence angle and decreases with increasing ion mass. For $\textrm{N}_2^+$ ions, the sputter yield is largely insensitive to ion energy and exceeds values expected for pure physical sputtering. The concurrent formation of volatile C 2 N 2 molecules during $\textrm{N}_2^+$ ion bombardment suggests that the process is dominated by chemical sputtering. The absence of a peak in the dependence of the sputter yield on the beam incidence angle expected at oblique incidence angles is attributed to the effects of surface roughness.

Materials science↗
DOE OSTI2026

Epitaxial registry and transport properties in V4O7 thin films

We report the growth of V4O7 thin films deposited simultaneously on amorphous SiO2 and crystalline c-cut Al2O3 substrates. x-ray diffraction shows that films grown on SiO2 are polycrystalline with no preferred orientation, while films grown on Al2O3 exhibit a well-defined out-of-plane orientation and an in-plane registry consistent with epitaxial growth. Transmission electron microscopy confirms the crystallographic relationship between the V4O7 film and the Al2O3 substrate. Atomic force microscopy indicates substantially lower surface roughness for films grown on Al2O3 (~6 nm) compared to those grown on SiO2 (~22 nm), and Raman spectroscopy confirms stabilization of the same V4O7 phase on both substrates. Electrical transport measurements reveal a metal–insulator transition near 244 K for both substrates, with thermal hysteresis not exceeding ~1 K. Although the transition temperature remains essentially unchanged, films grown on Al2O3 exhibit higher conductivity over the entire temperature range, exceeding that of films grown on SiO2 by approximately two orders of magnitude at 100 K and by a factor of five at 300 K. These results indicate that the conductivity differences are consistent with variations in microstructural connectivity associated with crystallographic order.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase-shift speckle-shearing interferometry

We have numerically simulated the process of measuring stress – strain states by the method of speckle-shearing interferometry using the phase-shift technique. A computer model with the possibility of setting its strain and roughness is developed, which includes a model of a diffusely reflecting test object corresponding to the characteristics of a real membrane, as well as a speckle interferometer model that allows speckle interferograms to be obtained for different speckle sizes and angles between interfering beams. The process of reconstructing the object surface topogram from model speckle interferograms by the phase-shift technique is implemented. Using the developed models, a two-dimensional shearogram are obtained, which is a derivative of the strain field of a circular membrane. Comparison of the results of numerical simulation with experimental data shows that the differences (rms deviations) do not exceed 0.02 μm. It is also shown that the error of interferogram reconstruction by the phase-shift technique increases significantly when the test object strains exceed 12 μm. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Radiation damage and dose limits in serial synchrotron crystallography at cryo- and room temperatures

Radiation damage limits the accuracy of macromolecular structures in X-ray crystallography. Cryogenic (cryo-) cooling reduces the global radiation damage rate and, therefore, became the method of choice over the past decades. The recent advent of serial crystallography, which spreads the absorbed energy over many crystals, thereby reducing damage, has rendered room temperature (RT) data collection more practical and also extendable to microcrystals, both enabling and requiring the study of specific and global radiation damage at RT. Here, we performed sequential serial raster-scanning crystallography using a microfocused synchrotron beam that allowed for the collection of two series of 40 and 90 full datasets at 2- and 1.9-Å resolution at a dose rate of 40.3 MGy/s on hen egg white lysozyme (HEWL) crystals at RT and cryotemperature, respectively. The diffraction intensity halved its initial value at average doses ( D 1/2 ) of 0.57 and 15.3 MGy at RT and 100 K, respectively. Specific radiation damage at RT was observed at disulfide bonds but not at acidic residues, increasing and then apparently reversing, a peculiar behavior that can be modeled by accounting for differential diffraction intensity decay due to the nonuniform illumination by the X-ray beam. Specific damage to disulfide bonds is evident early on at RT and proceeds at a fivefold higher rate than global damage. The decay modeling suggests it is advisable not to exceed a dose of 0.38 MGy per dataset in static and time-resolved synchrotron crystallography experiments at RT. This rough yardstick might change for proteins other than HEWL and at resolutions other than 2 Å.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Selective reduction of CO to acetaldehyde with CuAg electrocatalysts

Electrochemical CO reduction can serve as a sequential step in the transformation of CO 2 into multicarbon fuels and chemicals. In this study, we provide insights on how to steer selectivity for CO reduction almost exclusively toward a single multicarbon oxygenate by carefully controlling the catalyst composition and its surrounding reaction conditions. Under alkaline reaction conditions, we demonstrate that planar CuAg electrodes can reduce CO to acetaldehyde with over 50% Faradaic efficiency and over 90% selectivity on a carbon basis at a modest electrode potential of −0.536 V vs. the reversible hydrogen electrode. The Faradaic efficiency to acetaldehyde was further enhanced to 70% by increasing the roughness factor of the CuAg electrode. Density functional theory calculations indicate that Ag ad-atoms on Cu weaken the binding energy of the reduced acetaldehyde intermediate and inhibit its further reduction to ethanol, demonstrating that the improved selectivity to acetaldehyde is due to the electronic effect from Ag incorporation. These findings will aid in the design of catalysts that are able to guide complex reaction networks and achieve high selectivity for the desired product.

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

How a raindrop gets shattered on biological surfaces

Many biological surfaces of animals and plants (e.g., bird feathers, insect wings, plant leaves, etc.) are superhydrophobic with rough surfaces at different length scales. Previous studies have focused on a simple drop-bouncing behavior on biological surfaces with low-speed impacts. However, we observed that an impacting drop at high speeds exhibits more complicated dynamics with unexpected shock-like patterns: Hundreds of shock-like waves are formed on the spreading drop, and the drop is then abruptly fragmented along with multiple nucleating holes. Such drop dynamics result in the rapid retraction of the spreading drop and thereby a more than twofold decrease in contact time. Our results may shed light on potential biological advantages of hypothermia risk reduction for endothermic animals and spore spreading enhancement for fungi via wave-induced drop fragmentation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Large ecosystem-scale effects of restoration fail to mitigate impacts of land-use legacies in longleaf pine savannas

Ecological restoration is a global priority, with potential to reverse biodiversity declines and promote ecosystem functioning. Yet, successful restoration is challenged by lingering legacies of past land-use activities, which are pervasive on lands available for restoration. Although legacies can persist for centuries following cessation of human land uses such as agriculture, we currently lack understanding of how land-use legacies affect entire ecosystems, how they influence restoration outcomes, or whether restoration can mitigate legacy effects. Using a large-scale experiment, we evaluated how restoration by tree thinning and land-use legacies from prior cultivation and subsequent conversion to pine plantations affect fire-suppressed longleaf pine savannas. In this work, we evaluated 45 ecological properties across four categories: 1) abiotic attributes, 2) organism abundances, 3) species diversity, and 4) species interactions. The effects of restoration and land-use legacies were pervasive, shaping all categories of properties, with restoration effects roughly twice the magnitude of legacy effects. Restoration effects were of comparable magnitude in savannas with and without a history of intensive human land use; however, restoration did not mitigate numerous legacy effects present prior to restoration. As a result, savannas with a history of intensive human land use supported altered properties, especially related to soils, even after restoration. The signature of past human land-use activities can be remarkably persistent in the face of intensive restoration, influencing the outcome of restoration across diverse ecological properties. Understanding and mitigating land-use legacies will maximize the potential to restore degraded ecosystems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

The mesoscale order of nacreous pearls

A pearl’s distinguished beauty and toughness are attributable to the periodic stacking of aragonite tablets known as nacre. Nacre has naturally occurring mesoscale periodicity that remarkably arises in the absence of discrete translational symmetry. Gleaning the inspiring biomineral design of a pearl requires quantifying its structural coherence and understanding the stochastic processes that influence formation. By characterizing the entire structure of pearls (~3 mm) in a cross-section at high resolution, we show that nacre has medium-range mesoscale periodicity. Self-correcting growth mechanisms actively remedy disorder and topological defects of the tablets and act as a countervailing process to long-range disorder. Nacre has a correlation length of roughly 16 tablets (~5.5 µm) despite persistent fluctuations and topological defects. For longer distances (>25 tablets ~8.5 µm), the frequency spectrum of nacre tablets follows f –1.5 behavior, suggesting that growth is coupled to external stochastic processes—a universality found across disparate natural phenomena, which now includes pearls.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ground-state phase diagram of the t-t ′ -J model

Significance Theoretical efforts to understand the superconducting cuprates have focused on a few simplified but challenging models: notably the Hubbard and the t - t ′ - J models. Using density matrix renormalization group (DMRG) simulations on width-8 cylinders, we study the ground-state phase diagram of the t - t ′- J model. The next nearest-neighbor hopping t ′ allows one to distinguish between hole ( t ′ < 0 ) and electron doping ( t ′ > 0 ). We find for the magnetic, charge, and single-particle properties that this model roughly captures the differences between the electron- and hole-doped cuprates. However, the superconducting properties differ. For hole doping ( t ′ < 0 ), there is an absence of superconductivity, while for electron doping ( t ′ > 0), we find d -wave superconductivity coexisting with antiferromagnetism and a secondary triplet p -wave pairing component.

Jiang, Shengtao↗
DOE OSTI2023

Exceptional stratospheric contribution to human fingerprints on atmospheric temperature

In 1967, scientists used a simple climate model to predict that human-caused increases in atmospheric CO 2 should warm Earth’s troposphere and cool the stratosphere. This important signature of anthropogenic climate change has been documented in weather balloon and satellite temperature measurements extending from near-surface to the lower stratosphere. Stratospheric cooling has also been confirmed in the mid to upper stratosphere, a layer extending from roughly 25 to 50 km above the Earth’s surface (S 25 − 50 ). To date, however, S 25 − 50 temperatures have not been used in pattern-based attribution studies of anthropogenic climate change. Here, we perform such a “fingerprint” study with satellite-derived patterns of temperature change that extend from the lower troposphere to the upper stratosphere. Including S 25 − 50 information increases signal-to-noise ratios by a factor of five, markedly enhancing fingerprint detectability. Key features of this global-scale human fingerprint include stratospheric cooling and tropospheric warming at all latitudes, with stratospheric cooling amplifying with height. In contrast, the dominant modes of internal variability in S 25 − 50 have smaller-scale temperature changes and lack uniform sign. These pronounced spatial differences between S 25 − 50 signal and noise patterns are accompanied by large cooling of S 25 − 50 (1 to 2 ° C over 1986 to 2022) and low S 25 − 50 noise levels. Our results explain why extending “vertical fingerprinting” to the mid to upper stratosphere yields incontrovertible evidence of human effects on the thermal structure of Earth’s atmosphere.

climate modeling↗
DOE OSTI2024

Quantum computation of stopping power for inertial fusion target design

Stopping power is the rate at which a material absorbs the kinetic energy of a charged particle passing through it—one of many properties needed over a wide range of thermodynamic conditions in modeling inertial fusion implosions. First-principles stopping calculations are classically challenging because they involve the dynamics of large electronic systems far from equilibrium, with accuracies that are particularly difficult to constrain and assess in the warm-dense conditions preceding ignition. Here, we describe a protocol for using a fault-tolerant quantum computer to calculate stopping power from a first-quantized representation of the electrons and projectile. Our approach builds upon the electronic structure block encodings of Su et al. [ PRX Quant. 2 , 040332 (2021)], adapting and optimizing those algorithms to estimate observables of interest from the non-Born–Oppenheimer dynamics of multiple particle species at finite temperature. We also work out the constant factors associated with an implementation of a high-order Trotter approach to simulating a grid representation of these systems. Ultimately, we report logical qubit requirements and leading-order Toffoli costs for computing the stopping power of various projectile/target combinations relevant to interpreting and designing inertial fusion experiments. We estimate that scientifically interesting and classically intractable stopping power calculations can be quantum simulated with roughly the same number of logical qubits and about one hundred times more Toffoli gates than is required for state-of-the-art quantum simulations of industrially relevant molecules such as FeMoco or P450.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Observational evidence for groundwater influence on crop yields in the United States

As climate change shifts crop exposure to dry and wet extremes, a better understanding of factors governing crop response is needed. Recent studies identified shallow groundwater—groundwater within or near the crop rooting zone—as influential, yet existing evidence is largely based on theoretical crop model simulations, indirect or static groundwater data, or small-scale field studies. Here, we use observational satellite yield data and dynamic water table simulations from 1999 to 2018 to provide field-scale evidence for shallow groundwater effects on maize yields across the United States Corn Belt. We identify three lines of evidence supporting groundwater influence: 1) crop model simulations better match observed yields after improvements in groundwater representation; 2) machine learning analysis of observed yields and modeled groundwater levels reveals a subsidy zone between 1.1 and 2.5 m depths, with yield penalties at shallower depths and no effect at deeper depths; and 3) locations with groundwater typically in the subsidy zone display higher yield stability across time. We estimate an average 3.4% yield increase when groundwater levels are at optimum depth, and this effect roughly doubles in dry conditions. Groundwater yield subsidies occur ~35% of years on average across locations, with 75% of the region benefitting in at least 10% of years. Overall, we estimate that groundwater-yield interactions had a net monetary contribution of approximately $10 billion from 1999 to 2018. This study provides empirical evidence for region-wide groundwater yield impacts and further underlines the need for better quantification of groundwater levels and their dynamic responses to short- and long-term weather conditions.

60 APPLIED LIFE SCIENCES↗