Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessment of metadynamic recrystallization in single copper particle impacts by focused ion beam tomography

We study single Cu-on-Cu impacts relevant to cold spray deposition and quantitatively analyze the metadynamic recrystallization (mDRX) that takes place after the impact by virtue of lingering impact adiabatic heat. Unlike prior studies, the current full 3D tomographic analysis of the mDRX volume shows that mDRX is extremely common in such impacts, although it is often missed when examining 2D sections. We also report an unexpected trend: there is a “sour spot” for mDRX at velocities about 20–40 % above the velocity for particle adhesion. This non-monotonic trend is contrary to the expectations based on increasing adiabatic heating with velocity. With a schematic model, we show that the trend can be explained on the basis of heat transfer: cooling of the heat-affected region is limited by transport through the bonded regions at the particle-substrate interface. Thus, bonding has a prominent role in the heat dissipation process and the best bonded particles most rapidly bulk quench, avoiding mDRX. Here, the developed semi-empirical model aligns with the experimental findings and may help inform microstructural evolution during cold spray and post-spray processing.

FIB-SEM tomography↗
DOE OSTI2025

Temperature and dose effects on dislocation loops in self-ion irradiated high-purity iron

Body-centered cubic (BCC) Fe-based alloys are promising candidate materials for advanced nuclear reactors. However, a detailed understanding of irradiation induced dislocation loop microstructure development remains unresolved. It is a widespread belief that 〈001〉 loops become increasingly favorable over ½〈111〉 loops as irradiation temperature rises above ∼300 C. Unfortunately, the temperature effects on 〈001〉 loop have been primarily examined in in-situ irradiation on TEM thin foils but poorly explored on bulk Fe due to exceedingly limited experimental studies on bulk specimens, raising concerns about the potential influence of TEM thin foil artifacts on observed results. Here, in this study, we conducted experiments on ultra-high purity BCC Fe specimens irradiated with 6.7–8 MeV Fe ions over a wide temperature range on bulk samples. We investigated the effects of temperature (T irr = 250–500 °C), dose rate (10⁻⁵ to 10⁻³ dpa/s), and dose (0.35 to 3.5 dpa) on the formation and evolution of 〈001〉 and ½〈111〉 loops as well as cavity (void) formation. Post-irradiation Burgers vector analysis via g•b method on dislocation segments and loops revealed that 〈001〉 loop fraction does not show a monotonic positive correlation with irradiation temperature. Combined with previous and current theoretical as well as experimental findings, we explore the temperature effects on all existing models of 〈001〉 loop formation. We conclude that the prevailing reports regarding the dominance of 〈001〉 loops in Fe at elevated temperatures are mainly attributable to the loss of glissile ½〈111〉 clusters in TEM thin foil experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Estimating source-sink distributions and fluxes of reactive nitrogen and sulfur within a mixed forest canopy

The vertical source-sink distribution of air pollutants within and above forested canopies is necessary for describing the biological, physical, and chemical processes influencing the soil-vegetation-atmosphere exchange. Here, this study implemented inverse modeling methods to estimate the source-sink and flux profiles of reactive nitrogen (N) and sulfur (S) compounds from measurements of the mean concentration profiles of ammonia (NH 3 ), nitric acid (HNO 3 ), sulfur dioxide (SO 2 ), and particulate ammonium (NH 4 + ), nitrate (NO 3 − ), and sulfate (SO 4 2− ) at a forest site in the southern Appalachian Mountains. Three inverse approaches utilizing different approximations to scalar transport within the canopy were developed and evaluated against sensible heat flux measurements. The Eulerian model (EUL), which incorporates vertical velocity skewness, performed well in reproducing the turbulent heat fluxes and was subsequently used to calculate the chemical source-sink and flux profiles. Above-canopy fluxes of NH 3 were downward, indicating that the forest was a net sink of NH 3 . The soil/litter layer was both a source and a sink for NH 3 but the exchange rate at the forest floor was small. Fluxes of HNO 3 , SO 2 , NO 3 − , NH 4 + , and SO 4 2- were uni-directional (deposition only) between the air and the canopy/ground and increased monotonically from the forest floor to the canopy top. Crown foliage dominated the uptake of reactive N and S during the growing season, accounting for 80–90% of the total canopy-scale flux. Fluxes and canopy-ground partitioning estimated using the resistance-based Surface Tiled Aerosol and Gas Exchange (STAGE) model were generally comparable to EUL. The comparison highlights the need for improved parameterizations of litter exchange and NH 3 compensation points in resistance models for forest ecosystems. The findings here benefit the application of critical loads in forest ecosystems and guide further development of resistance-based exchange models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗
DOE OSTI2022

“Self-trapping” in solar cell hybrid inorganic-organic perovskite absorbers

In the simplest picture, a “self-trapped” polaron forms when an excess electron or hole deforms a crystal lattice, creating a potential well with bound states. Properties of self-trapped polarons in methylammonium lead iodide perovskite (MAPbI 3 ), which is widely used as solar cell absorber, are of great interest, and are a subject of ongoing investigations and debates concerning the existence of large polarons with the co-presence of metastable self-trapping. Herein, we employ a self-interaction-free density functional theory method to investigate the stability of small polarons in tetragonal MAPbI 3 phase. The electron small polaron is found to be unstable, while the hole small polaron is found to be metastable at realistic operation temperatures of solar cells. Further, the hole polaron is found to have a hole band close to the conduction band, which in conjunction with its metastability suggests that small polarons will have an appreciable effect on charge-carrier recombinations in MAPbI 3 . Further, in this study, we posit that the existence of the metastable polarons in addition to the large polarons may explain the experimentally observed non-monotonic temperature dependence of bimolecular charge-carrier recombination rate in tetragonal MAPbI 3 phase.

14 SOLAR ENERGY↗
DOE OSTI2021

Effect of substrate tilt on sputter-deposited AuTa films

Gold-tantalum alloy films are of interest for biomedical and magnetically-assisted inertial confinement fusion applications. However, growth mechanisms of such Au-Ta alloy films are not well understood. Here, we systematically study the effect of substrate tilt on properties of films deposited by direct current magnetron sputtering from a AuTa alloy target. Experimental data is correlated with distributions of landing energies and incident angles of sputtered and backscattered species calculated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. Results reveal that the deposition rate and film density monotonically decrease with increasing substrate tilt. Properties are similar for films with a tilt of 40° and below: films are amorphous, with an average compressive stress of 0.5 GPa, a density of 16 g/cm 3 , and electrical resistivity of 200 μΩ cm. The critical substrate tilt angle for the onset of porosity, a transition to a reduced stress state, and an increase in electrical resistivity lies between 40 and 60°. Based on Monte Carlo simulations, these changes in film properties are correlated with an increase in the average impact angle of the ballistic component of atomic flux. Implications of these results to deposition onto non-planar substrates are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mixed Tin-Titanium Oxides by Atomic Layer Deposition on Planar Substrates: Physical and Electronic Structure

Ternary atomic layer deposition of Sn-doped TiO x (Sn:TiO x ) and Ti-doped SnO x (Ti:SnO x ) was performed using tetrakisdimethylamidotitanium (IV) and tin(IV) with water at 150 °C on (1 0 0) Si and quartz substrates. The physical structure of the films was investigated using scanning electron, atomic force, and Kelvin probe force microscopies as well as spectroscopic ellipsometry, Raman spectroscopy, and X-ray diffraction. All as-deposited films are amorphous and possessed flat and conformal topographies with low surface roughness (<1 nm). Annealed Sn:TiO x and Ti:SnO x crystallize into the tetragonal anatase TiO 2 and cassiterite SnO 2 structures, respectively, without large variations to the film topography. The chemical composition of the as-deposited and annealed films was investigated through X-ray photoelectron spectroscopy. While dopants are homogenously mixed, the ternary ALD to form Sn:TiO x and Ti:SnO x surprisingly does not follow the rule of mixtures, despite similarities in both the growth rates of binary SnO x and TiO x and the precursors used. This finding is also supported by spectroscopic ellipsometry. Here, we postulate the deviation from the rule of mixtures is due to non-ideal surface chemistry of the binary ALD of SnO x . The optical bandgap of all films was quantified by the Tauc analysis and the valence band position was investigated by valence band X-ray photoelectron spectroscopy. Intra-bandgap states above the valence band maximum were observed for as-deposited Ti:SnO x and Sn:TiO x and ascribed to 5 s states from Sn 2+ ions. The presence of these intra-bandgap states induces enhanced visible absorption and bathochromic shifts to the apparent bandgap for all as-deposited films. Annealing the films removes the intra-bandgap states and increases the optical bandgap. Simplified band diagrams created using the data available demonstrated that the conduction and valence band positions did not vary between the positions of pure TiO 2 or SnO 2 with monotonic increases to the dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scanning tunneling microscope tip-induced formation of Bi bilayers on Bi 2 Te 3

We report the formation of Bi(111) bilayer islands and crater structures on Bi 2 Te 3 (111) surfaces induced by voltage pulses from a scanning tunneling microscope tip. Pulses above a threshold voltage (+3 V) produce craters ∼0.5μm in diameter, similar to the size of the tip. Redeposited material self-assembles into a network of atomically ordered islands with a lattice constant identical to the underlying Bi2Te3 surface. The island size monotonically decreases over several micrometers from the pulse site, until the pristine Bi 2 Te 3 surface is recovered. We assign these islands to Bi bilayer based on atomic resolution images, analysis of step heights, and tunneling spectroscopy. Here, the dependence of bilayer formation on bias polarity and the evidence for defect diffusion together suggest a mechanism driven by the interplay of field evaporation and tunneling-current-induced Joule heating.

Bi bilayer↗
DOE OSTI2023

Characteristics of upward-connecting-leader current leading to attachment in downward negative cloud-to-ground lightning strokes

For this work, we analyzed currents associated with upward connecting leaders (UCLs) initiated from the Kennedy Space Center Industrial Area Tower in two negative cloud-to-ground strokes that struck the tower. One stroke was also recorded using a high-speed video camera located 760 m from the KSC IAT. The NLDN-reported peak currents for the two strokes were – 31.7 and – 98.5 kA. During the UCL development phase the current waveforms exhibited a monotonically (quasi-exponentially) increasing “background” current overlaid with 10-μs scale pulses with a median amplitude of 51.1 A. The UCL current durations for the two strokes were 1039 and 449 μs, respectively. During the pre-attachment processes (UCL and slow front) the total negative charge effectively transferred to ground were 70.2 and 55 mC, respectively. For the stroke captured on high-speed video, the average line-charge-density for the 109-m long UCL was found to be 0.5 mC/m. The average UCL 2-D speed was 2.4 × 10 5 m/s, and it was observed to accelerate toward the downward leader prior to attachment. We observed that UCL-pulse amplitudes are larger, background currents are higher, and interpulse intervals are shorter at later times during UCLs, which can be attributed to the intensification of the local electric field due to the approaching downward negative leader. The median positive charge injected into the UCL by a pulse was 297 μC. The UCL associated with the higher peak-current stroke produced the highest injected pulse-charge values about three times sooner during its development, likely due to the 2–2.7 times faster average downward leader vertical speeds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Nematic cell alignment directs calcium waves in an epithelial monolayer

Tissues rely on supracellular signals to coordinate their cells over a long range. Two such tissuescale cues are calcium waves and patterns of cell-cell alignment or nematic order. During wound healing, for example, calcium waves propagate across a tissue to guide directed cell migration and reepithelialization. Defects in long-range cell-cell alignment, or nematic orientation, can act to localize morphogenetic events in a tissue. Although these two cues have been considered in isolation, we demonstrate a relationship in epithelial tissue between long-range calcium signaling and the cell’s nematic order: The speed of a wound-induced calcium wave depends monotonically on the angle between the wave vector and cell axis, with maximal wave speed occurring perpendicular to the tissue’s orientation. Including anisotropic di↵usive coupling between cells in a canonical reactiondi↵usion model recapitulates our measured calcium wave dynamics. Our model demonstrates how orientation defects can desynchronize information propagation across a tissue. A calcium wave front is bent around nematic defects, therefore cells the same distance from a wound can receive the calcium signal at di↵erent times. Our work elucidates how spatial patterns in global cell alignment can control collective communication via calcium signaling during development, wound healing, and disease.

Winterstrain, Annemarie C. [Brandeis University, W↗
DOE OSTI2024

Optimization-based, property-preserving algorithm for passive tracer transport

Here we present a new optimization-based property-preserving algorithm for passive tracer transport. The algorithm utilizes a semi-Lagrangian approach based on incremental remapping of the mass and the total tracer. However, unlike traditional semi-Lagrangian schemes, which remap the density and the tracer mixing ratio through monotone reconstruction or flux correction, we utilize an optimization-based remapping that enforces conservation and local bounds as optimization constraints. In so doing we separate accuracy considerations from preservation of physical properties to obtain a conservative, second-order accurate transport scheme that also has a notion of optimality. Moreover, we prove that the optimization-based algorithm preserves linear relationships between tracer mixing ratios. We illustrate the properties of the new algorithm using a series of standard tracer transport test problems in a plane and on a sphere.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

O 2 -Oxidation of Individual Graphite and Graphene Nanoparticles in the 1200 – 2200 K Range: Particle-to-Particle Variations and the Evolution of the Reaction Rates and Optical Properties

The kinetics for O 2 oxidation of individual graphite and graphene platelet nanoparticles (NPs) were studied as a function of temperature (1200 to 2200 K) at varying oxygen partial pressures, using a single nanoparticle mass spectrometry method. NP temperature (TNP) was measured by measuring the NP thermal emission spectra during the kinetics studies. The initial oxidation efficiency is found to peak in the 1200 to 1500 K range, dropping by an order of magnitude as TNP was increased above 2000 K. There were large NP-to-NP variations in the oxidation rates, attributed to variations in the NP surface structure. In addition, the oxidation efficiencies evolved, non-monotonically, as the NPs reacted, decreasing by factors of between 10 and 300. This evolution of reactivity is attributed to changes in the NP surface structure due to the combination of oxidation and annealing. In conclusion, the optical properties, including wavelength dependence of the emissivity and the absorption cross section for the 532 nm heating laser, also tended to evolve as the NPs oxidized, but differently for each individual NP, presumably reflecting differences in the initial structures and their evolution.

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

Chemical looping air separation with Sr 0.8 Ca0.2Fe 0.9 Co 0.1 O 3-δ perovskite sorbent: Packed bed modeling, verification, and optimization

Chemical looping air separation (CLAS) represents a promising approach for efficient O 2 production from the air. This present study aims at optimizing the absorber/desorber operations and the separation process with extensive experimental validation. Specifically, a one-dimensional packed bed model was developed to investigate the CLAS operation with a Sr 0.8 Ca 0.2 Fe 0.9 Co 0.1 O 3-δ perovskite sorbent. The redox thermodynamics of perovskite sorbent was measured by TGA and then incorporated into a linear driving force model to describe the O 2 absorption and desorption rates. Both 4-step and 5-step air separation cycle configurations, with various cyclic structures, were performed in a subpilot-scale packed bed. The model predicted O2 purity and productivity were consistent with experimental results, supporting its accuracy and applicability. Parametric analysis and multi-objective optimization were further carried out to assess the performance of CLAS. Both O 2 purity and recovery increased monotonically with the cycle time, airflow rate, steam flow rate, and absorption pressure. Meanwhile, optimal O 2 productivity and power consumption can only be achieved by specific combinations of these parameters. The optimized results showed that CLAS can be highly competitive when compared to conventional pressure swing adsorption (PSA) or cryogenic distillation. The 5-step cycle configuration achieved a minimum power consumption of 118 kW·h for producing 1 ton O 2 with ≥ 95% purity. The maximum O 2 productivity reached 0.0932 g O2 /(g sorbent ·h) with 390 kW·h/ton O 2 of energy consumption (95% pure). The optimization results also indicate that CLAS can potentially be more efficient than cryogenic distillation even when the required O 2 purity is above 99%.

42 ENGINEERING↗
DOE OSTI2025

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗
DOE OSTI2021

Hinderance of C-S-H sheet piling during first drying using a shrinkage reducing agent: A SAXS study

The role of a shrinkage reducing agent (SRA) in the hardened cement paste during a first desorption process was investigated by measuring the small-angle X-ray scattering (SAXS). Fractal disc analysis and Guinier approximation were adopted to interpret the SAXS profiles. We experimentally confirmed that the presence of SRA molecules mitigated the piling of calcium-silicate hydrate (C-S-H) sheets, which could be monitored by the increase in the thickness of the disc, representing the agglomeration of C-S-H sheets. In addition, the basal spacing of C-S-H did not change during the drying. By contrast, the diameter of the discs monotonically decreased and showed a similar trend to that of samples lacking SRA. Based on these findings, we concluded that SRA molecules were on the surface of the C-S-H sheets and mitigated the piling, facilitating water removal from the C-S-H agglomeration in the RH range of 80–40%.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗