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 865 records · Page 48

Diffusive crack growth at a bimaterial interface

The diffusional microcrack growth behavior in a bimaterial system is investigated with an aim at estimating service life of advanced ceramic composites under creep-rupture conditions. The crack is assumed to grow via a coupled surface and grain-boundary diffusion under steady state conditions. The tensile stress distribution along the interface ahead of the moving crack tip is solved, and it is found that a new length parameter exists as a scaling factor for which the solution becomes identical to the single phase case when plotted on the nondimensional physical plane. In contrast to the elastic stress solution which shows singularity at the tip, together with oscillatory character away from the tip, the creep stresses have a peak value away from the tip due to a wedging effect and interfacial sliding eliminates stress oscillation resulting in a decoupling between mode I and mode II loading. This solution ties the far-field loading parameter to the crack tip conditions in terms of the unknown crack velocity to give a specific V-K(1) relationship. It is shown that an exponent of 12 in the conventional crack growth power law emerges at the higher applied stress range.

Chuang, Tze-Jer↗

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

Shapes of Nonbuoyant Round Luminous Hydrocarbon/Air Laminar Jet Diffusion Flames

The shapes (luminous flame boundaries) of round luminous nonbuoyant soot-containing hydrocarbon/air laminar jet diffusion flames at microgravity were found from color video images obtained on orbit in the Space Shuttle Columbia. Test conditions included ethylene- and propane-fueled flames burning in still air at an ambient temperature of 300 K, ambient pressures of 35-130 kPa, initial jet diameters of 1.6 and 2.7 mm, and jet exit Reynolds numbers of 45-170. Present test times were 100-200 s and yielded steady axisymmetric flames that were close to the laminar smoke point (including flames both emitting and not emitting soot) with luminous flame lengths of 15-63 mm. The present soot-containing flames had larger luminous flame lengths than earlier ground-based observations having similar burner configurations: 40% larger than the luminous flame lengths of soot-containing low gravity flames observed using an aircraft (KC-135) facility due to reduced effects of accelerative disturbances and unsteadiness; roughly twice as large as the luminous flame lengths of soot-containing normal gravity flames due to the absence of effects of buoyant mixing and roughly twice as large as the luminous flame lengths of soot-free low gravity flames observed using drop tower facilities due to the presence of soot luminosity and possible reduced effects of unsteadiness. Simplified expressions to estimate the luminous flame boundaries of round nonbuoyant laminar jet diffusion flames were obtained from the classical analysis of Spalding (1979); this approach provided Successful Correlations of flame shapes for both soot-free and soot-containing flames, except when the soot-containing flames were in the opened-tip configuration that is reached at fuel flow rates near and greater than the laminar smoke point fuel flow rate.

Lin, K.-C.↗

Shapes of Nonbuoyant Round Luminous Hydrocarbon/Air Laminar Jet Diffusion Flames

The shapes (luminous flame boundaries) of round luminous nonbuoyant soot-containing hydrocarbon/air laminar jet diffusion flames at microgravity were found from color video images obtained on orbit in the Space Shuttle Columbia. Test conditions included ethylene- and propane-fueled flames burning in still air at an ambient temperature of 300 K ambient pressures of 35-130 kPa, initial jet diameters of 1.6 and 2.7 mm, and jet exit Reynolds numbers of 45-170. Present test times were 100-200 s and yielded steady axisymmetric flames that were close to the laminar smoke point (including flames both emitting and not emitting soot) with luminous flame lengths of 15-63 mm. The present soot-containing flames had larger luminous flame lengths than earlier ground-based observations having similar burner configurations: 40% larger than the luminous flame lengths of soot-containing low gravity flames observed using an aircraft (KC-135) facility due to reduced effects of accelerative disturbances and unsteadiness; roughly twice as large as the luminous flame lengths of soot-containing normal gravity flames due to the absence of effects of buoyant mixing and roughly twice as large as the luminous flame lengths of soot-free low gravity flames observed using drop tower facilities due to the presence of soot luminosity and possible reduced effects of unsteadiness, Simplified expressions to estimate the luminous flame boundaries of round nonbuoyant laminar jet diffusion flames were obtained from the classical analysis of Spalding; this approach provided successful correlations of flame shapes for both soot-free and soot-containing flames, except when the soot-containing flames were in the opened-tip configuration that is reached at fuel flow rates near and greater than the laminar smoke point fuel flow rate.

Lin, K.-C.↗

Measurement of turbulent velocity and bounds for thermal diffusivity in laser shock compressed foams by x-ray photon correlation spectroscopy

Experimental benchmarking of transport coefficients under extreme conditions is required for validation of differing theoretical models. To date, measurement of transport properties of dynamically compressed samples remains a challenge with only a limited number of studies able to quantify transport in high pressure and temperature matter. x-ray photon correlation spectroscopy utilizes coherent x-ray sources to measure time correlations of density fluctuations, thus providing measurements of length and timescale-dependent transport properties. Here, we present a first-of-a-kind experiment to conduct x-ray photon correlation spectroscopy in laser shock compression experiments. We report measurement of the turbulent velocity in the wake of a laser driven supersonic shock and place an upper bound on thermal diffusivity in a solid density plasma on nanosecond timescales.

High-energy-density plasmas↗

Buffered Coordination Modulation as a Means of Controlling Crystal Morphology and Molecular Diffusion in an Anisotropic Metal-Organic Framework

We know significant advances have been made in the synthesis of chemically selective environments within metal-organic frameworks, yet materials development and industrial implementation have been hindered by the inability to predictively control crystallite size and shape. One common strategy to control crystal growth is the inclusion of coordination modulators, which are molecular species designed to compete with the linker for metal coordination during synthesis. However, these modulators can simultaneously alter the pH of the reaction solution, an effect that can also significantly influence crystal morphology. Herein, noncoordinating buffers are used to independently control reaction pH during metal-organic framework synthesis, enabling direct interrogation of the role of the coordinating species on crystal growth. We demonstrate the efficacy of this strategy in the synthesis of low-dispersity single-crystals of the framework Co 2 (dobdc) (dobdc 4- = 2,5-dioxido-1,4-benzenedicarboxylate) in a pH 7-buffered solution using cobalt(II) acetate as the metal source. Density functional theory calculations reveal that acetate competitively binds to Co during crystallization, and by using a series of cobalt(II) salts with carboxylate anions of varying coordination strength, it is possible to control crystal growth along the c -direction. Finally, we use zero length column chromatography to show that crystal morphology has a direct impact on guest diffusional path length for the industrially important hydrocarbon m -xylene. Together, these results provide molecular-level insight into the use of modulators in governing crystallite morphology and a powerful strategy for the control of molecular diffusion rates within metal-organic frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating Multiscale Spatial Distribution of Sulfur in a CNT Scaffold and Its Impact on Li–S Cell Performance

Ever increasing demand on high energy density batteries positions sulfur as a very promising cathode material for next-generation energy storage due to its high theoretical capacity of 1675 mAh/g. Electronically sulfur is highly insulating, and therefore integration of a conducting framework such as a carbon nanotube (CNT) scaffold with sulfur is a key aspect of the cathode design. Despite numerous efforts dedicated to S-CNT cathode development, increasing sulfur loading to above 1 mg/cm 2 while maintaining the cycling stability of the Li–S cell remains challenging. This could be partly due to the lack of understanding of the spatial distribution of sulfur in the CNT matrix and its location with respect to the morphology of the CNT scaffold. We demonstrate herein that the sulfur has a hierarchical distribution in the CNT cathode at high sulfur loading (>5 mg/cm 2 ) spanning multiple length scales (from nanometer to submillimeter). Sulfur infiltration into the CNT rather than the sulfur loading plays a key role in determining the redox reaction kinetics, Li + ion diffusion, and the galvanostatic cycling capacity and stability of Li–S cells. This study provides new insights for the design and fabrication of high loading, binder-free sulfur–carbon-based cathode architectures for next-generation high energy Li–S batteries.

25 ENERGY STORAGE↗

Exploring the Structural Origins of Optically Efficient One-Dimensional Lead Halide Perovskite Nanostructures

Metal halide perovskites have excellent optoelectronic properties. This study aims to determine how the optoelectronic properties of a model perovskite, cesium lead bromide (CsPbBr3), change with length and thickness in one dimension (1D). By examining the photophysics of CsPbBr3 quantum dots (QDs), nanowires (NWs), and nanorods (NRs), we observe the influence of confinement, exciton diffusion, and trapping on their optical properties. Our findings reveal that exciton diffusion to trap states limits the photoluminescence quantum yield (PLQY) of 1D CsPbBr3 in the weakly confined regime (8-14 nm) and explains their long-lived exciton dynamics, while enhanced radiative rates contribute to achieving near-unity PLQY in the strongly confined regime (<7 nm). Consequently, blue-emitting, 2.4 nm-thick CsPbBr3 NRs were 3.6X more emissive than the conventional CsPbBr3 QDs. This study underscores how structural optimization can improve the optoelectronic performance of CsPbBr3 and provides insight into the complex interplay of radiative and nonradiative processes in 1D ionic semiconductors.

Oddo, Alexander M↗

Revisiting metal fluorides as lithium-ion battery cathodes

Metal fluorides, promising lithium-ion battery cathode materials, have been classified as conversion materials due to the reconstructive phase transitions widely presumed to occur upon lithiation. We challenge this view by studying FeF 3 using X-ray total scattering and electron diffraction techniques that measure structure over multiple length scales coupled with density functional theory calculations, and by revisiting prior experimental studies of FeF 2 and CuF 2 . Metal fluoride lithiation is instead dominated by diffusion-controlled displacement mechanisms, and a clear topological relationship between the metal fluoride F- sublattices and that of LiF is established. Initial lithiation of FeF3 forms FeF2 on the particle's surface, along with a cation-ordered and stacking-disordered phase, A-Li x Fe y F 3 , which is structurally related to alpha-/beta-LiMn 2+ Fe 3+ F 6 and which topotactically transforms to B- and then C-Li x Fe y F 3 , before forming LiF and Fe. Lithiation of FeF 2 and CuF 2 results in a buffer phase between FeF 2 /CuF 2 and LiF. Overall, the resulting principles will aid future developments of a wider range of isomorphic metal fluorides.

25 ENERGY STORAGE↗

Influence of controlled disorder on the dipolar spin-ice state of Ho-based pyrochlores

Pyrochlore magnets of the form 𝑅 2 ⁢𝐵 2 ⁢O 7 , in which rare-earth ions on the 𝑅 site form a three-dimensional network of corner-sharing tetrahedra, provide a canonical setting for geometrical frustration. Ho-based pyrochlores host a dipolar spin-ice ground state, characterized by Ising moments constrained by the ice rules and elementary excitations analogous to magnetic monopoles. Here, in this work, we examine how controlled chemical disorder influences this state by introducing site mixing on the nonmagnetic 𝐵 site in two compounds. Ho 2⁢ GaSbO 7 contains only Ga 3+ /Sb 5+ charge disorder, whereas Ho 2 ⁢ScSbO 7 exhibits both charge and substantial size disorder arising from the large ionic-radius mismatch between Sc 3+ and Sb 5+ . Although both materials retain the pyrochlore structure, neutron-scattering measurements reveal a reduced correlation length for the 𝑅/𝐵-site cation ordering and enhanced local structural distortions in Ho 2 ⁢ScSbO 7 . Despite these structural differences, bulk thermodynamic measurements and magnetic diffuse scattering demonstrate that both systems exhibit the defining signatures of a dipolar spin-ice state. Low-energy inelastic neutron spectroscopy further uncovers broad magnetic excitations that develop within the dipolar spin-ice regime, a feature absent in pristine Ho pyrochlores and indicative of disorder-induced splitting of the non-Kramers ground-state doublet. Together, these results show that controlled disorder generates tunable transverse-field-driven quantum fluctuations in Ho-based pyrochlores, although the dipolar spin-ice state is remarkably robust to this disorder.

magnetic anisotropy↗

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design and performance of controlled-diffusion stator compared with original double-circular-arc stator

The capabilities of two stators, one with controlled-diffusion (CD) blade sections and one with double-circular-arc (DCA) blade sections, were compared. A CD stator was designed and tested that had the same chord length but half the blades of the DCA stator. The same fan rotor (tip speed, 429 m/sec; pressure ratio, 1.65) was used with each stator row. The design and analysis system is briefly described. The overall stage and rotor performances with each stator are compared, as are selected blade element data. The minimum overall efficiency decrement across the stator was approximately 1 percentage point greater with the CD balde sections than with the DCA blade sections.

Gelder, Thomas F.↗

Design and performance of controlled-diffusion stator compared with original double-circular-arc stator

The capabilities of two stators, one with controlled-diffusion (CD) blade sections and one with double-circular-arc (DCA) blade sections, were compared. A CD stator was designed and tested that had the same chord length but half the blades of the DCA stator. The same fan rotor (tip speed, 429 m/sec; pressure ratio, 1.65) was used with each stator row. The design and analysis system is briefly described. The overall stage and rotor performances with each stator are compared, as are selected blade element data. The minimum overall efficiency decrement across the stator was approximately 1 percentage point greater with the CD blade sections than with the DCA blade sections.

Gelder, Thomas F.↗

Ballistic transport in planetary ring systems due to particle erosion mechanisms. I - Theory, numerical methods, and illustrative examples

Ballistic transport, defined as the net radial transport of mass and angular momentum due to exchanges of meteoroid hypersonic-impact ejecta by neighboring planetary ring regions on time-scales orders-of-magnitude shorter than the age of the solar system, is presently considered as a problem in mathematical physics. The preliminary results of a numerical scheme for following the combined effects of ballistic transport and viscous diffusion demonstrate that ballistic transport generates structure near sharp edges already present in the ring-mass distribution; the entire ring system ultimately develops an undulatory structure whose length scale is typically of the order of the radial excursion of the impact ejecta.

Durisen, Richard H.↗

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗

Mesospheric Odd Nitrogen Enhancements During Relativistic Electron Precipitation Events

The behavior of mesospheric odd nitrogen species during and following relativistic and diffuse auroral precipitation events is simulated, Below 75 km nitric oxide is enhanced in proportion to the ion pair production function associated with the electron precipitation and the length of the event. Nitrogen dioxide and nitric acid are also enhanced. At 65 km the percentage of odd nitrogen for N is 0.1%, HNO3 is 1.6%, NO2 is 15%, and NO is 83.3%. Between 75 and 85 km NO is depleted during particle events due to the faster destruction of NO by N relative to the production of NO by N reacting with O2. Recovery of NO depends on transport from the lower thermosphere, where NO is produced in abundant amounts during particle events.

Aikin, A. C.↗

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

42 ENGINEERING↗