Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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,585 records · Page 88

WaterTransport in PEM Fuel Cells: Advanced Modeling, Material Selection, Testing and Design Optimization

Water management in Proton Exchange Membrane, PEM, Fuel Cells is challenging because of the inherent conflicts between the requirements for efficient low and high power operation. Particularly at low powers, adequate water must be supplied to sufficiently humidify the membrane or protons will not move through it adequately and resistance losses will decrease the cell efficiency. At high power density operation, more water is produced at the cathode than is necessary for membrane hydration. This excess water must be removed effectively or it will accumulate in the Gas Diffusion Layers, GDLs, between the gas channels and catalysts, blocking diffusion paths for reactants to reach the catalysts and potentially flooding the electrode. As power density of the cells is increased, the challenges arising from water management are expected to become more difficult to overcome simply due to the increased rate of liquid water generation relative to fuel cell volume. Thus, effectively addressing water management based issues is a key challenge in successful application of PEMFC systems. In this project, CFDRC and our partners used a combination of experimental characterization, controlled experimental studies of important processes governing how water moves through the fuel cell materials, and detailed models and simulations to improve understanding of water management in operating hydrogen PEM fuel cells. The characterization studies provided key data that is used as inputs to all state-of-the-art models for commercially important GDL materials. Experimental studies and microscopic scale models of how water moves through the GDLs showed that the water follows preferential paths, not branching like a river, as it moves toward the surface of the material. Experimental studies and detailed models of water and airflow in fuel cells channels demonstrated that such models can be used as an effective design tool to reduce operating pressure drop in the channels and the associated costs and weight of blowers and pumps to force air and hydrogen gas through the fuel cell. Promising improvements to materials structure and surface treatments that can potentially aid in managing the distribution and removal of liquid water were developed; and improved steady-state and freeze-thaw performance was demonstrated for a fuel cell stack under the self-humidified operating conditions that are promising for stationary power generation with reduced operating costs.

30 DIRECT ENERGY CONVERSION↗

Sandia's Current Energy Conversion module for Flexible-Mesh Delft3D flow solver v. 1.0

The DOE has funded Sandia National Labs (SNL) to develop an open-source modeling tool to guide the design and layout of marine hydrokinetic (MHK) arrays to maximize power production while minimizing environmental effects. This modeling framework simulates flows through and around a MHK arrays while quantifying environmental responses. As an augmented version of the Dutch company, Deltares’s, environmental hydrodynamics code, Delft3D, SNL-Delft3D-CEC-FM includes a new module that simulates energy conversion (momentum withdrawal) by MHK current energy conversion devices with commensurate changes in the turbulent kinetic energy and its dissipation rate. SNL-Delft3D-CEC-FM modified the Delft3D flexible mesh flow solver, DFlowFM. The Following is a description of Deltares’s open-source code Delft3D from which Delft3D-CEC is built upon. “Delft3D is a world leading 3D modeling suite to investigate hydrodynamics, sediment transport and morphology and water quality for fluvial, estuarine and coastal environments. As per 1 January 2011, the Delft3D flow (FLOW), morphology (MOR) and waves (WAVE) modules are available in open source. The software is used and has proven his capabilities on many places around the world, like the Netherlands, USA, Hong Kong, Singapore, Australia, Venice, etc. The software is continuously improved and developed with innovating advanced modelling techniques as consequence of the research work of our institute and to stay world leading. The FLOW module is the heart of Delft3D and is a multi-dimensional (2D or 3D) hydrodynamic (and transport) simulation program which calculates non-steady flow and transport phenomena resulting from tidal and meteorological forcing on a curvilinear, boundary fitted grid or spherical coordinates. In 3D simulations, the vertical grid is defined following the so-called sigma coordinate approach or Z-layer approach. The MOR module computes sediment transport (both suspended and bed total load) and morphological changes for an arbitrary number of cohesive and non-cohesive fractions. Both currents and waves act as driving forces and a wide variety of transport formulae have been incorporated. For the suspended load this module connects to the 2D or 3D advection-diffusion solver of the FLOW module; density effects may be taken into account. An essential feature of the MOR module is the dynamic feedback with the FLOW and WAVE modules, which allow the flows and waves to adjust themselves to the local bathymetry and allows for simulations on any time scale from days (storm impact) to centuries (system dynamics). It can keep track of the bed composition to build up a stratigraphic record. The MOR module may be extended to include extensive features to simulate dredging and dumping scenarios. For over 30 years Deltares has been in the forefront of these types of combined morphological simulation techniques.”

Chartand, Chris↗

Investigation of the effect of local atomic polarization on field-enhanced ion transport in insulating binary oxides: The case of CeO 2

Over the course of more than 70 years, several theoretical models for electric field-enhanced ion diffusion in crystalline insulators have been presented. However, there is no assessment of the validity of these models despite the urge to have a correct description of field-enhanced ion diffusion in several emerging technologies. Herein, an assessment of five models was carried out by computing the field-dependent migration enthalpy ΔH mig (E) for the oxide ion in CeO 2 . The input to these models is the zero-field and zero-temperature migration pathway and/or the activation barrier obtained from an interatomic potential. ΔH mig (E) from these models was compared with reference values obtained from a set of classical molecular-dynamics simulations in the temperature range of 1000 ≤ T ≤ 1600 K and the electric field range of 0 ≤ E ≤ 30 MV/cm. It is revealed that the most successful theoretical models are those that consider local polarization effects induced in the diffusing ion itself and its vicinity. However, all the models did not account for the effect of the electric field and finite temperature on the local polarization, leading to discrepancies between the predicted and the reference ΔH mig (E), particularly at high fields. By comparing results from a rigid-ion potential and a polarizable core-shell potential, it is concluded that an intraionic polarization degree of freedom in the polarizable potential is an important factor in predicting ΔH mig (E), regardless of the field-enhanced diffusion model particularly for low to moderate field. Finally, future more accurate treatments should consider the field- and temperature-dependent interionic and intraionic local polarization effects.

36 MATERIALS SCIENCE↗

Temperature-dependent resistance switching in SrTiO{sub 3}

Resistance switching phenomena were studied by varying temperature in SrTiO{sub 3} single crystal. The resistance hysteresis loops appear at a certain temperature ranging from 340 K to 520 K. With the assistance of 375 nm ultraviolet continuous laser, the sample resistance is greatly reduced, leading to a stable effect than that in dark. These resistance switching phenomena only exist in samples with enough oxygen vacancies, which is confirmed by spherical aberration-corrected scanning transmission electron microscopy measurements, demonstrating an important role played by oxygen vacancies. At temperatures above 340 K, positively charged oxygen vacancies become mobile triggered by external electric field, and the resistance switching effect emerges. Our theoretical results based on drift-diffusion model reveal that the built-in field caused by oxygen vacancies can be altered under external electric field. Therefore, two resistance states are produced under the cooperative effect of built-in field and external field. However, the increasing mobility of oxygen vacancies caused by higher temperature promotes internal electric field to reach equilibrium states quickly, and suppresses the hysteresis loops above 420 K.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Communication: Modeling of concentration dependent water diffusivity in ionic solutions: Role of intermolecular charge transfer

The translational diffusivity of water in solutions of alkali halide salts depends on the identity of ions, exhibiting dramatically different behavior even in solutions of similar salts of NaCl and KCl. The water diffusion coefficient decreases as the salt concentration increases in NaCl. Yet, in KCl solution, it slightly increases and remains above bulk value as salt concentration increases. Previous classical molecular dynamics simulations have failed to describe this important behavior even when polarizable models were used. Here, we show that inclusion of dynamical charge transfer among water molecules produces results in a quantitative agreement with experiments. Our results indicate that the concentration-dependent diffusivity reflects the importance of many-body effects among the water molecules in aqueous ionic solutions. Comparison with quantum mechanical calculations shows that a heterogeneous and extended distribution of charges on water molecules around the ions due to ion-water and also water-water charge transfer plays a very important role in controlling water diffusivity. Explicit inclusion of the charge transfer allows us to model accurately the difference in the concentration-dependent water diffusivity between Na{sup +} and K{sup +} ions in simulations, and it is likely to impact modeling of a wide range of systems for medical and technological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of the multigrid amplitude function method for time-dependent transport equation using MOC

An efficient numerical method for time-dependent transport equation, the mutigrid amplitude function (MAF) method, is proposed. The method of characteristics (MOC) is being widely used for reactor analysis thanks to the advances of numerical algorithms and computer hardware. However, efficient kinetic calculation method for MOC is still desirable since it requires significant computation time. Various efficient numerical methods for solving the space-dependent kinetic equation, e.g., the improved quasi-static (IQS) and the frequency transform methods, have been developed so far mainly for diffusion calculation. These calculation methods are known as effective numerical methods and they offer a way for faster computation. However, they have not been applied to the kinetic calculation method using MOC as the authors' knowledge. Thus, the MAF method is applied to the kinetic calculation using MOC aiming to reduce computation time. The MAF method is a unified numerical framework of conventional kinetic calculation methods, e.g., the IQS, the frequency transform, and the theta methods. Although the MAF method is originally developed for the space-dependent kinetic calculation based on the diffusion theory, it is extended to transport theory in the present study. The accuracy and computational time are evaluated though the TWIGL benchmark problem. The calculation results show the effectiveness of the MAF method. (authors)

97 MATHEMATICS AND COMPUTING↗

Effects of dry bulk density and particle size fraction on gas transport parameters in variably saturated landfill cover soil

Highlights: > The effects of soil physical properties on gas transport parameters were investigated. > Higher values of D{sub p} and k{sub a} exhibited in the '+gravel' than the '-gravel' fraction at same soil-air content ({epsilon}). > Recent power law models for D{sub p} (WLR) and k{sub a} (RPL) were modified. > Model parameters were linearly related to easily measurable dry bulk density ({rho}{sub b}). - Abstract: Landfill sites are emerging in climate change scenarios as a significant source of greenhouse gases. The compacted final soil cover at landfill sites plays a vital role for the emission, fate and transport of landfill gases. This study investigated the effects of dry bulk density, {rho}{sub b}, and particle size fraction on the main soil-gas transport parameters - soil-gas diffusivity (D{sub p}/D{sub o}, ratio of gas diffusion coefficients in soil and free air) and air permeability (k{sub a}) - under variably-saturated moisture conditions. Soil samples were prepared by three different compaction methods (Standard and Modified Proctor compaction, and hand compaction) with resulting {rho}{sub b} values ranging from 1.40 to 2.10 g cm{sup -3}. Results showed that D{sub p} and k{sub a} values for the '+gravel' fraction (<35 mm) became larger than for the '-gravel' fraction (<2 mm) under variably-saturated conditions for a given soil-air content ({epsilon}), likely due to enhanced gas diffusion and advection through less tortuous, large-pore networks. The effect of dry bulk density on D{sub p} and k{sub a} was most pronounced for the '+gravel' fraction. Normalized ratios were introduced for all soil-gas parameters: (i) for gas diffusivity D{sub p}/D{sub f}, the ratio of measured D{sub p} to D{sub p} in total porosity (f), (ii) for air permeability k{sub a}/k{sub a,pF4.1}, the ratio of measured k{sub a} to k{sub a} at 1235 kPa matric potential (=pF 4.1), and (iii) for soil-air content, the ratio of soil-air content ({epsilon}) to total porosity (f) (air saturation). Based on the normalized parameters, predictive power-law models for D{sub p}({epsilon}/f) and k{sub a}({epsilon}/f) models were developed based on a single parameter (water blockage factor M for D{sub p} and P for k{sub a}). The water blockage factors, M and P, were found to be linearly correlated to {rho}{sub b} values, and the effects of dry bulk density on D{sub p} and k{sub a} for both '+gravel' and '-gravel' fractions were well accounted for by the new models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Long lifetime of thermally excited magnons in bulk yttrium iron garnet

Spin currents are generated within the bulk of magnetic materials due to heat flow, an effect called intrinsic spin Seebeck. This bulk bosonic spin current consists of a diffusing thermal magnon cloud, parametrized by the magnon chemical potential (μ m ), with a diffusion length of several microns in yttrium iron garnet (YIG). Transient optothermal measurements of the spin-Seebeck effect (SSE) as a function of temperature reveal the time evolution of μ m due to intrinsic SSE in YIG. The interface SSE develops at times <2 ns while the intrinsic SSE signal continues to evolve at times >500 μs, dominating the temperature dependence of SSE in bulk YIG. Time-dependent SSE data are fit to a multitemperature model of coupled spin/heat transport using the finite-element method (FEM), where the magnon spin lifetime (τ) and magnon-phonon thermalization time (τ mp ) are used as fit parameters. From 300 to 4 K, τ mp ) varies from 1 to 10 ns, whereas tau varies from 2 to 60 μ s with the spin lifetime peaking at 90 K. At low temperature, a reduction in τ is observed consistent with impurity relaxation reported in ferromagnetic resonance measurements. Finally, these results demonstrate that the thermal magnon cloud in YIG contains extremely low-frequency magnons (similar to 10 GHz), providing spectral insight to the microscopic scattering processes involved in magnon spin/heat diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

THE CIRCUMGALACTIC MEDIUM OF MASSIVE GALAXIES AT z {approx} 3: A TEST FOR STELLAR FEEDBACK, GALACTIC OUTFLOWS, AND COLD STREAMS

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift {approx}3, using the Eris suite of cosmological hydrodynamic ''zoom-in'' simulations. The reference run adopts a blastwave scheme for supernova feedback that produces large-scale galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. The effect of the local UV radiation field is added in post-processing. The CGM (defined as all gas at R > 0.2 R{sub vir} = 10 kpc, where R{sub vir} is the virial radius) contains multiple phases having a wide range of physical conditions, with more than half of its heavy elements locked in a warm-hot component at T > 10{sup 5} K. Synthetic spectra, generated by drawing sightlines through the CGM, produce interstellar absorption-line strengths of Ly{alpha}, C II, C IV, Si II, and Si IV as a function of the galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high redshift by Steidel et al. The covering factor of absorbing material declines less rapidly with impact parameter for Ly{alpha} and C IV compared to C II, Si IV, and Si II, with Ly{alpha} remaining strong (W{sub Ly{alpha}} > 300 mA) to {approx}> 5 R{sub vir} = 250 kpc. Only about one third of all the gas within R{sub vir} is outflowing. The fraction of sightlines within one virial radius that intercept optically thick, N{sub H{sub I}}>10{sup 17.2} cm{sup -2} material is 27%, in agreement with recent observations by Rudie et al. Such optically thick absorption is shown to trace inflowing ''cold'' streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar by previous episodes of star formation in the main host and in nearby dwarfs, are the origin of strong (N{sub C{sub II}}>10{sup 13} cm{sup -2}) C II absorption with a covering factor of 22% within R{sub vir} and 10% within 2 R{sub vir}. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large O VI halo, with a typical column density N{sub O{sub VI}} {approx}> 10{sup 14} cm{sup -2} and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. Our zoom-in simulations of this single system appear then to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗

On the study of threshold intensity dependence on the gain and loss processes in laser induced spark ignition of molecular hydrogen

In the present work, a numerical analysis is performed to investigate the comparative contribution of the mechanisms responsible for electron gain and losses in laser spark ignition and plasma formation of H{sub 2}. The analysis considered H{sub 2} over pressure range 150 -3000 torr irradiated by a Nd:YAG laser radiation at wavelengths 1064 and 532 nm with pulse length 5.5 ns. The study based on a modified electron cascade model by one of the authors which solves numerically the time dependent Boltzmann equation as well as a set of rate equations that describe the rate of change of the excited states population. The model includes most of the physical processes that might take place during the interaction. Computations of The threshold intensity are performed for the combined and separate contribution of each of the gain and loss processes. Reasonable agreement with the measured values over the tested pressure range is obtained only for the case of the combined contribution. Basing on the calculation of the electron energy distribution function, the determined relations of the time evolution of the electrons density for selected values of the tested gas pressure region revealed that photo-ionization of the excited states could determine the time of electron generation and hence spark ignition. Collisional ionization contributes to this phenomenon only at the high pressure regime. Loss processes due to electron diffusion, vibrational excitation are found to have significant effect over examined pressure values for the two applied laser wavelengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A weighted state redistribution algorithm for embedded boundary grids

State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions. Second, we discuss several algorithmic changes and improvements motivated by the more complicated cut cell geometries that can occur in higher dimensions. In particular, we introduce a weighted version with less dissipation in an easily generalizable framework. Third, we demonstrate that state redistribution can also stabilize a solution update that includes both advective and diffusive contributions. Notably, the stabilization algorithm is shown to be effective for incompressible as well as compressible reacting flows. Finally, we discuss the implementation of the algorithm for several exascale-ready simulation codes based on AMReX, demonstrating ease of use in combination with domain decomposition, hybrid parallelism and complex physics.

97 MATHEMATICS AND COMPUTING↗

Effects of temperature and dose rate on ion-irradiated γ-LiAlO 2 pellets

Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by competing processes of defect production, driven by the dose rate, and defect recovery, controlled by diffusion, interaction, and annihilation. Here, this study investigates the effects of irradiation temperature and the dose rate on microstructural evolution, deuterium retention, and lithium volatilization in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation. Experiments were performed to a total fluence of 3 × 10 17 (He + + D + )/cm 2 at 623, 673, 723, and 773 K with an average He + dose rate of 7.7 × 10 −4 dpa/s, and to 2 × 10 17 (He + + D + )/cm 2 at 773 K with dose rates of 6.8 × 10 −5 , 2.9 × 10 −4 , and 7.3 × 10 −4 dpa/s. At 623 K, the microstructure was dominated by cavities and fractures with no observable precipitate formation, while small precipitates emerged at 673 K. Increasing the irradiation temperature to 723–773 K promoted the formation of larger, faceted LiAl 5 O 8 precipitates, and surface amorphization, accompanied by pronounced lithium depletion and H–D isotopic exchange. At 773 K, medium and high dose rates produced an amorphized surface layer over a crystalline subsurface containing LiAl 5 O 8 precipitates and blisters at the crystalline–amorphous interface, whereas low-dose-rate irradiation preserved surface crystallinity with cavities distributed in the matrix, around precipitates, and along grain boundaries. Precipitate morphology was anisotropic with limited size dependence on the dose rate. These results elucidate the coupled effects of temperature and the dose rate and demonstrate that sequential He + and D 2 + irradiation at 773 K reproduces key microstructural features and H isotope behavior observed in neutron-irradiated γ-LiAlO 2 at 573 K.

dose rate effects↗

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

Coupled social and infrastructure approaches for enhancing solar energy adoption. Final Report

The goal of this project was to work with rural electric cooperatives to facilitate the diffusion of solar energy adoption in households located in the rural and semi-urban areas of Virginia by identifying social and behavioral factors that might be unique to rural regions; and develop a model to calculate the solar adoption propensity score for household based on their demographics, social and behavioral characteristics which would provide an objective metric to cooperatives that can be further used to do targeted marketing of rooftop solar panels. This was achieved through the following tasks: (1) Conducted a survey of the members of Virginia electric cooperatives to identify demographic, social, financial and behavioral attributes of individuals who are likely to adopt rooftop solar panels. (2) Developed a highly detailed, data-driven, agent-based model of the population of Virginia, focusing on the rural regions. (3) Developed diffusion models that use social, behavioral, and demographic factors, and peer effects to study their impact on solar adoption in rural areas. (4) Built a prototype tool based on the diffusion model to help study market segmentation in rural areas and made it available to National Rural Electric Cooperative Association (NRECA). (5) Results and recommendations derived from the model were provided to NRECA to be shared with participating cooperatives. (6) Results were published in peer reviewed journals, conference proceedings and book chapters, and ideas disseminated through presentations and newsletters. There were several important methodological contributions made under this project which are detailed in the published papers, including: (1) Built a decision-adjusted model for predicting adoptors with imbalanced training data; (2) designed seeding strategies to maximize adoption given a fixed budget; (3) built a methodology to compare different agent based models; (4) created models to identify important factors that influence decision to adopt solar panels; and (5) built a methodology for building household profiles of solar generation to study the duck curve phenomenon. The team included members from the University of Virginia (lead), National Rural Electric Cooperative Association (NRECA), Arizona State University, Virginia Tech and Sandia National Laboratory. Note that no individual entity or stakeholder has incentive to promote solar in rural regions. Most of the research and work focuses around urban regions where the potential for growth in solar adoption is higher due to higher population density. This puts rural areas at a disadvantage. By improving the diffusion of solar adoption in rural parts of the country, we can not only provide clean energy to rural areas but also promote job growth and improves energy independence.

14 SOLAR ENERGY↗

Multiphonon diffuse scattering in solids from first principles: Application to layered crystals and two-dimensional materials

Time-resolved diffuse scattering experiments have gained increasing attention due to their potential to reveal nonequilibrium dynamics of crystal lattice vibrations with full momentum resolution. Although progress has been made in interpreting experimental data on the basis of one-phonon scattering, understanding the role of individual phonons can be sometimes hindered by multiphonon excitations. In Ref. [M. Zacharias, H. Seiler, F. Caruso, D. Zahn, F. Giustino, P. C. Kelires, and R. Ernstorfer, Phys. Rev. Lett. 127, 207401 (2021)], we have introduced a rigorous approach for the calculation of the all-phonon inelastic scattering intensity of solids from first-principles. In the present work, we describe our implementation in detail and show that multiphonon interactions are captured efficiently by exploiting translational and time-reversal symmetries of the crystal. We demonstrate its predictive power by calculating the scattering patterns of monolayer molybdenum disulfide (MoS 2 ), bulk MoS 2 , and black phosphorus (bP), and we obtain excellent agreement with our measurements of thermal electron diffuse scattering. Remarkably, our results show that multiphonon excitations dominate in bP across multiple Brillouin zones, while in MoS 2 they play a less pronounced role. We expand our analysis for each system and examine the effect of individual atomic and interatomic vibrational motion on the diffuse scattering signals. We further demonstrate the high-throughput capability of our approach by reporting all-phonon scattering maps of two-dimensional MoSe 2 , WSe 2 , WS 2 , graphene, and CdI 2 , rationalizing in each case the effect of multiphonon processes. As a side point, we show that the special displacement method reproduces the thermally distorted configuration that generates precisely the all-phonon diffuse pattern. In conclusion, the present methodology opens the way for systematic calculations of the scattering intensity in crystals and the accurate interpretation of static and time-resolved inelastic scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Assessment of the Idaho National Laboratory Remote-Handled Low-Level Waste Disposal Facility Vault Concrete Data

This document provides a summary of the concrete performance data collected during vault fabrication and installation. The report includes the process by which quality data was collected during vault fabrication and installation, the inspection requirements, and a summary of resulting test data with the potential to impact vault system durability. This report also includes quantitative data collected on cured concrete samples poured for random concrete batches throughout the concrete component fabrication process. Quantitative data includes total porosity, bulk density, absorption, effective porosity, and gas-phase permeability. Chloride diffusivity data collected during the concrete mix selection process are also included for completeness. Data shows that the vault quality assurance program resulted in components received insignificant defects and damaged areas during the fabrication stage and during the vault installation stage. Data collected on the cured concrete samples show the gas-phase permeability and effective porosity are slightly higher than the data used to select the concrete mix designs.

99 GENERAL AND MISCELLANEOUS↗

Intermixing in Cu/Ni multilayers induced by cold rolling

Repeated cold rolling was performed on multilayers of Cu60/Ni40 and Cu40/Ni60 foil arrays to study the details of driven atomic scale interfacial mixing. With increasing deformation, there is a significant layer refinement down to the nm level that leads to the formation of a solid solution phase from the elemental end members. Intriguingly, the composition of the solid solution is revealed by an oscillation in the composition profile across the multilayers, which is different from the smoothly varying profile due to thermally activated diffusion. During the reaction, Cu mixed into Ni preferentially compared to Ni mixing into Cu, which is also in contrast to the thermal diffusion behavior. This is confirmed by observations from X-ray diffraction, electron energy loss spectrum and atom probe tomography. The diffusion coefficient induced by cold rolling is estimated as 1.7 × 10{sup −17} m{sup 2}/s, which cannot be attributed to any thermal effect. The effective temperature due to the deformation induced mixing is estimated as 1093 K and an intrinsic diffusivity d{sub b}, which quantifies the tendency towards equilibrium in the absence of thermal diffusion, is estimated as 6.38 × 10{sup −18} m{sup 2}/s. The fraction of the solid solution phase formed is illustrated by examining the layer thickness distribution and is described by using an error function representation. The evolution of mixing in the solid solution phase is described by a simplified sinusoid model, in which the amplitude decays with increased deformation level. The promoted diffusion coefficient could be related to the effective temperature concept, but the establishment of an oscillation in the composition profile is a characteristic behavior that develops due to deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗