Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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

Simulating the mechanisms of serrated flow in interstitial alloys with atomic resolution over diffusive timescales

Abstract The Portevin-Le Chatelier (PLC) effect is a phenomenon by which plastic slip in metallic materials becomes unstable, resulting in jerky flow and the onset of inhomogeneous deformation. The PLC effect is thought to be fundamentally caused by the dynamic interplay between dislocations and solute atoms. However, this interplay is almost always inaccessible experimentally due to the extremely fine length and time scales over which it occurs. In this paper, simulations of jerky flow in W-O interstitial solid solutions reveal three dynamic regimes emerging from the simulated strain rate-temperature space: one resembling standard solid solution strengthening, another one mimicking solute cloud formation, and a third one where dislocation/solute coevolution leads to jerky flow as a precursor of dynamic strain aging. The simulations are carried out in a stochastic framework that naturally captures rare events in a rigorous manner, providing atomistic resolution over diffusive time scales using no adjustable parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Beyond the universal Dyson singularity for 1-D chains with hopping disorder

We study a simple non-interacting nearest neighbor tight-binding model in one dimension with disorder, where the hopping terms are chosen randomly. This model exhibits a well-known singularity at the band center both in the density of states and localization length. If the probability distribution of the hopping terms is well-behaved, then the singularities exhibit universal behavior, the functional form of which was first discovered by Freeman Dyson in the context of a chain of classical harmonic oscillators. We show here that this universal form can be violated in a tunable manner if the hopping elements are chosen from a divergent probability distribution. We also demonstrate a connection between a breakdown of universality in this quantum problem and an analogous scenario in the classical domain — that of random walks and diffusion with anomalous exponents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A reactive molecular dynamics model for uranium/hydrogen containing systems

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. Here, we show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Capacitive response of biological membranes

We present a minimal model to analyze the capacitive response of a biological membrane subjected to a step voltage via blocking electrodes. Through a perturbative analysis of the underlying electrolyte transport equations, we show that the leading-order relaxation of the transmembrane potential is governed by a capacitive timescale, τ_{C}=λ_{D}L/D(2+Γδ^{M}/L/4+Γδ^{M}/λ_{D}), where λ_{D} is the Debye screening length, L is the electrolyte width, Γ is the ratio of the permittivity of the electrolyte to the membrane, δ^{M} is the membrane thickness, and D is the ionic diffusivity. This timescale is considerably shorter than the traditional RC timescale λ_{D}L/D for a bare electrolyte due to the membrane's low permittivity and finite thickness. Beyond the linear regime, however, salt diffusion in the bulk electrolyte drives a secondary, nonlinear relaxation process of the transmembrane potential over a longer timescale τ_{L}=L^{2}/4π^{2}D. A simple equivalent-circuit model accurately captures the linear behavior, and the perturbation expansion remains applicable across the entire range of observed physiological transmembrane potentials. Together, these findings underscore the importance of the faster capacitive timescale and nonlinear effects on the bulk diffusion timescale in determining transmembrane potential dynamics for a range of biological systems.

Farhadi, Jafar↗
DOE OSTI2025

Analysis of Liquids Using the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER)

We developed a sensor called the Submerged Plasma for Isotopic Detection and Elemental Resolution (SPIDER) probe, which uses an atmospheric pressure glow discharge below the surface of liquids to excite species in the liquid. Through emission spectroscopy of molten salts, liquid metals, and heavy water, we demonstrated the SPIDER probe’s high resolution, accuracy, and versatility. We successfully identified trace concentrations of transition and rare-earth metals in molten salts and detected the isotopic shift of the H β → D β emission line. Our analysis revealed unconventional spectral alkali line shapes, indicating two competing excitation modes: film explosion and droplet vaporization. The film explosion mode, characterized by dense plasma, exhibited self-reversal and broadband continuum emission, while the droplet vaporization mode, associated with diffusive plasma, produced narrow-line emissions. Furthermore, by analyzing circuit transients alongside individual plasma events, we observed that film explosions generate higher currents, likely due to a shorter plasma length as the current preferentially flows through the thin liquid layer. Altogether, our results highlight the SPIDER probe’s efficacy and flexibility, making it well-suited for online material quantification of liquids in extreme environments.

AES↗
DOE OSTI2025

Molecular aspect ratio effect on axial thermal transport in solution-spun carbon nanotube fibers

Neat, densely packed, and highly aligned carbon nanotube fibers (CNTFs) have appealing room-temperature axial thermal conductivity (⁠k⁠) and thermal diffusivity (⁠α⁠) for applications in lightweight heat spreading, flexible thermal connections, and thermoelectric active cooling. Although CNTFs are regularly produced from different input carbon nanotubes (CNTs), prior work has not quantified how the CNT molecular aspect ratio r (i.e., molecular length-to-diameter ratio) influences k and in well-aligned, packed CNTFs. Here, we perform self-heated steady-state and three-omega thermal measurements at room temperature on CNTF suspended in vacuum. Our results show that k increases from 150 to 380 W / mK for viscosity-averaged molecular aspect ratios increasing from r = 960 to 5600 and nanotube diameters of ~2 nm, which we attribute to the effects of thermal resistances between CNT bundles. CNTFs made with varying volume fraction Φ of constituent high-r and low-r CNT have properties that fall within or below the typical macroscopic rule-of-mixtures bounds. The thermal diffusivity α scales with k, leading to a sample-averaged volumetric heat capacity of 1.5 ± 0.3 MJ / m 3 K⁠. This work's findings that fibers made from longer CNT have larger k and α at room temperature motivate further investigation into thermal transport in solution-spun CNTF.

42 ENGINEERING↗
DOE OSTI2025

Nonlinear programming optimization of a single-stack electrodialysis desalination system for cost efficiency

Electrodialysis (ED) presents a competitive method for desalinating brackish waters. In this work, we perform cost optimization of a single-stack ED system across a range of feed salinities and water recoveries while optimizing operating voltage, number of cell pairs, and cell length. The results of our optimization show that the levelized cost of water (LCOW) increases with an increase in feed salinity. The outcomes of our optimization show that cost-optimal design generally increases cell length while decreasing cell pair number and operating voltage with an increase in salinity. These trends are nonlinear, with the number of cell pairs and applied voltage exhibiting local maxima when operating at low salinity and high recovery. We discuss the underlying mechanism for cell length becoming a leveraging design parameter by inspecting the length-dependent profiles of key electrochemical properties of the ED cell. Finally, we present how increasing performance metrics and decreasing costs impact LCOW, demonstrating that innovations that decrease counter-current diffusion have resulted in the highest decrease of LCOW.

42 ENGINEERING↗
DOE OSTI2022

The nano- and meso-scale structure of amorphous calcium carbonate

Abstract Understanding the underlying processes of biomineralization is crucial to a range of disciplines allowing us to quantify the effects of climate change on marine organisms, decipher the details of paleoclimate records and advance the development of biomimetic materials. Many biological minerals form via intermediate amorphous phases, which are hard to characterize due to their transient nature and a lack of long-range order. Here, using Monte Carlo simulations constrained by X-ray and neutron scattering data together with model building, we demonstrate a method for determining the structure of these intermediates with a study of amorphous calcium carbonate (ACC) which is a precursor in the bio-formation of crystalline calcium carbonates. We find that ACC consists of highly ordered anhydrous nano-domains of approx. 2 nm that can be described as nanocrystalline. These nano-domains are held together by an interstitial net-like matrix of water molecules which generate, on the mesoscale, a heterogeneous and gel-like structure of ACC. We probed the structural stability and dynamics of our model on the nanosecond timescale by molecular dynamics simulations. These simulations revealed a gel-like and glassy nature of ACC due to the water molecules and carbonate ions in the interstitial matrix featuring pronounced orientational and translational flexibility. This allows for viscous mobility with diffusion constants four to five orders of magnitude lower than those observed in solutions. Small and ultra-small angle neutron scattering indicates a hierarchically-ordered organization of ACC across length scales that allow us, based on our nano-domain model, to build a comprehensive picture of ACC formation by cluster assembly from solution. This contribution provides a new atomic-scale understanding of ACC and provides a framework for the general exploration of biomineralization and biomimetic processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗
DOE OSTI2026

Generative modeling enables molecular structure retrieval from Coulomb explosion imaging

Capturing the structural changes that molecules undergo during chemical reactions in real space and time is a long-standing dream and an essential prerequisite for understanding and ultimately controlling femtochemistry. A key approach to tackle this challenging task is Coulomb explosion imaging, which benefited decisively from recently emerging high-repetition-rate X-ray free-electron laser sources. With this technique, information on the molecular structure is inferred from the momentum distributions of the ions produced by the rapid Coulomb explosion of molecules. Retrieving molecular structures from these distributions poses a highly non-linear inverse problem that remains unsolved for molecules consisting of more than a few atoms. Here, we address this challenge using a diffusion-based Transformer neural network. We show that the network reconstructs unknown molecular geometries from ion-momentum distributions with a mean absolute error below one Bohr radius, which is half the length of a typical chemical bond.

Artificial Intelligence (cs.AI)↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Weak localization and universal conductance fluctuations in large-area twisted bilayer graphene

We study diffusive magnetotransport in highly 𝑝-doped large-area twisted bilayer graphene in 1∘, 7∘, 9∘, and 20∘ samples. All samples exhibit weak localization, from which we extract the phase coherence length and intervalley scattering lengths, and from that determine that dephasing is caused by electron-electron scattering and intervalley scattering is caused by point defects. We observe signatures of universal conductance fluctuations in the 9∘ sample, which has high mobility and is near the van Hove singularity. Further improvements in sample quality and applications to large-area moiré materials will open new avenues to observe quantum interference effects.

Talkington, Spenser [University of Pennsylvania]↗
DOE OSTI2022

The discovery of a radio galaxy of at least 5 Mpc

Context. Giant radio galaxies (GRGs, or colloquially ‘giants’) are the Universe’s largest structures generated by individual galaxies. They comprise synchrotron-radiating active galactic nucleus ejecta and attain cosmological (megaparsec-scale) lengths. However, the main mechanisms that drive their exceptional growth remain poorly understood. Aims. To deduce the main mechanisms that drive a phenomenon, it is usually instructive to study extreme examples. If there exist host galaxy characteristics that are an important cause for GRG growth, then the hosts of the largest GRGs are likely to possess them. Similarly, if there exist particular large-scale environments that are highly conducive to GRG growth, then the largest GRGs are likely to reside in them. For these reasons, we aim to perform a case study of the largest GRG available. Methods. We reprocessed the LOFAR Two-Metre Sky Survey DR2 by subtracting compact sources and performing multi-scale CLEAN de-convolutions at 60" and 90" resolution. The resulting images constitute the most sensitive survey yet for radio galaxy lobes, whose diffuse nature and steep synchrotron spectra have allowed them to evade previous detection attempts at higher resolution and shorter wavelengths. We visually searched these images for GRGs. Results. We have discovered Alcyoneus, a low-excitation radio galaxy with a projected proper length l p = 4.99 ± 0.04 Mpc. Both its jets and lobes are detected at very high significance, and the SDSS-based identification of the host, at spectroscopic redshift z spec = 0.24674 ± 6 × 10 -5 , is unambiguous. The total luminosity density at ν = 144 MHz is L ν = 8 ± 1 × 10 25 W Hz -1 , which is below average, though near median (percentile 45 ± 3%) for GRGs. The host is an elliptical galaxy with a stellar mass M* = 2.4 ± 0.4 × 10 11 M ⊙ and a super-massive black hole mass M• = 4 ± 2 × 10 8 M ⊙ , both of which tend towards the lower end of their respective GRG distributions (percentiles 25 ± 9% and 23 ± 11%). The host resides in a filament of the Cosmic Web. Through a new Bayesian model for radio galaxy lobes in three dimensions, we estimate the pressures in the megaparsec-cubed-scale northern and southern lobes to be P min,1 = 4.8 ± 0.3 × 10 -16 Pa and P min,2 = 4.9 ± 0.6 × 10 -16 Pa, respectively. The corresponding magnetic field strengths are B min,1 = 46 ± 1 pT and B min,2 = 46 ± 3 pT. Conclusions. We have discovered what is in projection the largest known structure made by a single galaxy – a GRG with a projected proper length l p = 4.99 ± 0.04 Mpc. The true proper length is at least l min = 5.04 ± 0.05 Mpc. Beyond geometry, Alcyoneus and its host are suspiciously ordinary: the total low-frequency luminosity density, stellar mass, and super-massive black hole mass are all lower than, though similar to, those of the medial GRG. Thus, very massive galaxies or central black holes are not necessary to grow large giants, and, if the observed state is representative of the source over its lifetime, neither is high radio power. A low-density environment remains a possible explanation. The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. The pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents the most promising radio galaxy yet to probe the warm–hot inter-galactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nanoparticle Diffusion in Miscible Polymer Nanocomposite Melts

X-ray photon correlation spectroscopy measurements were used to quantify the dynamics of bare and bimodal grafted silica nanoparticles mixed with PEO melts of different molecular weights. In dilute polymer nanocomposite (PNC) samples, we find diffusive NP behavior as described by the Stokes–Einstein relationship so long as the adsorbed PEO polymer layer is taken into account in determining both the effective NP size and its role on composite viscosity. The size of this bound layer was found to be approximately 2R g , where R g is the chain radius of gyration. We also expanded our system to investigate how the dynamics of grafted NPs differ from bare NPs with an adsorbed layer. We showed that the dynamics again can be determined by an effective NP radius at a scale smaller than the effective interparticle spacing; however, at larger length scales, the morphology and grafting parameters play a major role in the system dynamics. These results allow us to quantify NP ordering driven by polymer crystallization. It has previously been speculated that behavior is controlled by the relative ratio of time scale of crystal growth and the diffusive time scale of the NPs, a Peclet number. When the former time scale is longer, then the NPs are expected to be segregated into the interlamellar amorphous zones, while the NPs are too slow to be reorganized in the opposite case. We show here that this conjecture is quantitatively correct and the demarcation in behavior occurs for Pe = 1. Thus, we provide a way to estimate a critical spherulite growth rate for any semicrystalline PNC, at which a given NP can be ordered. Together, the results of this study permit us to tunably design PNCs through directed dispersion of NPs in a polymer matrix.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Stereochemically‐Controlled Fluorinated Copolymers for Selectively Permeable Barrier Applications

Selective oxygen permeability coupled with low water vapor transmission is essential for biomedical and packaging applications requiring controlled oxygen flux under humid conditions. However, most high‐performance barrier polymers depend on perfluoroalkyl substances (PFAS), whose persistence and regulatory restrictions limit their long‐term applicability. We designed a series of stereocontrolled thiol‐yne‐based polyesters, including both fluorinated and non‐fluorinated variants, for selective oxygen permeability with considerable water barrier performance. Tailoring polymer crystallinity and morphology tuned both oxygen transport and mechanical properties. Fluorinated polymers demonstrated enhanced hydrophobicity and water resistance while maintaining oxygen diffusivity within a range relevant to oxygen‐sensing applications. Structure–property relationships were elucidated through small‐ and wide‐angle X‐ray scattering, revealing semi‐crystalline domains influenced by fluorine content and dithiol chain length. Barrier performance was rigorously evaluated via water vapor transmission rate and dynamic vapor sorption, showing reduced water uptake with increasing dithiol monomer length and crystallinity. In conclusion, this work introduces a PFAS‐free alternative to conventional barrier materials and establishes a tunable materials platform with potential relevance for biomedical devices and packaging systems requiring controlled oxygen permeability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Laser powder bed fusion of Inconel 718 on 316 stainless steel

The joining of Inconel 718 and 316 stainless steel has notable industrial importance. However, traditional arc-welding has many issues. In particular, the formation of the detrimental Nb-rich Laves phase in the heat affected zone can promote crack initiation and propagation, which leads to poor mechanical properties. The small heat input and the rapid solidification of laser powder bed fusion (PBF) has the potential to suppress this problem. Despite the infancy of the multi-material PBF, it is of interest to investigate the relationship among the process parameter, the microstructure at the materials joining point, and the mechanical performance. In this study, we used high-speed x-ray imaging to in-situ monitor the printing of Inconel 718 tracks on a 316H SS substrate via laser powder bed fusion, followed by post-build characterization and testing with SEM, XRD and nanoindentation. We observed that a higher linear energy density resulted in a more diffused interface and a smaller compositional gradient, but as a trade-off, generated more keyhole pores and thermal cracks. Furthermore, we also showed that the built layer exhibited chemical inhomogeneity at two length scales: a tens-of-mu m-scale compositional vortex due to the melt flow, and a sub-mu m-scale Nb-rich dendritic structure due to the non-equilibrium solidification that led to the formation of the Laves phase. It was inferred that the chemical inhomogeneity may benefit the mechanical properties at the interface by providing more interlocking between the two materials and also by limiting the effect of the brittle Laves phase to potentially short range.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of Three Transport Coefficients and the Thermodynamic Factor in Block Copolymer Electrolytes with Different Morphologies

The design and engineering of composite materials is one strategy to satisfy the materials needs of systems with multiple orthogonal property requirements. In the case of rechargeable batteries with lithium metal anodes, the system requires a separator with fast lithium ion transport and good mechanical strength. In this work, we focus on the system polystyrene-block-poly(ethylene oxide) (SEO) with bis(trifluoromethane)sulfonimide lithium salt (LiTFSI). Ion transport occurs in the salt-containing poly(ethylene oxide)-rich domains. Mechanical rigidity arises due to the glassy nature of polystyrene (PS). If we assume that the salt does not interact with the PS-rich domains, we can describe ion transport in the electrolyte by three transport parameters (ionic conductivity, κ, salt diffusion coefficient, D, and cation transference number, t + 0 ) and a thermodynamic factor, T f . By systematically varying the volume fraction of the conducting phase, Φ c between 0.29 and 1.0, and chain length, N between 80 and 8000, we elucidate the role of morphology on ion transport. We find that κ is the strongest function of morphology, varying by three full orders of magnitude, while D is a weaker function of morphology. To calculate t + 0 and T f , we measure the current fraction, ρ + , and the open circuit potential, U, of concentration cells. We find that ρ + and U follow universal trends as a function of salt concentration, regardless of chain length, morphology, or Φ c , allowing us to calculate t + 0 for any SEO/LiTFSI or PEO/LiTFSI mixture when κ and D are known. Finally, the framework developed in this paper enables predicting the performance of any block copolymer electrolyte in a rechargeable battery.

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

Compact Diffusion Bonded Heat Exchanger Fatigue Life Simulations (Final Report)

Compact diffusion bonded heat exchangers are essential for high pressure heat exchange, but they are subject to thermal fatigue and ramp rate limitations. Simulation of these geometries is challenging with a large range of length and time scales from thousands of mm-sized microchannels inside a m-sized heat exchanger. Multi-physics simulations including thermal, fluid, and solid mechanics components are being used to predict stress within the heat exchangers under these conditions. These predictions can then be used to understand thermal ramp rate limitations while keeping maximum stresses low as well as fatigue life predictions from well-known empirical models.

42 ENGINEERING↗