Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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,045 records · Page 58

Understanding coarsening of a post-corrosion microstructure in a molten salt by combining phase-field modeling and in situ tomography

Alloys corroding in molten salt have been observed to form bicontinuous, nanoporous microstructures via dealloying, which subsequently undergo coarsening due to facile transport in high-temperature conditions. In this work, we describe a methodology to elucidate the underlying transport mechanisms during coarsening of a bicontinuous microstructure via quantitative comparisons between phase-field simulations and four-dimensional in situ experiments, in this case X-ray nanotomography of the coarsening of a dealloyed 80 wt% Ni-20 wt% Cr microwire in molten KCl-MgCl 2 at 800°C. We conduct phase-field simulations initialized from experimental data to model coarsening via three different transport mechanisms: surface diffusion, solid bulk diffusion, and liquid bulk diffusion. These simulations reproduce key features of the experiment, such as the densification of the outer layer of the dealloyed wire and the reduction in radius over time. We quantitatively compare different microstructural characteristics between the simulations and experiment and extract temporal scaling factors that optimally match the time scales of the simulations to that of the experiment. This allows us to evaluate morphological similarity between the simulations and experiment and relate the experimental coarsening kinetics to fundamental material properties. We find that surface diffusion is most likely to be the dominant coarsening mechanism, and its kinetics imply a surface diffusivity of D S = 8.9 x 10 -20 m 3 /s, which is within the range of reported values for Ni-vacuum interfaces at 800°C. However, the difference between the experiment and the surface diffusion simulation increases substantially at late times, suggesting that other mechanisms, such as the dissolution of residual Cr, may be at play.

36 MATERIALS SCIENCE↗

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic and Mechanistic Evaluation of Inorganic Arsenic Species Adsorption onto Humic Acid Grafted Magnetite Nanoparticles

In this study, humic acid coated magnetic iron oxide nanoparticles (HA-MNP) were synthesized, characterized, and studied for the removal of toxic inorganic arsenic species from aqueous media. The adsorption of As(III) and As(V) followed pseudo-second order kinetics and the observed data were accurately modeled employing Freundlich adsorption isotherm. Application of the Weber and Morris intraparticle diffusion model to the observed kinetic data suggest that the adsorption occurs in three distinct stages, fast, intermediate, and slow steps. We propose the initial stage is governed by surface association, followed by intraparticle diffusion of arsenic through the HA matrix and finally, chemical reaction or bonding between the arsenic species and HA functionality. The HAMNP nano-adsorbent absorbs > 95% of the inorganic arsenic species with an absorption capacity of 12.2-12.6 mg/g from aqueous media and is effective under a variety of conditions. Post arsenic adsorption characterization of the nanoparticles suggests that As(III) binds with the carboxylate group of HA through a proposed ester type linkage while electrophilic reactions can occur between the nucleophilic functional groups present in HA and the electrophilic arsenic atom in As(V). Finally, the results obtained demonstrated that HA-MNP are robust and have promise for effective As (III) and As(V) remediation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shifts in pore connectivity from precipitation versus groundwater rewetting increases soil carbon loss after drought

Droughts and other extreme precipitation events are predicted to increase in intensity, duration and extent, with uncertain implications for terrestrial carbon (C) sequestration. Soil wetting from above (precipitation) results in a characteristically different pattern of pore-filling than wetting from below (groundwater), with larger, well-connected pores filling before finer pore spaces, unlike groundwater rise in which capillary forces saturate the finest pores first. Here we demonstrate that pore-scale wetting patterns interact with antecedent soil moisture conditions to alter pore-, core- and field-scale C dynamics. Drought legacy and wetting direction are perhaps more important determinants of short-term C mineralization than current soil moisture content in these soils. Our results highlight that microbial access to C is not solely limited by physical protection, but also by drought or wetting-induced shifts in hydrologic connectivity. We argue that models should treat soil moisture within a three-dimensional framework emphasizing hydrologic conduits for C and resource diffusion.

Carbon↗

Sulphate resistance of low‐clinker engineered cementitious composites examined by MicroXRF imaging

Abstract Engineered cementitious composites (ECC) are a class of high‐performing fibre‐reinforced cementitious materials recognised for their increased ductility and durability compared to conventional cement‐based materials, owing to their autogenously controlled tight crack widths, even when subjected to high strains. To reduce ECC's environmental impact, this research examines the use of a low‐clinker binder − limestone‐calcined clay cement (LC3) − as an alternative to portland cement (PC), along with fly ash to further reduce the clinker proportion and the embodied CO 2 of the formulations. In conventional concrete, LC3 hydrates to a denser microstructure resulting from the synergistic reaction between limestone and calcined clay. At the lower water contents typical of ECC and with the presence of fly ash, the influence of the binder composition on the microstructure is difficult to anticipate. To examine the influence of these compositional variables on microstructure, permeability and durability, the sulphate resistance of LC3‐based ECC is explored. Specifically, the ECC‐LC3 blends are designed with high clinker replacement rate of 75% by mass of binder and contain either conventional fly ash or reclaimed fly ash at 50% by mass of binder. Expansion of ECC‐LC3 samples subjected to standard sodium sulphate test conditions was measured up to 12 months and the depth of penetration of sulphates into the ECC‐LC3 of varying compositions was quantified using micro‐X‐Ray Fluorescence (microXRF) imaging and modelling. The expansion results show that the ECC‐LC3 formulations performed better than the PC samples and can provide adequate resistance to external sulphate attack, even when reclaimed fly ashes are used in place of the conventional ash. In addition, the shallow penetration of sulphate into these cementitious composites demonstrates the low diffusion coefficients values that were determined using the quantitative data from MicroXRF imaging.

Microscopy↗

Moving beyond the Aerosol Climatology of WRF-Solar: A Case Study over the North China Plain

Numerical weather prediction (NWP), when accessible, is a crucial input to short-term solar power forecasting. WRF-Solar, the first NWP model specifically designed for solar energy applications, has shown promising predictive capability. Nevertheless, few attempts have been made to investigate its performance under high aerosol loading, which attenuates incoming radiation significantly. The North China Plain is a polluted region due to industrialization, which constitutes a proper testbed for such investigation. Here, in this paper, aerosol direct radiative effect (DRE) on three surface shortwave radiation components (i.e., global, beam, and diffuse) during five heavy pollution episodes is studied within the WRF-Solar framework. Results show that WRF-Solar overestimates instantaneous beam radiation up to 795.3 W m -2 when the aerosol DRE is not considered. Although such overestimation can be partially offset by an underestimation of the diffuse radiation of about 194.5 W m -2 , the overestimation of the global radiation still reaches 160.2 W m -2 . This undesirable bias can be reduced when WRF-Solar is powered by Copernicus Atmosphere Monitoring Service (CAMS) aerosol forecasts, which then translates to accuracy improvements in photovoltaic (PV) power forecasts. This work also compares the forecast performance of the CAMS-powered WRF-Solar with that of the European Centre for Medium-Range Weather Forecasts model. Under high aerosol loading conditions, the irradiance forecast accuracy generated by WRF-Solar increased by 53.2% and the PV power forecast accuracy increased by 6.8%.

54 ENVIRONMENTAL SCIENCES↗

Geothermal reservoir characterization using distributed temperature sensing at Brady Geothermal Field, Nevada

Distributed temperature sensing (DTS) systems provide near real-time data collection that captures borehole spatiotemporal temperature dynamics. For this study, temperature data was collected in an observation well at an active geothermal site for a period of eight days under geothermal production conditions. Collected temperature data showcase the ability of DTS systems to detect changes to the location of the steam-water interface, visualize borehole temperature recovery — following injection of a coldwater “slug” — and identify anomalously warm and/or cool zones. The high sampling rate and spatial resolution of DTS data also shows borehole temperature dynamics that are not captured by traditional pressure-temperature survey tools. Inversion of thermal recovery data using a finite-difference heat-transfer model produces a thermal-diffusivity profile that is consistent with laboratorymeasured values and correlates with identified lithologic changes within the borehole. Used alone or in conjunction with complementary data sets, DTS systems are useful tools for developing a better understanding of both reservoir rock thermal properties as well as within and near borehole fluid movement.

15 GEOTHERMAL ENERGY↗

Representative Atmospheric Plume Development for Elevated Releases

An atmospheric explosion of a low-yield nuclear device will produce a large number of radioactive isotopes, some of which can be measured with airborne detection systems. However, properly equipped aircraft may not arrive in the region where an explosion occurred for a number of hours after the event. Atmospheric conditions will have caused the radioactive plume to move and diffuse before the aircraft arrives. The science behind predicting atmospheric plume movement has advanced enough that the location of the maximum concentrations in the plume can be determined reasonably accurately in real time, or near real time. Given the assumption that an aircraft can follow a plume, this study addresses the amount of atmospheric dilution expected to occur in a representative plume as a function of time past the release event. The approach models atmospheric transport of hypothetical releases from a single location for every day in a year using the publically available HYSPLIT code. The effective dilution factors for the point of maximum concentration in an elevated plume based on a release of a non-decaying, non-depositing tracer can vary by orders of magnitude depending on the day of the release, even for the same number of hours after the release event. However, the median of the dilution factors based on releases for 365 consecutive days at one site follows a power law relationship in time, as shown in Figure S-1. The relationship is good enough to provide a general rule of thumb for estimating typical future dilution factors in a plume starting at the same point. However, the coefficients of the power law function may vary for different release point locations. Radioactive decay causes the effective dilution factors to decrease more quickly with the time past the release event than the dilution factors based on a non-decaying tracer. An analytical expression for the dilution factors of isotopes with different half-lives can be developed given the power law expression for the non-decaying tracer. If the power-law equation for the median dilution factor, Df, based on a non-decaying tracer has the general form Df=a(×t)^(-b) for time t after the release event, then the equation has the form Df=e^(-λt)×a×t^(-b) for a radioactive isotope, where λ is the decay constant for the isotope.

54 ENVIRONMENTAL SCIENCES↗

Characterization of Creep-Fatigue Crack Growth in Alloy 709 and Prediction of Service Lives in Nuclear Reactor Components (Final Report)

This final report summarizes the activities and results that were obtained by the research team working on the project “Characterization of Creep-Fatigue Crack Growth in Alloy 709 and Prediction of Service Lives in Nuclear Reactor Components”, DOE-NEUP Project 15-8623. This project was mainly concerned with the creep-fatigue crack growth characterization of austenitic stainless-steel Fe-20Cr-25Ni (Alloy 709). For this purpose, an array of experimental and computational activities was performed. The main goal of these activities was to describe and predict the crack growth behavior in Alloy 709 under different conditions of applied loading and temperature. Alloy 709 specimens were tested to describe both the uniaxial creep behavior and creep-fatigue crack growth rates at temperatures between 550°C and 700°C. Microstructure and fracture surface analysis was performed using several experimental techniques, and the defining features of crack growth propagation were identified and described. Thermodynamic and diffusion simulations were performed to understand the microstructure evolution of Alloy 709 at elevated temperatures, and to establish an aging schedule relevant for in-service conditions over extended periods. Experimental testing was performed on both as-received and aged Alloy 709 specimens. Two computational methods were developed to predict crack growth rates and remaining lives in specimens loaded under either fatigue of creep-fatigue. The first computational technique was the finite element method, while the second was strip-yield modeling. Both methods relied on the prediction of crack growth rates using the concept of plasticity-induced crack closure. Crack-tip opening loads were computed and used to perform crack growth rate predictions. The computational predictions matched well with the experimental results. The strip-yield model was shown to be an excellent predictive tool not only for crack growth in austenitic SS Alloy 709, but also in 316 SS, martensitic Cr-Mo steels, Ni-based alloys and Al-based alloys. The main conclusion of this study is that Alloy 709 exhibits excellent creep crack growth resistance at temperatures between 550°C and 700°C. In fact, the crack growth resistance of this alloy appears to be higher than that of 316 SS in similar testing conditions. Moreover, when creep loading with hold (dwell) times of 60s or 600s were intercalated every fatigue cycle, the crack growth rates did not change significantly. The conclusion is that the tests performed revealed that the crack growth is insignificant during the creep portion of a loading cycle, as compared with the crack growth during the fatigue portion of the cycle. In general, Alloy 709 shows excellent resistance to creep deformations and creep crack growth up to 700°C.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Compatibility Of U 3 Si 2 Fuel With Zr, FeCrAl and SiC/SiC Based Cladding

U 3 Si 2 fuel with FeCrAl or SiC/SiC composite cladding is being considered as an advanced concept and accident tolerant fuel (ATF) system for light water reactors. Understanding their chemical compatibility under operational and accident conditions is fundamental to predicting behavior and fuel qualification. In the present work, experimental determination of interactions and thermodynamic and first principle calculation methods were employed to assess U 3 Si 2 compatibility with novel cladding of FeCrAl-alloys and SiC composites, and were contrasted with compatibility with Zircaloy-4. Diffusion-couple experiments were performed to simulate fuelcladding conditions and the samples were subjected to characterization by scanning electron microscopy and x-ray diffraction. Thermodynamic calculations, based in part on newly determined thermodynamic properties, were used to predict the equilibrium state of the studied systems and to interpret experimental results. First principles calculations using density functional theory were employed to obtain incorporation energies for carbon, iron and uranium in U 3 Si 2 , α- Fe, and SiC structures. Together, the experimental and modeling results have allowed evaluation of the nature of interactions in the fuel-cladding systems.

36 MATERIALS SCIENCE↗

Development and validation of capabilities to measure thermal properties of layered monolithic U-Mo alloy plate-type fuel

The uranium-molybdenum (U-Mo) alloy in a monolithic form has been proposed as one fuel design capable of converting some of the world’s highest power research reactors from the use of high enriched uranium (HEU) to low enriched uranium (LEU). One aspect of the fuel development and qualification process is to demonstrate appropriate understanding of thermal conductivity behavior of the fuel system as a function of temperature and expected irradiation conditions. The purpose of this paper is to verify and validate the functionality of equipment methods installed in hot cells for eventual measurements on irradiated uranium-molybdenum (U-Mo) monolithic fuel specimens, procedures to operate the equipment, and models to extract the desired thermal properties. The results presented here demonstrate the adequacy of the equipment, procedures and models that have been developed for this purpose based on measurements conducted on surrogate depleted uranium-molybdenum (DU-Mo) alloy samples containing a zirconium diffusion barrier and clad in aluminum alloy 6061 (AA6061). The results are in excellent agreement with thermal property data reported in the literature for similar U-Mo alloys as a function of temperature.

multi-layer↗

Non-isothermal Crystallization Kinetics and Magnetic Properties of FeCoNiCrZr Metallic Glass

In the present paper, the non-isothermal crystallization kinetics, crystallization process, and magnetic properties of (Fe{sub 0.52}Co{sub 0.30}Ni{sub 0.18}){sub 73}Cr{sub 17}Zr{sub 10} metallic glass were studied by means of differential scanning calorimetry (DSC), Kissinger equation, Ozawa equation, Kempen model, and vibrating sample magnetometer (VSM), respectively. According to the classical nucleation theory, the value of Avrami exponent (n), initially being larger than 4 which implies that phase transformation mechanism is typical for interface-controlled three-dimensional (3-D) growth with increasing nucleation rates, along with the development of crystallization, decreases to 1.5 ∼ 2, which implies that the growth mechanism changes to diffusion-controlled 3-D growth with decreasing nucleation rates; finally, the nucleation rates decrease to zero. Through proper annealing conditions, the saturation magnetization(M{sub s}) values of the annealed sample can reach up to 106.8A m{sup 2} kg{sup − 1}, the crystallization process of this alloy is amorphous → α-Fe(Co) + amorphous′ → α-Fe(Co) + Cr{sub 2}Ni{sub 3} + Fe{sub 3}Ni{sub 2} + Cr{sub 2}Zr + unidentified phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Estimation and modeling of coal pore accessibility using small angle neutron scattering

Gas diffusion in coal is controlled by nano-structure of the pores. The interconnectivity of pores not only determines the dynamics of gas transport in the coal matrix but also influences the mechanical strength. In this study, small angle neutron scattering (SANS) was employed to quantify pore accessibility for two coal samples, one of sub-bituminous rank and the other of anthracite rank. Moreover, a theoretical pore accessibility model was proposed based on scattering intensities under both vacuum and zero average contrast (ZAC) conditions. Our results show that scattering intensity decreases with increasing gas pressure using deuterated methane (CD 4 ) at low Q values for both coals. Pores smaller than 40 nm in radius are less accessible for anthracite than sub-bituminous coal. On the contrary, when the pore radius is larger than 40 nm, the pore accessibility of anthracite becomes larger than that of sub-bituminous coal. Only 20% of pores are accessible to CD 4 for anthracite and 37% for sub-bituminous coal, where the pore radius is 16 nm. For these two coals, pore accessibility and pore radius follows a power-law relationship.

01 COAL, LIGNITE, AND PEAT↗

Reducing dielectric loss and enhancing electrical insulation for multilayer polymer films by nanoconfined ion transport under high poling electric fields

High temperature polar polymers have demonstrated potential for good thermal stability and high dielectric constant at the same time. However, polarization of contaminated impurity ions in polar polymers, even at the ppm level, can significantly increase the dielectric loss at high temperature and low frequencies. One effective strategy to mitigate this problem is to multilayer them with a high temperature nonpolar dielectric polymer to confine impurity ion transport at the nanometer scale. In this study, confined ion transport in high temperature polycarbonate (HTPC)/poly(vinylidene fluoride) (PVDF) multilayer films under high AC electric fields was studied using a direct analytical simulation method. Different from the ion transport under low fields, the ion diffusion model failed to describe the ion transport under high electric fields. An exponential ion distribution profile, which was observed for the DC poling situation, was employed to implement the direct analytical simulation. Confined impurity ion transport under high AC fields was quantitatively understood. As the AC field increased, the mobile ion concentration decreased whereas the diffusion coefficient increased. The decrease of mobile ion concentration was explained by the blockage of impurity ions by the HTPC layers. This knowledge helped in the determination of optimal conditions to polarize impurity ions from the PVDF layers into the HTPC layers. After cooling below the glass transition temperature of HTPC, polarized impurity ions were locked inside the HTPC layers. As a result, increased discharge efficiency and enhanced electrical insulation (i.e., increased dielectric breakdown strength) were achieved for the polarized multilayer films.

36 MATERIALS SCIENCE↗

Bayesian soft X-ray tomography using non-stationary Gaussian Processes

In this study, a Bayesian based non-stationary Gaussian Process (GP) method for the inference of soft X-ray emissivity distribution along with its associated uncertainties has been developed. For the investigation of equilibrium condition and fast magnetohydrodynamic behaviors in nuclear fusion plasmas, it is of importance to infer, especially in the plasma center, spatially resolved soft X-ray profiles from a limited number of noisy line integral measurements. For this ill-posed inversion problem, Bayesian probability theory can provide a posterior probability distribution over all possible solutions under given model assumptions. Specifically, the use of a non-stationary GP to model the emission allows the model to adapt to the varying length scales of the underlying diffusion process. In contrast to other conventional methods, the prior regularization is realized in a probability form which enhances the capability of uncertainty analysis, in consequence, scientists who concern the reliability of their results will benefit from it. Under the assumption of normally distributed noise, the posterior distribution evaluated at a discrete number of points becomes a multivariate normal distribution whose mean and covariance are analytically available, making inversions and calculation of uncertainty fast. Additionally, the hyper-parameters embedded in the model assumption can be optimized through a Bayesian Occam's Razor formalism and thereby automatically adjust the model complexity. This method is shown to produce convincing reconstructions and good agreements with independently calculated results from the Maximum Entropy and Equilibrium-Based Iterative Tomography Algorithm methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Enhanced photocurrent density in graphene/Si based solar cell (GSSC) by optimizing active layer thickness

Thickness dependent photocurrent density in active layer of graphene/Si based solar cell has been investigated via analytical – simulation study. This report is a preliminary comparison of experimental and analytical investigation of graphene/Si based solar cell. Graphene sheet was interfaced with Si thin film forming heterojunction solar cell that was treated as a device model for photocurrent generator. Such current can be enhanced by optimizing active layer thickness and involving metal oxide as supporting layer to shift photons absorption. In this case there are two type of devices model with and without TiO{sub 2} in which the silicon thickness varied at 20 – 100 nm. All of them have examined and also compared with each other to obtain an optimum value. From this calculation it found that generated currents almost linear with thickness but there are saturated conditions that no more enhancements will be achieved. Furthermore TiO{sub 2} layer is effectively increases photon absorption but reducing device stability, maximum current is fluctuates enough. This may caused by the disturbance of excitons diffusion and resistivity inside each layer. Finally by controlling active layer thickness, it is quite useful to estimate optimization in order to develop the next solar cell devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temperature Quenching and Fluorescence Depolarization of Carbon Nanodots Obtained via Paraffin Pyrolysis

A temperature effect on fluorescence intensity and polarization of a colloidal system of carbon nanodots in glycerol under linearly polarized pumping conditions is studied. Nanodots are obtained via pyrolysis of paraffin in nanopores of a mesoporous silica. An increase in temperature leads to quenching of nanodots fluorescence, and activation energy of the quenching process is assessed. An experimental relationship between the linear fluorescence polarization and temperature is described by the Levshin–Perrin equation, which takes into account the rotational diffusion of luminescent particles (fluorophores) in a liquid matrix. The size of fluorophores is noticeably smaller than that of carbon nanodots according to the Levshin–Perrin model. A difference between the dimensions of the fluorophore and the nanodot indicates that the small atomic groups responsible for luminescence of the nanodot possess high segmental mobility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multischeme transport with upwinding for Rattlesnake code

This report conclude our efforts at INL to make Rattlesnake, the MOOSE-based radiation multigroup transport application, fully support the multischeme transport capability. Specifically, it details all five types of interface condition with upwinding method for coupling eight available discretization schemes within Rattlesnake. These discretization schemes include SN (discrete ordinates method), PN (spherical harmonics expansion method) and diffusion for angular discretization, and continuous finite element method (CFEM) and discontinuous finite element method (DFEM) for spatial discretization. The other type of method based on Lagrange multipliers for coupling schemes based on CFEM has been studied in the past. Besides simple test problems, a more realistic application of multischeme transport for TREAT (transient reactor test facility) experiment modeling are presented. Several enhancements are reviewed at the end which could potentially be useful in the future.

97 MATHEMATICS AND COMPUTING↗