Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion processes”

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,225 records · Page 68

On the mechanism of effective chemical reactions with turbulent mixing of reactants and finite rate of molecular reactions

A generalization of the theory of chemical transformation processes under turbulent mixing of reactants and arbitrary values of the rate of molecular reactions is presented that was previously developed for the variant of an instantaneous reaction [13]. The use of the features of instantaneous reactions when considering the general case, namely, the introduction of the concept of effective reaction for the reactant volumes and writing a closing conservation equation for these volumes, became possible due to the partition of the whole amount of reactants into “active” and “passive” classes; the reactants of the first class are not mixed and react by the mechanism of instantaneous reactions, while the reactants of the second class approach each other only through molecular diffusion, and therefore their contribution to the reaction process can be neglected. The physical mechanism of reaction for the limit regime of an ideal mixing reactor (IMR) is revealed and described. Although formally the reaction rate in this regime depends on the concentration of passive fractions of the reactants, according to the theory presented, the true (hidden) mechanism of the reaction is associated only with the reaction of the active fractions of the reactants with vanishingly small concentration in the volume of the reactor. It is shown that the rate constant of fast chemical reactions can be evaluated when the mixing intensity of reactants is much less than that needed to reach the mixing conditions in an IMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Peak radial growth of diffuse-porous species occurs during periods of lower water availability than for ring-porous and coniferous trees

Climate models project warmer summer temperatures will increase the frequency and heat severity of droughts in temperate forests of Eastern North America. Hotter droughts are increasingly documented to affect tree growth and forest dynamics, with critical impacts on tree mortality, carbon sequestration and timber provision. The growing acknowledgement of the dominant role of drought timing on tree vulnerability to water deficit raises the issue of our limited understanding of radial growth phenology for most temperate tree species. Here, we use well-replicated dendrometer band data sampled frequently during the growing season to assess the growth phenology of 610 trees from 15 temperate species over 6 years. Patterns of diameter growth follow a typical logistic shape, with growth rates reaching a maximum in June, and then decreasing until process termination. On average, we find that diffuse-porous species take 16–18 days less than other wood-structure types to put on 50% of their annual diameter growth. However, their peak growth rate occurs almost a full month later than ring-porous and conifer species (ca. 24 ± 4 days; mean ± 95% credible interval). Unlike other species, the growth phenology of diffuse-porous species in our dataset is highly correlated with their spring foliar phenology. We also find that the later window of growth in diffuse-porous species, coinciding with peak evapotranspiration and lower water availability, exposes them to a higher water deficit of 88 ± 19 mm (mean ± SE) during their peak growth than ring-porous and coniferous species (15 ± 35 mm and 30 ± 30 mm, respectively). Given the high climatic sensitivity of wood formation, our findings highlight the importance of wood porosity as one predictor of species climatic sensitivity to the projected intensification of the drought regime in the coming decades.

59 BASIC BIOLOGICAL SCIENCES↗

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗

Kepler’S supernova: An OVERLUMINOUS type Ia event interacting with a massive circumstellar medium at a very late phase

We have analyzed XMM-Newton, Chandra, and Suzaku observations of three similarly aged, ejecta-dominated supernova remnants (SNRs), Kepler's SNR, Tycho's SNR, and SNR 0509-67.5, to investigate the properties of the SN ejecta and the circumstellar medium (CSM). By simply comparing the X-ray spectra, we find that line intensity ratios of iron-group elements (IGEs) to intermediate-mass elements (IMEs) for Kepler's SNR and SNR 0509-67.5 are much higher than those for Tycho's SNR. Given that Tycho's SNR and SNR 0509-67.5 are thought to be the remnants of a typical SN Ia and an overluminous SN Ia, respectively, we argue that Kepler is the product of an overlumious SN Ia. This inference is supported by our spectral modeling, which reveals the IGE and IME masses, respectively, to be 0.95{sub −0.37}{sup +0.34} M{sub ⊙} and 0.12{sub −0.05}{sup +0.19} M{sub ⊙} (Kepler's SNR), 0.75{sub −0.15}{sup +0.51} M{sub ⊙} and 0.34{sub −0.25}{sup +0.08} M{sub ⊙} (SNR 0509-67.5), and 0.35{sub −0.15}{sup +0.55} M{sub ⊙} and 0.70{sub −0.28}{sup +0.12} M{sub ⊙} (Tycho's SNR). We find that the CSM component in Kepler's SNR consists of tenuous diffuse gas (∼0.3 M{sub ⊙}) present throughout the entire remnant, plus dense knots (∼0.035 M{sub ⊙}). Since both of them show N overabundance, their origin would be CNO-processed material from the progenitor system. The mass of the diffuse CSM allows us to infer the pre-SN mass-loss rate to be ∼1.5×10{sup −5}(v{sub w}/10 km s{sup −1}) M{sub ⊙} yr{sup −1}. The dense knots have slow proper motions and relatively small ionization timescales; hence, they were likely located a few parsecs away from the progenitor at the explosion. We thus argue that Kepler's SN was an overluminous (91T-like) event that recently started to interact with the massive CSM. This supports the possible link between 91T-like SNe and “Ia-CSM” SNe.

79 ASTRONOMY AND ASTROPHYSICS↗

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗

COMPARISON OF IMPLICIT SCHEMES TO SOLVE EQUATIONS OF RADIATION HYDRODYNAMICS WITH A FLUX-LIMITED DIFFUSION APPROXIMATION: NEWTON–RAPHSON, OPERATOR SPLITTING, AND LINEARIZATION

Radiation is an important process of energy transport, a force, and a basis for synthetic observations, so radiation hydrodynamics (RHD) calculations have occupied an important place in astrophysics. However, although the progress in computational technology is remarkable, their high numerical cost is still a persistent problem. In this work, we compare the following schemes used to solve the nonlinear simultaneous equations of an RHD algorithm with the flux-limited diffusion approximation: the Newton–Raphson (NR) method, operator splitting, and linearization (LIN), from the perspective of the computational cost involved. For operator splitting, in addition to the traditional simple operator splitting (SOS) scheme, we examined the scheme developed by Douglas and Rachford (DROS). We solve three test problems (the thermal relaxation mode, the relaxation and the propagation of linear waves, and radiating shock) using these schemes and then compare their dependence on the time step size. As a result, we find the conditions of the time step size necessary for adopting each scheme. The LIN scheme is superior to other schemes if the ratio of radiation pressure to gas pressure is sufficiently low. On the other hand, DROS can be the most efficient scheme if the ratio is high. Although the NR scheme can be adopted independently of the regime, especially in a problem that involves optically thin regions, the convergence tends to be worse. In all cases, SOS is not practical.

79 ASTRONOMY AND ASTROPHYSICS↗

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗

Diffuser Braze Development

At the request of SRTE, through the PDRD project, SRNL is developing a vacuum brazing process to facilitate the fabrication of hydrogen purification systems. In particular, this effort is directed at tube-style hydrogen diffusers. The diffusers are comprised of numerous palladium-silver tubes that are brazed to a nickel fitting, which is subsequently welded to a stainless steel tube, fitting, or tube sheet. Historically, torch brazing has been used as a heat source; for this work vacuum furnace brazing was suggested to improve the uniformity of the process. This report describes the braze process development for braze selection and simple nickel fitting to tube brazing. Five braze alloys were selected and sessile drop testing was conducted over a range of temperatures from 25 to 75° C over the liquidus temperature. Samples were characterized for wetting angle, erosion, and brazability. The cup-to-tube braze joints were inspected using digital radiography, ultrasonic testing, and metallographic examination. Two of the five alloys exhibited excellent fabrication properties and are considered for further development.

36 MATERIALS SCIENCE↗

Cu Migration in Polycrystalline CdTe Solar Cells

An impurity reaction-diffusion model is applied to Cu defects and related intrinsic defects in polycrystalline CdTe for a better understanding of Cu’s role in the cell level reliability of CdTe PV devices. The simulation yields transient Cu distributions in polycrystalline CdTe during solar cell processing and stressing. Preliminary results for Cu migration using available diffusivity and solubility data show that Cu accumulates near the back contact, a phenomena that is commonly observed in devices after back-contact processing or stress conditions.

14 SOLAR ENERGY↗

Experimental and computational design tools for industrial drying processes: Challenges in process‐limit prediction

Abstract The coating and drying of inks and slurries are important steps to manufacture a plethora of products. Drying processes, particularly, comprise energy‐intensive steps that affect product cost and quality. Prior work has highlighted failures of various multicomponent diffusivity models to conserve mass in dryer modeling and challenges in predicting process limits given variability in published values of key thermodynamic parameters. Herein, we develop a computational model and benchtop drying experiments to investigate these concerns for drying polymer‐laden coatings. Model predictions of process limits in a single‐zone drying oven demonstrate that published variability in Flory–Huggins parameter yields large variations in predicted operating temperatures above which blistering occurs. This indicates that caution should be exercised when choosing approaches to obtain or predict the Flory–Huggins parameter, and that both benchtop drying experiments and a set of additional experiments, such as sorption experiments, are needed to fully characterize and optimize a given drying process.

Parrish, Chance↗

Performance Evaluation of Spectroscopic Detectors for LEU Hold-up Measurements

The hold-up measurement of low-enriched uranium materials may require use of alternate detector types relative to the measurement of highly enriched uranium. This is in part due to the difference in process scale (i.e., the components are generally larger for low-enriched uranium systems), but also because the characteristic gamma-ray lines from 235 U used for assay of highly enriched uranium will be present at a much reduced intensity (on a per gram of uranium basis) at lower enrichments. Researchers at Oak Ridge National Laboratory examined the performance of several standard detector types, e.g., NaI(Tl), LaBr3(Ce), and HPGe, to select a suitable candidate for measuring and quantifying low-enriched uranium hold-up in process pipes and equipment at the Portsmouth gaseous diffusion plant. Detector characteristics, such as energy resolution (full width at half maximum) and net peak count rates at gamma ray energies spanning a range of 60–1332 keV, were measured for the above-mentioned detector types using the same sources and in the same geometry. Uranium enrichment standards (Certified Reference Material no. 969 and Certified Reference Material no. 146) were measured using each of the detector candidates in the same geometry. The net count rates recorded by each detector at 186 keV and 1,001 keV were plotted as a function of enrichment (atom percentage). Background measurements were made in unshielded and shielded configurations under both ambient and elevated conditions of 238 U activity. The highly enriched uranium hold-up measurement campaign at the Portsmouth plant was performed on process equipment that had been cleaned out. Therefore, in most cases, the thickness of the uranium deposits was less than the “infinite thickness” for the 186 keV gamma rays to be completely self-attenuated. Because of this, in addition to measuring the 186 keV gamma, the 1,001 keV gamma ray from 234m Pa—a daughter of 238 U in secular equilibrium with its parent—will also need to be measured. Based on the performance criteria of detection efficiency, energy resolution, peak-to-continuum ratios, minimum detectable limits, and the weight of the shielded probe, a shielded (0.5 in. thick lead shield) 2 × 2 in. NaI(Tl) detector is recommended for use. The recommended approach is to carry out analysis using data from both 186 keV and 1,001 keV gamma rays, and select a best result based on propagated uncertainty estimates. It is also highly recommended that a two-point gain stabilization scheme based on an 241 Am seed embedded in the probe be implemented. Shielding configurations to reduce the impact of background interference on the measurement of 1,001 keV gamma-ray are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Stochastic modeling of catalytic processes in nanoporous materials: Beyond mean-field approach

Transport and reaction in zeolites and other porous materials, such as mesoporous silica particles, has been a focus of interest in recent years. This is in part due to the possibility of anomalous transport effects (e.g. single-file diffusion) and its impact in the reaction yield in catalytic processes. Computational simulations are often used to study these complex nonequilibrium systems. Computer simulations using Molecular Dynamics (MD) techniques are prohibitive, so instead coarse grained one-dimensional models with the aid of Kinetic Monte Carlo (KMC) simulations are used. Both techniques can be computationally expensive, both time and resource wise. These coarse-grained systems can be exactly described by a set of coupled stochastic master equations, that describe the reaction-diffusion kinetics of the system. The equations can be written exactly, however, coupling between the equations and terms within the equations make it impossible to solve them exactly; approximations must be made. One of the most common methods to obtain approximate solutions is to use Mean Field (MF) theory. MF treatments yield reasonable results at high ratios of reaction rate k to hop rate h of the particles, but fail completely at low k=h due to the over-estimation of fluxes of particles within the pore. We develop a method to estimate fluxes and intrapore diffusivity in simple one- dimensional reaction-diffusion models at high and low k=h, where the pores are coupled to an equilibrated three-dimensional fluid. We thus successfully describe analytically these simple reaction-diffusion one-dimensional systems. Extensions to models considering behavior with long range steric interactions and wider pores require determination of multiple boundary conditions. We give a prescription to estimate the required parameters for these simulations. For one dimensional systems, if single-file diffusion is relaxed, additional parameters to describe particle exchange have to be introduced. We use Langevin Molecular Dynamics (MD) simulations to assess these parameters.

36 MATERIALS SCIENCE↗

Nanoscale insight into the p-n junction of alkali-incorporated Cu(In,Ga)Se 2 solar cells

The effects of alkali diffusion and post-deposition treatment in three-stage processed Cu(In,Ga)Se 2 solar cells are examined by using atom probe tomography and electrical property measurements. Cells, for which the substrate was treated at 650 °C to induce alkali diffusion from the substrate prior to absorber deposition, exhibited high open-circuit voltage (758 mV) and efficiency (18.2%) and also exhibited a 50 to 100-nm-thick ordered vacancy compound layer at the metallurgical junction. Surprisingly, these high-temperature samples exhibited higher concentrations of K at the junction (1.8 at.%) than post-deposition treatment samples (0.4 at.%). A model that uses Ga/(Ga + In) and Cu/(Ga + In) profiles to predict bandgaps (+/-17.9 meV) of 22 Cu(In,Ga)Se2 solar cells reported in literature was discussed and ultimately used to predict band properties at the nanoscale by using atom probe tomography data. The high-temperature samples exhibited a greater drop in the valence band maximum (200 meV) due to a lower Cu/(Ga + In) ratio than the post-deposition treatment samples. There was an anticorrelation of K concentrations and Cu/(Ga + In) ratios for all samples, regardless of processing conditions. In conclusion, changes in elemental profiles at the active junctions correlate well with the electrical behaviour of these devices.

14 SOLAR ENERGY↗

Spatiotemporal Reaction Dynamics Control in Two‐Photon Polymerization for Enhancing Writing Characteristics

Since 2001, 3D microfabrication based on two‐photon polymerization (TPP) has drawn extensive attention and interest in biology, optics, photonics, material science, and high‐energy physics. The in‐volume fabrication capability due to the threshold behavior of two‐photon absorption enables TPP higher flexibility compared with other nanofabrication techniques. However, as determined by the in‐volume fabrication feature as well as various reaction dynamics, the writing characteristics of TPP, such as throughput, accuracy, surface quality, and fabrication capability, are still limited. Herein, a comprehensive study is performed on the spatiotemporal behavior of reaction dynamics during TPP fabrication, mainly focusing on spatiotemporal characteristics of radical diffusion, photothermal effect, microscale mechanics, and voxel stacking process. Based on the study, a nonsequential fabrication method is established to simultaneously improve key writing characteristics of TPP and realize sharp features, high speeds, large overhang structure, and smooth surfaces. The method established in this work can be applied to improve the performance of functional devices for various fields.

36 MATERIALS SCIENCE↗

Curvature effects in the reaction-diffusion system governing frontal polymerization

Frontal polymerization is a fast and energy-efficient curing process of thermoset polymers involving the self-sustained propagation of a sharp front that transforms a liquid monomer or gel into a polymer. The temperature, speed, and geometry of the polymerization front are closely intertwined quantities that can affect the curing process and the mechanical properties of the polymer. This work investigates the effect of the front curvature on the front propagation speed based on a reaction-diffusion model. The proposed analytical relations between front speed and curvature are verified by finite element simulations and validated through experiments. The models show that a concave front has a higher front temperature and speed than its convex counterpart. When the curvature of a convex front reaches a critical value, the reaction can no longer be sustained, and the front is quenched. The relation between the speed of a steadily propagating front and its geometry is derived, showing good agreement with experimental observations. This study also provides analytical estimates of the thickness of the polymerization front and insight into the relationship between the front curvature and the kinetics of the exothermic reaction.

42 ENGINEERING↗

A phase field study of the thermal migration of gas bubbles in UO 2 nuclear fuel under temperature gradient

Phase field models are developed to study the gas bubble migration in uranium dioxide nuclear fuel in which a large temperature gradient exists during the operation. In this work, thermal diffusion mechanism for nanosized gas bubbles and vapor transport process for micron-sized gas bubbles are considered, respectively. In both cases, gas bubbles migrate to the high-temperature area. Due to the velocity difference between leading and trailing edges of the gas bubbles, nanosized gas bubbles are elongated along the temperature gradient direction when thermal diffusion is dominated. Micron-sized gas bubbles are either compressed along temperature gradient direction to form lenticular shape bubbles or elongated along temperature gradient direction, depending on the location of the gas bubbles within the fuel pellet. Initial gas bubble radius has no significant effect on the gas bubble migration velocity for both thermal diffusion and vapor transport mechanisms. We notice that the shape change of the gas bubble due to vapor transport mechanism has no significant effect on the migration velocity. Furthermore, the center cavity formation is also captured by our model which is due to the migration and accumulation of lenticular gas bubbles at the center of the fuel pellet. The modeling results compare well with experimental observations and theoretical analysis in the literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Growth mechanism and magnetic and electrochemical properties of Na{sub 0.44}MnO{sub 2} nanorods as cathode material for Na-ion batteries

Nanorods of Na{sub 0.44}MnO{sub 2} are a promising cathode material for Na-ion batteries due to their large surface area and single crystalline structure. We report the growth mechanism of Na{sub 0.44}MnO{sub 2} nanorods via solid state synthesis and their physical properties. The structure and the morphology of the Na{sub 0.44}MnO{sub 2} nanorods are investigated by X-ray diffraction (XRD), scanning and tunneling electron microscopy (SEM and TEM), and energy-dispersive X-ray (EDX) techniques. The growth mechanism of the rods is investigated and the effects of vapor pressure and partial melting of Na-rich regions are discussed. The magnetic measurements show an antiferromagnetic phase transition at 25 K and the μ{sub eff} is determined as 3.41 and 3.24 μ{sub B} from the χ–T curve and theoretical calculation, respectively. The electronic configuration and spin state of Mn{sup 3+} and Mn{sup 4+} are discussed in detail. The electrochemical properties of the cell fabricated using the nanorods are investigated and the peaks in the voltammogram are attributed to the diffusion of Na ions from different sites. Na intercalation process is explained by one and two Margules and van Laar models. - Highlights: • We synthesized Na{sub 0.44}MnO{sub 2} nanorods via a simple solid state reaction technique. • Our studies show that excess Na plays a crucial role in the nanorod formation. • Magnetization measurements show that Mn{sup 3+} ions are in LS and HS states. • The electrochemical properties of the cell fabricated using the nanorods are investigated. • Na intercalation process is explained by one and two Margules and van Laar models.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Entropic and gradient flow formulations for nonlinear diffusion

Nonlinear diffusion ∂{sub t}ρ = Δ(Φ(ρ)) is considered for a class of nonlinearities Φ. It is shown that for suitable choices of Φ, an associated Lyapunov functional can be interpreted as thermodynamic entropy. This information is used to derive an associated metric, here called thermodynamic metric. The analysis is confined to nonlinear diffusion obtainable as hydrodynamic limit of a zero range process. The thermodynamic setting is linked to a large deviation principle for the underlying zero range process and the corresponding equation of fluctuating hydrodynamics. For the latter connections, the thermodynamic metric plays a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗